In This Package:

SimHitSplitAlgorithm.cc

Go to the documentation of this file.

#include "SimHitSplitSvc/SimHitSplitAlgorithm.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;

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

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

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

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

  vector<SimHitHeader*> shhs;
  SimHeader* sh = 0;
  sh = get<SimHeader>(evtSvc(), SimHeaderLocation::Default);
  StatusCode sc = m_simHitSplitSvc->getSimHitHeader(sh, shhs, m_hitEarliest, m_hitLatest);
  if(sc.isFailure()) return StatusCode::SUCCESS;

  for(vector<SimHitHeader*>::iterator vit = shhs.begin(); vit != shhs.end(); vit++)
  {
    if(!(*vit))
    {
      error()<< "Error when get one SimHitHeader!"<< endreq;
      return StatusCode::FAILURE;
    }
    info()<< "SimHitHeader's address: "<< (*vit)<< endreq;

    map<short int, SimHitCollection*> shc = (*vit)->hitCollection();
    for(map<short int, SimHitCollection*>::iterator mit = shc.begin(); mit != shc.end(); mit++)
    {
      vector<SimHit*> vsh = mit->second->collection();
      info()<< MSG::DEBUG<< "In Det"<< mit->first
            << " : SimHit's num = "<< vsh.size()<< endreq;
    }
    //delete (*vit);
  }     
  if(shhs.size() != 1 || shhs[0] != sh->hits())
  {
    for(vector<SimHitHeader*>::size_type i = 0; i < shhs.size(); i++)
    {
      SimHitHeader::hc_map& maphc = shhs[i]->hitCollection();
      for (SimHitHeader::hc_map::iterator it = maphc.begin(); it != maphc.end(); ++it)
      {
        it->second->collection().clear();
      }
      //shhs[i]->hitCollection().clear();
      delete shhs[i];
    }
  }
  info()<< "SimHitSplitAlgorithm execute() finalized!" << endreq;
  return StatusCode::SUCCESS;
}

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

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

---