#include <ZtoEEEventSelector.h>

Inheritance diagram for ZtoEEEventSelector:

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

	ZtoEEEventSelector (const edm::ParameterSet &)

Public Member Functions inherited from edm::stream::EDFilter<>
	EDFilter ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

	~EDFilterBase () 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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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 &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

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
const edm::InputTag	bsTag_

const edm::EDGetTokenT< reco::BeamSpot >	bsToken_

const edm::InputTag	electronTag_

const edm::EDGetTokenT< reco::GsfElectronCollection >	electronToken_

const double	maxD0_

const double	maxDeltaEtaInEB_

const double	maxDeltaEtaInEE_

const double	maxDeltaPhiInEB_

const double	maxDeltaPhiInEE_

const double	maxDz_

const double	maxEta_

const double	maxHOEEB_

const double	maxHOEEE_

const double	maxInvMass_

const double	maxIso_

const double	maxNormChi2_

const double	maxSigmaiEiEEB_

const double	maxSigmaiEiEEE_

const double	minInvMass_

const int	minPixelHits_

const double	minPt_

const double	minPtHighest_

const int	minStripHits_

Additional Inherited Members
Public Types inherited from edm::stream::EDFilter<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDFilterBase
typedef EDFilterAdaptorBase	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::stream::EDFilterBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)

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

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

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

Detailed Description

Definition at line 13 of file ZtoEEEventSelector.h.

Constructor & Destructor Documentation

◆ ZtoEEEventSelector()

ZtoEEEventSelector::ZtoEEEventSelector ( const edm::ParameterSet & ps )

explicit

Definition at line 22 of file ZtoEEEventSelector.cc.

     : electronTag_(ps.getUntrackedParameter<edm::InputTag>("electronInputTag", edm::InputTag("gedGsfElectrons"))),
       bsTag_(ps.getUntrackedParameter<edm::InputTag>("offlineBeamSpot", edm::InputTag("offlineBeamSpot"))),
       electronToken_(consumes<reco::GsfElectronCollection>(electronTag_)),
       bsToken_(consumes<reco::BeamSpot>(bsTag_)),
       maxEta_(ps.getUntrackedParameter<double>("maxEta", 2.4)),
       minPt_(ps.getUntrackedParameter<double>("minPt", 5)),
       maxDeltaPhiInEB_(ps.getUntrackedParameter<double>("maxDeltaPhiInEB", .15)),
       maxDeltaEtaInEB_(ps.getUntrackedParameter<double>("maxDeltaEtaInEB", .007)),
       maxHOEEB_(ps.getUntrackedParameter<double>("maxHOEEB", .12)),
       maxSigmaiEiEEB_(ps.getUntrackedParameter<double>("maxSigmaiEiEEB", .01)),
       maxDeltaPhiInEE_(ps.getUntrackedParameter<double>("maxDeltaPhiInEE", .1)),
       maxDeltaEtaInEE_(ps.getUntrackedParameter<double>("maxDeltaEtaInEE", .009)),
       maxHOEEE_(ps.getUntrackedParameter<double>("maxHOEEB_", .10)),
       maxSigmaiEiEEE_(ps.getUntrackedParameter<double>("maxSigmaiEiEEE", .03)),
       maxNormChi2_(ps.getUntrackedParameter<double>("maxNormChi2", 10)),
       maxD0_(ps.getUntrackedParameter<double>("maxD0", 0.02)),
       maxDz_(ps.getUntrackedParameter<double>("maxDz", 20.)),
       minPixelHits_(ps.getUntrackedParameter<uint32_t>("minPixelHits", 1)),
       minStripHits_(ps.getUntrackedParameter<uint32_t>("minStripHits", 8)),
       maxIso_(ps.getUntrackedParameter<double>("maxIso", 0.3)),
       minPtHighest_(ps.getUntrackedParameter<double>("minPtHighest", 24)),
       minInvMass_(ps.getUntrackedParameter<double>("minInvMass", 60)),
       maxInvMass_(ps.getUntrackedParameter<double>("maxInvMass", 120)) {}

Member Function Documentation

◆ filter()

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

overridevirtual

Implements edm::stream::EDFilterBase.

