OfflineEvtFilterAlg Class Reference

#include <OfflineEvtFilterAlg.h>

Inheritance diagram for OfflineEvtFilterAlg:

Public Member Functions
	OfflineEvtFilterAlg (const std::string &name, ISvcLocator *pSvcLocator)
	Standard constructor.
virtual	~OfflineEvtFilterAlg ()
	Destructor.
virtual StatusCode	initialize ()
	Algorithm initialization.
virtual StatusCode	execute ()
	Algorithm execution.
virtual StatusCode	finalize ()
	Algorithm finalization.

Detailed Description

Definition at line 9 of file OfflineEvtFilterAlg.h.

Constructor & Destructor Documentation

OfflineEvtFilterAlg()

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

Standard constructor.

Definition at line 24 of file OfflineEvtFilterAlg.cxx.

: GaudiAlgorithm ( name , pSvcLocator )
{
}

Referenced by OfflineEvtFilterAlg().

~OfflineEvtFilterAlg()

OfflineEvtFilterAlg::~OfflineEvtFilterAlg ( )

virtual

Destructor.

Definition at line 32 of file OfflineEvtFilterAlg.cxx.

{} 

Member Function Documentation

execute()

StatusCode OfflineEvtFilterAlg::execute ( )

virtual

Algorithm execution.

filter with ets-ist

filter with ets

Definition at line 58 of file OfflineEvtFilterAlg.cxx.

                                        {
 
    SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
 
    // whether to keep this event for further reconstruction
    bool keep = true;
    int npar = m_evtFilterSvc->getNpar();
 
    int iPar = 0;
    /// filter with ets-ist
    if ( eventHeader->etsFlag() != 21 ) {
        double dtEts = (long(eventHeader->etsT1()) - long(eventHeader->etsT2())) / double(2000.);  //in ms
        for ( ; iPar < npar && m_evtFilterSvc->getFlag(iPar) == 0; ++iPar ) {
            if ( dtEts < m_evtFilterSvc->getTBegin(iPar) ) {
                //keep = true by default
                break;
            }
            else if ( dtEts < m_evtFilterSvc->getTEnd(iPar) ) {
                keep = false;
                break;
            }
        }
    }
    /// filter with ets
    if ( keep && iPar < npar ) {
        double tEts  = eventHeader->etsT1() / double(2000000.);  //in second
        for ( ; iPar < npar; ++iPar ) {
            if ( m_evtFilterSvc->getFlag(iPar) == 0 ) {
                continue;
            }
            if ( tEts < m_evtFilterSvc->getTBegin(iPar) ) {
                //keep = true by default
                break;
            }
            else if ( tEts < m_evtFilterSvc->getTEnd(iPar) ) {
                keep = false;
                break;
            }
        }
    }
 
    if ( keep ) {
        // this is a good event, nothing to do
        return StatusCode::SUCCESS;
    }
 
    // Otherwise, this is a bad event ...
    // We register empty data objects without any hits to the Gaudi EDS
 
    {
        DigiEvent *digi = new DigiEvent();
        digi->initialize(false);  //not from MC
        m_dataSvc->registerObject(EventModel::Digi::Event, digi);
    }
 
    {
        DstEvent* dst = new DstEvent();
        m_dataSvc->registerObject(EventModel::Dst::Event, dst);
    }
 
    {
        TrigEvent * trig = new TrigEvent();
        trig->initialize(false);  //not from MC
        m_dataSvc->registerObject(EventModel::Trig::Event, trig);
        TrigData* trigData = new TrigData();
        m_dataSvc->registerObject(EventModel::Trig::TrigData, trigData);
    }
 
    {
        EvtRecObject* evtRecObject = new EvtRecObject();
        m_dataSvc->registerObject(EventModel::EvtRec::Event, evtRecObject);
    }
 
    {
        HltEvent * hlt = new HltEvent();
        hlt->initialize(false);  //not from MC
        m_dataSvc->registerObject(EventModel::Hlt::Event, hlt);
        DstHltInf* dstHltInf = new DstHltInf();
        m_dataSvc->registerObject(EventModel::Hlt::DstHltInf, dstHltInf);
    }
 
    {
        ReconEvent * recon = new ReconEvent();
        m_dataSvc->registerObject(EventModel::Recon::Event, recon);
    }
 
    debug() << "event " << eventHeader->eventNumber() << " is passed" << endmsg;
    m_incidentSvc->fireIncident(Incident(name(), IncidentType::AbortEvent));
 
    return StatusCode::SUCCESS;
}

finalize()

StatusCode OfflineEvtFilterAlg::finalize ( )

virtual

Algorithm finalization.

Definition at line 153 of file OfflineEvtFilterAlg.cxx.

                                         {
 
    debug() << "==> Finalize" << endmsg;
 
    return GaudiAlgorithm::finalize();  // must be called after all other actions
}

initialize()

StatusCode OfflineEvtFilterAlg::initialize ( )

virtual

Algorithm initialization.

Definition at line 37 of file OfflineEvtFilterAlg.cxx.

                                           {
    StatusCode sc = GaudiAlgorithm::initialize(); // must be executed first
    if ( sc.isFailure() ) return sc;  // error printed already by GaudiAlgorithm
 
    debug() << "==> Initialize" << endmsg;
 
    sc = service("OfflineEvtFilterSvc", m_evtFilterSvc);
    if ( sc != StatusCode::SUCCESS ) {
        fatal() << "can not use OfflineEvtFilterSvc" << endreq;
        return sc;
    }
 
    m_dataSvc = eventSvc();
    m_incidentSvc = svc<IIncidentSvc>("IncidentSvc",true);
 
    return StatusCode::SUCCESS;
}

---

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

• OfflineEvtFilterAlg.h
• OfflineEvtFilterAlg.cxx