#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	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 ()

Public Member Functions inherited from ISimplePIDSvc
virtual	~ISimplePIDSvc ()

virtual void	setdedxminchi (double x)=0

virtual void	settofminchi (double x)=0

virtual void	seteidratio (double x)=0

virtual void	preparePID (EvtRecTrack *track)=0

virtual bool	iselectron (bool emc=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

virtual double	getEmcLikelihood (int i)=0

virtual double	getTofDt (int m, int n)=0

virtual double	getTofPh (int n)=0

virtual double	getTofZR (int n)=0

virtual int	getTofCounter (int n)=0

virtual int	getTofBarrel ()=0

virtual double	getEmcE ()=0

virtual double	getEmcE13 ()=0

virtual double	getEmcE35 ()=0

virtual double	getEmcLatMoment ()=0

virtual double	getEmcSecMoment ()=0

virtual double	getElectronLikelihood ()=0

Friends
class	CnvFactory< SimplePIDSvc >

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

Detailed Description

Definition at line 24 of file SimplePIDSvc.h.

Constructor & Destructor Documentation

◆ SimplePIDSvc()

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

Definition at line 20 of file SimplePIDSvc.cxx.

                                                                     : Service(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 29 of file SimplePIDSvc.cxx.

29{;}

Member Function Documentation

◆ finalize()

StatusCode SimplePIDSvc::finalize ( )

virtual

Definition at line 45 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 )

virtual

Implements ISimplePIDSvc.

Definition at line 911 of file SimplePIDSvc.cxx.

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

◆ getdEdxChi()

double SimplePIDSvc::getdEdxChi ( int i )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 50 of file SimplePIDSvc.h.

50{ return m_dedx_chi[i]; }

◆ getElectronLikelihood()

double SimplePIDSvc::getElectronLikelihood ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 66 of file SimplePIDSvc.h.

66{ return m_lh_electron; }

◆ getEmcE()

double SimplePIDSvc::getEmcE ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 61 of file SimplePIDSvc.h.

61{ return m_emc_e; }

◆ getEmcE13()

double SimplePIDSvc::getEmcE13 ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 62 of file SimplePIDSvc.h.

62{ return m_emc_e13; }

◆ getEmcE35()

double SimplePIDSvc::getEmcE35 ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 63 of file SimplePIDSvc.h.

63{ return m_emc_e35; }

◆ getEmcLatMoment()

double SimplePIDSvc::getEmcLatMoment ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 64 of file SimplePIDSvc.h.

64{ return m_emc_lat; }

◆ getEmcLikelihood()

double SimplePIDSvc::getEmcLikelihood ( int i )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 52 of file SimplePIDSvc.h.

52{ return m_emc_likelihood[i]; }

◆ getEmcSecMoment()

double SimplePIDSvc::getEmcSecMoment ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 65 of file SimplePIDSvc.h.

65{ return m_emc_sec; }

◆ getTofBarrel()

int SimplePIDSvc::getTofBarrel ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 60 of file SimplePIDSvc.h.

60{ return m_tof_barrel; }

◆ getTOFChi()

double SimplePIDSvc::getTOFChi ( int i )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 51 of file SimplePIDSvc.h.

51{ return m_tof_chi[i]; }

◆ getTofCounter()

int SimplePIDSvc::getTofCounter ( int n )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 59 of file SimplePIDSvc.h.

59{ return m_tof_counter[n]; }

◆ getTofDt()

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

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 56 of file SimplePIDSvc.h.

56{ return m_tof_dt[m][n]; }

◆ getTofPh()

double SimplePIDSvc::getTofPh ( int n )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 57 of file SimplePIDSvc.h.

57{ return m_tof_ph[n]; }

◆ getTofZR()

double SimplePIDSvc::getTofZR ( int n )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 58 of file SimplePIDSvc.h.

58{ return m_tof_zr[n]; }

◆ initialize()

StatusCode SimplePIDSvc::initialize ( )

virtual

Definition at line 31 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();
 
    return sc;
}

◆ iselectron()

bool SimplePIDSvc::iselectron ( bool emc = true )

virtual

Implements ISimplePIDSvc.

Definition at line 1047 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 ( )

virtual

Implements ISimplePIDSvc.

Definition at line 1020 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 ( )

virtual

Implements ISimplePIDSvc.

Definition at line 1012 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 )

virtual

Implements ISimplePIDSvc.

Definition at line 91 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();
    }
    //TOF PID
    loadTOFInfo(track);
    if (m_tof_corr)
    {
        if (m_tof_barrel == 1)
        {
            tofBarrelCorrection();
        }
        else if (m_tof_barrel == 0)
            tofEndcapCorrection();
    }
    //EMC
    loadEMCInfo(track);
 
    calprob();
}

◆ probElectron()

double SimplePIDSvc::probElectron ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 44 of file SimplePIDSvc.h.

44{ return m_prob[0]; }

◆ probKaon()

double SimplePIDSvc::probKaon ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 47 of file SimplePIDSvc.h.

47{ return m_prob[3]; }

◆ probMuon()

double SimplePIDSvc::probMuon ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 45 of file SimplePIDSvc.h.

45{ return m_prob[1]; }

◆ probPion()

double SimplePIDSvc::probPion ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 46 of file SimplePIDSvc.h.

46{ return m_prob[2]; }

◆ probProton()

double SimplePIDSvc::probProton ( )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 48 of file SimplePIDSvc.h.

48{ return m_prob[4]; }

◆ queryInterface()

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

virtual

Definition at line 77 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 36 of file SimplePIDSvc.h.

36{ m_dedx_chi_cut = x; }

x

double x[1000]

Definition: draw_charge_space_coarse_II.cxx:8

◆ seteidratio()

void SimplePIDSvc::seteidratio ( double x )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 38 of file SimplePIDSvc.h.

38{ m_eid_ratio = x; }

◆ settofminchi()

void SimplePIDSvc::settofminchi ( double x )

inlinevirtual

Implements ISimplePIDSvc.

Definition at line 37 of file SimplePIDSvc.h.

37{ 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:

Public Member Functions

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SimplePIDSvc()

◆ ~SimplePIDSvc()

Member Function Documentation

◆ finalize()

◆ getChi2()

◆ getdEdxChi()

◆ getElectronLikelihood()

◆ getEmcE()

◆ getEmcE13()

◆ getEmcE35()

◆ getEmcLatMoment()

◆ getEmcLikelihood()

◆ getEmcSecMoment()

◆ getTofBarrel()

◆ getTOFChi()

◆ getTofCounter()

◆ getTofDt()

◆ getTofPh()

◆ getTofZR()

◆ initialize()

◆ iselectron()

◆ iskaon()

◆ ispion()

◆ preparePID()

◆ probElectron()

◆ probKaon()

◆ probMuon()

◆ probPion()

◆ probProton()

◆ queryInterface()

◆ setdedxminchi()

◆ seteidratio()

◆ settofminchi()

Friends And Related Function Documentation

◆ CnvFactory< SimplePIDSvc >