#include <EmcPID.h>

Inheritance diagram for EmcPID:

Public Member Functions
	~EmcPID ()

void	init ()

void	calculate ()

bool	IsPidInfoValid () const

double	chi (int n) const

double	prob (int n) const

int	ndof () const

double	energy () const

double	eseed () const

double	e3x3 () const

double	e5x5 () const

double	delta_theta () const

double	delta_phi () const

double	secondmoment () const

double	val_emc1 () const

int	neuronPID () const

double	pars (int n) 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 EmcPID *	instance ()

Protected Member Functions
int	particleIDCalculation ()

int	neuronPIDCalculation ()

int	LikelihoodCalculation ()

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 16 of file EmcPID.h.

Constructor & Destructor Documentation

◆ ~EmcPID()

EmcPID::~EmcPID ( )

inline

Definition at line 21 of file EmcPID.h.

21{;}

Member Function Documentation

◆ calculate()

void EmcPID::calculate ( )

virtual

Implements ParticleIDBase.

Definition at line 153 of file EmcPID.cxx.

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

◆ chi()

double EmcPID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 26 of file EmcPID.h.

26{return m_chi[n];}

◆ delta_phi()

double EmcPID::delta_phi ( ) const

inline

Definition at line 34 of file EmcPID.h.

34{return m_delta_phi;}

◆ delta_theta()

double EmcPID::delta_theta ( ) const

inline

Definition at line 33 of file EmcPID.h.

33{return m_delta_theta;}

◆ e3x3()

double EmcPID::e3x3 ( ) const

inline

Definition at line 31 of file EmcPID.h.

31{return m_e3x3;}

◆ e5x5()

double EmcPID::e5x5 ( ) const

inline

Definition at line 32 of file EmcPID.h.

32{return m_e5x5;}

◆ energy()

double EmcPID::energy ( ) const

inline

Definition at line 29 of file EmcPID.h.

29{return m_energy;}

◆ eseed()

double EmcPID::eseed ( ) const

inline

Definition at line 30 of file EmcPID.h.

30{return m_eseed;}

◆ init()

void EmcPID::init ( )

virtual

Implements ParticleIDBase.

Definition at line 95 of file EmcPID.cxx.

                  {
   for(int i = 0; i < 5; i++) {
      m_chi[i] = 99.0;
      m_prob[i] = -1.0;
   }
   m_chimin = 99.;
   m_ndof = 0;
   m_energy = -99;
   m_eseed = -99;
   m_e3x3 = -99;
   m_e5x5 = -99;
   m_delta_theta = -99;
   m_delta_phi = -99;
   m_secondmoment = -99;
   m_val_emc = -99;
   //   std::string emc =  path + "/share/emc_epimu.txt";
   /*   if(!m_trainTree_emconly){
     m_trainTree_emconly = new TTree("m_trainTree_emconly","m_trainTree_emconly");
     m_trainTree_emconly->Branch("ptrk",&m_ptrk,"ptrk/D");
     m_trainTree_emconly->Branch("pt",&m_pt,"pt/D");
     m_trainTree_emconly->Branch("type",&m_type,"type/D");
     m_trainTree_emconly->Branch("energy",&m_energy,"energy/D");
     m_trainTree_emconly->Branch("eseed",&m_eseed,"eseed/D");
     m_trainTree_emconly->Branch("e3x3",&m_e3x3,"e3x3/D");
     m_trainTree_emconly->Branch("e5x5",&m_e5x5,"e5x5/D");
     m_trainTree_emconly->Branch("secondmoment",&m_secondmoment,"secondmoment/D");
     m_trainTree_emconly->Branch("delta_phi",&m_delta_phi,"delta_phi/D");
     m_trainTree_emconly->Branch("delta_theta",&m_delta_theta,"delta_theta/D");
      }
     if(!m_mlp_emc){
    m_mlp_emc = new TMultiLayerPerceptron("ptrk,pt,energy,eseed,e3x3,e5x5,secondmoment,delta_theta,delta_phi:18:type",m_trainTree_emconly);
     m_mlp_emc->LoadWeights(emc.c_str());
   }
   */
 
   /*if(!m_trainTree_emconly){
   m_trainTree_emconly = new TTree("m_trainTree_emconly","m_trainTree_emconly");
   m_trainTree_emconly->Branch("ptrk",&m_ptrk,"ptrk/D");
   m_trainTree_emconly->Branch("pt",&m_pt,"pt/D");
   m_trainTree_emconly->Branch("type",&m_type,"type/D");
   m_trainTree_emconly->Branch("energy",&m_energy,"energy/D");
   m_trainTree_emconly->Branch("eseed",&m_eseed,"eseed/D");
   m_trainTree_emconly->Branch("e3x3",&m_e3x3,"e3x3/D");
   m_trainTree_emconly->Branch("e5x5",&m_e5x5,"e5x5/D");
   m_trainTree_emconly->Branch("secondmoment",&m_secondmoment,"secondmoment/D");
   m_trainTree_emconly->Branch("delta_phi",&m_delta_phi,"delta_phi/D");
   m_trainTree_emconly->Branch("delta_theta",&m_delta_theta,"delta_theta/D");
   }
   if(!m_mlp_emc){
   m_mlp_emc = new TMultiLayerPerceptron("ptrk,pt,energy,eseed,e3x3,e5x5,secondmoment,delta_theta,delta_phi:18:type",m_trainTree_emconly);
   m_mlp_emc->LoadWeights(emc.c_str());
   }
   */
 
 
 
}

