ReadRawData.cxx

Go to the documentation of this file.

// Include files.
#include "ReadRawData.h"
#include "GaudiKernel/IDataManagerSvc.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/PropertyMgr.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "RawDataCnv/EventManagement/RAWEVENT.h"
#include <vector>
 
// Event Model related classes
#include "EmcRawEvent/EmcDigi.h"
#include "EventModel/Event.h"
#include "MdcRawEvent/MdcDigi.h"
#include "MucRawEvent/MucDigi.h"
#include "RawEvent/DigiEvent.h"
#include "TofRawEvent/TofDigi.h"
 
#include "EventModel/EventModel.h"
 
using namespace std;
 
// Constructor.
ReadRawData::ReadRawData( const string& name, ISvcLocator* pSvcLocator )
    : Algorithm( name, pSvcLocator ) {}
 
// ReadRawData::initialize().
StatusCode ReadRawData::initialize() {
 
  MsgStream log( msgSvc(), name() );
  m_count = 0;
 
  log << MSG::INFO << "in initialize()" << endmsg;
 
  return StatusCode::SUCCESS;
}
 
// ReadRawData::execute().
StatusCode ReadRawData::execute() {
  StatusCode sc;
  MsgStream  log( msgSvc(), name() );
 
  log << MSG::INFO << "================================================" << endmsg;
  log << MSG::INFO << "In execute()" << endmsg;
  log << MSG::INFO << "++++++++++++++++++++++++++++++++++++++++++++++++" << endmsg;
 
  SmartDataPtr<Event::EventHeader> evt( eventSvc(), EventModel::EventHeader );
  if ( !evt )
  {
    log << MSG::ERROR << "Did not retrieve event" << endmsg;
    return StatusCode::FAILURE;
  }
  // Check event header
 
  int eventNo = evt->eventNumber();
  int runNo   = evt->runNumber();
 
  log << MSG::INFO << "ReadRawData: retrieved event: " << eventNo << "  run: " << runNo
      << endmsg;
 
  log << MSG::INFO << "ReadRawData: start to read MdcDigiData" << endmsg;
  sc = readMdcDigiData();
 
  log << MSG::INFO << "ReadRawData: start to read EmcDigiData" << endmsg;
  sc = readEmcDigiData();
 
  log << MSG::INFO << "ReadRawData: start to read TofDigiData" << endmsg;
  sc = readTofDigiData();
 
  log << MSG::INFO << "ReadRawData: start to read MucDigiData" << endmsg;
  sc = readMucDigiData();
  if ( sc.isFailure() ) { return sc; }
 
  m_count++;
 
  return StatusCode::SUCCESS;
}
 
StatusCode ReadRawData::readMdcDigiData() {
  MsgStream  log( msgSvc(), name() );
  StatusCode sc = StatusCode::SUCCESS;
 
  SmartDataPtr<DigiEvent> digiEvt( eventSvc(), EventModel::Digi::Event );
  if ( !digiEvt )
  {
    log << MSG::ERROR << "Did not retrieve DigiEvent" << endmsg;
    return StatusCode::FAILURE;
  }
 
  SmartDataPtr<MdcDigiCol> digiCol( eventSvc(), EventModel::Digi::MdcDigiCol );
 
  if ( digiCol == 0 )
  {
    log << "no MdcDigiCol found" << endmsg;
    return sc;
  }
  else
  {
    log << "Event No. " << m_count << " --->   " << digiCol->size() << " Mdc digis found "
        << endmsg;
 
    log << MSG::INFO << "Detailed dump of this event: " << endmsg;
    int                        ndigi    = 0;
    MdcDigiCol::const_iterator pMdcDigi = digiCol->begin();
    for ( ; pMdcDigi != digiCol->end(); pMdcDigi++ )
    {
      log << MSG::INFO << "Digi " << ndigi++ << " ";
      // output the digi proper:  " << digi " no longer works because
      // "digi" now returns an integer sort order
      ( **pMdcDigi ).fillStream( log.stream() );
      log << endmsg;
    }
  }
 
  return StatusCode::SUCCESS;
}
 
StatusCode ReadRawData::readEmcDigiData() {
  MsgStream  log( msgSvc(), name() );
  StatusCode sc = StatusCode::SUCCESS;
 
  SmartDataPtr<EmcDigiCol> digiCol( eventSvc(), EventModel::Digi::EmcDigiCol );
 
  if ( digiCol == 0 )
  {
    log << "no EmcDigiCol found" << endmsg;
    return sc;
  }
  else
  {
    log << "Event No. " << m_count << " --->  " << digiCol->size() << " Emc digis found "
        << endmsg;
    log << MSG::INFO << "Detailed dump of this event: " << endmsg;
    int                        ndigi    = 0;
    EmcDigiCol::const_iterator pEmcDigi = digiCol->begin();
    for ( ; pEmcDigi != digiCol->end(); pEmcDigi++ )
    {
      log << MSG::INFO << "Digi " << ndigi++ << " ";
      ( **pEmcDigi ).fillStream( log.stream() );
      log << endmsg;
    }
  }
 
  return StatusCode::SUCCESS;
}
 
StatusCode ReadRawData::readTofDigiData() {
  MsgStream  log( msgSvc(), name() );
  StatusCode sc = StatusCode::SUCCESS;
 
  SmartDataPtr<TofDigiCol> digiCol( eventSvc(), EventModel::Digi::TofDigiCol );
 
  if ( digiCol == 0 )
  {
    log << "no TofDigiCol found" << endmsg;
    return sc;
  }
  else
  {
    log << "Event No. " << m_count << " --->  " << digiCol->size() << " Tof digis found "
        << endmsg;
    log << MSG::INFO << "Detailed dump of this event: " << endmsg;
    int                        ndigi    = 0;
    TofDigiCol::const_iterator pTofDigi = digiCol->begin();
    for ( ; pTofDigi != digiCol->end(); pTofDigi++ )
    {
      log << MSG::INFO << "Digi " << ndigi++ << " ";
      ( **pTofDigi ).fillStream( log.stream() );
      log << endmsg;
    }
  }
 
  return StatusCode::SUCCESS;
}
 
StatusCode ReadRawData::readMucDigiData() {
  MsgStream  log( msgSvc(), name() );
  StatusCode sc = StatusCode::SUCCESS;
 
  SmartDataPtr<MucDigiCol> digiCol( eventSvc(), EventModel::Digi::MucDigiCol );
 
  if ( digiCol == 0 )
  {
    log << "no MucDigiCol found" << endmsg;
    return sc;
  }
  else
  {
    log << "Event No. " << m_count << " --->  " << digiCol->size() << " Muc digis found "
        << endmsg;
    log << MSG::INFO << "Detailed dump of this event: " << endmsg;
    int                        ndigi    = 0;
    MucDigiCol::const_iterator pMucDigi = digiCol->begin();
    for ( ; pMucDigi != digiCol->end(); pMucDigi++ )
    {
      log << MSG::INFO << "Digi " << ndigi++ << " ";
      ( **pMucDigi ).fillStream( log.stream() );
      log << endmsg;
    }
  }
 
  return StatusCode::SUCCESS;
}
 
// ReadRawData::finalize().
StatusCode ReadRawData::finalize() {
 
  MsgStream log( msgSvc(), name() );
 
  log << MSG::INFO << "In finalize()" << endmsg;
 
  return StatusCode::SUCCESS;
}