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 25 of file SimplePIDSvc.h.

Constructor & Destructor Documentation

SimplePIDSvc()

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

Definition at line 21 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);
}

~SimplePIDSvc()

SimplePIDSvc::~SimplePIDSvc ( )

virtual

Definition at line 30 of file SimplePIDSvc.cxx.

{;}

Member Function Documentation

finalize()

StatusCode SimplePIDSvc::finalize ( )

virtual

Definition at line 47 of file SimplePIDSvc.cxx.

{
    MsgStream log(messageService(), name());
    log << MSG::INFO << "in SimplePIDSvc finalize()" << endreq;
 
    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 1233 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 51 of file SimplePIDSvc.h.

{ return m_dedx_chi[i]; }

getElectronLikelihood()

double SimplePIDSvc::getElectronLikelihood ( )

inline

Definition at line 67 of file SimplePIDSvc.h.

{ return m_lh_electron; }

getEmcE()

double SimplePIDSvc::getEmcE ( )

inline

Definition at line 62 of file SimplePIDSvc.h.

{ return m_emc_e; }

getEmcE13()

double SimplePIDSvc::getEmcE13 ( )

inline

Definition at line 63 of file SimplePIDSvc.h.

{ return m_emc_e13; }

getEmcE35()

double SimplePIDSvc::getEmcE35 ( )

inline

Definition at line 64 of file SimplePIDSvc.h.

{ return m_emc_e35; }

getEmcLatMoment()

double SimplePIDSvc::getEmcLatMoment ( )

inline

Definition at line 65 of file SimplePIDSvc.h.

{ return m_emc_lat; }

getEmcLikelihood()

double SimplePIDSvc::getEmcLikelihood ( int  i )

inline

Definition at line 53 of file SimplePIDSvc.h.

{ return m_emc_likelihood[i]; }

getEmcSecMoment()

double SimplePIDSvc::getEmcSecMoment ( )

inline

Definition at line 66 of file SimplePIDSvc.h.

{ return m_emc_sec; }

getTofBarrel()

int SimplePIDSvc::getTofBarrel ( )

inline

Definition at line 61 of file SimplePIDSvc.h.

{ return m_tof_barrel; }

getTOFChi()

double SimplePIDSvc::getTOFChi ( int  i )

inline

Definition at line 52 of file SimplePIDSvc.h.

{ return m_tof_chi[i]; }

getTofCounter()

int SimplePIDSvc::getTofCounter ( int  n )

inline

Definition at line 60 of file SimplePIDSvc.h.

{ return m_tof_counter[n]; }

getTofDt()

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

inline

Definition at line 57 of file SimplePIDSvc.h.

{ return m_tof_dt[m][n]; }

getTofPh()

double SimplePIDSvc::getTofPh ( int  n )

inline

Definition at line 58 of file SimplePIDSvc.h.

{ return m_tof_ph[n]; }

getTofZR()

double SimplePIDSvc::getTofZR ( int  n )

inline

Definition at line 59 of file SimplePIDSvc.h.

{ return m_tof_zr[n]; }

initialize()

StatusCode SimplePIDSvc::initialize ( )

virtual

Definition at line 32 of file SimplePIDSvc.cxx.

{
    MsgStream log(messageService(), name());
    log << MSG::INFO << "in SimplePIDSvc initialize()" << endreq;
 
    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 1369 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 1342 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 1334 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 95 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 45 of file SimplePIDSvc.h.

{ return m_prob[0]; }

probKaon()

double SimplePIDSvc::probKaon ( )

inline

Definition at line 48 of file SimplePIDSvc.h.

{ return m_prob[3]; }

probMuon()

double SimplePIDSvc::probMuon ( )

inline

Definition at line 46 of file SimplePIDSvc.h.

{ return m_prob[1]; }

probPion()

double SimplePIDSvc::probPion ( )

inline

Definition at line 47 of file SimplePIDSvc.h.

{ return m_prob[2]; }

probProton()

double SimplePIDSvc::probProton ( )

inline

Definition at line 49 of file SimplePIDSvc.h.

{ return m_prob[4]; }

setdedxminchi()

void SimplePIDSvc::setdedxminchi ( double  x )

inline

Definition at line 37 of file SimplePIDSvc.h.

{ m_dedx_chi_cut = x; }

seteidratio()

void SimplePIDSvc::seteidratio ( double  x )

inline

Definition at line 39 of file SimplePIDSvc.h.

{ m_eid_ratio = x; }

settofminchi()

void SimplePIDSvc::settofminchi ( double  x )

inline

Definition at line 38 of file SimplePIDSvc.h.

{ m_tof_chi_cut = 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