#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

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	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	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)

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

int	ipmt ()

double	dt (int ipar, int ipmt)

double	dtCorr (int ipar, int ipmt)

double	sigCorr (int ipar, int ipmt)

double	chiCorr (int ipar, int ipmt)

double	offset (int ipar)

double	sigma (int ipar)

Public Member Functions inherited from ParticleIDBase
	ParticleIDBase ()

virtual	~ParticleIDBase ()

virtual void	init ()=0

virtual void	calculate ()=0

virtual bool	IsPidInfoValid () const =0

virtual double	chi (int n) const =0

virtual double	prob (int n) const =0

virtual int	ndof () const =0

EvtRecTrack *	PidTrk () const

void	setRecTrack (EvtRecTrack *trk)

double	chiMinCut () const

void	setChiMinCut (const double chi=4)

double	chiMaxCut () const

void	setChiMaxCut (const double chi=6)

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	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

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 ()

virtual int	particleIDCalculation ()=0

virtual int	LikelihoodCalculation ()=0

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

Detailed Description

Definition at line 36 of file ParticleID.h.

Constructor & Destructor Documentation

◆ ~ParticleID()

ParticleID::~ParticleID ( )

Definition at line 88 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 97 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();
   }
   // 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->setChiMaxCut(chiMaxCut());
     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();
   }
   // 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 Gam4pikp::execute(), Rhopi::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQAPi2p2::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), DTagSetAlg::execute(), LTagSetAlg::execute(), TagSetAlg::execute(), XYZTagSetAlg::execute(), BeamParams::execute(), LocalKaonSelector::operator()(), LocalPionSelector::operator()(), LocalLambdaSelector::operator()(), and LocalProtonSelector::operator()().

◆ chi()

double ParticleID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 113 of file ParticleID.h.

113{return m_chisq[n];}

◆ chiCorr()

double ParticleID::chiCorr	(	int	ipar,
		int	ipmt
	)

inline

Definition at line 165 of file ParticleID.h.

                                     {
    return m_tofcorrpid->chiCorr(ipar, ipmt);
  }

◆ chiDedx()

double ParticleID::chiDedx ( int n ) const

inline

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

◆ chiEmc()

double ParticleID::chiEmc ( 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 Gam4pikp::execute(), Rhopi::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), and DQASelHadron::execute().

◆ chiTof2()

double ParticleID::chiTof2 ( int n ) const

inline

Referenced by Gam4pikp::execute(), Rhopi::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), and DQASelHadron::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 126 of file ParticleID.h.

126{return m_discard;}

◆ dt()

double ParticleID::dt	(	int	ipar,
		int	ipmt
	)

inline

Definition at line 156 of file ParticleID.h.

                                {
    return m_tofcorrpid->dt(ipar, ipmt);
  }

◆ dtCorr()

double ParticleID::dtCorr	(	int	ipar,
		int	ipmt
	)

inline

Definition at line 159 of file ParticleID.h.

                                    {
    return m_tofcorrpid->dtCorr(ipar, ipmt);
  }

◆ getNhitCut()

int ParticleID::getNhitCut ( ) const

inline

Definition at line 88 of file ParticleID.h.

88{return m_nhitcut;}

◆ identify()

void ParticleID::identify ( const int pidcase )

inline

Definition at line 103 of file ParticleID.h.

103{m_pidcase = (m_pidcase | pidcase);}

Referenced by Gam4pikp::execute(), Rhopi::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQAPi2p2::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), DTagSetAlg::execute(), LTagSetAlg::execute(), TagSetAlg::execute(), XYZTagSetAlg::execute(), BeamParams::execute(), LocalKaonSelector::operator()(), LocalPionSelector::operator()(), LocalLambdaSelector::operator()(), and LocalProtonSelector::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();
   }
 
   // global info.
   m_pidsys = 0;
   m_pidcase = 0;
   m_method = 0;
   m_TotalLikelihood =0;
   m_discard =1;
   // probability
   m_ndof = 0;
   m_nhitcut=5;
   // chi cut
   setChiMinCut( 4.0 );
   setChiMaxCut( 6.0 );
   for(int i = 0; i < 5; i++) {
      m_chisq[i] = 9999.;
      m_prob[i] = -1.0;
   }
}

