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 ()
void	setdedxminchi (double x)
void	settofminchi (double x)
void	seteidratio (double x)
void	preparePID (EvtRecTrack *track)
bool	iselectron (bool emc=true)
bool	ispion ()
bool	iskaon ()
double	probElectron ()
double	probMuon ()
double	probPion ()
double	probKaon ()
double	probProton ()
double	getdEdxChi (int i)
double	getTOFChi (int i)
double	getEmcLikelihood (int i)
double	getChi2 (int i)
double	getTofDt (int m, int n)
double	getTofPh (int n)
double	getTofZR (int n)
int	getTofCounter (int n)
int	getTofBarrel ()
double	getEmcE ()
double	getEmcE13 ()
double	getEmcE35 ()
double	getEmcLatMoment ()
double	getEmcSecMoment ()
double	getElectronLikelihood ()

Friends
class	CnvFactory< SimplePIDSvc >

Detailed Description

Definition at line 27 of file SimplePIDSvc.h.

Constructor & Destructor Documentation

SimplePIDSvc()

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

Definition at line 19 of file SimplePIDSvc.cxx.

    : base_class( name, svcLoc ) {
  declareProperty( "DedxChiCut", m_dedx_chi_cut = 4 );
  declareProperty( "TofChiCut", m_tof_chi_cut = 4 );
  declareProperty( "IsTofCorr", m_tof_corr = true );
  declareProperty( "IsDedxCorr", m_dedx_corr = true );
  declareProperty( "EidRatio", m_eid_ratio = 0.80 );
}

Referenced by SimplePIDSvc().

~SimplePIDSvc()

SimplePIDSvc::~SimplePIDSvc ( )

virtual

Definition at line 28 of file SimplePIDSvc.cxx.

{ ; }

Member Function Documentation

finalize()

StatusCode SimplePIDSvc::finalize ( )

virtual

Definition at line 44 of file SimplePIDSvc.cxx.

                                  {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in SimplePIDSvc finalize()" << endmsg;
 
  StatusCode sc = Service::finalize();
 
  for ( unsigned int i = 0; i < 2; i++ )
  {
    for ( unsigned int j = 0; j < 4; j++ )
    {
      f_dedx[i][j]->Close();
      f_tof_q[i][j]->Close();
      f_tof_bgcost[i][j]->Close();
      f_tof_wgt[i][j]->Close();
      f_tof_final[i][j]->Close();
      f_tofec_q[i][j]->Close();
      f_tofec_bg[i][j]->Close();
      f_tofec_cost[i][j]->Close();
    }
  }
  for ( unsigned int i = 0; i < 3; i++ )
  {
    for ( unsigned int j = 0; j < 4; j++ ) { f_emc[i][j]->Close(); }
  }
 
  return sc;
}

getChi2()

double SimplePIDSvc::getChi2

(

int

i

)

Definition at line 1557 of file SimplePIDSvc.cxx.

                                    {
  return pow( m_dedx_chi[i], 2 ) + pow( m_tof_chi[i], 2 );
}

getdEdxChi()

double SimplePIDSvc::getdEdxChi ( int i )

inline

Definition at line 52 of file SimplePIDSvc.h.

{ return m_dedx_chi[i]; }

getElectronLikelihood()

double SimplePIDSvc::getElectronLikelihood ( )

inline

Definition at line 68 of file SimplePIDSvc.h.

{ return m_lh_electron; }

getEmcE()

double SimplePIDSvc::getEmcE ( )

inline

Definition at line 63 of file SimplePIDSvc.h.

{ return m_emc_e; }

getEmcE13()

double SimplePIDSvc::getEmcE13 ( )

inline

Definition at line 64 of file SimplePIDSvc.h.

{ return m_emc_e13; }

getEmcE35()

double SimplePIDSvc::getEmcE35 ( )

inline

Definition at line 65 of file SimplePIDSvc.h.

{ return m_emc_e35; }

getEmcLatMoment()

double SimplePIDSvc::getEmcLatMoment ( )

inline

Definition at line 66 of file SimplePIDSvc.h.

{ return m_emc_lat; }

getEmcLikelihood()

double SimplePIDSvc::getEmcLikelihood ( int i )

inline

Definition at line 54 of file SimplePIDSvc.h.

{ return m_emc_likelihood[i]; }

getEmcSecMoment()

double SimplePIDSvc::getEmcSecMoment ( )

inline

Definition at line 67 of file SimplePIDSvc.h.

{ return m_emc_sec; }

getTofBarrel()

int SimplePIDSvc::getTofBarrel ( )

inline

Definition at line 62 of file SimplePIDSvc.h.

{ return m_tof_barrel; }

getTOFChi()

double SimplePIDSvc::getTOFChi ( int i )

inline

Definition at line 53 of file SimplePIDSvc.h.

{ return m_tof_chi[i]; }

getTofCounter()

int SimplePIDSvc::getTofCounter ( int n )

inline

Definition at line 61 of file SimplePIDSvc.h.

{ return m_tof_counter[n]; }

getTofDt()

double SimplePIDSvc::getTofDt ( int m, int n )

inline

Definition at line 58 of file SimplePIDSvc.h.

{ return m_tof_dt[m][n]; }

getTofPh()

double SimplePIDSvc::getTofPh ( int n )

inline

Definition at line 59 of file SimplePIDSvc.h.