Referenced by ParticleID::init().

◆ instance()

EmcPID * EmcPID::instance ( )

static

Definition at line 23 of file EmcPID.cxx.

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

Referenced by ParticleID::init().

◆ IsPidInfoValid()

bool EmcPID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 25 of file EmcPID.h.

25{return (m_ndof> 0);}

◆ LikelihoodCalculation()

int EmcPID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 44 of file EmcPID.h.

44{return -1;}

◆ ndof()

int EmcPID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 28 of file EmcPID.h.

28{return m_ndof;}

◆ neuronPID()

int EmcPID::neuronPID ( ) const

inline

Definition at line 37 of file EmcPID.h.

37{return -1;}

◆ neuronPIDCalculation()

int EmcPID::neuronPIDCalculation ( )

inlineprotected

Definition at line 43 of file EmcPID.h.

43{return -1;}

◆ pars()

double EmcPID::pars ( int n ) const

inline

Definition at line 38 of file EmcPID.h.

38{return params_emc1[n];}

◆ particleIDCalculation()

int EmcPID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 157 of file EmcPID.cxx.

                                  {
   int irc = -1;
   EvtRecTrack* recTrk = PidTrk();
   if(!(recTrk->isMdcTrackValid())) return irc;
   RecMdcTrack* mdcTrk = recTrk->mdcTrack();
   if (!(recTrk)->isMdcKalTrackValid()) return irc;
   RecMdcKalTrack* mdcKalTrk = (recTrk)->mdcKalTrack();
   mdcKalTrk->setPidType(RecMdcKalTrack::electron);
   double ptrk=mdcKalTrk->p();
   // double ptrk = mdcTrk->p();
   m_ptrk = ptrk;
 
   m_pt=mdcKalTrk->pxy();
   m_pt = m_pt*mdcTrk->charge();
   //   double cost = cos(mdcTrk->theta());
 
   if(!(recTrk->isExtTrackValid())) return irc;
   RecExtTrack* extTrk = recTrk->extTrack();
   if(extTrk->emcVolumeNumber() == -1) return irc;
   if (!(recTrk->isEmcShowerValid())) return irc;
   RecEmcShower* emcTrk = recTrk->emcShower();
 
   m_energy = emcTrk->energy();
   m_eseed = emcTrk->eSeed();
   m_e3x3 = emcTrk->e3x3();
   m_e5x5 = emcTrk->e5x5();
 
   double m_emc_theta = emcTrk->theta();
   double m_emc_phi = emcTrk->phi();
 
   Hep3Vector mc = extTrk->emcPosition();
   double m_ext_theta = mc.theta();
   double m_ext_phi = mc.phi();
 
 
   m_delta_theta = m_emc_theta - m_ext_theta;
   m_delta_phi = m_emc_phi - m_ext_phi;
   if(m_delta_phi>1)   m_delta_phi=m_delta_phi-6.283;
   if(m_delta_phi<-1)  m_delta_phi=m_delta_phi+6.283;
 
 
   m_secondmoment = emcTrk->secondMoment()/1000.;
   m_a20moment=emcTrk->a20Moment();
   m_latmoment=emcTrk->latMoment();
   m_a42moment=emcTrk->a42Moment();
 
   if(emcTrk->energy() <= 0) return irc;
   //if(emcTrk->energy() > ptrk) return irc;
 
   /*  params_emc1[0] = m_ptrk;
     params_emc1[1] = m_pt;
     params_emc1[2] = m_energy;
     params_emc1[3] = m_eseed;
     params_emc1[4] = m_e3x3;
     params_emc1[5] = m_e5x5;
     params_emc1[6] = m_secondmoment;
     params_emc1[7] = m_delta_theta;
     params_emc1[8] = m_delta_phi;*/
   params_emc1[0] =m_ptrk;
   params_emc1[1] =m_pt;
   params_emc1[2] =m_energy;
   params_emc1[3] =m_eseed;
   params_emc1[4] =m_e3x3;
   params_emc1[5] =m_e5x5;
   params_emc1[6] =m_secondmoment;
   params_emc1[7] =m_latmoment;
   params_emc1[8] =m_a20moment;
   params_emc1[9] =m_a42moment;
   params_emc1[10] =m_delta_theta;
   params_emc1[11] =m_delta_phi;
 
   m_val_emc = m_mlp_emc->Evaluate(0,params_emc1);
   int pindex = int((m_ptrk-0.2)/0.1);
   int bindex = int((m_val_emc-0.5)/0.01);
   if(bindex>300||bindex<0) return irc;
   if(pindex>17) pindex=17;
   if(pindex<0) pindex=0;
   //   double bin_pos[3];
   m_prob[0] = m_e_h[pindex][bindex];
   m_prob[1] = m_m_h[pindex][bindex];
   m_prob[2] = m_p_h[pindex][bindex];
   m_prob[3] = m_p_h[pindex][bindex];
   m_prob[4] = m_p_h[pindex][bindex];
   for(int i =0; i<5; i++) {
      if(m_prob[i]==0) m_prob[i] = 0.001;
   }
   //calculate the chisq value using GAUSIN
   float ppp[5];
   for(int i=0; i<5; i++) {
      ppp[i]=0;
   }
   for(int j=0; j<=bindex; j++) {
      ppp[0]+= m_e_h[pindex][j];
      ppp[1]+= m_m_h[pindex][j];
      ppp[2]+= m_p_h[pindex][j];
   }
   for(int i=0; i<3; i++) {
      ppp[i]=ppp[i]*0.01;
      if(ppp[i]>0&&ppp[i]<1) {
         CALG(ppp[i],m_chi[i]);
      }
      if(ppp[i]<=0||ppp[i]>=1) m_chi[i]=-99;
   }
   // if(fabs(m_chi[2])==-99)
   m_chi[3]=m_chi[2];
   m_chi[4]=m_chi[2];
 
   m_ndof = 1;
   irc = 0;
   return irc;
}

Referenced by calculate().

◆ prob()

double EmcPID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 27 of file EmcPID.h.

27{return m_prob[n];}

◆ secondmoment()

double EmcPID::secondmoment ( ) const

inline

Definition at line 35 of file EmcPID.h.

35{return m_secondmoment;}

◆ val_emc1()

double EmcPID::val_emc1 ( ) const

inline

Definition at line 36 of file EmcPID.h.

36{return m_val_emc;}

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

◆ ~EmcPID()

Member Function Documentation

◆ calculate()

◆ chi()

◆ delta_phi()

◆ delta_theta()

◆ e3x3()

◆ e5x5()

◆ energy()

◆ eseed()

◆ init()

◆ instance()

◆ IsPidInfoValid()

◆ LikelihoodCalculation()

◆ ndof()

◆ neuronPID()

◆ neuronPIDCalculation()

◆ pars()

◆ particleIDCalculation()

◆ prob()

◆ secondmoment()

◆ val_emc1()