Referenced by Gam4pikp::execute(), Rhopi::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQAPi2p2::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), DTagSetAlg::execute(), LTagSetAlg::execute(), TagSetAlg::execute(), XYZTagSetAlg::execute(), BeamParams::execute(), LocalKaonSelector::operator()(), LocalPionSelector::operator()(), LocalLambdaSelector::operator()(), and LocalProtonSelector::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 Gam4pikp::execute(), Rhopi::execute(), inclks::execute(), inclkstar::execute(), incllambda::execute(), inclphi::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQAPi2p2::execute(), DQARhopi::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), DTagSetAlg::execute(), LTagSetAlg::execute(), TagSetAlg::execute(), XYZTagSetAlg::execute(), BeamParams::execute(), LocalKaonSelector::operator()(), LocalPionSelector::operator()(), LocalLambdaSelector::operator()(), and LocalProtonSelector::operator()().

◆ ipmt()

int ParticleID::ipmt ( )

inline

Definition at line 153 of file ParticleID.h.

                {
    return m_tofcorrpid->ipmt();
  }

Referenced by chiCorr(), dt(), dtCorr(), and sigCorr().

◆ IsDedxInfoValid()

bool ParticleID::IsDedxInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsEmcInfoValid()

bool ParticleID::IsEmcInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsMucInfoValid()

bool ParticleID::IsMucInfoValid ( ) const

inline

Referenced by LikelihoodCalculation().

◆ IsPidInfoValid()

bool ParticleID::IsPidInfoValid ( ) const

inlinevirtual

◆ 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 119 of file ParticleID.h.

119{return m_likelihoodfraction[n];}

Referenced by neuronVal().

◆ LikelihoodCalculation()

