#include <ParticleID.h>

Inheritance diagram for ParticleID:

Public Member Functions
	~ParticleID ()

void	init ()

bool	IsDedxInfoValid () const

bool	IsTofInfoValid () const

bool	IsTofEInfoValid () const

bool	IsTofQInfoValid () const

bool	IsTofCInfoValid () const

bool	IsTofCorrInfoValid () const

bool	IsEmcInfoValid () const

bool	IsMucInfoValid () const

bool	IsMrpcInfoValid () const

double	chiDedx (int n) const

double	chiTof (int n) const

double	chiTof1 (int n) const

double	chiTof2 (int n) const

double	chiTofE (int n) const

double	chiTofQ (int n) const

double	chiTofC (int n) const

double	chiTofCorr (int n) const

double	chiEmc (int n) const

double	chiMuc (int n) const

double	chiMrpc (int n) const

double	probDedx (int n) const

double	probTof (int n) const

double	probTofE (int n) const

double	probTofQ (int n) const

double	probTofC (int n) const

double	probTofCorr (int n) const

double	probEmc (int n) const

double	probMuc (int n) const

double	probMrpc (int n) const

double	pdfDedx (int n)

double	pdfTof (int n)

double	pdfTofE (int n)

double	pdfTofQ (int n)

double	pdfTofC (int n)

double	pdfTofCorr (int n)

double	pdfEmc (int n)

double	pdfMuc (int n)

double	pdfMrpc (int n)

int	getNhitCut () const

void	setNhitCut (const int nhitcuth=5)

void	setMethod (const int method)

void	usePidSys (const int pidsys)

void	identify (const int pidcase)

void	calculate ()

bool	IsPidInfoValid () const

double	chi (int n) const

double	prob (int n) const

int	ndof () const

double	pdf (int n) const

double	lhf (int n) const

double	probElectron () const

double	probMuon () const

double	probPion () const

double	probKaon () const

double	probProton () const

int	discard () const

double	pdfElectron () const

double	pdfMuon () const

double	pdfPion () const

double	pdfKaon () const

double	pdfProton () const

int	neuronStat () const

double	neuronVal () const

double	offset (int i)

double	sigma (int i)

Public Member Functions inherited from ParticleIDBase
	ParticleIDBase ()

virtual	~ParticleIDBase ()

EvtRecTrack *	PidTrk () const

void	setRecTrack (EvtRecTrack *trk)

double	chiMinCut () const

void	setChiMinCut (const double chi=4)

double	pdfMinSigmaCut () const

void	setPdfMinSigmaCut (const double pdf=4)

double	getRunNo () const

void	setRunNo (const double runh=8093)

double	p ()

double	pt ()

double	charge ()

double	xmass (int n)

double	velc ()

double	probCalculate (double chi2, int n)

double	pdfCalculate (double offset, double sigma)

double	interpolation (double x, double y, double x1)

double	pol2 (double x, double *par)

double	pol3 (double x, double *par)

double	pol4 (double x, double *par)

int	useDedx () const

int	useTof () const

int	useTof1 () const

int	useTof2 () const

int	useTofE () const

int	useTofQ () const

int	useTofC () const

int	useTofCorr () const

int	useEmc () const

int	useMuc () const

int	useMrpc () const

int	onlyPionKaon () const

int	onlyPionKaonProton () const

int	onlyPionKaonElectron () const

int	all () const

int	onlyElectron () const

int	onlyMuon () const

int	onlyPion () const

int	onlyKaon () const

int	onlyProton () const

int	methodLikelihood () const

int	methodProbability () const

int	methodNeuronNetwork () const

int	dedxValid () const

int	tofValid () const

int	tofeValid () const

int	tofqValid () const

int	tofcValid () const

int	tofcorrValid () const

int	emcValid () const

int	mucValid () const

int	MrpcValid () const

void	set_path (const char *s_path=0)

void	set_path (std::string s_path)

Static Public Member Functions
static ParticleID *	instance ()

Protected Member Functions
int	particleIDCalculation ()

int	LikelihoodCalculation ()

int	neuronPIDCalculation ()

Protected Member Functions inherited from ParticleIDBase

Additional Inherited Members
Static Protected Attributes inherited from ParticleIDBase
static std::string	path = ""

Detailed Description

Definition at line 37 of file ParticleID.h.

Constructor & Destructor Documentation

◆ ~ParticleID()

ParticleID::~ParticleID ( )

Definition at line 92 of file ParticleID.cxx.

                        {
   //  if(m_dedxpid) delete m_dedxpid;
   //  if(m_tof1pid) delete m_tof1pid;
   //  if(m_tof2pid) delete m_tof2pid;
   //  if(m_tofepid) delete m_tofepid;
   //  if(m_tofqpid) delete m_tofqpid;
   //  if(m_emcpid) delete m_emcpid;
}

