BFieldCorr Class Reference

#include <BFieldCorr.h>

Inheritance diagram for BFieldCorr:

Public Member Functions
	BFieldCorr (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()

Detailed Description

Definition at line 19 of file BFieldCorr.h.

Constructor & Destructor Documentation

BFieldCorr()

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

Definition at line 45 of file BFieldCorr.cxx.

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

Member Function Documentation

execute()

StatusCode BFieldCorr::execute

(

)

Definition at line 68 of file BFieldCorr.cxx.

{
    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;
}

finalize()

StatusCode BFieldCorr::finalize

(

)

Definition at line 151 of file BFieldCorr.cxx.

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

initialize()

StatusCode BFieldCorr::initialize

(

)

Definition at line 50 of file BFieldCorr.cxx.

{
    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;
}

---

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

• BFieldCorr.h
• BFieldCorr.cxx