#include <Tof1PID.h>

Inheritance diagram for Tof1PID:

Public Member Functions
	~Tof1PID ()

void	init ()

void	calculate ()

bool	IsPidInfoValid () const

double	chi (int n) const

double	prob (int n) const

double	offset (int n) const

double	sigma (int n) const

int	ndof () const

double	mass2 () const

double	ph1 () const

double	zhit1 () const

double	path1 () const

double	tof1 () 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 Tof1PID *	instance ()

Protected Member Functions
int	particleIDCalculation ()

int	neuronPIDCalculation ()

int	LikelihoodCalculation ()

double	offsetTof1 (int n, int cntr, double ptrk, double ztof, double m_ph1, double charge)

double	sigmaTof1 (int n, int cntr, double ptrk, double ztof, double m_ph1, double charge)

double	sampleQ0 (double betagamma, double beta)

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 14 of file Tof1PID.h.

Constructor & Destructor Documentation

◆ ~Tof1PID()

Tof1PID::~Tof1PID ( )

inline

Definition at line 19 of file Tof1PID.h.

19{;}

Member Function Documentation

◆ calculate()

void Tof1PID::calculate ( )

virtual

Implements ParticleIDBase.

Definition at line 55 of file Tof1PID.cxx.

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

Referenced by TofPID::particleIDCalculation().

◆ chi()

double Tof1PID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 24 of file Tof1PID.h.

24{return m_chi[n];}

Referenced by TofPID::particleIDCalculation().

◆ init()

void Tof1PID::init ( )

virtual

Implements ParticleIDBase.

Definition at line 40 of file Tof1PID.cxx.

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

Referenced by TofPID::particleIDCalculation().

◆ instance()

Tof1PID * Tof1PID::instance ( )

static

Definition at line 16 of file Tof1PID.cxx.

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

Referenced by TofPID::particleIDCalculation().

◆ IsPidInfoValid()

bool Tof1PID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 23 of file Tof1PID.h.

23{return (m_ndof > 0); }

Referenced by TofPID::particleIDCalculation().

◆ LikelihoodCalculation()

int Tof1PID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 40 of file Tof1PID.h.

40{return -1;}

◆ mass2()

double Tof1PID::mass2 ( ) const

inline

Definition at line 29 of file Tof1PID.h.

29{return m_mass2;}

◆ ndof()

int Tof1PID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 28 of file Tof1PID.h.

28{return m_ndof;}

◆ neuronPID()

int Tof1PID::neuronPID ( ) const

inline

Definition at line 34 of file Tof1PID.h.

34{return -1;}

◆ neuronPIDCalculation()

int Tof1PID::neuronPIDCalculation ( )

inlineprotected

Definition at line 39 of file Tof1PID.h.

39{ return -1;}

◆ offset()

double Tof1PID::offset ( int n ) const

inline

Definition at line 26 of file Tof1PID.h.

26{return m_offset[n];}

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

◆ offsetTof1()

double Tof1PID::offsetTof1	(	int	n,
		int	cntr,
		double	ptrk,
		double	ztof,
		double	m_ph1,
		double	charge
	)

protected