Definition at line 47 of file ZtoEEEventSelector.cc.

                                                                              {
   // Read Electron Collection
   edm::Handle<reco::GsfElectronCollection> electronColl;
   iEvent.getByToken(electronToken_, electronColl);
  
   edm::Handle<reco::BeamSpot> beamSpot;
   iEvent.getByToken(bsToken_, beamSpot);
  
   std::vector<TLorentzVector> list;
   std::vector<int> chrgeList;
  
   if (electronColl.isValid()) {
     for (auto const& ele : *electronColl) {
       if (!ele.ecalDriven())
         continue;
       if (ele.pt() < minPt_)
         continue;
       // set a max Eta cut
       if (!(ele.isEB() || ele.isEE()))
         continue;
  
       double hOverE = ele.hadronicOverEm();
       double sigmaee = ele.sigmaIetaIeta();
       double deltaPhiIn = ele.deltaPhiSuperClusterTrackAtVtx();
       double deltaEtaIn = ele.deltaEtaSuperClusterTrackAtVtx();
  
       // separate cut for barrel and endcap
       if (ele.isEB()) {
         if (fabs(deltaPhiIn) >= maxDeltaPhiInEB_ && fabs(deltaEtaIn) >= maxDeltaEtaInEB_ && hOverE >= maxHOEEB_ &&
             sigmaee >= maxSigmaiEiEEB_)
           continue;
       } else if (ele.isEE()) {
         if (fabs(deltaPhiIn) >= maxDeltaPhiInEE_ && fabs(deltaEtaIn) >= maxDeltaEtaInEE_ && hOverE >= maxHOEEE_ &&
             sigmaee >= maxSigmaiEiEEE_)
           continue;
       }
  
       reco::GsfTrackRef trk = ele.gsfTrack();
       if (!trk.isNonnull())
         continue;  // only electrons with tracks
       double chi2 = trk->chi2();
       double ndof = trk->ndof();
       double chbyndof = (ndof > 0) ? chi2 / ndof : 0;
       if (chbyndof >= maxNormChi2_)
         continue;
  
       double trkd0 = trk->d0();
       if (beamSpot.isValid()) {
         trkd0 = -(trk->dxy(beamSpot->position()));
       } else {
         edm::LogError("ZtoEEEventSelector") << "Error >> Failed to get BeamSpot for label: " << bsTag_;
       }
       if (std::fabs(trkd0) >= maxD0_)
         continue;
  
       const reco::HitPattern& hitp = trk->hitPattern();
       int nPixelHits = hitp.numberOfValidPixelHits();
       if (nPixelHits < minPixelHits_)
         continue;
  
       int nStripHits = hitp.numberOfValidStripHits();
       if (nStripHits < minStripHits_)
         continue;
  
       // PF Isolation
       reco::GsfElectron::PflowIsolationVariables pfIso = ele.pfIsolationVariables();
       float absiso =
           pfIso.sumChargedHadronPt + std::max(0.0, pfIso.sumNeutralHadronEt + pfIso.sumPhotonEt - 0.5 * pfIso.sumPUPt);
       float eiso = absiso / (ele.pt());
       if (eiso > maxIso_)
         continue;
  
       TLorentzVector lv;
       lv.SetPtEtaPhiE(ele.pt(), ele.eta(), ele.phi(), ele.energy());
       list.push_back(lv);
       chrgeList.push_back(ele.charge());
     }
   } else {
     edm::LogError("ZtoEEEventSelector") << "Error >> Failed to get ElectronCollection for label: " << electronTag_;
   }
   if (list.size() < 2)
     return false;
   if (chrgeList[0] + chrgeList[1] != 0)
     return false;
  
   if (list[0].Pt() < minPtHighest_)
     return false;
   TLorentzVector zv = list[0] + list[1];
   if (zv.M() < minInvMass_ || zv.M() > maxInvMass_)
     return false;
  
   return true;
 }