Member Function Documentation

◆ calculate()

void ParticleID::calculate ( )

virtual

Implements ParticleIDBase.

Definition at line 101 of file ParticleID.cxx.

                           {
#ifdef BEAN
   cout << " please use ParticleID::calculate(run) ! " << endl;
   exit(1);
}
 
 
void ParticleID::calculate(int run)
{
#endif
   int nhitcutpid=getNhitCut();
 
#ifndef BEAN
   IDataProviderSvc* m_eventSvc;
   Gaudi::svcLocator()->service("EventDataSvc", m_eventSvc, true);
 
   SmartDataPtr<Event::EventHeader> eventHeaderpid(m_eventSvc,"/Event/EventHeader");
   int runpid=eventHeaderpid->runNumber();
   int eventpid=eventHeaderpid->eventNumber();
   //   cout<<"runpid="<<runpid<<endl;
   //   cout<<"eventpid="<<eventpid<<endl;
#else
   int runpid=run;
#endif
 
   EvtRecTrack* recTrk = PidTrk();
   //  int runnum=getRunNo();
   //  cout<<"runnum="<<runnum<<endl;
   // if user did not specify sub sys, sepcify the default value
   if(m_pidsys == 0) {
     m_pidsys = useDedx() | useTof() | useTofE() | useEmc() | useMuc() | useTofQ() | useTofC() | useTofCorr() | useMrpc();
   }
   // if user did not set the seperate case, set the default value
 
   if(m_pidcase == 0 ) {
      m_pidcase = all();
   }
   //dedx sys
   if(IsDedxInfoUsed()) {
      if(!m_dedxpid) m_dedxpid = DedxPID::instance();
      m_dedxpid->init();
      m_dedxpid->setRunNo(runpid);
      m_dedxpid->setNhitCutDx(nhitcutpid);
      m_dedxpid->setRecTrack(recTrk);
      m_dedxpid->setChiMinCut(chiMinCut());
      m_dedxpid->setPdfMinSigmaCut(pdfMinSigmaCut());
      m_dedxpid->calculate();
   }
 
   // tof1 and tof2 sys
   if(IsTofInfoUsed()|IsTof1InfoUsed()|IsTof2InfoUsed()|IsTofCInfoUsed())
   {
      if(IsTofCInfoUsed())
      {
         if(!m_tofcpid) m_tofcpid = TofCPID::instance();
         m_tofcpid->init();
         m_tofcpid->setRunNo(runpid);
         m_tofcpid->setRecTrack(recTrk);
         m_tofcpid->setChiMinCut(chiMinCut());
         m_tofcpid->setPdfMinSigmaCut(pdfMinSigmaCut());
         m_tofcpid->calculate();
      }
      else
      {
         if(!m_tofpid) m_tofpid = TofPID::instance();
         m_tofpid->init();
         m_tofpid->setRecTrack(recTrk);
         m_tofpid->setChiMinCut(chiMinCut());
         m_tofpid->setPdfMinSigmaCut(pdfMinSigmaCut());
         m_tofpid->calculate();
      }
 
   }
   // tof secondary correction sys
   if(IsTofCorrInfoUsed()) {
     if(!m_tofcorrpid) m_tofcorrpid = TofCorrPID::instance();
     m_tofcorrpid->setRunNo(runpid);
     m_tofcorrpid->init();
     m_tofcorrpid->setRecTrack(recTrk);
     m_tofcorrpid->setChiMinCut(chiMinCut());
     m_tofcorrpid->setPdfMinSigmaCut(pdfMinSigmaCut());
     m_tofcorrpid->calculate();
   }
 
   /*
     // tof1 sys
     if(IsTof1InfoUsed()){
       if(!m_tof1pid) m_tof1pid = Tof1PID::instance();
       m_tof1pid->init();
       m_tof1pid->setRecTrack(recTrk);
       m_tof1pid->setChiMinCut(4);
       m_tof1pid->setPdfMinSigmaCut(4);
       m_tof1pid->calculate();
     }
 
     // tof2 sys
     if(IsTof2InfoUsed()){
       if(!m_tof2pid) m_tof2pid = Tof2PID::instance();
       m_tof2pid->init();
       m_tof2pid->setRecTrack(recTrk);
       m_tof2pid->setChiMinCut(4);
       m_tof2pid->setPdfMinSigmaCut(4);
       m_tof2pid->calculate();
     }
 
   */
   // tofq sys
   if(IsTofQInfoUsed()) {
      if(!m_tofqpid) m_tofqpid = TofQPID::instance();
      m_tofqpid->init();
      m_tofqpid->setRecTrack(recTrk);
      m_tofqpid->setChiMinCut(chiMinCut());
      m_tofqpid->calculate();
   }
 
   // endcap tof sys
   if(IsTofEInfoUsed()&&(!IsTofCorrInfoUsed())) {
      if(!m_tofepid) m_tofepid = TofEPID::instance();
      m_tofepid->init();
      m_tofepid->setRecTrack(recTrk);
      m_tofepid->setChiMinCut(chiMinCut());
      m_tofepid->setPdfMinSigmaCut(pdfMinSigmaCut());
      m_tofepid->calculate();
   }
   // emc sys
   if(IsEmcInfoUsed()) {
      if(!m_emcpid) m_emcpid = EmcPID::instance();
      m_emcpid->init();
      m_emcpid->setRecTrack(recTrk);
      m_emcpid->setChiMinCut(chiMinCut());
      m_emcpid->calculate();
   }
 
   // muc sys
   if(IsMucInfoUsed()) {
      if(!m_mucpid) m_mucpid = MucPID::instance();
      m_mucpid->init();
      m_mucpid->setRecTrack(recTrk);
      m_mucpid->setChiMinCut(chiMinCut());
      m_mucpid->calculate();
   }
   //mrpc sys
   if(IsMrpcInfoUsed()) {
      if(!m_mrpcpid) m_mrpcpid = MrpcPID::instance();
      m_mrpcpid->init();
      m_mrpcpid->setRecTrack(recTrk);
      m_mrpcpid->setChiMinCut(chiMinCut());
      m_mrpcpid->setPdfMinSigmaCut(pdfMinSigmaCut());
      m_mrpcpid->calculate();
   }
 
 
 
   // probability method
   if((m_method & methodProbability()) == methodProbability())
      if(particleIDCalculation() < 0) m_ndof = 0;
   //    std::cout<<"m_ndof="<<m_ndof<<std::endl;
 
   //likelihood method
   if((m_method & methodLikelihood()) == methodLikelihood())
      if(LikelihoodCalculation() < 0) m_discard =0;
   // neuron network
   if((m_method & methodNeuronNetwork()) == methodNeuronNetwork())
      //    m_neuronPid = neuronPIDCalculation();
      if(LikelihoodCalculation() < 0) m_discard =0;
 
}

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), Rhopi::execute(), LocalKaonSelector::operator()(), and LocalPionSelector::operator()().