int ParticleID::LikelihoodCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 432 of file ParticleID.cxx.

                                      {
   int irc = -1;
 
   bool valid = IsDedxInfoValid() || IsTofInfoValid() || IsTofEInfoValid() || IsTofQInfoValid() || IsTofCInfoValid() || IsTofCorrInfoValid() || IsEmcInfoValid()||IsMucInfoValid();
   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(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 115 of file ParticleID.h.

115{return m_ndof;}

◆ neuronPIDCalculation()

int ParticleID::neuronPIDCalculation ( )

protected

◆ neuronStat()

int ParticleID::neuronStat ( ) const

inline

Definition at line 144 of file ParticleID.h.

144{return IsPidInfoValid();}

ParticleID::IsPidInfoValid

bool IsPidInfoValid() const

◆ neuronVal()

double ParticleID::neuronVal ( ) const

inline

Definition at line 145 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 ipar )

inline

Definition at line 168 of file ParticleID.h.

                          {
    return m_tofcorrpid->offset(ipar);
  }

◆ particleIDCalculation()

int ParticleID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 256 of file ParticleID.cxx.

                                       {
   int irc = -1;
   bool valid = IsDedxInfoValid() || IsTofInfoValid()||IsTofEInfoValid()
                || IsTofQInfoValid() || IsEmcInfoValid() || IsMucInfoValid()
                || IsTofCInfoValid() || IsTofCorrInfoValid();
 
   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] && ( chiTofCorr(i) < m_tofcorrpid->chiMaxCut() ) && ( chiTofCorr(i) > ( 0.0 - 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
   }
 
 
   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 118 of file ParticleID.h.

118{ 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 129 of file ParticleID.h.

129{return m_pdf[0];}

◆ pdfEmc()

double ParticleID::pdfEmc ( int n )

inline

Referenced by LikelihoodCalculation().

◆ pdfKaon()

double ParticleID::pdfKaon ( ) const

inline

Definition at line 132 of file ParticleID.h.

132{return m_pdf[3];}

◆ pdfMuc()

double ParticleID::pdfMuc ( int n )

inline

Referenced by LikelihoodCalculation().

◆ pdfMuon()

double ParticleID::pdfMuon ( ) const

inline

Definition at line 130 of file ParticleID.h.

130{return m_pdf[1];}

◆ pdfPion()

double ParticleID::pdfPion ( ) const

inline

Definition at line 131 of file ParticleID.h.

131{return m_pdf[2];}

◆ pdfProton()

double ParticleID::pdfProton ( ) const

inline

Definition at line 133 of file ParticleID.h.

133{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 114 of file ParticleID.h.

114{ return m_prob[n];}

Referenced by DQAKsKpi::execute(), DQAKsKpiDEDX::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 121 of file ParticleID.h.

121{return m_prob[0];}

Referenced by Gam4pikp::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), DTagSetAlg::execute(), LTagSetAlg::execute(), TagSetAlg::execute(), XYZTagSetAlg::execute(), and BeamParams::execute().

◆ probEmc()

double ParticleID::probEmc ( int n ) const

inline

◆ probKaon()

double ParticleID::probKaon ( ) const

inline

◆ probMuc()

double ParticleID::probMuc ( int n ) const

inline

◆ probMuon()

double ParticleID::probMuon ( ) const

inline

Definition at line 122 of file ParticleID.h.

122{return m_prob[1];}

Referenced by Gam4pikp::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), DTagSetAlg::execute(), LTagSetAlg::execute(), TagSetAlg::execute(), XYZTagSetAlg::execute(), and BeamParams::execute().

◆ probPion()

double ParticleID::probPion ( ) const

inline

◆ probProton()

double ParticleID::probProton ( ) const

inline

Definition at line 125 of file ParticleID.h.

125{return m_prob[4];}

Referenced by Gam4pikp::execute(), incllambda::execute(), DQAJpsi2PPbarAlg::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQAPi2p2::execute(), DQASelBhabha::execute(), DQASelDimu::execute(), DQASelHadron::execute(), DTagSetAlg::execute(), LTagSetAlg::execute(), TagSetAlg::execute(), XYZTagSetAlg::execute(), BeamParams::execute(), LocalLambdaSelector::operator()(), and LocalProtonSelector::operator()().

◆ 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

◆ setNhitCut()

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

inline

Definition at line 89 of file ParticleID.h.

89{m_nhitcut = nhitcuth;}

◆ sigCorr()

double ParticleID::sigCorr	(	int	ipar,
		int	ipmt
	)

inline

Definition at line 162 of file ParticleID.h.

                                     {
    return m_tofcorrpid->sigCorr(ipar, ipmt);
  }

◆ sigma()

double ParticleID::sigma ( int ipar )

inline

Definition at line 171 of file ParticleID.h.

                         {
    return m_tofcorrpid->sigma(ipar);
  }

◆ usePidSys()

void ParticleID::usePidSys ( const int pidsys )

inline

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()

◆ chiCorr()

◆ chiDedx()

◆ chiEmc()

◆ chiMuc()

◆ chiTof()

◆ chiTof1()

◆ chiTof2()

◆ chiTofC()

◆ chiTofCorr()

◆ chiTofE()

◆ chiTofQ()

◆ discard()

◆ dt()

◆ dtCorr()

◆ getNhitCut()

◆ identify()

◆ init()

◆ instance()

◆ ipmt()

◆ IsDedxInfoValid()

◆ IsEmcInfoValid()

◆ IsMucInfoValid()

◆ IsPidInfoValid()

◆ IsTofCInfoValid()

◆ IsTofCorrInfoValid()

◆ IsTofEInfoValid()

◆ IsTofInfoValid()

◆ IsTofQInfoValid()

◆ lhf()

◆ LikelihoodCalculation()

◆ ndof()

◆ neuronPIDCalculation()

◆ neuronStat()

◆ neuronVal()

◆ offset()

◆ particleIDCalculation()

◆ pdf()

◆ pdfDedx()

◆ pdfElectron()

◆ pdfEmc()

◆ pdfKaon()

◆ pdfMuc()

◆ pdfMuon()

◆ pdfPion()

◆ pdfProton()

◆ pdfTof()

◆ pdfTofC()

◆ pdfTofCorr()

◆ pdfTofE()

◆ pdfTofQ()

◆ prob()

◆ probDedx()

◆ probElectron()

◆ probEmc()

◆ probKaon()

◆ probMuc()

◆ probMuon()

◆ probPion()

◆ probProton()

◆ probTof()

◆ probTofC()

◆ probTofCorr()

◆ probTofE()

◆ probTofQ()

◆ setMethod()

◆ setNhitCut()

◆ sigCorr()

◆ sigma()

◆ usePidSys()