In This Package:

MixInputAlgorithm.cc

Go to the documentation of this file.

#include "MixInputSvc/MixInputAlgorithm.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;

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

StatusCode MixInputAlgorithm::initialize()
{
  StatusCode status= this->GaudiAlgorithm::initialize();
  if(status.isFailure()) 
  { 
    error()<< "Error when initializing  GaudiAlgorithm"<< endreq;
    return status;
  }
   
  info()<< "MixInputAlgorithm ininalize()"<< endreq; 
  IService* isvc = 0;
  StatusCode sc = serviceLocator()->service("RootIOCnvSvc", isvc, false);
  if(sc.isFailure()) 
  {
    error()<< "Conversion service RootIOCnvSvc"<< " could not be retrieved"<< endreq;
    return sc;
  }
  isvc->addRef();
  sc = isvc->queryInterface(IMixInputSvc::interfaceID(), (void**)&m_mixInputSvc);
  if(sc.isFailure()) 
  {
    error()<< "Conversion service RootIOCnvSvc"<< " does not implement IRootIOCnvSvc"<< endreq;
  }
  return sc;
}

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

  SimHeader* sh = new SimHeader();
  GenHeader* gh = new GenHeader();
  m_mixInputSvc->prepareStream();
  m_mixInputSvc->getSimHeader(gh, sh);

  const SimHitHeader* shh = sh->hits();
  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;
  }

  info()<< "GenHeader: "<< gh<< "--"<< "SimHeader: "<< sh<< endreq;
                              
  if(0 != sh && 0 != gh)
  {
    // Using some functions in RegistrationSequence
    RegistrationSequence* seq = get<RegistrationSequence>(evtSvc(), RegistrationSequenceLocation::Default);
    debug()<< "ATTENTION!===========> Size = "<< seq->size()<< endreq;
    seq->setEarliest(sh->earliest());
    seq->setLatest(sh->latest());

    //If the sample coming from MDC, the codes below needed
    /*IRegistrationSequence::Registrations& regs = const_cast<IRegistrationSequence::Registrations&>(seq->registrations());
    IRegistrationSequence::Registrations::iterator regit;
    for(regit = regs.begin(); regit != regs.end(); regit++)
    {
      info()<< regit->path()<< endreq;
      if(regit->path() == "/Event/Readout/ReadoutHeader")
      {
        DataObject*    child = const_cast<DataObject*>(regit->object());   
        ReadoutHeader* reoh  = dynamic_cast<ReadoutHeader*>(child);
        info()<< reoh->readout()->triggerTime().GetSeconds()<< endreq;
        regit->setStore(false);
        info()<< "Not Storing: "<< regit->path()<< endreq;
      }
      }*/

    // Put SimHeader(GenHeader) into TES
    put(evtSvc(), gh, GenHeaderLocation::Default);
    put(evtSvc(), sh, SimHeaderLocation::Default);
  }

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

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

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

---