◆ chi()

double ParticleID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 120 of file ParticleID.h.

120{return m_chisq[n];}

n

const Int_t n

Definition DataBase/tau_mode.c:65

◆ chiDedx()

double ParticleID::chiDedx ( int n ) const

inline

Referenced by DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), Rhopi::execute(), and LikelihoodCalculation().

◆ chiEmc()

double ParticleID::chiEmc ( int n ) const

inline

◆ chiMrpc()

double ParticleID::chiMrpc ( int n ) const

inline

◆ chiMuc()

double ParticleID::chiMuc ( int n ) const

inline

◆ chiTof()

double ParticleID::chiTof ( int n ) const

inline

Referenced by LikelihoodCalculation().

◆ chiTof1()

double ParticleID::chiTof1 ( int n ) const

inline

Referenced by DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), and Rhopi::execute().

◆ chiTof2()

double ParticleID::chiTof2 ( int n ) const

inline

Referenced by DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), and Rhopi::execute().

◆ chiTofC()

double ParticleID::chiTofC ( int n ) const

inline

◆ chiTofCorr()

double ParticleID::chiTofCorr ( int n ) const

inline

Referenced by LikelihoodCalculation().

◆ chiTofE()

double ParticleID::chiTofE ( int n ) const

inline

Referenced by LikelihoodCalculation().

◆ chiTofQ()

double ParticleID::chiTofQ ( int n ) const

inline

◆ discard()

int ParticleID::discard ( ) const

inline

Definition at line 133 of file ParticleID.h.

133{return m_discard;}

◆ getNhitCut()

int ParticleID::getNhitCut ( ) const

inline

Definition at line 95 of file ParticleID.h.

95{return m_nhitcut;}

◆ identify()

void ParticleID::identify ( const int pidcase )

inline

Definition at line 110 of file ParticleID.h.

110{m_pidcase = (m_pidcase | pidcase);}

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), Rhopi::execute(), LocalKaonSelector::operator()(), and LocalPionSelector::operator()().

◆ init()

void ParticleID::init ( )

virtual

Implements ParticleIDBase.

