MdcDummy.cxx

Go to the documentation of this file.

#include "MdcDummyAlg/MdcDummy.h"
#include "EmcRawEvent/EmcDigi.h"
#include "EventModel/EventHeader.h"
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/PropertyMgr.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "McTruth/DecayMode.h"
#include "MdcRawEvent/MdcDigi.h"
#include "TofRawEvent/TofDigi.h"
 
/*#include "McTruth/McParticle.h"
 
#include "MdcGeomSvc/IMdcGeomSvc.h"
#include "MdcGeomSvc/MdcGeoLayer.h"
#include "MdcGeomSvc/MdcGeoWire.h"
*/
#include <iostream>
#include <vector>
 
using namespace std;
using namespace Event;
 
/////////////////////////////////////////////////////////////////////////////
DECLARE_COMPONENT( MdcDummy )
MdcDummy::MdcDummy( const std::string& name, ISvcLocator* pSvcLocator )
    : Algorithm( name, pSvcLocator ) {
  // Declare the properties
  declareProperty( "FittingMethod", m_fittingMethod = 2 );
  declareProperty( "ConfigFile", m_configFile = "MDCConfig.xml" );
}
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode MdcDummy::initialize() {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in initialize()" << endmsg;
 
  return StatusCode::SUCCESS;
}
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode MdcDummy::execute() {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in execute()" << endmsg;
 
  // Part 1: Get the event header, print out event and run number
  int event, run;
 
  SmartDataPtr<Event::EventHeader> eventHeader( eventSvc(), "/Event/EventHeader" );
  if ( !eventHeader )
  {
    log << MSG::FATAL << "Could not find Event Header" << endmsg;
    return ( StatusCode::FAILURE );
  }
  log << MSG::INFO << "MdcDummy: retrieved event: " << eventHeader->eventNumber()
      << "  run: " << eventHeader->runNumber() << endmsg;
 
  int digiId;
 
  // Part 2: Retrieve MC truth
  SmartDataPtr<DecayMode> decayMode( eventSvc(), "/Event/MC/DecayMode" );
  if ( !decayMode )
  {
    log << MSG::FATAL << "Could not find DecayMode" << endmsg;
    return ( StatusCode::FAILURE );
  }
  unsigned int size = 10;
  int          decayId[10];
  decayMode->getData( decayId, size );
 
  for ( unsigned int index = 0; index < 10; index++ )
  { log << MSG::INFO << " Decay Mode  = " << decayId[index] << endmsg; }
 
  // Part 3: Retrieve MDC digi
  SmartDataPtr<MdcDigiCol> mdcDigiCol( eventSvc(), "/Event/Digi/MdcDigiCol" );
  if ( !mdcDigiCol )
  {
    log << MSG::FATAL << "Could not find MDC digi" << endmsg;
    return ( StatusCode::FAILURE );
  }
 
  MdcDigiCol::iterator iter1 = mdcDigiCol->begin();
  digiId                     = 0;
  for ( ; iter1 != mdcDigiCol->end(); iter1++, digiId++ )
  {
    log << MSG::INFO << "MDC digit No: " << digiId << endmsg;
 
    log << MSG::INFO << " time_channel = " << ( *iter1 )->getTimeChannel()
        << " charge_channel = " << ( *iter1 )->getChargeChannel() << endmsg;
  }
 
  // Part 4: Retrieve TOF digi
  SmartDataPtr<TofDigiCol> tofDigiCol( eventSvc(), "/Event/Digi/TofDigiCol" );
  if ( !tofDigiCol )
  {
    log << MSG::FATAL << "Could not find TOF digi" << endmsg;
    return ( StatusCode::FAILURE );
  }
 
  TofDigiCol::iterator iter2 = tofDigiCol->begin();
  digiId                     = 0;
  for ( ; iter2 != tofDigiCol->end(); iter2++, digiId++ )
  { log << MSG::INFO << "TOF digit No: " << digiId << endmsg; }
 
  // Part 5: Retrieve EMC digi
  SmartDataPtr<EmcDigiCol> emcDigiCol( eventSvc(), "/Event/Digi/EmcDigiCol" );
  if ( !emcDigiCol )
  {
    log << MSG::FATAL << "Could not find EMC digi" << endmsg;
    return ( StatusCode::FAILURE );
  }
 
  EmcDigiCol::iterator iter3 = emcDigiCol->begin();
  digiId                     = 0;
  for ( ; iter3 != emcDigiCol->end(); iter3++, digiId++ )
  {
    log << MSG::INFO << "Emc digit No: " << digiId << endmsg;
 
    log << MSG::INFO << " time_channel = " << ( *iter3 )->getTimeChannel()
        << " charge_channel = " << ( *iter3 )->getChargeChannel() << endmsg;
  }
 
  return StatusCode::SUCCESS;
}
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode MdcDummy::finalize() {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in finalize()" << endmsg;
 
  return StatusCode::SUCCESS;
}