#include <HiggsToZZ4LeptonsPreFilter.h>

Inheritance diagram for HiggsToZZ4LeptonsPreFilter:

Public Member Functions
bool	filter (edm::Event &, const edm::EventSetup &) override
	Get event properties to send to builder to fill seed collection. More...

	HiggsToZZ4LeptonsPreFilter (const edm::ParameterSet &)

	~HiggsToZZ4LeptonsPreFilter () 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 Attributes
bool	debug

int	evt

edm::EDGetTokenT< reco::GenParticleCollection >	genToken

int	ikept

int	leptonFlavour

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 28 of file HiggsToZZ4LeptonsPreFilter.h.

Constructor & Destructor Documentation

HiggsToZZ4LeptonsPreFilter::HiggsToZZ4LeptonsPreFilter ( const edm::ParameterSet & pset )

explicit

Definition at line 32 of file HiggsToZZ4LeptonsPreFilter.cc.

References debug, edm::ParameterSet::getParameter(), and StEvtSolProducer_cfi::leptonFlavour.

                                                                                   {
 
 // LeptonFlavour
 // 0 = no tau
 // 1 = 4 mu
 // 2 = 4 e
 // 3 = 2e 2mu
 
   // Local Debug flag
   debug              = pset.getParameter<bool>("DebugHiggsToZZ4LeptonsPreFilter");
   leptonFlavour      = pset.getParameter<int>("HiggsToZZ4LeptonsPreFilterLeptonFlavour");
   genToken           = consumes<GenParticleCollection>(edm::InputTag("genParticles"));
 
   ikept = 0;
   evt = 0;
 }

HiggsToZZ4LeptonsPreFilter::~HiggsToZZ4LeptonsPreFilter ( )

override

Definition at line 51 of file HiggsToZZ4LeptonsPreFilter.cc.

References gather_cfg::cout.

                                                         {
 
   std::cout << "number of events processed: " << evt << std::endl;
   std::cout << "number of events kept: " << ikept << std::endl;
 
 }

Member Function Documentation

bool HiggsToZZ4LeptonsPreFilter::filter	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

override

Get event properties to send to builder to fill seed collection.

Definition at line 60 of file HiggsToZZ4LeptonsPreFilter.cc.

References GenHFHadronMatcher_cfi::genParticles, edm::HandleBase::isValid(), and StEvtSolProducer_cfi::leptonFlavour.

                                                                                     {
 
   bool keepEvent   = false;
   evt++;
 
   bool FourL    = false;
   bool FourE    = false;
   bool FourM    = false;
   bool TwoETwoM = false;
 
   // get gen particle candidates
   Handle<GenParticleCollection> genParticles;
   event.getByToken(genToken, genParticles);
 
   if ( genParticles.isValid() ) {
 
   int nElec = 0;
   int nMuon = 0;
 
   for (GenParticleCollection::const_iterator mcIter=genParticles->begin(); mcIter!=genParticles->end(); ++mcIter) {
     // Muons:
     if ( mcIter->pdgId() == 13 || mcIter->pdgId() == -13) {
       // Mother is a Z
       if ( mcIter->mother()->pdgId() == 23 ) {
        // In fiducial volume:
         if ( mcIter->pt() < 3 ) continue;
         if ( mcIter->eta() > -2.4 && mcIter->eta() < 2.4 ) nMuon++;
       }
     }
     // Electrons:
     if ( mcIter->pdgId() == 11 || mcIter->pdgId() == -11)
       // Mother is a Z
       if ( mcIter->mother()->pdgId() == 23 ) {
         // In fiducial volume:
         if ( mcIter->pt() < 3 ) continue;
         if ( mcIter->eta() > -2.5 && mcIter->eta() < 2.5 ) nElec++;
       }
     }
 
 
     if (nElec > 3) FourE = true;
     if (nMuon > 3) FourM = true;
     if (nMuon > 1 && nElec > 1) TwoETwoM = true;
     if ( FourE || FourM || TwoETwoM ) FourL = true;
 
     if ( leptonFlavour == 0 && FourL    ) keepEvent = true;
     if ( leptonFlavour == 1 && FourM    ) keepEvent = true;
     if ( leptonFlavour == 2 && FourE    ) keepEvent = true;
     if ( leptonFlavour == 3 && TwoETwoM ) keepEvent = true;
 
   }
 
   if (keepEvent ) ikept++;
 
   return keepEvent;
 
 }