BFieldCorr Class Reference

#include <BFieldCorr.h>

List of all members.

Public Member Functions
	BFieldCorr (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
Private Attributes
double	m_factor
HepSymMatrix	m_Ew
int	RUN_BEGIN_10
int	RUN_END_10
int	RUN_BEGIN_11
int	RUN_END_11

---

Detailed Description

Definition at line 19 of file BFieldCorr.h.

---

Constructor & Destructor Documentation

BFieldCorr::BFieldCorr	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 45 of file BFieldCorr.cxx.

References m_factor.

                                                                     :Algorithm(name,pSvcLocator)
{
        declareProperty("CorrectionFactor",  m_factor = 1.0000);
}

---

Member Function Documentation

StatusCode BFieldCorr::execute

(

)

Definition at line 68 of file BFieldCorr.cxx.

References DstMdcKalTrack::electron, DstMdcKalTrack::err(), DstMdcTrack::err(), EventModel::EvtRec::EvtRecEvent, EventModel::EvtRec::EvtRecTrackCol, Bes_Common::FATAL, DstMdcKalTrack::helix(), DstMdcTrack::helix(), genRecEmupikp::i, Bes_Common::INFO, ganga-rec::j, DstMdcKalTrack::kaon, m_Ew, m_factor, msgSvc(), DstMdcKalTrack::muon, DstMdcKalTrack::p(), DstMdcTrack::p(), DstMdcKalTrack::pion, DstMdcKalTrack::proton, DstMdcKalTrack::px(), DstMdcTrack::px(), DstMdcKalTrack::pxy(), DstMdcTrack::pxy(), DstMdcKalTrack::py(), DstMdcTrack::py(), DstMdcKalTrack::pz(), DstMdcTrack::pz(), RUN_BEGIN_10, RUN_BEGIN_11, RUN_END_10, RUN_END_11, DstMdcTrack::setError(), DstMdcTrack::setHelix(), DstMdcKalTrack::setP(), DstMdcTrack::setP(), DstMdcKalTrack::setPidType(), DstMdcKalTrack::setPx(), DstMdcTrack::setPx(), DstMdcKalTrack::setPxy(), DstMdcTrack::setPxy(), DstMdcKalTrack::setPy(), DstMdcTrack::setPy(), DstMdcKalTrack::setPz(), DstMdcTrack::setPz(), DstMdcKalTrack::setZError(), and DstMdcKalTrack::setZHelix().

{
        MsgStream log(msgSvc(), name());
        log << MSG::INFO << "in BFieldCorr::execute()" << endreq;
        SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
        int run_no = eventHeader->runNumber();
        if (run_no < 0) return StatusCode::SUCCESS;

        SmartDataPtr<EvtRecEvent> evtRecEvent(eventSvc(), EventModel::EvtRec::EvtRecEvent);
        if (!evtRecEvent)
        {
                log << MSG::FATAL << "Could not find EvtRecEvent" << endreq;
                return StatusCode::SUCCESS;
        }

        SmartDataPtr<EvtRecTrackCol> evtRecTrkCol(eventSvc(),  EventModel::EvtRec::EvtRecTrackCol);
        if (!evtRecTrkCol)
        {
                log << MSG::FATAL << "Could not find EvtRecTrackCol" << endreq;
                return StatusCode::SUCCESS;
        }

        if ((unsigned int)evtRecEvent->totalTracks() > evtRecTrkCol->size()) return StatusCode::SUCCESS;
        //if (evtRecEvent->totalTracks() > 50) return StatusCode::SUCCESS;

        //set Ew if by default
        double factor;
        if (fabs(m_factor - 1.000) < 1e-6)
        {
                if (run_no >= RUN_BEGIN_10 && run_no <= RUN_END_10)
                        factor = 1.0004;
                else if (run_no >= RUN_BEGIN_11 && run_no <= RUN_END_11)
                        factor = 1.0002;
                else
                        factor = 1.0000;
        }       
        else
        {
                factor = m_factor;
        }
        m_Ew[2][2] = fabs(factor) < 1e-6 ? 1.0 : 1.0 / factor; 

        for (int i = 0; i < evtRecEvent->totalCharged(); i++)
        {
                EvtRecTrackIterator itTrk = evtRecTrkCol->begin() + i;
                if (!(*itTrk)->isMdcTrackValid()) continue;
                if (!(*itTrk)->isMdcKalTrackValid()) continue;

                //RecMdcTrack correction 
                RecMdcTrack *mdc_trk = (*itTrk)->mdcTrack();
                mdc_trk->setHelix( m_Ew * mdc_trk->helix() );
                mdc_trk->setError( mdc_trk->err().similarity(m_Ew) );
                mdc_trk->setP(   factor * mdc_trk->p()   );
                mdc_trk->setPxy( factor * mdc_trk->pxy() );
                mdc_trk->setPx(  factor * mdc_trk->px()  );
                mdc_trk->setPy(  factor * mdc_trk->py()  );
                mdc_trk->setPz(  factor * mdc_trk->pz()  );

                //RecMdcKalTrack correction 
                RecMdcKalTrack *kal_trk = (*itTrk)->mdcKalTrack();
                RecMdcKalTrack::PidType trk_type[] = {
                        RecMdcKalTrack::electron,
                        RecMdcKalTrack::muon,
                        RecMdcKalTrack::pion,
                        RecMdcKalTrack::kaon,
                        RecMdcKalTrack::proton,
                };
                
                for (int j = 0; j < 5; j++)
                {
                        kal_trk->setPidType( trk_type[j] );
                        kal_trk->setZHelix( m_Ew * kal_trk->helix(), j );
                        kal_trk->setZError( kal_trk->err().similarity(m_Ew), j );
                        kal_trk->setP(   factor * kal_trk->p(),   j );
                        kal_trk->setPxy( factor * kal_trk->pxy(), j );
                        kal_trk->setPx(  factor * kal_trk->px(),  j );
                        kal_trk->setPy(  factor * kal_trk->py(),  j );
                        kal_trk->setPz(  factor * kal_trk->pz(),  j );
                }
        }
        return StatusCode::SUCCESS;
}

StatusCode BFieldCorr::finalize

(

)

Definition at line 151 of file BFieldCorr.cxx.

References Bes_Common::INFO, and msgSvc().

{
        MsgStream log(msgSvc(), name());
        log << MSG::INFO << "in BFieldCorr::finalize()" << endreq;
        return StatusCode::SUCCESS;
}

StatusCode BFieldCorr::initialize

(

)

Definition at line 50 of file BFieldCorr.cxx.

References Bes_Common::INFO, m_Ew, m_factor, msgSvc(), RUN_BEGIN_10, RUN_BEGIN_11, RUN_END_10, and RUN_END_11.

{
        m_Ew = HepSymMatrix(5, 0);
        m_Ew[0][0] = 1.0; 
        m_Ew[1][1] = 1.0; 
        m_Ew[2][2] = fabs(m_factor) < 1e-6 ? 1.0 : 1.0 / m_factor; 
        m_Ew[3][3] = 1.0; 
        m_Ew[4][4] = 1.0; 
        RUN_BEGIN_10 = 11414;
        RUN_END_10   = 14604;
        RUN_BEGIN_11 = 20448;
        RUN_END_11   = 23454;

        MsgStream log(msgSvc(), name());
        log << MSG::INFO << "in MagCorr::initialize()" <<endmsg;
        return StatusCode::SUCCESS;
}

---

Member Data Documentation

HepSymMatrix BFieldCorr::m_Ew [private]

Definition at line 30 of file BFieldCorr.h.

Referenced by execute(), and initialize().

double BFieldCorr::m_factor [private]

Definition at line 28 of file BFieldCorr.h.

Referenced by BFieldCorr(), execute(), and initialize().

int BFieldCorr::RUN_BEGIN_10 [private]

Definition at line 31 of file BFieldCorr.h.

Referenced by execute(), and initialize().

int BFieldCorr::RUN_BEGIN_11 [private]

Definition at line 33 of file BFieldCorr.h.

Referenced by execute(), and initialize().

int BFieldCorr::RUN_END_10 [private]

Definition at line 32 of file BFieldCorr.h.

Referenced by execute(), and initialize().

int BFieldCorr::RUN_END_11 [private]

Definition at line 34 of file BFieldCorr.h.

Referenced by execute(), and initialize().

---