#include <EGamma/EGammaCutBasedEleIdAnalyzer/src/EGammaCutBasedEleIdAnalyzer.cc>

Inheritance diagram for EGammaCutBasedEleIdAnalyzer:

Public Types
typedef std::vector< edm::Handle< edm::ValueMap< reco::IsoDeposit > > >	IsoDepositMaps

typedef std::vector< edm::Handle< edm::ValueMap< double > > >	IsoDepositVals

Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Public Member Functions
	EGammaCutBasedEleIdAnalyzer (const edm::ParameterSet &)

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Public Attributes
ElectronEffectiveArea::ElectronEffectiveAreaTarget	EAtarget

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_

edm::EDGetTokenT< reco::ConversionCollection >	conversionsToken_

std::string	EAtargetToken_

edm::EDGetTokenT< reco::GsfElectronCollection >	electronsToken_

TH1F *	h1_pt_

TH1F *	h1_pt_fbremeopin_

TH1F *	h1_pt_loose_

TH1F *	h1_pt_medium_

TH1F *	h1_pt_tight_

TH1F *	h1_pt_trig_

TH1F *	h1_pt_veto_

std::vector< edm::EDGetTokenT< edm::ValueMap< double > > >	isoValTokens_

edm::EDGetTokenT< reco::VertexCollection >	primaryVertexToken_

bool	printDebug_

edm::EDGetTokenT< double >	rhoIsoToken_

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 48 of file EGammaCutBasedEleIdAnalyzer.cc.

Member Typedef Documentation

◆ IsoDepositMaps

typedef std::vector<edm::Handle<edm::ValueMap<reco::IsoDeposit> > > EGammaCutBasedEleIdAnalyzer::IsoDepositMaps

Definition at line 50 of file EGammaCutBasedEleIdAnalyzer.cc.

◆ IsoDepositVals

typedef std::vector<edm::Handle<edm::ValueMap<double> > > EGammaCutBasedEleIdAnalyzer::IsoDepositVals

Definition at line 51 of file EGammaCutBasedEleIdAnalyzer.cc.

Constructor & Destructor Documentation

◆ EGammaCutBasedEleIdAnalyzer()

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

explicit

Definition at line 93 of file EGammaCutBasedEleIdAnalyzer.cc.

References beamSpotToken_, conversionsToken_, EAtargetToken_, electronsToken_, compareTotals::fs, edm::ParameterSet::getParameter(), h1_pt_, h1_pt_fbremeopin_, h1_pt_loose_, h1_pt_medium_, h1_pt_tight_, h1_pt_trig_, h1_pt_veto_, isoValTokens_, primaryVertexToken_, printDebug_, rhoIsoToken_, AlCaHLTBitMon_QueryRunRegistry::string, makeGlobalPositionRcd_cfg::tag, and edm::vector_transform().

                                                                                        {
   // get input parameters
   electronsToken_ = consumes<reco::GsfElectronCollection>(iConfig.getParameter<edm::InputTag>("electronsInputTag"));
   conversionsToken_ = consumes<reco::ConversionCollection>(iConfig.getParameter<edm::InputTag>("conversionsInputTag"));
   beamSpotToken_ = consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpotInputTag"));
   rhoIsoToken_ = consumes<double>(iConfig.getParameter<edm::InputTag>("rhoIsoInputTag"));
   primaryVertexToken_ = consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("primaryVertexInputTag"));
   isoValTokens_ =
       edm::vector_transform(iConfig.getParameter<std::vector<edm::InputTag> >("isoValInputTags"),
                             [this](edm::InputTag const &tag) { return consumes<edm::ValueMap<double> >(tag); });
   EAtargetToken_ = iConfig.getParameter<std::string>(
       "EAtarget");  //EleEANoCorr, EleEAData2011, EleEASummer11MC,EleEAFall11MC, EleEAData2012/
 
   // debug
   printDebug_ = iConfig.getParameter<bool>("printDebug");
 
   // output histograms
   edm::Service<TFileService> fs;
 
   h1_pt_ = fs->make<TH1F>("h1_pt", "pt", 100, 0.0, 100.0);
   h1_pt_veto_ = fs->make<TH1F>("h1_pt_veto", "pt (veto)", 100, 0.0, 100.0);
   h1_pt_loose_ = fs->make<TH1F>("h1_pt_loose", "pt (loose)", 100, 0.0, 100.0);
   h1_pt_medium_ = fs->make<TH1F>("h1_pt_medium", "pt (medium)", 100, 0.0, 100.0);
   h1_pt_tight_ = fs->make<TH1F>("h1_pt_tight", "pt (tight)", 100, 0.0, 100.0);
   h1_pt_trig_ = fs->make<TH1F>("h1_pt_trig", "pt (trig)", 100, 0.0, 100.0);
   h1_pt_fbremeopin_ = fs->make<TH1F>("h1_pt_fbremeopin", "pt (fbremeopin)", 100, 0.0, 100.0);
 }

Member Function Documentation

◆ analyze()

void EGammaCutBasedEleIdAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 126 of file EGammaCutBasedEleIdAnalyzer.cc.

                                                                                              {
   // electrons
   edm::Handle<reco::GsfElectronCollection> els_h;
   iEvent.getByToken(electronsToken_, els_h);
 
   // conversions
   edm::Handle<reco::ConversionCollection> conversions_h;
   iEvent.getByToken(conversionsToken_, conversions_h);
 
   // iso deposits
   IsoDepositVals isoVals(isoValTokens_.size());
   for (size_t j = 0; j < isoValTokens_.size(); ++j) {
     iEvent.getByToken(isoValTokens_[j], isoVals[j]);
   }
 
   // beam spot
   edm::Handle<reco::BeamSpot> beamspot_h;
   iEvent.getByToken(beamSpotToken_, beamspot_h);
   const reco::BeamSpot &beamSpot = *(beamspot_h.product());
 
   // vertices
   edm::Handle<reco::VertexCollection> vtx_h;
   iEvent.getByToken(primaryVertexToken_, vtx_h);
 
   // rho for isolation
   edm::Handle<double> rhoIso_h;
   iEvent.getByToken(rhoIsoToken_, rhoIso_h);
   double rhoIso = *(rhoIso_h.product());
 
   // loop on electrons
   unsigned int n = els_h->size();
   for (unsigned int i = 0; i < n; ++i) {
     // get reference to electron
     reco::GsfElectronRef ele(els_h, i);
 
     //
     // get particle flow isolation
     //
 
     double iso_ch = (*(isoVals)[0])[ele];
     double iso_em = (*(isoVals)[1])[ele];
     double iso_nh = (*(isoVals)[2])[ele];
 
     //
     // test ID
     //
 
     // working points
     bool veto = EgammaCutBasedEleId::PassWP(
         EgammaCutBasedEleId::VETO, ele, conversions_h, beamSpot, vtx_h, iso_ch, iso_em, iso_nh, rhoIso, EAtarget);
     bool loose = EgammaCutBasedEleId::PassWP(
         EgammaCutBasedEleId::LOOSE, ele, conversions_h, beamSpot, vtx_h, iso_ch, iso_em, iso_nh, rhoIso, EAtarget);
     bool medium = EgammaCutBasedEleId::PassWP(
         EgammaCutBasedEleId::MEDIUM, ele, conversions_h, beamSpot, vtx_h, iso_ch, iso_em, iso_nh, rhoIso, EAtarget);
     bool tight = EgammaCutBasedEleId::PassWP(
         EgammaCutBasedEleId::TIGHT, ele, conversions_h, beamSpot, vtx_h, iso_ch, iso_em, iso_nh, rhoIso, EAtarget);
 
     // eop/fbrem cuts for extra tight ID
     bool fbremeopin = EgammaCutBasedEleId::PassEoverPCuts(ele);
 
     // cuts to match tight trigger requirements
     bool trigtight = EgammaCutBasedEleId::PassTriggerCuts(EgammaCutBasedEleId::TRIGGERTIGHT, ele);
 
     // for 2011 WP70 trigger
     bool trigwp70 = EgammaCutBasedEleId::PassTriggerCuts(EgammaCutBasedEleId::TRIGGERWP70, ele);
 
     //
     // fill histograms
     //
 
     h1_pt_->Fill(ele->pt());
     if (veto)
       h1_pt_veto_->Fill(ele->pt());
     if (loose)
       h1_pt_loose_->Fill(ele->pt());
     if (medium)
       h1_pt_medium_->Fill(ele->pt());
     if (tight)
       h1_pt_tight_->Fill(ele->pt());
     if (trigtight)
       h1_pt_trig_->Fill(ele->pt());
     if (fbremeopin)
       h1_pt_fbremeopin_->Fill(ele->pt());
 
     //
     // print decisions
     //
 
     if (printDebug_) {
       printf("%u %u %llu : ", iEvent.id().run(), iEvent.luminosityBlock(), iEvent.id().event());
       printf("veto(%i), ", veto);
       printf("loose(%i), ", loose);
       printf("medium(%i), ", medium);
       printf("tight(%i), ", tight);
       printf("trigtight(%i), ", trigtight);
       printf("trigwp70(%i), ", trigwp70);
       printf("fbremeopin(%i)\n", fbremeopin);
     }
   }
 }

◆ beginJob()

void EGammaCutBasedEleIdAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 228 of file EGammaCutBasedEleIdAnalyzer.cc.

References EAtarget, EAtargetToken_, ElectronEffectiveArea::kEleEAData2011, ElectronEffectiveArea::kEleEAData2012, ElectronEffectiveArea::kEleEAFall11MC, ElectronEffectiveArea::kEleEANoCorr, and ElectronEffectiveArea::kEleEASummer11MC.

                                            {
   if (EAtargetToken_ == "EleEANoCorr")
     EAtarget = ElectronEffectiveArea::kEleEANoCorr;
   else if (EAtargetToken_ == "EleEAData2011")
     EAtarget = ElectronEffectiveArea::kEleEAData2011;
   else if (EAtargetToken_ == "EleEASummer11MC")
     EAtarget = ElectronEffectiveArea::kEleEASummer11MC;
   else if (EAtargetToken_ == "EleEAFall11MC")
     EAtarget = ElectronEffectiveArea::kEleEAFall11MC;
   else if (EAtargetToken_ == "EleEAData2012")
     EAtarget = ElectronEffectiveArea::kEleEAData2012;
   else
     EAtarget = ElectronEffectiveArea::kEleEAData2012;
 }

◆ fillDescriptions()

void EGammaCutBasedEleIdAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 244 of file EGammaCutBasedEleIdAnalyzer.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                                              {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }