SimplePIDSvc Class Reference

#include <SimplePIDSvc.h>

Inheritance diagram for SimplePIDSvc:

Public Member Functions
	SimplePIDSvc (const std::string &name, ISvcLocator *svcLoc)
virtual	~SimplePIDSvc ()
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvIF)
void	setdedxminchi (double x)
void	settofminchi (double x)
void	preparePID (EvtRecTrack *track)
void	getConst (int run)
void	combineTOF ()
void	calprob ()
bool	iselectron (bool eop=false)
bool	ispion ()
bool	iskaon ()
double	probElectron ()
double	probMuon ()
double	probPion ()
double	probKaon ()
double	probProton ()
double	getdEdxChi (int i)
double	getTOFChi (int i)
double	getChi2 (int i)
Public Member Functions inherited from ISimplePIDSvc
virtual	~ISimplePIDSvc ()
virtual void	setdedxminchi (double x)=0
virtual void	settofminchi (double x)=0
virtual void	preparePID (EvtRecTrack *track)=0
virtual bool	iselectron (bool eop=false)=0
virtual bool	ispion ()=0
virtual bool	iskaon ()=0
virtual double	probElectron ()=0
virtual double	probMuon ()=0
virtual double	probPion ()=0
virtual double	probKaon ()=0
virtual double	probProton ()=0
virtual double	getdEdxChi (int i)=0
virtual double	getTOFChi (int i)=0
virtual double	getChi2 (int i)=0

Friends
class	CnvFactory< SimplePIDSvc >

Additional Inherited Members
Static Public Member Functions inherited from ISimplePIDSvc
static const InterfaceID &	interfaceID ()

Detailed Description

Definition at line 12 of file SimplePIDSvc.h.

Constructor & Destructor Documentation

SimplePIDSvc()

SimplePIDSvc::SimplePIDSvc	(	const std::string &	name,
		ISvcLocator *	svcLoc
	)

Definition at line 18 of file SimplePIDSvc.cxx.

   : Service(name, svcLoc)
{
  declareProperty("mindedxchi", m_dedxminchi=4);
  declareProperty("mintofchi", m_tofminchi=2);
  declareProperty("tofcorrection", m_tofcorrec=true);
}

~SimplePIDSvc()

SimplePIDSvc::~SimplePIDSvc ( )

virtual

Definition at line 26 of file SimplePIDSvc.cxx.

{
}

Member Function Documentation

calprob()

void SimplePIDSvc::calprob

(

)

Definition at line 430 of file SimplePIDSvc.cxx.

                          {
  
  bool usededx=false;
  bool usetof=false;
  
  for(int i=0; i<5 ;i++){
    
    if(!usededx && fabs(m_dedxchi[i])<m_dedxminchi)
      usededx=true;
    if(!usetof && fabs(m_tofchi[i])<m_tofminchi)
      usetof=true;
  
    m_dedxonly[i]=false;
  }
  
  if(!usededx){
 
    for(int i=0; i<5; i++)
      m_dedxchi[i]=-99;
  }
  
  if(!usetof){
 
    for(int i=0; i<5; i++)
      m_tofchi[i]=-99;
  }
 
  /*
  if( fabs(m_p[2])<0.3)
    usetof=false;
  */
 
  for(int i=0; i<5; ++i){
    m_prob[i]=-99;
    double chi2=0;
    int ndf=0;
 
 
    if( usededx && usetof ){
      chi2= pow(m_dedxchi[i],2)+pow(m_tofchi[i],2);
      ndf=2;
    }
    else if( usededx && !usetof ){
      chi2= pow(m_dedxchi[i],2);
      ndf=1;
      m_dedxonly[i]=true;
    }
    else if( !usededx && usetof ){
      chi2= pow(m_tofchi[i],2);
      ndf=1;
    }
    /*
    if( fabs(m_dedxchi[i])<10 && fabs(m_tofchi[i])<10 ){
      chi2= pow(m_dedxchi[i],2)+pow(m_tofchi[i],2);
      ndf=2;
    }
    else if( fabs(m_dedxchi[i])<10 && fabs(m_tofchi[i])>10 ){
      chi2= pow(m_dedxchi[i],2);
      ndf=1;
      m_dedxonly[i]=true;
    }
    else if( fabs(m_dedxchi[i])>10 && fabs(m_tofchi[i])<10 ){
      chi2= pow(m_tofchi[i],2);
      ndf=1;
    }
    
    */
 
    if(ndf>0 && chi2>0)
      m_prob[i]=TMath::Prob(chi2, ndf);
    
  }
  
  
 
}