{ return m_tof_ph[n]; }

getTofZR()

double SimplePIDSvc::getTofZR ( int n )

inline

Definition at line 60 of file SimplePIDSvc.h.

{ return m_tof_zr[n]; }

initialize()

StatusCode SimplePIDSvc::initialize ( )

virtual

Definition at line 30 of file SimplePIDSvc.cxx.

                                    {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in SimplePIDSvc initialize()" << endmsg;
 
  StatusCode sc = Service::initialize();
 
  sc = serviceLocator()->service( "EventDataSvc", eventSvc_, true );
 
  loadHistogram();
  loadSecondPar(); // the second correct factor
 
  return sc;
}

iselectron()

bool SimplePIDSvc::iselectron

(

bool

emc = true

)

Definition at line 1675 of file SimplePIDSvc.cxx.

                                        {
  if ( !emc )
  {
    if ( m_prob[0] > 0 && m_prob[0] > m_prob[2] && m_prob[0] > m_prob[3] ) return true;
    else return false;
  }
  else
  {
    if ( calEMCLikelihood() )
    {
      if ( m_lh_electron > m_eid_ratio ) return true;
      else return false;
    }
    else
    {
      if ( m_prob[0] > 0 && m_prob[0] > m_prob[2] && m_prob[0] > m_prob[3] ) return true;
      else return false;
    }
  }
}

iskaon()

bool SimplePIDSvc::iskaon

(

)

Definition at line 1655 of file SimplePIDSvc.cxx.

                          {
  if ( m_prob[3] > 0.00 && m_prob[3] > m_prob[2] ) return true;
  else return false;
}

ispion()

bool SimplePIDSvc::ispion

(

)

Definition at line 1650 of file SimplePIDSvc.cxx.

                          {
  if ( m_prob[2] > 0.00 && m_prob[2] > m_prob[3] ) return true;
  else return false;
}

preparePID()

void SimplePIDSvc::preparePID

(

EvtRecTrack *

track

)

Definition at line 87 of file SimplePIDSvc.cxx.

                                                  {
 
  SmartDataPtr<Event::EventHeader> eventHeaderpid( eventSvc_, "/Event/EventHeader" );
  m_run = eventHeaderpid->runNumber();
 
  if ( track->isMdcKalTrackValid() )
  {
    RecMdcKalTrack*         mdckalTrk   = track->mdcKalTrack();
    RecMdcKalTrack::PidType trk_type[5] = {
        RecMdcKalTrack::electron, RecMdcKalTrack::muon,   RecMdcKalTrack::pion,
        RecMdcKalTrack::kaon,     RecMdcKalTrack::proton,
    };
    double mass[5] = {
        0.000511, 0.105658, 0.13957, 0.493677, 0.938272,
    };
    for ( unsigned int pid = 0; pid < 5; pid++ )
    {
      mdckalTrk->setPidType( trk_type[pid] );
      m_p[pid]         = mdckalTrk->p();
      m_betagamma[pid] = m_p[pid] / mass[pid];
      m_charge[pid]    = mdckalTrk->charge();
      m_cost[pid]      = cos( mdckalTrk->theta() );
    }
  }
  else
  {
    for ( unsigned int i = 0; i < 5; i++ )
    {
      m_p[i]         = -99;
      m_betagamma[i] = -99;
      m_cost[i]      = -99;
      m_charge[i]    = 0;
    }
  }
 
  // dE/dx PID
  loadDedxInfo( track );
  if ( m_dedx_corr )
  {
    dedxCorrection();
    dedxSecondCorrection();
  }
  // TOF PID
  loadTOFInfo( track );
  if ( m_tof_corr )
  {
    if ( m_tof_barrel == 1 )
    {
      tofBarrelCorrection();
      tofBarrelSecondCorrection();
    }
    else if ( m_tof_barrel == 0 )
    {
      tofEndcapCorrection();
      tofEndcapSecondCorrection();
    }
  }
  // EMC
  loadEMCInfo( track );
 
  calprob();
}

probElectron()

double SimplePIDSvc::probElectron ( )

inline

Definition at line 46 of file SimplePIDSvc.h.

{ return m_prob[0]; }

probKaon()

double SimplePIDSvc::probKaon ( )

inline

Definition at line 49 of file SimplePIDSvc.h.

{ return m_prob[3]; }

probMuon()

double SimplePIDSvc::probMuon ( )

inline

Definition at line 47 of file SimplePIDSvc.h.

{ return m_prob[1]; }

probPion()

double SimplePIDSvc::probPion ( )

inline

Definition at line 48 of file SimplePIDSvc.h.

{ return m_prob[2]; }

probProton()

double SimplePIDSvc::probProton ( )

inline

Definition at line 50 of file SimplePIDSvc.h.

{ return m_prob[4]; }

setdedxminchi()

void SimplePIDSvc::setdedxminchi ( double x )

inline

Definition at line 38 of file SimplePIDSvc.h.

{ m_dedx_chi_cut = x; }

seteidratio()

void SimplePIDSvc::seteidratio ( double x )

inline

Definition at line 40 of file SimplePIDSvc.h.

{ m_eid_ratio = x; }

settofminchi()

void SimplePIDSvc::settofminchi ( double x )

inline

Definition at line 39 of file SimplePIDSvc.h.

{ m_tof_chi_cut = x; }

Friends And Related Symbol 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