#include <Tof2PID.h>

Inheritance diagram for Tof2PID:

Public Member Functions
	~Tof2PID ()

void	init ()

void	calculate ()

bool	IsPidInfoValid () const

double	chi (int n) const

double	sigma (int n) const

double	offset (int n) const

double	prob (int n) const

int	ndof () const

double	mass2 () const

double	ph2 () const

double	zhit2 () const

double	path2 () const

double	tof2 () const

int	neuronPID () const

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

Protected Member Functions
int	particleIDCalculation ()

int	neuronPIDCalculation ()

int	LikelihoodCalculation ()

double	offsetTof2 (int n, int cntr, double ptrk, double ztof, double m_ph2, double charge)

double	sigmaTof2 (int n, int cntr, double ptrk, double ztof, double m_ph2, double charge)

double	sampleQ0 (double betagamma, double beta)

Protected Member Functions inherited from ParticleIDBase

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

Detailed Description

Definition at line 14 of file Tof2PID.h.

Constructor & Destructor Documentation

◆ ~Tof2PID()

Tof2PID::~Tof2PID ( )

inline

Definition at line 18 of file Tof2PID.h.

18{;}

Member Function Documentation

◆ calculate()

void Tof2PID::calculate ( )

virtual

Implements ParticleIDBase.

Definition at line 58 of file Tof2PID.cxx.

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

Referenced by TofPID::particleIDCalculation().

◆ chi()

double Tof2PID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 23 of file Tof2PID.h.

23{return m_chi[n];}

n

const Int_t n

Definition DataBase/tau_mode.c:65

Referenced by TofPID::particleIDCalculation().

◆ init()

void Tof2PID::init ( )

virtual

Implements ParticleIDBase.

Definition at line 42 of file Tof2PID.cxx.

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

Referenced by TofPID::particleIDCalculation().

◆ instance()

Tof2PID * Tof2PID::instance ( )

static

Definition at line 17 of file Tof2PID.cxx.

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

Referenced by TofPID::particleIDCalculation().

◆ IsPidInfoValid()

bool Tof2PID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 22 of file Tof2PID.h.

22{return (m_ndof > 0); }

Referenced by TofPID::particleIDCalculation().

◆ LikelihoodCalculation()

int Tof2PID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 38 of file Tof2PID.h.

38{return -1;}

◆ mass2()

double Tof2PID::mass2 ( ) const

inline

Definition at line 28 of file Tof2PID.h.

28{return m_mass2;}

◆ ndof()

int Tof2PID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 27 of file Tof2PID.h.

27{return m_ndof;}

◆ neuronPID()

int Tof2PID::neuronPID ( ) const

inline

Definition at line 33 of file Tof2PID.h.

33{return -1;}

◆ neuronPIDCalculation()

int Tof2PID::neuronPIDCalculation ( )

inlineprotected

Definition at line 37 of file Tof2PID.h.

37{return -1;}

◆ offset()

double Tof2PID::offset ( int n ) const

inline

Definition at line 25 of file Tof2PID.h.

25{return m_offset[n];}

Referenced by offsetTof2(), and TofPID::particleIDCalculation().

◆ offsetTof2()

double Tof2PID::offsetTof2	(	int	n,
		int	cntr,
		double	ptrk,
		double	ztof,
		double	m_ph2,
		double	charge )

protected

