#include <MCPdgIndexFilter.h>

Inheritance diagram for MCPdgIndexFilter:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)

	MCPdgIndexFilter (const edm::ParameterSet &)

	~MCPdgIndexFilter ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

ModuleDescription const &	moduleDescription () const

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
bool	pass (const edm::Event &)

Private Attributes
const std::vector< unsigned >	index

const unsigned	maxIndex

const std::vector< int >	pdgID

const std::string	tag

const bool	taggingMode

const edm::EDGetTokenT < edm::HepMCProduct >	token_

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 24 of file MCPdgIndexFilter.h.

Constructor & Destructor Documentation

MCPdgIndexFilter::MCPdgIndexFilter ( const edm::ParameterSet & cfg )

explicit

Definition at line 5 of file MCPdgIndexFilter.cc.

References index, pdgID, tag, and taggingMode.

                                                              :
   token_(consumes<edm::HepMCProduct>(cfg.getUntrackedParameter("moduleLabel",std::string("generator")))),
   pdgID(cfg.getParameter<std::vector<int> >("PdgId")),
   index(cfg.getParameter<std::vector<unsigned> >("Index")),
   maxIndex(*std::max_element(index.begin(),index.end())),
   taggingMode(cfg.getUntrackedParameter<bool>("TagMode",false)),
   tag(cfg.getUntrackedParameter<std::string>("Tag",""))
 {
   if (pdgID.size() != index.size())
     edm::LogWarning("MCPdgIndexFilter")
       << "Configuration Error :"
       << "Sizes of array parameters 'PdgId' and 'Index' differ.";
 
   if (taggingMode) {
     produces<bool>(tag);
     edm::LogInfo("TagMode") << "Filter result in '" << tag << "', filtering disabled.";
   }
 }

MCPdgIndexFilter::~MCPdgIndexFilter ( )

inline

Definition at line 27 of file MCPdgIndexFilter.h.

27 {};

Member Function Documentation

bool MCPdgIndexFilter::filter	(	edm::Event &	evt,
		const edm::EventSetup &
	)

virtual

Implements edm::EDFilter.

Definition at line 25 of file MCPdgIndexFilter.cc.

References LogDebug, pass(), edm::Event::put(), query::result, tag, and taggingMode.

                                                                  {
   bool result = pass(evt);
   LogDebug("FilterResult") << (result?"Pass":"Fail");
   if (!taggingMode) return result;
   evt.put( std::auto_ptr<bool>(new bool(result)), tag);
   return true;
 }

bool MCPdgIndexFilter::pass ( const edm::Event & evt )

private

Definition at line 34 of file MCPdgIndexFilter.cc.

References MCTruth::genEvent, edm::Event::getByToken(), i, index, j, LogDebug, maxIndex, AlCaHLTBitMon_ParallelJobs::p, pdgID, and token_.

Referenced by filter().

                                                {
   edm::Handle<edm::HepMCProduct> hepmc;
   evt.getByToken(token_, hepmc);
 
   const HepMC::GenEvent * genEvent = hepmc->GetEvent();
 
   HepMC::GenEvent::particle_const_iterator
     p(genEvent->particles_begin()),
     p_end(genEvent->particles_end());
 
   for ( unsigned i=0; p!=p_end && i<=maxIndex; ++p, i++ ) {
     LogDebug("Particle") << "index: " << i << "   pdgID: " << (*p)->pdg_id();
     for (unsigned j = 0; j < pdgID.size(); j++) {
       if (i==index[j] && pdgID[j] != (*p)->pdg_id())
     return false;
     }
   }
   return true;
 }