#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::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

	~EGammaCutBasedEleIdAnalyzer () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

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

std::vector< ESProxyIndex > 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 &	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)

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Public Attributes
ElectronEffectiveArea::ElectronEffectiveAreaTarget	EAtarget

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

void	beginJob () override

void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

void	endJob () override

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endRun (edm::Run const &, edm::EventSetup const &) 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
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 = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

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

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

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

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 100 of file EGammaCutBasedEleIdAnalyzer.cc.

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

References beamSpotToken_, conversionsToken_, EAtargetToken_, electronsToken_, 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, GlobalPosition_Frontier_DevDB_cff::tag, and edm::vector_transform().

◆ ~EGammaCutBasedEleIdAnalyzer()

EGammaCutBasedEleIdAnalyzer::~EGammaCutBasedEleIdAnalyzer ( )

override

Definition at line 128 of file EGammaCutBasedEleIdAnalyzer.cc.

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

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 138 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::EDAnalyzer.

Definition at line 240 of file EGammaCutBasedEleIdAnalyzer.cc.

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

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

◆ beginLuminosityBlock()

void EGammaCutBasedEleIdAnalyzer::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 265 of file EGammaCutBasedEleIdAnalyzer.cc.

265 {}

◆ beginRun()

void EGammaCutBasedEleIdAnalyzer::beginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 259 of file EGammaCutBasedEleIdAnalyzer.cc.

259 {}

◆ endJob()

void EGammaCutBasedEleIdAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 256 of file EGammaCutBasedEleIdAnalyzer.cc.

256 {}

◆ endLuminosityBlock()

void EGammaCutBasedEleIdAnalyzer::endLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 268 of file EGammaCutBasedEleIdAnalyzer.cc.

268 {}

◆ endRun()

void EGammaCutBasedEleIdAnalyzer::endRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 262 of file EGammaCutBasedEleIdAnalyzer.cc.

262 {}

◆ fillDescriptions()

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

static

Definition at line 271 of file EGammaCutBasedEleIdAnalyzer.cc.

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