PreElecAlgorithm.cc

Go to the documentation of this file.

#include "PreElecSimSvc/PreElecAlgorithm.h"

#include "Event/GenHeader.h"
#include "Event/SimHeader.h"
#include "Event/RegistrationSequence.h"

#include "RootIOSvc/IRootIOSvc.h"
#include "RootIOSvc/RootIOAddress.h"
#include "PerGenEvent/PerGenHeader.h"
#include "PerSimEvent/PerSimHeader.h"

#include "GaudiKernel/IDataManagerSvc.h"

using namespace std;
using namespace DayaBay;

PreElecAlgorithm::PreElecAlgorithm(const string& name, ISvcLocator* pSvcLocator) 
                 :DybAlgorithm<SimHeader>(name,pSvcLocator)
{
  info()<< "Construct:PreElecAlgorithm()!"<< endreq;
}

StatusCode PreElecAlgorithm::initialize()
{
  StatusCode status= this->GaudiAlgorithm::initialize();
  if(status.isFailure()) 
  { 
    error()<< "Error when initializing  GaudiAlgorithm"<< endreq;
    return status;
  }
   
  info()<< "PreElecAlgorithm ininalize()"<< endreq;

  //Initialize the PreElecSimSvc
  IService* isvc = 0;
  StatusCode sc = serviceLocator()->service("PreElecSimSvc", isvc, false);
  if(sc.isFailure()) 
  {
    error()<< "Service : PreElecSimSvc could not be retrieved"<< endreq;
    return sc;
  }
  isvc->addRef();
  sc = isvc->queryInterface(IPreElecSimSvc::interfaceID(), (void**)&m_preElecSimSvc);
  if(sc.isFailure()) 
  {
    error()<< "Service PreElecSvc does not implement IPreElecSimSvc"<< endreq;
  }
  return sc;
}

StatusCode PreElecAlgorithm::execute() 
{
  info()<< "Now into PreElecAlgorithm executation!"<< endreq;

  SimHitHeader* shh = 0;
  vector<const SimHeader*> shs;
  shh = m_preElecSimSvc->getSimHitHeader(shs, m_hitEarliest, m_hitLatest);

  if(!shh)
  {
    error()<< "Error when get one SimHitHeader!"<< endreq;
    return StatusCode::FAILURE;
  }
  info()<< "SimHitHeader's address: "<< shh<< endreq;

  map<short int, SimHitCollection*> shc = shh->hitCollection();
  for(map<short int, SimHitCollection*>::iterator mit = shc.begin(); mit != shc.end(); mit++)
  {
    vector<SimHit*> vsh = mit->second->collection();
    info()<< MSG::INFO<< "In Det"<< mit->first
          << " : SimHit's num = "<< vsh.size()<< endreq;
  }

  //delete shh;
  for(vector<const SimHeader*>::size_type i = 0; i < shs.size(); i++)
  {
    delete shs[i];
  }

  info()<< "PreElecAlgorithm execute() finalized!" << endreq;
  return StatusCode::SUCCESS;
}

StatusCode PreElecAlgorithm::finalize() 
{
  info()<< "PreElecAlgorithm finalize()"<< endreq;
  return StatusCode::SUCCESS;
}

PreElecAlgorithm::~PreElecAlgorithm()
{
  info()<< "Destruct:~PreElecAlgorithm()!"<< endreq;
}