Definition at line 27 of file ParticleID.cxx.

                      {
 
   if(IsDedxInfoUsed()) {
      if(!m_dedxpid) m_dedxpid = DedxPID::instance();
      m_dedxpid->init();
   }
 
   if(IsTofInfoUsed()|IsTof1InfoUsed()|IsTof2InfoUsed()) {
      if(!m_tofpid) m_tofpid = TofPID::instance();
      m_tofpid->init();
   }
 
   if(IsTofCorrInfoUsed()) {
      if(!m_tofcorrpid) m_tofcorrpid = TofCorrPID::instance();
      //      m_tofcorrpid->init();
   }
 
   if(IsEmcInfoUsed()) {
      if(!m_emcpid) m_emcpid = EmcPID::instance();
      m_emcpid->init();
   }
 
   if(IsMucInfoUsed()) {
      if(!m_mucpid) m_mucpid = MucPID::instance();
      m_mucpid->init();
   }
 
 
   if(IsMrpcInfoUsed()) {
      if(!m_mrpcpid) m_mrpcpid = MrpcPID::instance();
      m_mrpcpid->init();
   }
 
 
   // global info.
   m_pidsys = 0;
   m_pidcase = 0;
   m_method = 0;
   m_TotalLikelihood =0;
   m_discard =1;
   // probability
   m_ndof = 0;
   m_nhitcut=5;
   for(int i = 0; i < 4; i++) {
      m_chisq[i] = 9999.;
      m_prob[i] = -1.0;
   }
}

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), Rhopi::execute(), LocalKaonSelector::operator()(), and LocalPionSelector::operator()().

◆ instance()

ParticleID * ParticleID::instance ( )

static

Definition at line 22 of file ParticleID.cxx.

                                  {
   if(!m_pointer) m_pointer = new ParticleID();
   return m_pointer;
}

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), Rhopi::execute(), LocalKaonSelector::operator()(), and LocalPionSelector::operator()().

◆ IsDedxInfoValid()

bool ParticleID::IsDedxInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsEmcInfoValid()

bool ParticleID::IsEmcInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsMrpcInfoValid()

bool ParticleID::IsMrpcInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsMucInfoValid()

bool ParticleID::IsMucInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsPidInfoValid()

bool ParticleID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), Rhopi::execute(), and neuronStat().

◆ IsTofCInfoValid()

bool ParticleID::IsTofCInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsTofCorrInfoValid()

bool ParticleID::IsTofCorrInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsTofEInfoValid()

bool ParticleID::IsTofEInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsTofInfoValid()

bool ParticleID::IsTofInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsTofQInfoValid()

bool ParticleID::IsTofQInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ lhf()

double ParticleID::lhf ( int n ) const

inline

Definition at line 126 of file ParticleID.h.

126{return m_likelihoodfraction[n];}

Referenced by neuronVal().

◆ LikelihoodCalculation()

