TofQPID.cxx

Go to the documentation of this file.

#include <cmath>
 
#include "ParticleID/TofQPID.h"
 
#ifndef BEAN
#  include "EvtRecEvent/EvtRecTrack.h"
#  include "MdcRecEvent/RecMdcTrack.h"
#  include "TofRecEvent/RecTofTrack.h"
#endif
 
TofQPID* TofQPID::m_pointer = 0;
 
TofQPID* TofQPID::instance() {
  if ( !m_pointer ) m_pointer = new TofQPID();
  return m_pointer;
}
 
TofQPID::TofQPID() : ParticleIDBase() { ; }
 
void TofQPID::init() {
  for ( int i = 0; i < 5; i++ )
  {
    m_chi[i]  = 99.0;
    m_prob[i] = -1.0;
  }
  m_chimin    = 99.;
  m_ndof      = 0;
  m_hitstatus = -2;
  //  m_q1 = -1;
  //  m_q2 = -1;
}
 
void TofQPID::calculate() {
  if ( particleIDCalculation() == 0 ) m_ndof = 1;
}
int TofQPID::particleIDCalculation() {
  int          irc    = -1;
  EvtRecTrack* recTrk = PidTrk();
  if ( !( recTrk->isMdcTrackValid() ) ) return irc;
  //   RecMdcTrack* mdcTrk = recTrk->mdcTrack();
 
  //   double ptrk = mdcTrk->p();
  //   double cost = cos(mdcTrk->theta());
 
  if ( !( recTrk->isTofTrackValid() ) ) return irc;
  // RecTofTrack* tofTrk = recTrk->tofTrack();
 
  //
  // two layer hit in barrel TOF are required
  //
  //  if(tofTrk->getPart() != 1) return irc;
  // double ph1 = tofTrk->getPh1();
  // double ph2 = tofTrk->getPh2();
  // if(ph1 <= 0 || ph2 <= 0) return irc;
  //
  //  Tof Q calibration is needed
  //
  //  m_hitstatus = tofTrk->getHitStatus();
  // m_hitstatus = 1;
  // if(m_hitstatus == 0 || m_hitstatus == 1) m_q1 = tofTrk->getQ1();
  // if(m_hitstatus == 0 || m_hitstatus == 2) m_q2 = tofTrk->getQ2();
 
  // m_q1 = 1.0;
  // m_q2 = 1.0;
  m_ndof = 1;
  irc    = 0;
  return irc;
}