Referenced by preparePID().

combineTOF()

void SimplePIDSvc::combineTOF

(

)

Definition at line 372 of file SimplePIDSvc.cxx.

                             {
  
  
  for(int i=0; i<5 ; ++i){
    
    double dt1= m_tofdt1[i];
    double dt2= m_tofdt2[i];
    double sigma1= m_sigma1;
    double sigma2= m_sigma2;
    
    double sqrtcov=0.041;
    if(m_tofcorrec){
      sqrtcov=tofsqrcov(m_p[i]/m_mass[i],m_datatof);
      if(m_run<0)
    sqrtcov=mctofsqrcov(m_p[i]/m_mass[i],m_mctof);
    }
        
    double weight1=(sigma2*sigma2-sqrtcov*sqrtcov)/(sigma1*sigma1+sigma2*sigma2-2*sqrtcov*sqrtcov);
    double weight2= (sigma1*sigma1-sqrtcov*sqrtcov)/(sigma1*sigma1+sigma2*sigma2-2*sqrtcov*sqrtcov);
    
    double dt_comb=-99;
    double sigma_comb=-99;
    
    //  if(sigma1==0 || sigma2==0) cout<<"zero sigma"<<endl;
 
    if( dt1!=-99 && dt2!=-99){
            
      dt_comb= weight1*dt1+weight2*dt2;
      sigma_comb= sqrt((sigma1*sigma1*sigma2*sigma2 - sqrtcov*sqrtcov*sqrtcov*sqrtcov)/(sigma1*sigma1+sigma2*sigma2-2*sqrtcov*sqrtcov));
    }
    else if(dt1!=-99 && dt2==-99){
      dt_comb = dt1;
      sigma_comb=sigma1;
    }
    else if ( dt1==-99 && dt2!=-99){
      dt_comb = dt2;
      sigma_comb=sigma2;
    }
    
    if(sigma_comb>0){
      double sigmascale=1.0;
      if(m_tofcorrec){
    sigmascale=tofsigmascale(m_p[i]/m_mass[i],m_datatof);
    if(m_run<0)
      sigmascale=mctofsigmascale(m_p[i]/m_mass[i],m_mctof);
      }
      
      m_tofchi[i]= dt_comb/sigma_comb/sigmascale;
 
    }
    else
      m_tofchi[i]=-99;
  }
 
  
  
}

Referenced by preparePID().

finalize()

StatusCode SimplePIDSvc::finalize ( )

virtual

Definition at line 59 of file SimplePIDSvc.cxx.

{
   MsgStream log(messageService(), name());
   log << MSG::INFO << "@initialize()" << endreq;
 
   StatusCode sc = Service::finalize();
 
   delete h_ebarlp;
   delete h_ebarhp;
   delete h_eendlp;
   delete h_eendhp;
   delete h_kbar;
   delete h_kend;
   delete h_pibar;
   delete h_piend;
 
 
   return sc;
}

getChi2()

double SimplePIDSvc::getChi2 ( int  i )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 45 of file SimplePIDSvc.h.

{return pow(m_dedxchi[i],2)+pow(m_tofchi[i],2);}

getConst()

void SimplePIDSvc::getConst

(

int

run

)

Definition at line 91 of file SimplePIDSvc.cxx.

                                  {
  string filename;
  //filename = string(getenv("SIMPLEPIDSVCROOT"))+"/share/constant/"+string(getenv("BES_RELEASE"))+"_const.txt";
  filename = string(getenv("SIMPLEPIDSVCROOT"))+"/share/constant/const.txt";
  //cout<<"const file is:"<<filename<<endl;
  
  m_datatof.clear();
  m_mctof.clear();
  run = fabs(run);
  std::ifstream fs(filename.c_str());
  double tag, num;
  int run_begin, run_end, method;
  bool isdata=false, ismc=false;
  while(fs.good() && !fs.eof() ){
    
    fs >> num;
    if(num==999){
      isdata=false;
      ismc=false;
      fs >> tag;
      fs >> run_begin;
      fs >> run_end;
      fs >> method;
      if(tag==99 && run>=run_begin && run<=run_end){
    isdata=true;
    ismc=false;
      }
      if(tag==-99 && run>=run_begin && run<=run_end){
    isdata=false;
    ismc=true;
      }
      fs >> num;
    }
    
    if(isdata)
      m_datatof.push_back(num);
    if(ismc)
      m_mctof.push_back(num);
    
  }//end of while
 
}

