#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

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

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)

virtual	~ProducerBase () noexcept(false)

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

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

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::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 12 of file ZtoEEEventSelector.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 23 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

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

overridevirtual

Implements edm::stream::EDFilterBase.

Definition at line 50 of file ZtoEEEventSelector.cc.

References ecalDrivenElectronSeedsParameters_cff::beamSpot, bsTag_, bsToken_, vertices_cff::chi2, DEFINE_FWK_MODULE, electronIdCutBased_cfi::deltaEtaIn, electronIdCutBased_cfi::deltaPhiIn, electronTag_, electronToken_, edm::Event::getByToken(), electronIdCutBased_cfi::hOverE, 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::BeamSpot::position(), reco::tau::disc::Pt(), reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt, reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt, reco::GsfElectron::PflowIsolationVariables::sumPhotonEt, and reco::GsfElectron::PflowIsolationVariables::sumPUPt.

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

                                                                              {
   // 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;
 }

Member Data Documentation

const edm::InputTag ZtoEEEventSelector::bsTag_

private

Definition at line 21 of file ZtoEEEventSelector.h.

Referenced by filter().

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

private

Definition at line 23 of file ZtoEEEventSelector.h.

Referenced by filter().

const edm::InputTag ZtoEEEventSelector::electronTag_

private

Definition at line 20 of file ZtoEEEventSelector.h.

Referenced by filter().

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

private

Definition at line 22 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxD0_

private

Definition at line 36 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxDeltaEtaInEB_

private

Definition at line 28 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxDeltaEtaInEE_

private

Definition at line 32 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxDeltaPhiInEB_

private

Definition at line 27 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxDeltaPhiInEE_

private

Definition at line 31 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxDz_

private

Definition at line 37 of file ZtoEEEventSelector.h.

const double ZtoEEEventSelector::maxEta_

private

Definition at line 25 of file ZtoEEEventSelector.h.

const double ZtoEEEventSelector::maxHOEEB_

private

Definition at line 29 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxHOEEE_

private

Definition at line 33 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxInvMass_

private

Definition at line 43 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxIso_

private

Definition at line 40 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxNormChi2_

private

Definition at line 35 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxSigmaiEiEEB_

private

Definition at line 30 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::maxSigmaiEiEEE_

private

Definition at line 34 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::minInvMass_

private

Definition at line 42 of file ZtoEEEventSelector.h.

Referenced by filter().

const int ZtoEEEventSelector::minPixelHits_

private

Definition at line 38 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::minPt_

private

Definition at line 26 of file ZtoEEEventSelector.h.

Referenced by filter().

const double ZtoEEEventSelector::minPtHighest_

private

Definition at line 41 of file ZtoEEEventSelector.h.

Referenced by filter().

const int ZtoEEEventSelector::minStripHits_

private

Definition at line 39 of file ZtoEEEventSelector.h.

Referenced by filter().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation