LumTau Class Reference

#include <LumTau.h>

Inheritance diagram for LumTau:

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

Detailed Description

Definition at line 10 of file LumTau.h.

Constructor & Destructor Documentation

LumTau()

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

Definition at line 55 of file LumTau.cxx.

                                                              :
    Algorithm(name, pSvcLocator) 
{  
}

Member Function Documentation

execute()

StatusCode LumTau::execute

(

)

Definition at line 105 of file LumTau.cxx.

{
    StatusCode sc=StatusCode::SUCCESS;
 
    MsgStream log(msgSvc(),name());
    log<<MSG::INFO<<"in execute()"<<endreq;
 
 
    SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
    SmartDataPtr<EvtRecEvent> evtRecEvent(eventSvc(), EventModel::EvtRec::EvtRecEvent);
 
    log<<MSG::DEBUG<<"ncharg, nneu, tottks = "<<evtRecEvent->totalCharged()<<" , "<<evtRecEvent->totalNeutral()<<" , "<<evtRecEvent->totalTracks()<<endreq;
    SmartDataPtr<EvtRecTrackCol> evtRecTrkCol(eventSvc(),EventModel::EvtRec::EvtRecTrackCol);
 
    log<<MSG::DEBUG <<"ncharg, nneu, tottks = "<<evtRecEvent->totalCharged()<<" , "<<evtRecEvent->totalNeutral()<<" , "<<evtRecEvent->totalTracks()<<endreq;
    SmartDataPtr<RecEsTimeCol> evTimeCol(eventSvc(),"/Event/Recon/RecEsTimeCol");
 
        int interval = 200;
    interval = m_intervalSvc->getTInterval();
    log<<MSG::DEBUG << "interval = " << interval << endreq;
    m_topup = false;
    if (interval <= 0 ) { 
       log << MSG::FATAL << "error in reading the interval!" << endreq;
       return StatusCode::FAILURE;
    }
    else if ( interval < 100 ) {  // top-up mode
        m_topup = true;
    }
 
    double time = eventHeader->time();
    log << MSG::DEBUG << "time = " << time << endreq;
    if (time <=0 ) return StatusCode::SUCCESS;
    m_time = time;
    m_etsT1 = eventHeader->etsT1()/2000000.; // convert to second
    log << MSG::DEBUG << "ets time = " << m_etsT1 << endreq;
 
    m_run = eventHeader->runNumber();
    m_rec = eventHeader->eventNumber();
 
    if(m_rec%1000==0) log << MSG::INFO <<"Run   "<<m_run<<"     Event   "<<m_rec<<endreq;
 
    double Emax1 = -1;
    double Emax2 = -1;
    int Imax[2];
 
    if((evtRecEvent->totalTracks() >= 2)){
        double etot = 0.;
        for(int i = 0;i < evtRecEvent->totalTracks(); i++)
        {
            if(i>=evtRecTrkCol->size()) break;
            EvtRecTrackIterator  itTrk=evtRecTrkCol->begin() + i;
            if(!(*itTrk)->isEmcShowerValid()) continue;
            RecEmcShower *emcTrk = (*itTrk)->emcShower();
            double Ener=emcTrk->energy();
            if(Ener>Emax2)
            {
                Emax2=Ener;
                Imax[1]=i;
            }
            if(Ener>Emax1)
            {
                Emax2=Emax1;
                Imax[1]=Imax[0];
                Emax1=Ener;
                Imax[0]=i;
            }
            etot += Ener;
        }
 
        m_etot = etot;
        m_e1 = Emax1;
        m_e2 = Emax2;
 
        log << MSG::INFO << "Emax1 = " << Emax1 <<"Emax2= "<<Emax2<< endreq;
        if(Emax1 > 0 && Emax2 > 0){
            double emcphi[2],emctht[2];
            for(int i = 0;i < 2; i++)
            {
                if(i>=evtRecTrkCol->size()) break;
                EvtRecTrackIterator  itTrk=evtRecTrkCol->begin() + Imax[i];
                if(!(*itTrk)->isEmcShowerValid()) continue;
                RecEmcShower *emcTrk = (*itTrk)->emcShower();
                emcphi[i]=emcTrk->phi();
                emctht[i]=emcTrk->theta();
            }
 
            double dltphi=(fabs(emcphi[0]-emcphi[1])-pai)*180./pai;
            double dlttht=(fabs(emctht[0]+emctht[1])-pai)*180./pai;
            m_costht1=cos(emctht[0]);
            m_costht2=cos(emctht[1]);
            m_phi1=emcphi[0];
            m_phi2=emcphi[1];
            m_dlttht=dlttht;
            m_dltphi=dltphi;
        }
        else{
            m_etot = -1;
            m_e1 = -1;
            m_e2 = -1;
            m_costht1 = -1;
            m_costht2 = -1;
            m_phi1 = -1;
            m_phi2 = -1;
            m_dlttht = -1;
            m_dltphi = -1;
        }
    }
    else{
        m_etot = -1;
        m_e1 = -1;
        m_e2 = -1;
        m_costht1 = -1;
        m_costht2 = -1;
        m_phi1 = -1;
        m_phi2 = -1;
        m_dlttht = -1;
        m_dltphi = -1;
    }
 
//  m_tuple2->write();
 
    int DiskWrite = m_tuple2->write();
    if(DiskWrite != 1){
        log<<MSG::FATAL<<"ERROR In LumTau DiskWrite!"<<endreq;
        exit(1);
    }
 
    return StatusCode::SUCCESS;
}

