#include <TofEPID.h>

Inheritance diagram for TofEPID:

Public Member Functions
	~TofEPID ()

void	init ()

void	calculate ()

bool	IsPidInfoValid () const

double	chi (int n) const

double	prob (int n) const

double	sigma (int n) const

double	offset (int n) const

int	ndof () const

double	mass2 () const

int	part () const

double	rhit () const

int	neuronPID () const

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 TofEPID *	instance ()

Protected Member Functions
int	neuronPIDCalculation ()

int	particleIDCalculation ()

int	LikelihoodCalculation ()

double	offsetTofE (int n, int cntr, double ptrk, double rtof, double ph, double charge)

double	sigmaTofE (int n, int cntr, double ptrk, double rtof, double ph, double charge)

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 13 of file TofEPID.h.

Constructor & Destructor Documentation

◆ ~TofEPID()

TofEPID::~TofEPID ( )

inline

Definition at line 17 of file TofEPID.h.

17{;}

Member Function Documentation

◆ calculate()

void TofEPID::calculate ( )

virtual

Implements ParticleIDBase.

Definition at line 39 of file TofEPID.cxx.

                        {
   if(particleIDCalculation()==0) m_ndof=1;
}

◆ chi()

double TofEPID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 22 of file TofEPID.h.

22{return m_chi[n];}

◆ init()

void TofEPID::init ( )

virtual

Implements ParticleIDBase.

Definition at line 25 of file TofEPID.cxx.

                   {
   for(int i = 0; i < 5; i++) {
      m_chi[i] = 99.0;
      m_prob[i] = -1.0;
      m_offset[i] = 99.0;
      m_sigma[i] = 1.0;
   }
   m_chimin = 99.;
   m_pdfmin =99.;
   m_ndof = 0;
   m_mass2 = -999;
   m_rhit = -99;
}

◆ instance()

TofEPID * TofEPID::instance ( )

static

Definition at line 16 of file TofEPID.cxx.

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

◆ IsPidInfoValid()

bool TofEPID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 21 of file TofEPID.h.

21{return (m_ndof > 0); }

◆ LikelihoodCalculation()

int TofEPID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 35 of file TofEPID.h.

35{return -1;}

◆ mass2()

double TofEPID::mass2 ( ) const

inline

Definition at line 27 of file TofEPID.h.

27{return m_mass2;}

◆ ndof()

int TofEPID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 26 of file TofEPID.h.

26{return m_ndof;}

◆ neuronPID()

int TofEPID::neuronPID ( ) const

inline

Definition at line 30 of file TofEPID.h.

30{return -1;}

◆ neuronPIDCalculation()

int TofEPID::neuronPIDCalculation ( )

inlineprotected

Definition at line 33 of file TofEPID.h.

33{ return -1;}

◆ offset()

double TofEPID::offset ( int n ) const

inline

Definition at line 25 of file TofEPID.h.

25{return m_offset[n];}

Referenced by offsetTofE().

◆ offsetTofE()

double TofEPID::offsetTofE	(	int	n,
		int	cntr,
		double	ptrk,
		double	rtof,
		double	ph,
		double	charge
	)

protected

Definition at line 158 of file TofEPID.cxx.

                                                                                             {
   double offset;
   //   double gb = ptrk/xmass(n);
   switch(n) {
   case 0: {  // Electron
      double ptemp = ptrk;
      if(ptrk<0.2) ptemp = 0.2;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      offset = 0.001*(-28.8481+138.159*plog-249.334*plog*plog);
      break;
   }
 
   case 1: { // Muon
      double ptemp = ptrk;
      if(ptrk<0.2) ptemp = 0.2;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      offset = 0.001*(-33.6966+1.91915*plog-0.592320*plog*plog);
      break;
   }
   case 2: { // Pion
      double ptemp = ptrk;
      if(ptrk<0.2) ptemp = 0.2;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      offset = 0.001*(-27.9965 + 1.213 * plog - 2.02740 * plog * plog);
      break;
   }
   case 3: { // Kaon
      double ptemp = ptrk;
      if(ptrk<0.3) ptemp = 0.3;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      offset = 0.001*(-23.4842 -28.7569 * plog + 78.21* plog *plog);
      break;
   }
 
   case 4: { // Proton
      double ptemp = ptrk;
      if(ptrk<0.4) ptemp = 0.4;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      if(charge>0)
         offset = 0.001*(-4.854-110.540*plog+99.8732*plog*plog);
      if(charge<0)
         offset = 0.001*(27.047-145.120*plog+167.014*plog*plog);
      break;
   }
 
   default:
      offset = 0.0;
      break;
   }
   //  offset = 0.0;
   return offset;
}