int ParticleID::LikelihoodCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 460 of file ParticleID.cxx.

                                      {
   int irc = -1;
 
   bool valid = IsDedxInfoValid() || IsTofInfoValid() || IsTofEInfoValid() || IsTofQInfoValid() || IsTofCInfoValid() || IsTofCorrInfoValid() || IsEmcInfoValid()||IsMucInfoValid() || IsMrpcInfoValid();
   if(!valid) return irc;
   double pdf[5];
   bool pidcase[5];
   for(int i = 0; i < 5; i++) {
      pdf[i] = 1;
      pidcase[i] = false;
   }
 
   if((m_pidcase & onlyElectron()) == onlyElectron()) pidcase[0] = true;
   if((m_pidcase & onlyMuon()) == onlyMuon()) pidcase[1] = true;
   if((m_pidcase & onlyPion()) == onlyPion()) pidcase[2] = true;
   if((m_pidcase & onlyKaon()) == onlyKaon()) pidcase[3] = true;
   if((m_pidcase & onlyProton()) == onlyProton()) pidcase[4] = true;
 
   for(int i = 0; i < 5; i++) {
      if(pidcase[i]==0)
         pdf[i]=0;
   }
 
   //
   //  dEdx PID
   //
   if(IsDedxInfoUsed()) {
      if(IsDedxInfoValid()) {
         bool okpid = false;
         for(int i = 0; i < 5; i++) {
            if(pidcase[i] && pdfCalculate(chiDedx(i),1) > pdfCalculate(m_dedxpid->pdfMinSigmaCut(),1.5))
               if(!okpid) okpid = true;
         }
         if(okpid) {
            m_ndof++;
            for(int i = 0; i < 5; i++) {
               pdf[i] *= pdfDedx(i);
            }
         }
      } // dE/dx
   }
 
 
   //
   //  Barrel TOF
   //
   if(IsTofInfoUsed()|IsTof1InfoUsed()|IsTof2InfoUsed()|IsTofCInfoUsed())
   {  if(IsTofCInfoUsed())
      {
 
         if(IsTofCInfoValid()) {
            bool okpid = false;
            for(int i = 0; i < 5; i++) {
               if(pidcase[i] && pdfCalculate(chiTof(i),1) > pdfCalculate(m_tofcpid->pdfMinSigmaCut(),1.5))
                  if(!okpid) okpid = true;
            }
            if(okpid) {
               m_ndof++;
               for(int i = 0; i < 5; i++) {
                  pdf[i] *= pdfTofC(i);
               }
            }
         }  // TOF
      }
 
      else {
         if(IsTofInfoValid()) {
            bool okpid = false;
            for(int i = 0; i < 5; i++) {
               if(pidcase[i] && pdfCalculate(chiTof(i),1) > pdfCalculate(m_tofpid->pdfMinSigmaCut(),1.5))
                  if(!okpid) okpid = true;
            }
            if(okpid) {
               m_ndof++;
               for(int i = 0; i < 5; i++) {
                  pdf[i] *= pdfTof(i);
               }
            }
         }  // TOF
 
 
 
         //
         //  EndCap Tof
         //
 
         if(IsTofEInfoUsed()) {
            if(IsTofEInfoValid()) {
               bool okpid = false;
               for(int i = 0; i < 5; i++) {
                  if(pidcase[i]&& pdfCalculate(chiTofE(i),1) > pdfCalculate(m_tofepid->pdfMinSigmaCut(),1.5))
                     if(!okpid) okpid = true;
               }
               if(okpid) {
                  //      m_ndof++;
                  //      for(int i = 0; i < 5; i++) pdf[i] *= pdfTofE(i);
               }
            } // EndCap TOF
         }
      }
 
   }
 
   // Secondary TOF correction
   if(IsTofCorrInfoUsed()) {
     if(IsTofCorrInfoValid()) {
       bool okpid = false;
       for(int i = 0; i < 5; i++) {
     if(pidcase[i] && pdfCalculate(chiTofCorr(i),1) > pdfCalculate(m_tofcorrpid->pdfMinSigmaCut(),1.5))
       if(!okpid) okpid = true;
       }
       if(okpid) {
     m_ndof++;
     for(int i = 0; i < 5; i++) {
       pdf[i] *= pdfTofCorr(i);
     }
       }
     }
   }
 
 
   //
   //  Barrel TOF Q
   //
 
   if(IsTofQInfoValid()) {
      bool okpid = false;
      for(int i = 0; i < 5; i++) {
         if(pidcase[i])
            if(!okpid) okpid = true;
      }
      if(okpid) {
         //    m_ndof++;
         for(int i = 0; i < 5; i++) pdf[i] *= pdfTofQ(i);
      }
   } // TofQ
 
   //
   //  Emc PID
   //
   if(IsEmcInfoUsed()) {
      if(IsEmcInfoValid()) {
         bool okpid = false;
         for(int i = 0; i < 5; i++) {
            if(pidcase[i]&&pdfEmc(i)>0)
               if(!okpid) okpid = true;
         }
         if(okpid) {
            m_ndof++;
            for(int i = 0; i < 5; i++) {
               pdf[i] *= pdfEmc(i);
            }
         } // Emc Pid
      }
   }
   if(IsMucInfoUsed()) {
      if(IsMucInfoValid()) {
         bool okpid = false;
         for(int i = 0; i < 5; i++) {
            if(pidcase[i]&&pdfMuc(i)>0)
               if(!okpid) okpid = true;
         }
         if(okpid) {
            m_ndof++;
            for(int i = 0; i < 5; i++) {
               pdf[i] *= pdfMuc(i);
            }
         }
      } // Emc Pid
   }
   if(IsMrpcInfoUsed()) {
      if(IsMrpcInfoValid()) {
         bool okpid = false;
         for(int i = 0; i < 5; i++) {
            if(pidcase[i]&&pdfMrpc(i)>0)
               if(!okpid) okpid = true;
         }
         if(okpid) {
            m_ndof++;
            for(int i = 0; i < 5; i++) {
               pdf[i] *= pdfMrpc(i);
            }
         }
      } // Emc Pid
   }
 
 
 
 
   if(m_ndof <= 0) return irc;
   for(int i = 0; i < 5; i++) {
      m_pdf[i] = pdf[i];
      m_TotalLikelihood += pdf[i];
   }
   for(int i = 0; i < 5; i++) {
      m_likelihoodfraction[i] = pdf[i]/m_TotalLikelihood;
   }
 
 
   irc = 0;
   return irc;
}

◆ ndof()

int ParticleID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 122 of file ParticleID.h.

122{return m_ndof;}

◆ neuronPIDCalculation()

int ParticleID::neuronPIDCalculation ( )

protected

◆ neuronStat()

int ParticleID::neuronStat ( ) const

inline

Definition at line 151 of file ParticleID.h.

151{return IsPidInfoValid();}

ParticleID::IsPidInfoValid

bool IsPidInfoValid() const

◆ neuronVal()

double ParticleID::neuronVal ( ) const

inline

Definition at line 152 of file ParticleID.h.

                           {
    if(lhf(0)>0.02) return 1.0;
    if(lhf(1)>0.02) return 8.5;
    if(lhf(2)>0.02) return 2.0;
    if(lhf(3)>0.02) return 4.0;
    if(lhf(4)>0.02) return 5.0;
    }

◆ offset()

double ParticleID::offset ( int i )

inline

Definition at line 160 of file ParticleID.h.

                       {
    return m_tofcorrpid->offset(i);
  }

◆ particleIDCalculation()

int ParticleID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 269 of file ParticleID.cxx.

                                       {
   int irc = -1;
   bool valid = IsDedxInfoValid() || IsTofInfoValid()||IsTofEInfoValid()
     || IsTofQInfoValid() || IsEmcInfoValid() || IsMucInfoValid()
     || IsTofCInfoValid() || IsTofCorrInfoValid() || IsMrpcInfoValid();
 
   if(!valid) return irc;
 
   double chisq[5];
   bool pidcase[5];
   for(int i = 0; i < 5; i++) {
      chisq[i] = 0;
      pidcase[i] = false;
   }
 
   if((m_pidcase & onlyElectron()) == onlyElectron()) pidcase[0] = true;
   if((m_pidcase & onlyMuon()) == onlyMuon()) pidcase[1] = true;
   if((m_pidcase & onlyPion()) == onlyPion()) pidcase[2] = true;
   if((m_pidcase & onlyKaon()) == onlyKaon()) pidcase[3] = true;
   if((m_pidcase & onlyProton()) == onlyProton()) pidcase[4] = true;
 
   //
   //  dEdx PID
   //
   if(IsDedxInfoUsed()) {
      if(IsDedxInfoValid()) {
         bool okpid = false;
         for(int i = 0; i < 5; i++) {
            if(pidcase[i] && (fabs(chiDedx(i)) < m_dedxpid->chiMinCut()))
               if(!okpid) okpid = true;
         }
         if(okpid) {
            m_ndof++;
            for(int i = 0; i < 5; i++) chisq[i] += chiDedx(i)*chiDedx(i);
 
 
         }
      } // dE/dx
   }
   //
   //  Barrel TOF
   //
 
   if(IsTofInfoUsed()|IsTof1InfoUsed()|IsTof2InfoUsed() | IsTofCInfoUsed())
   {  if(IsTofCInfoUsed())
      {
         if(IsTofCInfoValid()) {
            bool okpid = false;
            for(int i = 0; i < 5; i++) {
               if(pidcase[i] && (fabs(chiTofC(i)) < m_tofcpid->chiMinCut()))
                  if(!okpid) okpid = true;
            }
            if(okpid) {
               m_ndof++;
               for(int i = 0; i < 5; i++) chisq[i] += chiTofC(i)*chiTofC(i);
            }
         }  // TOF1
      }
      else {
         if(IsTofInfoValid()) {
            bool okpid = false;
            for(int i = 0; i < 5; i++) {
               if(pidcase[i] && (fabs(chiTof(i)) < m_tofpid->chiMinCut()))
                  if(!okpid) okpid = true;
            }
            if(okpid) {
               m_ndof++;
               for(int i = 0; i < 5; i++) chisq[i] += chiTof(i)*chiTof(i);
            }
         }  // TOF1
 
 
         //
         //  EndCap Tof
         //
 
         if(IsTofEInfoUsed()) {
            if(IsTofEInfoValid()) {
               bool okpid = false;
               for(int i = 0; i < 5; i++) {
                  if(pidcase[i] && (fabs(chiTofE(i)) < m_tofepid->chiMinCut()))
                     if(!okpid) okpid = true;
               }
               if(okpid) {
                  m_ndof++;
                  for(int i = 0; i < 5; i++) chisq[i] += chiTofE(i)*chiTofE(i);
               }
            } // EndCap TOF
         }
 
      }
   }
 
   // Secondary TOF correction
   if(IsTofCorrInfoUsed()) {
     if(IsTofCorrInfoValid()) {
       bool okpid = false;
       for(int i = 0; i < 5; i++) {
     //  if(pidcase[i] && (fabs(chiTofCorr(i)) < m_tofcorrpid->chiMinCut()))
     if(pidcase[i] && ( chiTofCorr(i) < 6.0 && chiTofCorr(i) > -4.0 ) )
       if(!okpid) okpid = true;
       }
       if(okpid) {
     m_ndof++;
     for(int i = 0; i < 5; i++) chisq[i] += chiTofCorr(i)*chiTofCorr(i);
       }
     }
   }
 
   //
   //  Barrel TOF Q
   //
 
   if(IsTofQInfoUsed()) {
      if(IsTofQInfoValid()) {
         bool okpid = false;
         for(int i = 0; i < 5; i++) {
            if(pidcase[i] && (fabs(chiTofQ(i)) < m_tofqpid->chiMinCut()))
               if(!okpid) okpid = true;
         }
         if(okpid) {
            m_ndof++;
            for(int i = 0; i < 5; i++) chisq[i] += chiTofQ(i)*chiTofQ(i);
         }
      } // TofQ
   }
 
   // Muc Pid
   if(IsMucInfoUsed()) {
      if(IsMucInfoValid()) {
         bool okpid = false;
         for(int i = 0; i < 5; i++) {
            if(pidcase[i] && (fabs(chiMuc(i)) < m_mucpid->chiMinCut()))
               if(!okpid) okpid = true;
         }
         if(okpid) {
            m_ndof++;
            for(int i = 0; i < 5; i++) chisq[i] += chiMuc(i)*chiMuc(i);
         }
      } // Muc Pid
   }
 
 
   //  Emc PID
   if(IsEmcInfoUsed()) {
      if(IsEmcInfoValid()) {
         bool okpid = false;
         for(int i = 0; i < 5; i++) {
            if(pidcase[i] && (fabs(chiEmc(i)) < m_emcpid->chiMinCut()))
               if(!okpid) okpid = true;
         }
         if(okpid) {
            m_ndof++;
            for(int i = 0; i < 5; i++) chisq[i] += chiEmc(i)*chiEmc(i);
         }
      } // Emc Pid
   }
 
   //  Mrpc PID
   if(IsMrpcInfoUsed()) {
      if(IsMrpcInfoValid()) {
         bool okpid = false;
         for(int i = 0; i < 5; i++) {
            if(pidcase[i] && (fabs(chiMrpc(i)) < m_mrpcpid->chiMinCut()))
               if(!okpid) okpid = true;
         }
         if(okpid) {
            m_ndof++;
            for(int i = 0; i < 5; i++) chisq[i] += chiMrpc(i)*chiMrpc(i);
         }
      } // Mrpc Pid
   }
 
 
   if(m_ndof <= 0) return irc;
 
 
   for(int i = 0; i < 5; i++) {
      m_chisq[i] = chisq[i];
      m_prob[i] = probCalculate(chisq[i], m_ndof);
   }
 
 
   irc = 0;
   return irc;
}