References pwdgSkimBPark_cfi::beamSpot, bsTag_, bsToken_, hltPixelTracks_cff::chi2, electronIdCutBasedExt_cfi::deltaEtaIn, electronIdCutBasedExt_cfi::deltaPhiIn, electronTag_, electronToken_, EgHLTOffHistBins_cfi::hOverE, iEvent, edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), list(), SiStripPI::max, maxD0_, maxDeltaEtaInEB_, maxDeltaEtaInEE_, maxDeltaPhiInEB_, maxDeltaPhiInEE_, maxHOEEB_, maxHOEEE_, maxInvMass_, maxIso_, maxNormChi2_, maxSigmaiEiEEB_, maxSigmaiEiEEE_, minInvMass_, minPixelHits_, minPt_, minPtHighest_, minStripHits_, ndof, reco::HitPattern::numberOfValidPixelHits(), reco::HitPattern::numberOfValidStripHits(), reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, reco::GsfElectron::PflowIsolationVariables::sumPhotonEt, and reco::GsfElectron::PflowIsolationVariables::sumPUPt.

Member Data Documentation

◆ bsTag_

const edm::InputTag ZtoEEEventSelector::bsTag_

private

Definition at line 22 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ bsToken_

const edm::EDGetTokenT<reco::BeamSpot> ZtoEEEventSelector::bsToken_

private

Definition at line 24 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ electronTag_

const edm::InputTag ZtoEEEventSelector::electronTag_

private

Definition at line 21 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ electronToken_

const edm::EDGetTokenT<reco::GsfElectronCollection> ZtoEEEventSelector::electronToken_

private

Definition at line 23 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxD0_

const double ZtoEEEventSelector::maxD0_

private

Definition at line 37 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxDeltaEtaInEB_

const double ZtoEEEventSelector::maxDeltaEtaInEB_

private

Definition at line 29 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxDeltaEtaInEE_

const double ZtoEEEventSelector::maxDeltaEtaInEE_

private

Definition at line 33 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxDeltaPhiInEB_

const double ZtoEEEventSelector::maxDeltaPhiInEB_

private

Definition at line 28 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxDeltaPhiInEE_

const double ZtoEEEventSelector::maxDeltaPhiInEE_

private

Definition at line 32 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxDz_

const double ZtoEEEventSelector::maxDz_

private

Definition at line 38 of file ZtoEEEventSelector.h.

◆ maxEta_

const double ZtoEEEventSelector::maxEta_

private

Definition at line 26 of file ZtoEEEventSelector.h.

◆ maxHOEEB_

const double ZtoEEEventSelector::maxHOEEB_

private

Definition at line 30 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxHOEEE_

const double ZtoEEEventSelector::maxHOEEE_

private

Definition at line 34 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxInvMass_

const double ZtoEEEventSelector::maxInvMass_

private

Definition at line 44 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxIso_

const double ZtoEEEventSelector::maxIso_

private

Definition at line 41 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxNormChi2_

const double ZtoEEEventSelector::maxNormChi2_

private

Definition at line 36 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxSigmaiEiEEB_

const double ZtoEEEventSelector::maxSigmaiEiEEB_

private

Definition at line 31 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ maxSigmaiEiEEE_

const double ZtoEEEventSelector::maxSigmaiEiEEE_

private

Definition at line 35 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ minInvMass_

const double ZtoEEEventSelector::minInvMass_

private

Definition at line 43 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ minPixelHits_

const int ZtoEEEventSelector::minPixelHits_

private

Definition at line 39 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ minPt_

const double ZtoEEEventSelector::minPt_

private

Definition at line 27 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ minPtHighest_

const double ZtoEEEventSelector::minPtHighest_

private

Definition at line 42 of file ZtoEEEventSelector.h.

Referenced by filter().

◆ minStripHits_

const int ZtoEEEventSelector::minStripHits_

private

Definition at line 40 of file ZtoEEEventSelector.h.

Referenced by filter().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ZtoEEEventSelector()

Member Function Documentation

◆ filter()

Member Data Documentation

◆ bsTag_

◆ bsToken_

◆ electronTag_

◆ electronToken_

◆ maxD0_

◆ maxDeltaEtaInEB_

◆ maxDeltaEtaInEE_

◆ maxDeltaPhiInEB_

◆ maxDeltaPhiInEE_

◆ maxDz_

◆ maxEta_

◆ maxHOEEB_

◆ maxHOEEE_

◆ maxInvMass_

◆ maxIso_

◆ maxNormChi2_

◆ maxSigmaiEiEEB_

◆ maxSigmaiEiEEE_

◆ minInvMass_

◆ minPixelHits_

◆ minPt_

◆ minPtHighest_

◆ minStripHits_