◆ part()

int TofEPID::part ( ) const

inline

Definition at line 28 of file TofEPID.h.

28{return m_part;}

◆ particleIDCalculation()

int TofEPID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 43 of file TofEPID.cxx.

                                   {
   int irc = -1;
   EvtRecTrack* recTrk = PidTrk();
   if(!(recTrk->isMdcTrackValid())) return irc;
   if(!(recTrk->isTofTrackValid())) return irc;
 
#ifndef BEAN
   SmartRefVector<RecTofTrack> tofTrackCol = recTrk->tofTrack();
   SmartRefVector<RecTofTrack>::iterator it;
#else
   const std::vector<TTofTrack* >& tofTrackCol = recTrk->tofTrack();
   std::vector<TTofTrack* >::const_iterator it;
#endif
 
   bool isMRPC = false;
   if( tofTrackCol.size()!=1 && tofTrackCol.size()!=3 ) { return irc; }
   else {
     TofHitStatus *hitst = new TofHitStatus;
     if( tofTrackCol.size()==1 ) {
       it = tofTrackCol.begin();
       hitst->setStatus( (*it)->status() );
       bool isRaw = hitst->is_raw();
       if( isRaw ) { return irc; }
       isMRPC = hitst->is_mrpc();
       bool isReadout = hitst->is_readout();
       bool isCounter = hitst->is_counter();
       bool isCluster = hitst->is_cluster();
       if( !isReadout || !isCounter || !isCluster ) { return irc; }
     }
     else if( tofTrackCol.size()==3 ) {
       unsigned int icounter = -1;
       unsigned int inumber  = 0;
       for( it = tofTrackCol.begin(); it != tofTrackCol.end(); it++, inumber++ ) {
     hitst->setStatus( (*it)->status() );
     bool isRaw = hitst->is_raw();
     if( isRaw ) continue;
     isMRPC = hitst->is_mrpc();
     if( !isMRPC ) { return irc; } // MRPC is possible to have 3 record.
     bool isReadout = hitst->is_readout();
     bool isCounter = hitst->is_counter();
     bool isCluster = hitst->is_cluster();
     if( !isReadout && isCounter && isCluster ) { icounter = inumber; }
       }
       if( icounter == -1 ) { return irc; }
       it = tofTrackCol.begin() + icounter;
     }
   }
 
   double tof = (*it)->tof();
   if( tof <= 0 ) { return irc; }
   double path = (*it)->path();
   m_rhit = (*it)->zrhit();
   double beta = path/velc()/tof;
   double beta2 = beta*beta;
   RecMdcTrack* mdcTrk = recTrk->mdcTrack();
   double ptrk = mdcTrk->p();
   m_mass2 = ptrk*ptrk*(1.0/beta2 -1.0);
 
   double chitemp = 99.;
   double pdftemp = 0;
   for( unsigned int i=0; i<5; i++ ) {
     double texp = (*it)->texp(i);
     m_offset[i] = tof - texp - (*it)->toffset(i);
     double sigma_tmp = (*it)->sigma(0);
     if( !isMRPC ) {
       if( sigma_tmp>0 ) {
     if( i<2 ) { m_sigma[i] = sigma_tmp; }
     else { m_sigma[i] = 1.2*sigma_tmp; }
       }
       else { m_sigma[i] = 0.12; }
     }
     else {
       //       m_sigma[i] = 0.065;
       int strip = (*it)->strip();
       if( strip<0 || strip>11 ) { m_sigma[i] = 0.065; }
       else {
     double p0, p1;
     if( strip==0 )       { p0=0.16;     p1=0.0;      }
     else if( strip==1 )  { p0=0.051094; p1=0.019467; }
     else if( strip==2 )  { p0=0.056019; p1=0.012964; }
     else if( strip==3 )  { p0=0.043901; p1=0.015933; }
     else if( strip==4 )  { p0=0.049750; p1=0.010473; }
     else if( strip==5 )  { p0=0.048345; p1=0.008545; }
     else if( strip==6 )  { p0=0.046518; p1=0.009038; }
     else if( strip==7 )  { p0=0.048803; p1=0.007251; }
     else if( strip==8 )  { p0=0.047523; p1=0.008434; }
     else if( strip==9 )  { p0=0.042187; p1=0.010307; }
     else if( strip==10 ) { p0=0.050337; p1=0.005951; }
     else if( strip==11 ) { p0=0.054740; p1=0.005629; }
     if( ptrk<0.05 ) { m_sigma[i] = p0+p1/0.05; }
     else { m_sigma[i] = p0+p1/ptrk; }
       }
     }
     m_chi[i] = m_offset[i]/m_sigma[i];
 
     if( fabs(m_chi[i]) < chitemp ) { chitemp = fabs(m_chi[i]); }
     double ppp = pdfCalculate(m_chi[i],1);
     if( fabs(ppp) > pdftemp ) { pdftemp = fabs(ppp); }
   }
   m_chimin = chitemp;
   //   if( pdftemp < pdfCalculate(pdfMinSigmaCut(),1) ) { return irc; }
   //   if( fabs(m_chimin) > chiMinCut() ) { return irc; }
 
   for( unsigned int i=0; i<5; i++ ) {
     m_prob[i] = probCalculate( m_chi[i]*m_chi[i], 1 );
   }
 
   irc = 0;
   return irc;
}