◆ pdf()

double ParticleID::pdf ( int n ) const

inline

Definition at line 125 of file ParticleID.h.

125{ return m_pdf[n];}

Referenced by LikelihoodCalculation().

◆ pdfDedx()

double ParticleID::pdfDedx ( int n )

inline

Referenced by LikelihoodCalculation().

◆ pdfElectron()

double ParticleID::pdfElectron ( ) const

inline

Definition at line 136 of file ParticleID.h.

136{return m_pdf[0];}

◆ pdfEmc()

double ParticleID::pdfEmc ( int n )

inline

Referenced by LikelihoodCalculation().

◆ pdfKaon()

double ParticleID::pdfKaon ( ) const

inline

Definition at line 139 of file ParticleID.h.

139{return m_pdf[3];}

◆ pdfMrpc()

double ParticleID::pdfMrpc ( int n )

inline

Referenced by LikelihoodCalculation().

◆ pdfMuc()

double ParticleID::pdfMuc ( int n )

inline

Referenced by LikelihoodCalculation().

◆ pdfMuon()

double ParticleID::pdfMuon ( ) const

inline

Definition at line 137 of file ParticleID.h.

137{return m_pdf[1];}

◆ pdfPion()

double ParticleID::pdfPion ( ) const

inline

Definition at line 138 of file ParticleID.h.

138{return m_pdf[2];}

◆ pdfProton()

double ParticleID::pdfProton ( ) const

inline

Definition at line 140 of file ParticleID.h.

140{return m_pdf[4];}

◆ pdfTof()

double ParticleID::pdfTof ( int n )

inline

Referenced by LikelihoodCalculation().

◆ pdfTofC()

double ParticleID::pdfTofC ( int n )

inline

Referenced by LikelihoodCalculation().

◆ pdfTofCorr()

double ParticleID::pdfTofCorr ( int n )

inline

Referenced by LikelihoodCalculation().

◆ pdfTofE()

double ParticleID::pdfTofE ( int n )

inline

◆ pdfTofQ()

double ParticleID::pdfTofQ ( int n )

inline

Referenced by LikelihoodCalculation().

◆ prob()

double ParticleID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 121 of file ParticleID.h.

121{ return m_prob[n];}

Referenced by DQAKsKpi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), and DQASelHadron::execute().

◆ probDedx()

double ParticleID::probDedx ( int n ) const

inline

◆ probElectron()

double ParticleID::probElectron ( ) const

inline

Definition at line 128 of file ParticleID.h.

128{return m_prob[0];}

Referenced by BeamParams::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), and Gam4pikp::execute().

◆ probEmc()

double ParticleID::probEmc ( int n ) const

inline

◆ probKaon()

double ParticleID::probKaon ( ) const

inline

Definition at line 131 of file ParticleID.h.

131{return m_prob[3];}

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), inclks::execute(), inclkstar::execute(), inclphi::execute(), LocalKaonSelector::operator()(), and LocalPionSelector::operator()().

◆ probMrpc()

double ParticleID::probMrpc ( int n ) const

inline

◆ probMuc()

double ParticleID::probMuc ( int n ) const

inline

◆ probMuon()

double ParticleID::probMuon ( ) const

inline

Definition at line 129 of file ParticleID.h.

129{return m_prob[1];}

Referenced by BeamParams::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), and Gam4pikp::execute().

◆ probPion()

double ParticleID::probPion ( ) const

inline

Definition at line 130 of file ParticleID.h.

130{return m_prob[2];}

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), Rhopi::execute(), LocalKaonSelector::operator()(), and LocalPionSelector::operator()().

◆ probProton()

double ParticleID::probProton ( ) const

inline

Definition at line 132 of file ParticleID.h.

132{return m_prob[4];}

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), and incllambda::execute().

◆ probTof()

double ParticleID::probTof ( int n ) const

inline

◆ probTofC()

double ParticleID::probTofC ( int n ) const

inline

◆ probTofCorr()

double ParticleID::probTofCorr ( int n ) const

inline

◆ probTofE()

double ParticleID::probTofE ( int n ) const

inline

◆ probTofQ()

double ParticleID::probTofQ ( int n ) const

inline

◆ setMethod()

void ParticleID::setMethod ( const int method )

inline

Definition at line 101 of file ParticleID.h.

101{m_method = (m_method | method);}

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), Rhopi::execute(), LocalKaonSelector::operator()(), and LocalPionSelector::operator()().

◆ setNhitCut()

void ParticleID::setNhitCut ( const int nhitcuth = 5 )

inline

Definition at line 96 of file ParticleID.h.

96{m_nhitcut = nhitcuth;}

◆ sigma()

double ParticleID::sigma ( int i )

inline

Definition at line 163 of file ParticleID.h.

                      {
    return m_tofcorrpid->sigma(i);
  }

◆ usePidSys()

void ParticleID::usePidSys ( const int pidsys )

inline

Definition at line 104 of file ParticleID.h.

104{m_pidsys = (m_pidsys | pidsys);}

Referenced by BeamParams::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), Gam4pikp::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), Rhopi::execute(), LocalKaonSelector::operator()(), and LocalPionSelector::operator()().

The documentation for this class was generated from the following files:

Public Member Functions

Static Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ~ParticleID()

Member Function Documentation

◆ calculate()

◆ chi()

◆ chiDedx()

◆ chiEmc()

◆ chiMrpc()

◆ chiMuc()

◆ chiTof()

◆ chiTof1()

◆ chiTof2()

◆ chiTofC()

◆ chiTofCorr()

◆ chiTofE()

◆ chiTofQ()

◆ discard()

◆ getNhitCut()

◆ identify()

◆ init()

◆ instance()

◆ IsDedxInfoValid()

◆ IsEmcInfoValid()

◆ IsMrpcInfoValid()

◆ IsMucInfoValid()

◆ IsPidInfoValid()

◆ IsTofCInfoValid()

◆ IsTofCorrInfoValid()

◆ IsTofEInfoValid()

◆ IsTofInfoValid()

◆ IsTofQInfoValid()

◆ lhf()

◆ LikelihoodCalculation()

◆ ndof()

◆ neuronPIDCalculation()

◆ neuronStat()

◆ neuronVal()

◆ offset()

◆ particleIDCalculation()

◆ pdf()

◆ pdfDedx()

◆ pdfElectron()

◆ pdfEmc()

◆ pdfKaon()

◆ pdfMrpc()

◆ pdfMuc()

◆ pdfMuon()

◆ pdfPion()

◆ pdfProton()

◆ pdfTof()

◆ pdfTofC()

◆ pdfTofCorr()

◆ pdfTofE()

◆ pdfTofQ()

◆ prob()

◆ probDedx()

◆ probElectron()

◆ probEmc()

◆ probKaon()

◆ probMrpc()

◆ probMuc()

◆ probMuon()

◆ probPion()

◆ probProton()

◆ probTof()

◆ probTofC()

◆ probTofCorr()

◆ probTofE()

◆ probTofQ()

◆ setMethod()

◆ setNhitCut()

◆ sigma()

◆ usePidSys()