Referenced by preparePID().

getdEdxChi()

double SimplePIDSvc::getdEdxChi ( int  i )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 43 of file SimplePIDSvc.h.

{return m_dedxchi[i];}

getTOFChi()

double SimplePIDSvc::getTOFChi ( int  i )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 44 of file SimplePIDSvc.h.

{return m_tofchi[i];}

initialize()

StatusCode SimplePIDSvc::initialize ( )

virtual

Definition at line 30 of file SimplePIDSvc.cxx.

{
   MsgStream log(messageService(), name());
   log << MSG::INFO << "@initialize()" << endreq;
 
   StatusCode sc = Service::initialize();
 
   sc = serviceLocator()->service("EventDataSvc", eventSvc_, true);
   //sc = serviceLocator()->service("VertexDbSvc", vtxSvc_, true);
 
   //m_dedxminchi=4;
   //m_tofminchi=2;
   
   //get eop files
   string filename;
   filename = string(getenv("SIMPLEPIDSVCROOT"))+"/share/eop.root";
   TFile* file=new TFile(filename.c_str(),"read");
   h_ebarlp=(TH1F*)file->Get("barlp");
   h_ebarhp=(TH1F*)file->Get("barhp");
   h_eendlp=(TH1F*)file->Get("endlp");
   h_eendhp=(TH1F*)file->Get("endhp");
   h_kbar=(TH1F*)file->Get("barkaon");
   h_kend=(TH1F*)file->Get("endkaon");
   h_pibar=(TH1F*)file->Get("barpion");
   h_piend=(TH1F*)file->Get("endpion");
   
   return sc;
}

iselectron()

bool SimplePIDSvc::iselectron ( bool  eop = false )

virtual

Implements ISimplePIDSvc.

Definition at line 508 of file SimplePIDSvc.cxx.

                                     {
 
 
  if(m_prob[0]>0 && m_prob[0]>m_prob[2] && m_prob[0]>m_prob[3]){
   
    if(!eop)
      return true;
    else{
      
      if( fabs(m_costh[0])<0.7 && m_eop>0 && m_eop<0.8)
    return false;
      
      if( fabs(m_costh[0])>=0.7 &&   fabs(m_costh[0])<0.8 && m_eop>0 && m_eop<-7.5*fabs(m_costh[0])+6.05)
    return false;
      
      if( fabs(m_costh[0])>0.85 && m_eop>0 && m_eop<0.6)
    return false;
     
      return true;
      
    }
 
  }
  return false;
    
  
  /*
  double pe=-99, pk=-99, ppi=-99;
  
  if( fabs(m_dedxchi[0])>5 )
    return false;
 
  if(m_eop>0.5 && m_eop<1.5){
    if( fabs(m_costh[0])<0.81 ){
      pk=h_kbar->GetBinContent(h_kbar->FindBin(m_eop));
      ppi=h_pibar->GetBinContent(h_pibar->FindBin(m_eop));
      if(fabs(m_p[0])<1.0)
    pe=h_ebarlp->GetBinContent(h_ebarlp->FindBin(m_eop));
      else if(fabs(m_p[0])>=1.0)
    pe=h_ebarhp->GetBinContent(h_ebarhp->FindBin(m_eop));
        
    }
    else if(fabs(m_costh[0])>0.85 ){
      pk=h_kend->GetBinContent(h_kend->FindBin(m_eop));
      ppi=h_piend->GetBinContent(h_piend->FindBin(m_eop));
      if(fabs(m_p[0])<1.0)
    pe=h_eendlp->GetBinContent(h_eendlp->FindBin(m_eop));
      else if(fabs(m_p[0])>=1.0)
    pe=h_eendhp->GetBinContent(h_eendhp->FindBin(m_eop));
    }
  }
 
  double lr=0;
 
  if( pe>0 ){
    
    if( fabs(m_dedxchi[0])<5 && fabs(m_tofchi[0])<5){
      
      lr= pe*TMath::Gaus(m_dedxchi[0])*TMath::Gaus(m_tofchi[0])/(pe*TMath::Gaus(m_dedxchi[0])*TMath::Gaus(m_tofchi[0])+
                                 0.9*ppi*TMath::Gaus(m_dedxchi[2])*TMath::Gaus(m_tofchi[2])+
                                 0.1*pk*TMath::Gaus(m_dedxchi[3])*TMath::Gaus(m_tofchi[3]));
    }
    
    else if( fabs(m_dedxchi[0])<5 && fabs(m_tofchi[0])>=5){
      
      lr= pe*TMath::Gaus(m_dedxchi[0])/(pe*TMath::Gaus(m_dedxchi[0])+
                     0.9*ppi*TMath::Gaus(m_dedxchi[2])+
                     0.1*pk*TMath::Gaus(m_dedxchi[3]));
    }
    
    else if( fabs(m_tofchi[0])<5 && fabs(m_dedxchi[0])>=5){
          
      lr= pe*TMath::Gaus(m_tofchi[0])/(pe*TMath::Gaus(m_tofchi[0])+
                    0.9*ppi*TMath::Gaus(m_tofchi[2])+
                    0.1*pk*TMath::Gaus(m_tofchi[3]));
    }
       
 
  }
  else{
    if( fabs(m_dedxchi[0])<5 && fabs(m_tofchi[0])<5){
          
      lr= TMath::Gaus(m_dedxchi[0])*TMath::Gaus(m_tofchi[0])/(TMath::Gaus(m_dedxchi[0])*TMath::Gaus(m_tofchi[0])+
                                   0.9*TMath::Gaus(m_dedxchi[2])*TMath::Gaus(m_tofchi[2])+
                                   0.1*TMath::Gaus(m_dedxchi[3])*TMath::Gaus(m_tofchi[3]));
    }
    
    else if( fabs(m_dedxchi[0])<5 && fabs(m_tofchi[0])>=5){
          
      lr= TMath::Gaus(m_dedxchi[0])/(TMath::Gaus(m_dedxchi[0])+
                     0.9*TMath::Gaus(m_dedxchi[2])+
                     0.1*TMath::Gaus(m_dedxchi[3]));
    }
    
    else if( fabs(m_tofchi[0])<5 && fabs(m_dedxchi[0])>=5){
          
      lr= TMath::Gaus(m_tofchi[0])/(TMath::Gaus(m_tofchi[0])+
                    0.9*TMath::Gaus(m_tofchi[2])+
                    0.1*TMath::Gaus(m_tofchi[3]));
    }
 
  }
  
  if( lr>0.8 )
    return true;
  return false;
  */
 
}

