GenRdmTrgIdxAlg.cxx

Go to the documentation of this file.

#include "GenRdmTrgIdxAlg.h"
#include "EventModel/EventHeader.h"
#include "EventModel/EventModel.h"
#include "GaudiKernel/SmartDataPtr.h"
 
// #include "GaudiKernel/AlgFactory.h"
#include "TrigEvent/TrigData.h"
 
DECLARE_COMPONENT( GenRdmTrgIdxAlg )
 
GenRdmTrgIdxAlg::GenRdmTrgIdxAlg( const std::string& name, ISvcLocator* pSvcLocator )
    : Algorithm( name, pSvcLocator ), m_fcount( 0 ) {
  declareProperty( "RandomTriggerFiles", m_rdmFiles );
}
 
GenRdmTrgIdxAlg::~GenRdmTrgIdxAlg() {}
 
StatusCode GenRdmTrgIdxAlg::initialize() {
  // StatusCode sc = GaudiAlgorithm::initialize(); // must be executed first
  // StatusCode sc;
 
  if ( m_rdmFiles.empty() )
  {
    error() << "RandomTriggerFiles is empty" << endmsg;
    return StatusCode::FAILURE;
  }
 
  m_idxhandler = new EvtIdxHandler;
  m_freader    = new SimpleRawReader( m_rdmFiles[0] );
 
  return StatusCode::SUCCESS;
}
 
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 StatusCode::SUCCESS;
}