PurgeUnstableParticles.cxx

Go to the documentation of this file.

#include "GeneratorModule/PurgeUnstableParticles.h"
#include "GeneratorObjects/McEventCollection.h"
 
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/ISvcLocator.h"
 
#include "GaudiKernel/IDataManagerSvc.h"
 
#include "GaudiKernel/DataSvc.h"
#include "GaudiKernel/SmartDataPtr.h"
 
static const AlgFactory<PurgeUnstableParticles> Factory;
const IAlgFactory&                              PurgeUnstableParticlesFactory = Factory;
 
PurgeUnstableParticles::PurgeUnstableParticles( const std::string& name,
                                                ISvcLocator*       pSvcLocator )
    : Algorithm( name, pSvcLocator ) {}
 
StatusCode PurgeUnstableParticles::initialize() {
 
  IDataManagerSvc* m_DMS;
  StatusCode       sc =
      serviceLocator()->getService( "EventDataSvc", IID_IDataManagerSvc, (IInterface*&)m_DMS );
  if ( sc.isFailure() ) { return sc; }
 
  std::string rootEventName = m_DMS->rootEventName();
  m_mcEventColl             = rootEventName + "/McEventCollection";
 
  return StatusCode::SUCCESS;
}
 
StatusCode PurgeUnstableParticles::execute() {
 
  SmartDataPtr<McEventCollection> mcCollptr( eventDataService(), m_mcEventColl );
  if ( mcCollptr )
  {
    McEventCollection::iterator it;
    for ( it = mcCollptr->begin(); it != mcCollptr->end(); it++ )
    { ( *it )->purgeUnstableParticles(); }
  }
 
  return StatusCode::SUCCESS;
}
 
StatusCode PurgeUnstableParticles::finalize() { return StatusCode::SUCCESS; }