iskaon()

bool SimplePIDSvc::iskaon ( )

virtual

Implements ISimplePIDSvc.

Definition at line 631 of file SimplePIDSvc.cxx.

                         {
  
  if(m_prob[3]>=0.00 && m_prob[3]>=m_prob[2]){
    return true;
  }
  /*
  if(m_dedxonly[3] && m_prob[3]>0.10){
    return true;
  }
  */
  
  return false;
}

ispion()

bool SimplePIDSvc::ispion ( )

virtual

Implements ISimplePIDSvc.

Definition at line 618 of file SimplePIDSvc.cxx.

                         {
   
  if(m_prob[2]>=0.00 && m_prob[2]>=m_prob[3])
    return true;
 
  /*
  if(m_dedxonly[2] && m_prob[2]>0.10){
    return true;
  }
  */
  return false;
}

preparePID()

void SimplePIDSvc::preparePID ( EvtRecTrack *  track )

virtual

Implements ISimplePIDSvc.

Definition at line 133 of file SimplePIDSvc.cxx.

                                               {
 
 
  SmartDataPtr<Event::EventHeader> eventHeaderpid(eventSvc_,"/Event/EventHeader");
  m_run=eventHeaderpid->runNumber();
 
  getConst(m_run);
  //cout<<"size of constant:"<<m_mctof.size()<<endl;
  //for(int i=0; i<m_mctof.size();i++)
  //  cout<<m_mctof[i]<<endl;
  if(m_datatof.size()==0 || m_mctof.size()==0)
    m_tofcorrec=false;
 
  //cout<<"track id :"<<track->trackId()<<endl;
  //track quantity
  if( track->isMdcKalTrackValid()) {    
    RecMdcKalTrack *mdckalTrk = track->mdcKalTrack();
    
    for(int pid=0; pid<5; ++pid){
      HepVector zhelix;
      if(pid==0){
    mdckalTrk->setPidType(RecMdcKalTrack::electron);    
    zhelix=mdckalTrk->getZHelixE();
      }
      else if(pid==1){
    mdckalTrk->setPidType(RecMdcKalTrack::muon);
    zhelix=mdckalTrk->getZHelixMu();
      }
      else if(pid==2){
    mdckalTrk->setPidType(RecMdcKalTrack::pion);
    zhelix=mdckalTrk->getZHelix();
      }
      else if(pid==3){
    mdckalTrk->setPidType(RecMdcKalTrack::kaon);
    zhelix=mdckalTrk->getZHelixK();
      }
      else{
    mdckalTrk->setPidType(RecMdcKalTrack::proton);
    zhelix=mdckalTrk->getZHelixP();
      }
 
      double kappa=zhelix[2];
      double theta=CLHEP::pi/2.0-atan(zhelix[4]);
      m_p[pid]=1.0/fabs(kappa)/sin(theta);
      m_costh[pid]=cos(theta);
    }
  }
  else{
    for(int i=0; i<5; ++i){
      m_p[i]=-99;
      m_costh[i]=-99;
    }
 
  }
  
  
  //mass
  m_mass[0]=0.000511;
  m_mass[1]=0.105658;
  m_mass[2]=0.14;
  m_mass[3]=0.494;
  m_mass[4]=0.94;
 
  //tof correction factors
  for(int i=0; i<5; ++i){
    if(m_tofcorrec){
      if(m_run>0){
    m_tofscale1[i]=toflayer1scale(m_p[i]/m_mass[i],fabs(m_costh[i]),m_datatof);
    m_tofscale2[i]=toflayer2scale(m_p[i]/m_mass[i],fabs(m_costh[i]),m_datatof);
      }
      else{
    m_tofscale1[i]=mctoflayer1scale(m_p[i]/m_mass[i],fabs(m_costh[i]),m_mctof);
    m_tofscale2[i]=mctoflayer2scale(m_p[i]/m_mass[i],fabs(m_costh[i]),m_mctof);
      }
    }
    else{
      m_tofscale1[i]=1.0;
      m_tofscale2[i]=1.0;
    }
  }
    
  
  //eop for electrons
  m_eop=-99;
  if( track->isEmcShowerValid() ){
    RecEmcShower* emcTrk = track->emcShower();
    double energy= emcTrk->energy();
    if(fabs(m_p[0])>0)
      m_eop=energy/fabs(m_p[0]);
  }
 
  
  //dE/dx PID
  
  if( track->isMdcDedxValid()){
    RecMdcDedx* dedxTrk = track->mdcDedx();
    m_dedxchi[0]= dedxTrk->chi(0);
    m_dedxchi[1]= dedxTrk->chi(1);
    m_dedxchi[2]= dedxTrk->chi(2);
    m_dedxchi[3]= dedxTrk->chi(3);
    m_dedxchi[4]= dedxTrk->chi(4);    
  }
  else{
    m_dedxchi[0]=-99;
    m_dedxchi[1]=-99;
    m_dedxchi[2]=-99;
    m_dedxchi[3]=-99;
    m_dedxchi[4]=-99;
  }
 
  
  //TOF PID
  
  m_tofchi[0]=-99;
  m_tofchi[1]=-99;
  m_tofchi[2]=-99;
  m_tofchi[3]=-99;
  m_tofchi[4]=-99;
 
  m_tofdt[0]=-99;
  m_tofdt[1]=-99;
  m_tofdt[2]=-99;
  m_tofdt[3]=-99;
  m_tofdt[4]=-99;
 
 
  m_tofdt1[0]=-99;
  m_tofdt1[1]=-99;
  m_tofdt1[2]=-99;
  m_tofdt1[3]=-99;
  m_tofdt1[4]=-99;
 
  m_tofdt2[0]=-99;
  m_tofdt2[1]=-99;
  m_tofdt2[2]=-99;
  m_tofdt2[3]=-99;
  m_tofdt2[4]=-99;
  
  int layer1size=0;
  int layer2size=0;
  
 
  if(track->isTofTrackValid() ){
    SmartRefVector<RecTofTrack> tofTrkCol=track->tofTrack();
    SmartRefVector<RecTofTrack>::iterator iter_tof=tofTrkCol.begin();
 
    for(;iter_tof!=tofTrkCol.end();iter_tof++){
      TofHitStatus* status =new TofHitStatus;
      status->setStatus( (*iter_tof)->status() );
      
      /*
      if( status->is_cluster()){
        double  time=(*iter_tof)->tof();
    double sigma=(*iter_tof)->sigma(0);
    double expE_tof=(*iter_tof)->texpElectron();
    double expMu_tof=(*iter_tof)->texpMuon();
        double expPi_tof=(*iter_tof)->texpPion();
        double expK_tof=(*iter_tof)->texpKaon();
    double expP_tof=(*iter_tof)->texpProton();
    
    m_tofchi[0]=  (time- expE_tof)/sigma  ;
    m_tofchi[1]=  (time- expMu_tof)/sigma  ;
    m_tofchi[2]=  (time- expPi_tof)/sigma  ;
    m_tofchi[3]=  (time- expK_tof)/sigma  ;
    m_tofchi[4]=  (time- expP_tof)/sigma  ;
 
    m_tofdt[0]=  (time- expE_tof);
    m_tofdt[1]=  (time- expMu_tof);
    m_tofdt[2]=  (time- expPi_tof);
    m_tofdt[3]=  (time- expK_tof);
    m_tofdt[4]=  (time- expP_tof);
    
      }
      */
      
      //two layers
      
      if(!(status->is_barrel())){//endcap
    
        if( status->is_counter() && status->layer()==1 ) {
      layer1size++;
      for(int i=0; i<5; i++)
        m_tofdt1[i]=(*iter_tof)->tof() - m_tofscale1[i]*(*iter_tof)->texp(i);  
      m_sigma1=(*iter_tof)->sigma(0);;
    }
      }
      else if( status->is_barrel()) { //barrel
        if( status->is_counter() ){
      if(status->layer()==1){ //layer1
        layer1size++;
 
        for(int i=0; i<5; i++){
          m_tofdt1[i]=(*iter_tof)->tof() - m_tofscale1[i]*(*iter_tof)->texp(i);
        }
        m_sigma1=(*iter_tof)->sigma(0);;
      }
      else if(status->layer()==2){//layer2
        layer2size++;
        
        for(int i=0; i<5; i++){
          m_tofdt2[i]=(*iter_tof)->tof() - m_tofscale2[i]*(*iter_tof)->texp(i);
        }
        m_sigma2=(*iter_tof)->sigma(0);;
      }
      
    }
    
      }//end of barrel
      
      
 
 
      delete status;
    }
 
  }//end of tof
  
  if(layer1size>1){
    m_tofdt1[0]=-99;
    m_tofdt1[1]=-99;
    m_tofdt1[2]=-99;
    m_tofdt1[3]=-99;
    m_tofdt1[4]=-99;
  }
  if(layer2size>1){
    m_tofdt2[0]=-99;
    m_tofdt2[1]=-99;
    m_tofdt2[2]=-99;
    m_tofdt2[3]=-99;
    m_tofdt2[4]=-99;
  }
  
  //combine two layers
  
  combineTOF();
  calprob();
  
}

