EGammaCutBasedEleIdAnalyzer Class Reference

#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::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)

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

References pwdgSkimBPark_cfi::beamSpot, beamSpotToken_, conversionsToken_, EAtarget, electronsToken_, h1_pt_, h1_pt_fbremeopin_, h1_pt_loose_, h1_pt_medium_, h1_pt_tight_, h1_pt_trig_, h1_pt_veto_, mps_fire::i, iEvent, isoValTokens_, dqmiolumiharvest::j, EgammaCutBasedEleId::LOOSE, EgammaCutBasedEleId::MEDIUM, dqmiodumpmetadata::n, EgammaCutBasedEleId::PassEoverPCuts(), EgammaCutBasedEleId::PassTriggerCuts(), EgammaCutBasedEleId::PassWP(), primaryVertexToken_, printDebug_, edm::Handle< T >::product(), electronIsolatorFromEffectiveArea_cfi::rhoIso, rhoIsoToken_, EgammaCutBasedEleId::TIGHT, EgammaCutBasedEleId::TRIGGERTIGHT, EgammaCutBasedEleId::TRIGGERWP70, EgammaCutBasedEleId::VETO, and PbPb_ZMuSkimMuonDPG_cff::veto.

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.

{}

beginRun()

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

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 259 of file EGammaCutBasedEleIdAnalyzer.cc.

{}

endJob()

void EGammaCutBasedEleIdAnalyzer::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 256 of file EGammaCutBasedEleIdAnalyzer.cc.

{}

endLuminosityBlock()

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

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 268 of file EGammaCutBasedEleIdAnalyzer.cc.

{}

endRun()

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

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 262 of file EGammaCutBasedEleIdAnalyzer.cc.

{}

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

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

Member Data Documentation

beamSpotToken_

edm::EDGetTokenT<reco::BeamSpot> EGammaCutBasedEleIdAnalyzer::beamSpotToken_

private

Definition at line 74 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

conversionsToken_

edm::EDGetTokenT<reco::ConversionCollection> EGammaCutBasedEleIdAnalyzer::conversionsToken_

private

Definition at line 73 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

EAtarget

ElectronEffectiveArea::ElectronEffectiveAreaTarget EGammaCutBasedEleIdAnalyzer::EAtarget

Definition at line 57 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and beginJob().

EAtargetToken_

std::string EGammaCutBasedEleIdAnalyzer::EAtargetToken_

private

Definition at line 78 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by beginJob(), and EGammaCutBasedEleIdAnalyzer().

electronsToken_

edm::EDGetTokenT<reco::GsfElectronCollection> EGammaCutBasedEleIdAnalyzer::electronsToken_

private

Definition at line 72 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

h1_pt_

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_

private

Definition at line 84 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

h1_pt_fbremeopin_

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_fbremeopin_

private

Definition at line 90 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

h1_pt_loose_

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_loose_

private

Definition at line 86 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

h1_pt_medium_

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_medium_

private

Definition at line 87 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

h1_pt_tight_

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_tight_

private

Definition at line 88 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

h1_pt_trig_

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_trig_

private

Definition at line 89 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

h1_pt_veto_

TH1F* EGammaCutBasedEleIdAnalyzer::h1_pt_veto_

private

Definition at line 85 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

isoValTokens_

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

private

Definition at line 77 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

primaryVertexToken_

edm::EDGetTokenT<reco::VertexCollection> EGammaCutBasedEleIdAnalyzer::primaryVertexToken_

private

Definition at line 76 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

printDebug_

bool EGammaCutBasedEleIdAnalyzer::printDebug_

private

Definition at line 81 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().

rhoIsoToken_

edm::EDGetTokenT<double> EGammaCutBasedEleIdAnalyzer::rhoIsoToken_

private

Definition at line 75 of file EGammaCutBasedEleIdAnalyzer.cc.

Referenced by analyze(), and EGammaCutBasedEleIdAnalyzer().