LHEGenericFilter Class Reference

#include <LHEGenericFilter.h>

Inheritance diagram for LHEGenericFilter:

Public Member Functions
	LHEGenericFilter (const edm::ParameterSet &)
	~LHEGenericFilter () override
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDFilter () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Types
enum	logic_ { LT, GT, EQ, NE }

Private Member Functions
void	endJob () override
bool	filter (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::string	acceptLogic_
int	numRequired_
std::vector< int >	particleID_
int	passedEvents_
edm::EDGetTokenT< LHEEventProduct >	src_
int	totalEvents_
logic_	whichlogic

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
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 &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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 39 of file LHEGenericFilter.h.

Member Enumeration Documentation

enum LHEGenericFilter::logic_

private

Enumerator
LT
GT
EQ
NE

Definition at line 59 of file LHEGenericFilter.h.

{ LT, GT, EQ, NE};

Constructor & Destructor Documentation

LHEGenericFilter::LHEGenericFilter ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 6 of file LHEGenericFilter.cc.

References acceptLogic_, EQ, edm::ParameterSet::getParameter(), GT, LT, NE, src_, and whichlogic.

                                                                 :
numRequired_(iConfig.getParameter<int>("NumRequired")),
acceptLogic_(iConfig.getParameter<std::string>("AcceptLogic")),
particleID_(iConfig.getParameter< std::vector<int> >("ParticleID")),
totalEvents_(0), passedEvents_(0)
{
  //here do whatever other initialization is needed
  src_ = consumes<LHEEventProduct>(iConfig.getParameter<edm::InputTag>("src"));
  
  if(acceptLogic_=="LT") whichlogic = LT;
  else if(acceptLogic_=="GT") whichlogic = GT;
  else if(acceptLogic_=="EQ") whichlogic = EQ;
  else if(acceptLogic_=="NE") whichlogic = NE;
  else edm::LogError  ("cat_A")   << "wrong input for AcceptLogic string";

  
}

LHEGenericFilter::~LHEGenericFilter ( )

override

Definition at line 24 of file LHEGenericFilter.cc.

{

  // do anything here that needs to be done at destruction time
  // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void LHEGenericFilter::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 71 of file LHEGenericFilter.cc.

References DEFINE_FWK_MODULE, passedEvents_, and totalEvents_.

                              {
  edm::LogInfo("LHEGenericFilter") << "=== Results of LHEGenericFilter: passed "
  << passedEvents_ << "/" << totalEvents_ << " events" << std::endl;
}

bool LHEGenericFilter::filter ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 34 of file LHEGenericFilter.cc.

References funct::abs(), EQ, edm::Event::getByToken(), GT, LHEEventProduct::hepeup(), mps_fire::i, lhef::HEPEUP::IDUP, lhef::HEPEUP::ISTUP, LT, NE, numRequired_, lhef::HEPEUP::NUP, particleID_, passedEvents_, src_, totalEvents_, and whichlogic.

{
  edm::Handle< LHEEventProduct > EvtHandle ;
  iEvent.getByToken( src_ , EvtHandle ) ;
  
  totalEvents_++;
  int nFound = 0;
  
  for (int i = 0; i < EvtHandle->hepeup().NUP; ++i) {
    if (EvtHandle->hepeup().ISTUP[i] != 1) { // keep only outgoing particles
      continue;
    }
    for (unsigned int j = 0; j < particleID_.size(); ++j) {
      if (particleID_[j] == 0 || abs(particleID_[j]) == abs(EvtHandle->hepeup().IDUP[i]) ) {
        nFound++;
        break; // only match a given particle once!
      }
    } // loop over targets
    
  } // loop over particles
  
  // event accept/reject logic
  if (
         (whichlogic==LT && nFound < numRequired_)
      || (whichlogic==GT && nFound > numRequired_)
      || (whichlogic==EQ && nFound == numRequired_)
      || (whichlogic==NE && nFound != numRequired_)
      ) {
    passedEvents_++;
    return true;
  } else {
    return false;
  }
  
}

Member Data Documentation

std::string LHEGenericFilter::acceptLogic_

private

Definition at line 52 of file LHEGenericFilter.h.

Referenced by LHEGenericFilter().

int LHEGenericFilter::numRequired_

private

Definition at line 51 of file LHEGenericFilter.h.

Referenced by filter().

std::vector<int> LHEGenericFilter::particleID_

private

Definition at line 56 of file LHEGenericFilter.h.

Referenced by filter().

int LHEGenericFilter::passedEvents_

private

Definition at line 58 of file LHEGenericFilter.h.

Referenced by endJob(), and filter().

edm::EDGetTokenT<LHEEventProduct> LHEGenericFilter::src_

private

Definition at line 50 of file LHEGenericFilter.h.

Referenced by filter(), and LHEGenericFilter().

int LHEGenericFilter::totalEvents_

private

Definition at line 57 of file LHEGenericFilter.h.

Referenced by endJob(), and filter().

logic_ LHEGenericFilter::whichlogic

private

Definition at line 60 of file LHEGenericFilter.h.

Referenced by filter(), and LHEGenericFilter().