probElectron()

double SimplePIDSvc::probElectron ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 37 of file SimplePIDSvc.h.

{ return m_prob[0];};

probKaon()

double SimplePIDSvc::probKaon ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 40 of file SimplePIDSvc.h.

{ return m_prob[3];};

probMuon()

double SimplePIDSvc::probMuon ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 38 of file SimplePIDSvc.h.

{ return m_prob[1];};

probPion()

double SimplePIDSvc::probPion ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 39 of file SimplePIDSvc.h.

{ return m_prob[2];};

probProton()

double SimplePIDSvc::probProton ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 41 of file SimplePIDSvc.h.

{ return m_prob[4];};

queryInterface()

StatusCode SimplePIDSvc::queryInterface ( const InterfaceID &  riid, void **  ppvIF  )

virtual

Definition at line 79 of file SimplePIDSvc.cxx.

{
   if ( ISimplePIDSvc::interfaceID().versionMatch(riid) ) {
      *ppvIF = dynamic_cast<ISimplePIDSvc*>(this);
   }
   else {
      return Service::queryInterface(riid, ppvIF);
   }
   addRef();
   return StatusCode::SUCCESS;
}

setdedxminchi()

void SimplePIDSvc::setdedxminchi ( double  x )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 26 of file SimplePIDSvc.h.

{m_dedxminchi=x;}

settofminchi()

void SimplePIDSvc::settofminchi ( double  x )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 27 of file SimplePIDSvc.h.

{m_tofminchi=x;}

Friends And Related Function Documentation

CnvFactory< SimplePIDSvc >

friend class CnvFactory< SimplePIDSvc >

friend

Definition at line 1 of file SimplePIDSvc.h.

---

The documentation for this class was generated from the following files:

• SimplePIDSvc.h
• SimplePIDSvc.cxx