finalize()

StatusCode LumTau::finalize

(

)

Definition at line 235 of file LumTau.cxx.

{
    MsgStream log(msgSvc(),name());
    log << MSG::INFO <<"Event Finalize"<<endreq;
    return StatusCode::SUCCESS;
}

hello()

StatusCode LumTau::hello

(

)

initialize()

StatusCode LumTau::initialize

(

)

Definition at line 61 of file LumTau.cxx.

{
    MsgStream log(msgSvc(), name());
 
    log << MSG::INFO << "in initialize()" << endmsg;
 
    StatusCode status;
 
    status = service("InjSigIntervalSvc", m_intervalSvc);
    if ( status != StatusCode::SUCCESS ){
             log << MSG::FATAL << "can not use InjSigIntervalSvc" << endreq;
        }
 
    NTuplePtr nt2(ntupleSvc(), "LumTau/event");
    if ( nt2 ) m_tuple2 = nt2;
    else 
    {
        m_tuple2 = ntupleSvc()->book ("LumTau/event", CLID_ColumnWiseTuple, "Bhabha N-Tuple signal");
        if ( m_tuple2 )
        {
            status = m_tuple2->addItem ("topup",    m_topup);
            status = m_tuple2->addItem ("run",    m_run);
            status = m_tuple2->addItem ("rec",    m_rec);
            status = m_tuple2->addItem ("time",  m_time);  
            status = m_tuple2->addItem ("etsT1",  m_etsT1);  
            status = m_tuple2->addItem ("etot",  m_etot);
            status = m_tuple2->addItem ("e1",    m_e1);
            status = m_tuple2->addItem ("e2",    m_e2);
            status = m_tuple2->addItem ("costht1",    m_costht1);
            status = m_tuple2->addItem ("costht2",    m_costht2);
            status = m_tuple2->addItem ("dltphi",    m_dltphi);
            status = m_tuple2->addItem ("dlttht",    m_dlttht);
            status = m_tuple2->addItem ("phi1",    m_phi1);
            status = m_tuple2->addItem ("phi2",    m_phi2);
        }
        else
        {
            log << MSG::ERROR << "Cannot book N-tuple2:"<<long(m_tuple2)<<endmsg;
            return StatusCode::FAILURE;
        }
    }
    return StatusCode::SUCCESS;
}

---

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

• LumTau.h
• LumTau.cxx