Definition at line 153 of file Tof2PID.cxx.

                                                                                          {
   double offset;
   //   double gb = ptrk/xmass(n);
   double betagamma;
   switch(n) {
 
   case 0: {  // Electron
      offset = 0.0;
      return offset;
   }
 
   case 1: { // Muon
      offset = 0.0;
      return offset;
   }
 
   case 2: { // Pion
      betagamma=ptrk/139.57;
      break;
   }
   case 3: { // Kaon
      betagamma=ptrk/493.68;
      break;
   }
 
   case 4: { // Proton
      betagamma=ptrk/938.27;
      break;
   }
 
   default: {
      offset = 0.0;
      return offset;
   }
   break;
   }
   z=z/1000.0;
   double Q = ph;
   double beta = betagamma / TMath::Sqrt(1 + betagamma*betagamma);
   double Q0 = sampleQ0(betagamma,beta);
   double func[15]= {0.};
   func[0]=1.;
   func[1]=1./sqrt(Q);
   func[2]=z/sqrt(Q);
   func[3]=z*z/sqrt(Q);
   func[4]=Q;
   func[5]=Q*Q;
   func[6]=Q*Q*Q;
   func[7]=1./(0.89*0.89+z*z);
   func[8]=z;
   func[9]=z*z;
   func[10]=z*z*z;
   func[11]=1./sqrt(Q0);
   func[12]=z/sqrt(Q0);
   func[13]=z*z/sqrt(Q0);
   func[14]=z*z*z/sqrt(Q0);
   offset=0.;
   for(int i=0; i<15; i++) {
      offset+= m_pars[i]*func[i];
   }
 
   return offset;
}

◆ particleIDCalculation()

int Tof2PID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 62 of file Tof2PID.cxx.

                                   {
   int irc = -1;
   EvtRecTrack* recTrk = PidTrk();
   if(!(recTrk->isMdcTrackValid())) return irc;
   RecMdcTrack* mdcTrk = recTrk->mdcTrack();
 
   double ptrk = mdcTrk->p();
   //   double cost = cos(mdcTrk->theta());
   double charge = mdcTrk->charge();
 
   if(!(recTrk->isTofTrackValid())) return irc;
 
#ifndef BEAN
   SmartRefVector<RecTofTrack> tofTrk = recTrk->tofTrack();
   SmartRefVector<RecTofTrack>::iterator it;//=tofTrk.begin();
#else
   const std::vector<TTofTrack* >& tofTrk = recTrk->tofTrack();
   std::vector<TTofTrack* >::const_iterator it;//=tofTrk.begin();
#endif
 
   TofHitStatus *hitst = new TofHitStatus;
   std::vector<int> tof2count;
   int goodtof2trk=0;
   for(it = tofTrk.begin(); it!=tofTrk.end(); it++,goodtof2trk++) {
      unsigned int st = (*it)->status();
      hitst->setStatus(st);
 
      if( !(hitst->is_barrel()) ) continue;
      if( !(hitst->is_counter()) ) continue;
      if( hitst->layer()==2 ) tof2count.push_back(goodtof2trk);
   }
   delete hitst;
 
   if(tof2count.size()!=1) return irc;//not tof2 track or more than 1 tracks
   it = tofTrk.begin()+tof2count[0];
   double tof  = (*it)->tof();
   m_tof2 = tof;
 
   if(tof <=0 ) return irc;
   //   int qual = (*it)->quality();
   //  if(qual != 1) return irc;
   int cntr = (*it)->tofID();
   double path = ((*it)->path())*10.0;//the unit from mm to cm
 
   m_path2 = path;
   //   m_path2 = path;
   // fix the bugs the 6.0.0, He K.L.
   // double path = tofTrk->getPath1();
   //
   m_ph2   = (*it)->ph();
   m_zhit2 = ((*it)->zrhit())*10; //the unit from mm to cm
 
   double beta2 = path*path/velc()/velc()/tof/tof;
   m_mass2 = ptrk * ptrk * (1/beta2 -1);
   if ((m_mass2>20)||(m_mass2<-1)) return irc;
 
   double chitemp = 99.;
   double pdftemp = 0;
   double sig_tmp = (*it)->sigma(0);
   for(int i = 0; i < 5; i++) {
      m_offset[i] = tof - (*it)->texp(i)-(*it)->toffset(i);//- offsetTof1(i, cntr, ptrk, m_zhit1, m_ph1,charge);
      if(sig_tmp!=0)  {
         m_sigma[i] = 1.2*sig_tmp;
         if(i<2) m_sigma[i]=sig_tmp;
      }
      else m_sigma[i]= sigmaTof2(i, cntr,ptrk,m_zhit2, m_ph2,charge);
      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);
   }
   //  if(ptrk<0.65) m_chi[2]=m_chi[3];
   m_chimin = chitemp;
   m_pdfmin = pdftemp;
   if(pdftemp < pdfCalculate(pdfMinSigmaCut(),1)) return irc;
   if(fabs(m_chimin) > chiMinCut()) return irc;
 
   // calculate prob
 
   for(int i = 0; i < 5; i++) {
      m_prob[i] = probCalculate(m_chi[i]*m_chi[i], 1);
   }
   irc = 0;
   return irc;
}