Definition at line 137 of file Tof1PID.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;
   }
   double Q = ph1;
   z=z/1000.0;
   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.84*0.84+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 Tof1PID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 59 of file Tof1PID.cxx.

                                   {
   int irc = -1;
 
   EvtRecTrack* recTrk = PidTrk();
   if(!(recTrk->isMdcTrackValid())) return irc;
   RecMdcTrack* mdcTrk = recTrk->mdcTrack();
   double ptrk = mdcTrk->p();
   double charge = mdcTrk->charge();
   //   double cost = cos(mdcTrk->theta());
   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> tof1count;
   int goodtof1trk=0;
   for(it = tofTrk.begin(); it!=tofTrk.end(); it++,goodtof1trk++) {
      unsigned int st = (*it)->status();
      hitst->setStatus(st);
      if( !(hitst->is_barrel()) ) continue;
      if( !(hitst->is_counter()) ) continue;
      if( hitst->layer()==1 )  tof1count.push_back(goodtof1trk);
   }
   delete hitst;
   if(tof1count.size()!=1) return irc;//not tof2 track or more than 1 tracks
   it = tofTrk.begin()+tof1count[0];
   double tof  = (*it)->tof();
   m_tof1 = tof;
   //   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_path1 = path;
   m_ph1  = (*it)->ph(); //no change
   m_zhit1 = ((*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;
   if(tof <=0 ) return irc;
   double chitemp = 99.;
   double pdftemp = 0;
   double sigma_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(sigma_tmp!=0) {
         m_sigma[i] = 1.2*sigma_tmp;
         if(i<2) m_sigma[i]=sigma_tmp;
      }
      else m_sigma[i]=sigmaTof1(i, cntr,ptrk,m_zhit1, m_ph1,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);
   }
   m_chimin = chitemp;
   m_pdfmin = pdftemp;
   if(pdftemp < pdfCalculate(pdfMinSigmaCut(),1)) return irc;
   if(fabs(m_chimin) > chiMinCut()) return irc;
   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().

◆ path1()

double Tof1PID::path1 ( ) const

inline

Definition at line 32 of file Tof1PID.h.

32{return m_path1;}

◆ ph1()

double Tof1PID::ph1 ( ) const

inline

Definition at line 30 of file Tof1PID.h.

30{return m_ph1;}

Referenced by offsetTof1().

◆ prob()

double Tof1PID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 25 of file Tof1PID.h.

25{return m_prob[n];}

◆ sampleQ0()

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

protected

Definition at line 242 of file Tof1PID.cxx.

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

Referenced by offsetTof1().

◆ sigma()

double Tof1PID::sigma ( int n ) const

inline

Definition at line 27 of file Tof1PID.h.

27{return m_sigma[n];}

Referenced by TofPID::particleIDCalculation().

◆ sigmaTof1()

double Tof1PID::sigmaTof1	(	int	n,
		int	cntr,
		double	ptrk,
		double	ztof,
		double	m_ph1,
		double	charge
	)

protected

Definition at line 201 of file Tof1PID.cxx.

                                                                                             {
   double sigma;
   //   double gb = ptrk/xmass(n);
   double cosz = cos(ztof/1000.0);
 
   switch(n) {
 
   case 0: { // Electron
      sigma = 0.132405-0.135638*cosz+0.105528*cosz*cosz;
      break;
   }
 
   case 1: { // Muon
      sigma = 0.164057 -0.199648*cosz+ 0.139481*cosz*cosz;
      break;
   }
 
   case 2: { // Pion
      sigma = 0.132943-0.0996678*cosz+0.0785578*cosz*cosz;
      break;
   }
   case 3: { // Kaon
      sigma = 0.146572 -0.124672*cosz+0.0920656*cosz*cosz;
      break;
   }
 
   case 4: { // Proton
      sigma = 0.1513 -0.0806081*cosz+0.0607409*cosz*cosz;
      break;
   }
 
   default:
      sigma = 0.100;
      break;
   }
   sigma = 0.110;
   //sigma = 1.0;
   return sigma;
}

Referenced by particleIDCalculation().

◆ tof1()

double Tof1PID::tof1 ( ) const

inline

Definition at line 33 of file Tof1PID.h.

33{return m_tof1;}

◆ zhit1()

double Tof1PID::zhit1 ( ) const

inline

Definition at line 31 of file Tof1PID.h.

31{return m_zhit1;}

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

◆ ~Tof1PID()

Member Function Documentation

◆ calculate()

◆ chi()

◆ init()

◆ instance()

◆ IsPidInfoValid()

◆ LikelihoodCalculation()

◆ mass2()

◆ ndof()

◆ neuronPID()

◆ neuronPIDCalculation()

◆ offset()

◆ offsetTof1()

◆ particleIDCalculation()

◆ path1()

◆ ph1()

◆ prob()

◆ sampleQ0()

◆ sigma()

◆ sigmaTof1()

◆ tof1()

◆ zhit1()