Referenced by calculate().

◆ prob()

double TofEPID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 23 of file TofEPID.h.

23{return m_prob[n];}

◆ rhit()

double TofEPID::rhit ( ) const

inline

Definition at line 29 of file TofEPID.h.

29{return m_rhit;}

◆ sigma()

double TofEPID::sigma ( int n ) const

inline

Definition at line 24 of file TofEPID.h.

24{return m_sigma[n];}

◆ sigmaTofE()

double TofEPID::sigmaTofE	(	int	n,
		int	cntr,
		double	ptrk,
		double	rtof,
		double	ph,
		double	charge
	)

protected

Definition at line 216 of file TofEPID.cxx.

                                                                                             {
 
   double sigma;
   //   double gb = ptrk/xmass(n);
   switch(n) {
 
   case 0: { // Electron
      double ptemp = ptrk;
      if(ptrk < 0.2) ptemp = 0.2;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      sigma = 0.001 * (109.974 +15.2457 * plog + 36.8139 * plog * plog);
 
      break;
   }
 
   case 1: { // Muon
      double ptemp = ptrk;
      if(ptrk < 0.2) ptemp = 0.2;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      sigma = 0.001 * (96.5077 -2.96232 * plog + 3.12910 * plog * plog);
      break;
   }
 
   case 2: { // pion
      double ptemp = ptrk;
      if(ptrk < 0.2) ptemp = 0.2;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      sigma = 0.001 * (105.447 - 2.08044 * plog + 3.44846 * plog * plog);
      break;
   }
 
   case 3: { // Kaon
      double ptemp = ptrk;
      if(ptrk < 0.3) ptemp = 0.3;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      sigma = 0.001*(88.8806 - 26.8464 * plog + 113.672 * plog * plog);
      break;
   }
   case 4: { // Proton
      double ptemp = ptrk;
      if(ptrk < 0.5) ptemp = 0.5;
      if(ptrk > 2.1) ptemp = 2.1;
      double plog = log(ptemp);
      if(charge>0)
         sigma = 0.001 * (96.3534 -44.1139 * plog + 53.9805 * plog * plog);
      if(charge<0)
         sigma = 0.001 * (157.345 -98.7357 * plog + 55.1145 * plog * plog);
      break;
   }
 
   default:
      sigma = 0.100;
 
      break;
   }
   // sigma =1;
   return sigma;
}

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

◆ ~TofEPID()

Member Function Documentation

◆ calculate()

◆ chi()

◆ init()

◆ instance()

◆ IsPidInfoValid()

◆ LikelihoodCalculation()

◆ mass2()

◆ ndof()

◆ neuronPID()

◆ neuronPIDCalculation()

◆ offset()

◆ offsetTofE()

◆ part()

◆ particleIDCalculation()

◆ prob()

◆ rhit()

◆ sigma()

◆ sigmaTofE()