GenRdmTrgIdxAlg.cxx

Go to the documentation of this file.

#include "GenRdmTrgIdxAlg.h"
#include "EventModel/EventModel.h"
#include "EventModel/EventHeader.h"
#include "TrigEvent/TrigData.h"
#include "GaudiKernel/AlgFactory.h"
 
DECLARE_ALGORITHM_FACTORY( GenRdmTrgIdxAlg )
 
GenRdmTrgIdxAlg::GenRdmTrgIdxAlg(const std::string& name, ISvcLocator* pSvcLocator)
    : GaudiAlgorithm(name, pSvcLocator),
      m_fcount(0)
{
    declareProperty("RandomTriggerFiles",  m_rdmFiles);
}
 
GenRdmTrgIdxAlg::~GenRdmTrgIdxAlg()
{
}
 
StatusCode GenRdmTrgIdxAlg::initialize()
{
    StatusCode sc = GaudiAlgorithm::initialize(); // must be executed first
 
    if ( m_rdmFiles.empty() ) {
        error() << "RandomTriggerFiles is empty" << endmsg;
        return StatusCode::FAILURE;
    }
 
    m_idxhandler = new EvtIdxHandler;
    m_freader = new SimpleRawReader(m_rdmFiles[0]);
 
    return sc;
}
 
StatusCode GenRdmTrgIdxAlg::execute()
{
    SmartDataPtr<TrigData> trigData(eventSvc(), EventModel::Trig::TrigData);
    if ( trigData->getTrigChannel(9) != 1 ) {
        // this is not a random trigger event, pass it
        return StatusCode::SUCCESS;
    }
 
    SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(), "/Event/EventHeader");
    uint32_t evtNo = eventHeader->eventNumber();
    uint32_t runNo = eventHeader->runNumber();
 
    uint32_t thePos = 0;
    const uint32_t* data = NULL;
    uint32_t rdmEvtNo = 0xFFFFFFFF;
 
    //std::cout << "BEGIN " << evtNo << std::endl;
    do {
        thePos = m_freader->raw_stream()->tellg();
        try {
            data = m_freader->nextEvent();
        }
        catch ( ReachEndOfFile& e ) {
            m_idxhandler->write( RawFileTools::fname2idxname(m_freader->raw_stream()->currentFile()) );
            delete m_freader;
            m_freader = NULL;
            ++m_fcount;
            if ( m_fcount < m_rdmFiles.size() ) {
                m_freader = new SimpleRawReader(m_rdmFiles[m_fcount]);
                rdmEvtNo = 0xFFFFFFFF;
                continue;
            }
            else {
                fatal() << "can not find event " << evtNo << " in random trigger data files" << endmsg;
                return StatusCode::FAILURE;
            }
        }
        catch ( RawFileException& e ) {
            e.print();
            throw e;
        }
 
        //std::cout << "run " << data[9+data[5]] << ",  event " << data[8+data[5]] << std::endl;
 
        if ( data[9+data[5]] != runNo ) {
            fatal() << "bad runNo in random trigger data" << endmsg;
            return StatusCode::FAILURE;
        }
        rdmEvtNo = data[8+data[5]];
    } while ( rdmEvtNo != evtNo );
 
    m_idxhandler->addPos(data[8+data[5]], thePos);
 
    return StatusCode::SUCCESS;
}
 
StatusCode GenRdmTrgIdxAlg::finalize()
{
    if ( m_freader != NULL ) {
        m_idxhandler->write( RawFileTools::fname2idxname(m_freader->raw_stream()->currentFile()) );
        delete m_freader;
    }
    delete m_idxhandler;
 
    return GaudiAlgorithm::finalize();  // must be called after all other actions
}