Referenced by calculate().

◆ path2()

double Tof2PID::path2 ( ) const

inline

Definition at line 31 of file Tof2PID.h.

31{return m_path2;}

◆ ph2()

double Tof2PID::ph2 ( ) const

inline

Definition at line 29 of file Tof2PID.h.

29{return m_ph2;}

◆ prob()

double Tof2PID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 26 of file Tof2PID.h.

26{return m_prob[n];}

◆ sampleQ0()

double Tof2PID::sampleQ0	(	double	betagamma,
		double	beta )

protected

Definition at line 258 of file Tof2PID.cxx.

                                                     {
   double p1 = 0.261;
   double p2 = 5.43;
   double p3 = 1.12;
   double p4 = 1.96;
   double p5 = 0.386;
   double val = p1 / TMath::Power(beta, p4) *
                ( p2- TMath::Power(beta, p4) -
                  TMath::Log(p3 + TMath::Power(1/betagamma, p5)));
 
   return val*100.;
}

Referenced by offsetTof2().

◆ sigma()

double Tof2PID::sigma ( int n ) const

inline

Definition at line 24 of file Tof2PID.h.

24{return m_sigma[n];}

Referenced by TofPID::particleIDCalculation().

◆ sigmaTof2()

double Tof2PID::sigmaTof2	(	int	n,
		int	cntr,
		double	ptrk,
		double	ztof,
		double	m_ph2,
		double	charge )

protected

Definition at line 217 of file Tof2PID.cxx.

                                                                                             {
   double sigma;
   //   double gb = ptrk/xmass(n);
   double cosz = cos(ztof/1000.0);
   //  double phtemp=ph;
   switch(n) {
 
   case 0: { // Electron
      sigma = 0.115816 -0.0933215*cosz+ 0.0788252*cosz*cosz;
      break;
   }
 
   case 1: { // Muon
      sigma = 0.121536 -0.0935898*cosz+0.0748771*cosz*cosz;
      break;
   }
 
   case 2: { // Pion
      sigma =  0.106882-0.0147375*cosz+0.027491*cosz*cosz;
      break;
   }
   case 3: { // Kaon
      sigma =0.106882 -0.0147375*cosz +0.027491 *cosz*cosz;
      break;
   }
 
   case 4: { // Proton
      sigma = 0.168489 -0.131762*cosz+ 0.105708*cosz*cosz;
      break;
   }
 
   default:
      sigma = 0.100;
      break;
   }
   sigma = 0.110;
   //sigma = 1.0;
   return sigma;
}

Referenced by particleIDCalculation().

◆ tof2()

double Tof2PID::tof2 ( ) const

inline

Definition at line 32 of file Tof2PID.h.

32{return m_tof2;}

◆ zhit2()

double Tof2PID::zhit2 ( ) const

inline

Definition at line 30 of file Tof2PID.h.

30{return m_zhit2;}

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

◆ ~Tof2PID()

Member Function Documentation

◆ calculate()

◆ chi()

◆ init()

◆ instance()

◆ IsPidInfoValid()

◆ LikelihoodCalculation()

◆ mass2()

◆ ndof()

◆ neuronPID()

◆ neuronPIDCalculation()

◆ offset()

◆ offsetTof2()

◆ particleIDCalculation()

◆ path2()

◆ ph2()

◆ prob()

◆ sampleQ0()

◆ sigma()

◆ sigmaTof2()

◆ tof2()

◆ zhit2()