CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Member Functions | Private Attributes
EfficiencyPlotter Class Reference

#include <EfficiencyPlotter.h>

Inheritance diagram for EfficiencyPlotter:
DQMEDHarvester edm::one::EDProducer< edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::EndProcessBlockProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources, edm::Accumulator > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 EfficiencyPlotter (const edm::ParameterSet &ps)
 Constructor. More...
 
 ~EfficiencyPlotter () override
 Destructor. More...
 
- Public Member Functions inherited from DQMEDHarvester
void accumulate (edm::Event const &ev, edm::EventSetup const &es) final
 
void beginJob () override
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &) final
 
void beginRun (edm::Run const &run, edm::EventSetup const &) override
 
virtual void dqmAnalyze (DQMStore::IBooker &, DQMStore::IGetter &, edm::Event const &, edm::EventSetup const &)
 
 DQMEDHarvester ()
 
 DQMEDHarvester (edm::ParameterSet const &iConfig)
 
virtual void dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void dqmEndRun (DQMStore::IBooker &, DQMStore::IGetter &, edm::Run const &, edm::EventSetup const &)
 
void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
 
void endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &es) final
 
void endProcessBlockProduce (edm::ProcessBlock &) final
 
void endRun (edm::Run const &, edm::EventSetup const &) override
 
void endRunProduce (edm::Run &run, edm::EventSetup const &es) final
 
 ~DQMEDHarvester () override=default
 
- Public Member Functions inherited from edm::one::EDProducer< edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::EndProcessBlockProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns, edm::one::SharedResources, edm::Accumulator >
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProduceroperator= (const EDProducer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDProducerBase () 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< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Protected Member Functions

void dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
 DQM Client Diagnostic. More...
 
- 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 = 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)
 

Private Attributes

int etaBin
 
double etaMax
 
double etaMin
 
MonitorElementh_eff_eta_ID
 
MonitorElementh_eff_hp_eta_ID
 
MonitorElementh_eff_inner_eta_ID
 
MonitorElementh_eff_inner_phi_ID
 
MonitorElementh_eff_inner_pt_ID
 
MonitorElementh_eff_phi_ID
 
MonitorElementh_eff_pt_detIsoID
 
MonitorElementh_eff_pt_EB_detIsoID
 
MonitorElementh_eff_pt_EB_ID
 
MonitorElementh_eff_pt_EB_pfIsodBID
 
MonitorElementh_eff_pt_EB_pfIsoID
 
MonitorElementh_eff_pt_EE_detIsoID
 
MonitorElementh_eff_pt_EE_ID
 
MonitorElementh_eff_pt_EE_pfIsodBID
 
MonitorElementh_eff_pt_EE_pfIsoID
 
MonitorElementh_eff_pt_ID
 
MonitorElementh_eff_pt_pfIsodBID
 
MonitorElementh_eff_pt_pfIsoID
 
MonitorElementh_eff_vtx_detIsoID
 
MonitorElementh_eff_vtx_EB_detIsoID
 
MonitorElementh_eff_vtx_EB_pfIsodBID
 
MonitorElementh_eff_vtx_EB_pfIsoID
 
MonitorElementh_eff_vtx_EE_detIsoID
 
MonitorElementh_eff_vtx_EE_pfIsodBID
 
MonitorElementh_eff_vtx_EE_pfIsoID
 
MonitorElementh_eff_vtx_pfIsodBID
 
MonitorElementh_eff_vtx_pfIsoID
 
std::string ID_
 
std::string metname
 
int nevents
 
unsigned int nLumiSegs
 
edm::ParameterSet parameters
 
int phiBin
 
double phiMax
 
double phiMin
 
int prescaleFactor
 
int ptBin
 
double ptMax
 
double ptMin
 
int run
 
std::string theFolder
 
int vtxBin
 
double vtxMax
 
double vtxMin
 

Additional Inherited Members

- Public Types inherited from DQMEDHarvester
typedef dqm::harvesting::DQMStore DQMStore
 
typedef dqm::harvesting::MonitorElement MonitorElement
 
- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase 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::one::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Attributes inherited from DQMEDHarvester
DQMStoredqmstore_
 
edm::GetterOfProducts< DQMTokenjobmegetter_
 
edm::EDPutTokenT< DQMTokenjobToken_
 
edm::GetterOfProducts< DQMTokenlumimegetter_
 
edm::EDPutTokenT< DQMTokenlumiToken_
 
edm::GetterOfProducts< DQMTokenrunmegetter_
 
edm::EDPutTokenT< DQMTokenrunToken_
 

Detailed Description

Definition at line 24 of file EfficiencyPlotter.h.

Constructor & Destructor Documentation

◆ EfficiencyPlotter()

EfficiencyPlotter::EfficiencyPlotter ( const edm::ParameterSet ps)

Constructor.

Definition at line 27 of file EfficiencyPlotter.cc.

27  {
28 #ifdef DEBUG
29  cout << "EfficiencyPlotter(): Constructor " << endl;
30 #endif
31  parameters = ps;
32 
33  etaBin = parameters.getParameter<int>("etaBin");
34  etaMin = parameters.getParameter<double>("etaMin");
35  etaMax = parameters.getParameter<double>("etaMax");
36 
37  phiBin = parameters.getParameter<int>("phiBin");
38  phiMin = parameters.getParameter<double>("phiMin");
39  phiMax = parameters.getParameter<double>("phiMax");
40 
41  ptBin = parameters.getParameter<int>("ptBin");
42  ptMin = parameters.getParameter<double>("ptMin");
43  ptMax = parameters.getParameter<double>("ptMax");
44 
45  vtxBin = parameters.getParameter<int>("vtxBin");
46  vtxMin = parameters.getParameter<double>("vtxMin");
47  vtxMax = parameters.getParameter<double>("vtxMax");
48 
49  ID_ = parameters.getParameter<string>("MuonID");
50  theFolder = parameters.getParameter<string>("folder");
51 }

References gather_cfg::cout, etaBin(), ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, BeamMonitor_cff::phiBin, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, muonRecoAnalyzer_cfi::ptBin, AlignmentTrackSelector_cfi::ptMax, ptMin, EfficencyPlotter_cfi::vtxBin, EfficencyPlotter_cfi::vtxMax, and EfficencyPlotter_cfi::vtxMin.

◆ ~EfficiencyPlotter()

EfficiencyPlotter::~EfficiencyPlotter ( )
override

Destructor.

Definition at line 52 of file EfficiencyPlotter.cc.

52 {}

Member Function Documentation

◆ dqmEndJob()

void EfficiencyPlotter::dqmEndJob ( DQMStore::IBooker ibooker,
DQMStore::IGetter igetter 
)
overrideprotectedvirtual

DQM Client Diagnostic.

— Tight Muon efficiency vs muon Pt

— Tight Muon efficiency vs muon inner Pt

— Tight Muon efficiency vs muon inner eta

— Tight Muon efficiency vs muon inner phi

— Tight Muon efficiency vs muon Pt [EB]

— Tight Muon efficiency vs muon Pt [EE]

— Tight Muon efficiency vs muon Eta

— Tight Muon efficiency vs muon Eta [ pt > 20 ]

— Tight Muon efficiency vs muon Phi

— Tight Muon + Detector Isolation efficiency vs muon Pt

— Tight Muon + Detector Isolation efficiency vs muon Pt [EB]

— Tight Muon + Detector Isolation efficiency vs muon Pt [EE]

— Tight Muon + PF Isolation efficiency vs muon Pt

— Tight Muon + PF Isolation efficiency vs muon Pt [EB]

— Tight Muon + PF Isolation efficiency vs muon Pt [EE]

— Tight Muon + PF Isolation efficiency vs muon number of reco Vertex

— Tight Muon + detector-Based Isolation efficiency vs muon number of reco Vertex

— Tight Muon + detector-Based Isolation efficiency vs muon number of reco Vertex [EB]

— Tight Muon + detector-Based Isolation efficiency vs muon number of reco Vertex [EE]

— Tight Muon + PF-Based Isolation efficiency vs muon number of reco Vertex [EB]

— Tight Muon + PF-Based Isolation efficiency vs muon number of reco Vertex [EE]

— Tight Muon + PF IsodBlation corrected by deltaBeta efficiency vs muon Pt

— Tight Muon + PF IsodBlation corrected by deltaBeta efficiency vs muon Pt [EB]

— Tight Muon + PF IsodBlation corrected by deltaBeta efficiency vs muon Pt [EE]

— Tight Muon + PF Isolation corrected by deltaBeta efficiency vs muon number of reco Vertex

— Tight Muon + PF-Based Isolation corrected by deltaBeta efficiency vs muon number of reco Vertex [EB]

— Tight Muon + PF-Based Isolation corrected by deltaBeta efficiency vs muon number of reco Vertex [EE]

Implements DQMEDHarvester.

Definition at line 54 of file EfficiencyPlotter.cc.

54  {
55  ibooker.setCurrentFolder(theFolder);
56 
57  // efficiency plot
58  h_eff_eta_ID = ibooker.book1D("Eff_eta_" + ID_, ID_ + " Eff. vs #eta", etaBin, etaMin, etaMax);
59  h_eff_inner_pt_ID = ibooker.book1D("Eff_inner_pt_" + ID_, ID_ + " Eff. vs pt", ptBin, ptMin, ptMax);
60  h_eff_inner_eta_ID = ibooker.book1D("Eff_inner_eta_" + ID_, ID_ + " Eff. vs #eta", etaBin, etaMin, etaMax);
61  h_eff_inner_phi_ID = ibooker.book1D("Eff_inner_phi_" + ID_, ID_ + " Eff. vs #phi", phiBin, phiMin, phiMax);
62  h_eff_hp_eta_ID = ibooker.book1D("Eff_hp_eta_" + ID_, "High Pt " + ID_ + " Eff. vs #eta", etaBin, etaMin, etaMax);
63  h_eff_phi_ID = ibooker.book1D("Eff_phi_" + ID_, ID_ + " Eff. vs #phi", phiBin, phiMin, phiMax);
64  h_eff_pt_ID = ibooker.book1D("Eff_pt_" + ID_, ID_ + " Eff. vs Pt", ptBin, ptMin, ptMax);
65  h_eff_pt_EB_ID = ibooker.book1D("Eff_pt_EB_" + ID_, "Barrel: " + ID_ + " Eff. vs Pt", ptBin, ptMin, ptMax);
66  h_eff_pt_EE_ID = ibooker.book1D("Eff_pt_EE_" + ID_, "Endcap: " + ID_ + " Eff. vs Pt", ptBin, ptMin, ptMax);
67  h_eff_pt_detIsoID = ibooker.book1D("Eff_pt_detIso" + ID_, "detIso" + ID_ + " Efficiency vs Pt", ptBin, ptMin, ptMax);
69  ibooker.book1D("Eff_pt_EB_detIso" + ID_, "Barrel: detIso" + ID_ + " Eff. vs Pt", ptBin, ptMin, ptMax);
71  ibooker.book1D("Eff_pt_EE_detIso" + ID_, "Endcap: detIso" + ID_ + " Eff. vs Pt", ptBin, ptMin, ptMax);
72  h_eff_pt_pfIsoID = ibooker.book1D("Eff_pt_pfIso" + ID_, "pfIso" + ID_ + " Eff. vs Pt", ptBin, ptMin, ptMax);
74  ibooker.book1D("Eff_pt_EB_pfIso" + ID_, "Barrel: pfIso" + ID_ + " Eff. vs Pt", ptBin, ptMin, ptMax);
76  ibooker.book1D("Eff_pt_EE_pfIso" + ID_, "Endcap: pfIso" + ID_ + " Eff. vs Pt", ptBin, ptMin, ptMax);
77  h_eff_vtx_detIsoID = ibooker.book1D("Eff_vtx_detIso" + ID_, "detIso" + ID_ + " Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
78  h_eff_vtx_pfIsoID = ibooker.book1D("Eff_vtx_pfIso" + ID_, "pfIso" + ID_ + " Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
80  ibooker.book1D("Eff_vtx_EB_detIso" + ID_, "Barrel: detIso" + ID_ + " Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
82  ibooker.book1D("Eff_vtx_EB_pfIso" + ID_, "Barrel: pfIso" + ID_ + " Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
84  ibooker.book1D("Eff_vtx_EE_detIso" + ID_, "Endcap: detIso" + ID_ + " Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
86  ibooker.book1D("Eff_vtx_EE_pfIso" + ID_, "Endcap: pfIso" + ID_ + " Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
87 
89  ibooker.book1D("Eff_pt_pfIsodB" + ID_, "pfIso" + ID_ + " (deltaBeta) Eff. vs Pt", ptBin, ptMin, ptMax);
91  ibooker.book1D("Eff_pt_EB_pfIsodB" + ID_, "Barrel: pfIso" + ID_ + " (deltaBeta) Eff. vs Pt", ptBin, ptMin, ptMax);
93  ibooker.book1D("Eff_pt_EE_pfIsodB" + ID_, "Endcap: pfIso" + ID_ + " (deltaBeta) Eff. vs Pt", ptBin, ptMin, ptMax);
95  ibooker.book1D("Eff_vtx_pfIsodB" + ID_, "pfIso" + ID_ + " (deltaBeta) Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
97  "Eff_vtx_EB_pfIsodB" + ID_, "Barrel: pfIso" + ID_ + " (deltaBeta) Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
99  "Eff_vtx_EE_pfIsodB" + ID_, "Endcap: pfIso" + ID_ + " (deltaBeta) Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
100 
101  // This prevents this ME to be normalized when drawn into the GUI
123 
130 
131  // AXIS TITLES....
132  h_eff_hp_eta_ID->setAxisTitle("#eta", 1);
133  h_eff_inner_eta_ID->setAxisTitle("#eta", 1);
134  h_eff_inner_pt_ID->setAxisTitle("#eta", 1);
135  h_eff_inner_phi_ID->setAxisTitle("#eta", 1);
136  h_eff_eta_ID->setAxisTitle("#eta", 1);
137  h_eff_phi_ID->setAxisTitle("#phi", 1);
138  h_eff_pt_ID->setAxisTitle("p_{T} (GeV)", 1);
139  h_eff_pt_EB_ID->setAxisTitle("p_{T} (GeV)", 1);
140  h_eff_pt_EE_ID->setAxisTitle("p_{T} (GeV)", 1);
141  h_eff_pt_detIsoID->setAxisTitle("p_{T} (GeV)", 1);
142  h_eff_pt_EB_detIsoID->setAxisTitle("p_{T} (GeV)", 1);
143  h_eff_pt_EE_detIsoID->setAxisTitle("p_{T} (GeV)", 1);
144  h_eff_pt_pfIsoID->setAxisTitle("p_{T} (GeV)", 1);
145  h_eff_pt_EB_pfIsoID->setAxisTitle("p_{T} (GeV)", 1);
146  h_eff_pt_EE_pfIsoID->setAxisTitle("p_{T} (GeV)", 1);
147  h_eff_vtx_detIsoID->setAxisTitle("Number of PV", 1);
148  h_eff_vtx_pfIsoID->setAxisTitle("Number of PV", 1);
149  h_eff_vtx_EB_detIsoID->setAxisTitle("Number of PV", 1);
150  h_eff_vtx_EB_pfIsoID->setAxisTitle("Number of PV", 1);
151  h_eff_vtx_EE_detIsoID->setAxisTitle("Number of PV", 1);
152  h_eff_vtx_EE_pfIsoID->setAxisTitle("Number of PV", 1);
153 
154  h_eff_pt_pfIsodBID->setAxisTitle("p_{T} (GeV)", 1);
155  h_eff_pt_EB_pfIsodBID->setAxisTitle("p_{T} (GeV)", 1);
156  h_eff_pt_EE_pfIsodBID->setAxisTitle("p_{T} (GeV)", 1);
157  h_eff_vtx_pfIsodBID->setAxisTitle("Number of PV", 1);
158  h_eff_vtx_EB_pfIsodBID->setAxisTitle("Number of PV", 1);
159  h_eff_vtx_EE_pfIsodBID->setAxisTitle("Number of PV", 1);
160 
162  string inputdir = "Muons/EfficiencyAnalyzer/" + ID_;
163  string numpath_pt = inputdir + "/passProbes_ID_pt";
164  string denpath_pt = inputdir + "/allProbes_pt";
165 
166  MonitorElement *Numerator_pt = igetter.get(numpath_pt);
167  MonitorElement *Denominator_pt = igetter.get(denpath_pt);
168 
169  if (Numerator_pt && Denominator_pt) {
170  TH1F *h_numerator_pt = Numerator_pt->getTH1F();
171  TH1F *h_denominator_pt = Denominator_pt->getTH1F();
172  TH1F *h_eff_pt = h_eff_pt_ID->getTH1F();
173 
174  if (h_eff_pt->GetSumw2N() == 0)
175  h_eff_pt->Sumw2();
176  h_eff_pt->Divide(h_numerator_pt, h_denominator_pt, 1., 1., "B");
177  }
178 
180  string numpath_inner_pt = inputdir + "/passProbes_ID_inner_pt";
181  string denpath_inner_pt = inputdir + "/allProbes_inner_pt";
182 
183  MonitorElement *Numerator_inner_pt = igetter.get(numpath_inner_pt);
184  MonitorElement *Denominator_inner_pt = igetter.get(denpath_inner_pt);
185 
186  if (Numerator_inner_pt && Denominator_inner_pt) {
187  TH1F *h_numerator_inner_pt = Numerator_inner_pt->getTH1F();
188  TH1F *h_denominator_inner_pt = Denominator_inner_pt->getTH1F();
189  TH1F *h_eff_inner_pt = h_eff_inner_pt_ID->getTH1F();
190 
191  if (h_eff_inner_pt->GetSumw2N() == 0)
192  h_eff_inner_pt->Sumw2();
193  h_eff_inner_pt->Divide(h_numerator_inner_pt, h_denominator_inner_pt, 1., 1., "B");
194  }
195 
197  string numpath_inner_eta = inputdir + "/passProbes_ID_inner_eta";
198  string denpath_inner_eta = inputdir + "/allProbes_inner_eta";
199 
200  MonitorElement *Numerator_inner_eta = igetter.get(numpath_inner_eta);
201  MonitorElement *Denominator_inner_eta = igetter.get(denpath_inner_eta);
202 
203  if (Numerator_inner_eta && Denominator_inner_eta) {
204  TH1F *h_numerator_inner_eta = Numerator_inner_eta->getTH1F();
205  TH1F *h_denominator_inner_eta = Denominator_inner_eta->getTH1F();
206  TH1F *h_eff_inner_eta = h_eff_inner_eta_ID->getTH1F();
207 
208  if (h_eff_inner_eta->GetSumw2N() == 0)
209  h_eff_inner_eta->Sumw2();
210  h_eff_inner_eta->Divide(h_numerator_inner_eta, h_denominator_inner_eta, 1., 1., "B");
211  }
212 
214  string numpath_inner_phi = inputdir + "/passProbes_ID_inner_phi";
215  string denpath_inner_phi = inputdir + "/allProbes_inner_phi";
216 
217  MonitorElement *Numerator_inner_phi = igetter.get(numpath_inner_phi);
218  MonitorElement *Denominator_inner_phi = igetter.get(denpath_inner_phi);
219 
220  if (Numerator_inner_phi && Denominator_inner_phi) {
221  TH1F *h_numerator_inner_phi = Numerator_inner_phi->getTH1F();
222  TH1F *h_denominator_inner_phi = Denominator_inner_phi->getTH1F();
223  TH1F *h_eff_inner_phi = h_eff_inner_phi_ID->getTH1F();
224 
225  if (h_eff_inner_phi->GetSumw2N() == 0)
226  h_eff_inner_phi->Sumw2();
227  h_eff_inner_phi->Divide(h_numerator_inner_phi, h_denominator_inner_phi, 1., 1., "B");
228  }
229 
231  string numpath_EB_pt = inputdir + "/passProbes_ID_EB_pt";
232  string denpath_EB_pt = inputdir + "/allProbes_EB_pt";
233 
234  MonitorElement *Numerator_EB_pt = igetter.get(numpath_EB_pt);
235  MonitorElement *Denominator_EB_pt = igetter.get(denpath_EB_pt);
236 
237  if (Numerator_EB_pt && Denominator_EB_pt) {
238  TH1F *h_numerator_EB_pt = Numerator_EB_pt->getTH1F();
239  TH1F *h_denominator_EB_pt = Denominator_EB_pt->getTH1F();
240  TH1F *h_eff_EB_pt = h_eff_pt_EB_ID->getTH1F();
241 
242  if (h_eff_EB_pt->GetSumw2N() == 0)
243  h_eff_EB_pt->Sumw2();
244  h_eff_EB_pt->Divide(h_numerator_EB_pt, h_denominator_EB_pt, 1., 1., "B");
245  }
246 
248  string numpath_EE_pt = inputdir + "/passProbes_ID_EE_pt";
249  string denpath_EE_pt = inputdir + "/allProbes_EE_pt";
250 
251  MonitorElement *Numerator_EE_pt = igetter.get(numpath_EE_pt);
252  MonitorElement *Denominator_EE_pt = igetter.get(denpath_EE_pt);
253 
254  if (Numerator_EE_pt && Denominator_EE_pt) {
255  TH1F *h_numerator_EE_pt = Numerator_EE_pt->getTH1F();
256  TH1F *h_denominator_EE_pt = Denominator_EE_pt->getTH1F();
257  TH1F *h_eff_EE_pt = h_eff_pt_EE_ID->getTH1F();
258 
259  if (h_eff_EE_pt->GetSumw2N() == 0)
260  h_eff_EE_pt->Sumw2();
261  h_eff_EE_pt->Divide(h_numerator_EE_pt, h_denominator_EE_pt, 1., 1., "B");
262  }
263 
265  string numpath_eta = inputdir + "/passProbes_ID_eta";
266  string denpath_eta = inputdir + "/allProbes_eta";
267 
268  MonitorElement *Numerator_eta = igetter.get(numpath_eta);
269  MonitorElement *Denominator_eta = igetter.get(denpath_eta);
270 
271  if (Numerator_eta && Denominator_eta) {
272  TH1F *h_numerator_eta = Numerator_eta->getTH1F();
273  TH1F *h_denominator_eta = Denominator_eta->getTH1F();
274 
275  TH1F *h_eff_eta = h_eff_eta_ID->getTH1F();
276 
277  if (h_eff_eta->GetSumw2N() == 0)
278  h_eff_eta->Sumw2();
279 
280  h_eff_eta->Divide(h_numerator_eta, h_denominator_eta, 1., 1., "B");
281  }
282 
284 
285  string numpath_hp_eta = inputdir + "/passProbes_ID_hp_eta";
286  string denpath_hp_eta = inputdir + "/allProbes_hp_eta";
287 
288  MonitorElement *Numerator_hp_eta = igetter.get(numpath_hp_eta);
289  MonitorElement *Denominator_hp_eta = igetter.get(denpath_hp_eta);
290 
291  if (Numerator_hp_eta && Denominator_hp_eta) {
292  TH1F *h_numerator_hp_eta = Numerator_hp_eta->getTH1F();
293  TH1F *h_denominator_hp_eta = Denominator_hp_eta->getTH1F();
294 
295  TH1F *h_eff_hp_eta = h_eff_hp_eta_ID->getTH1F();
296 
297  if (h_eff_hp_eta->GetSumw2N() == 0)
298  h_eff_hp_eta->Sumw2();
299 
300  h_eff_hp_eta->Divide(h_numerator_hp_eta, h_denominator_hp_eta, 1., 1., "B");
301  }
302 
304 
305  string numpath_phi = inputdir + "/passProbes_ID_phi";
306  string denpath_phi = inputdir + "/allProbes_phi";
307 
308  MonitorElement *Numerator_phi = igetter.get(numpath_phi);
309  MonitorElement *Denominator_phi = igetter.get(denpath_phi);
310 
311  if (Numerator_phi && Denominator_phi) {
312  TH1F *h_numerator_phi = Numerator_phi->getTH1F();
313  TH1F *h_denominator_phi = Denominator_phi->getTH1F();
314 
315  TH1F *h_eff_phi = h_eff_phi_ID->getTH1F();
316 
317  if (h_eff_phi->GetSumw2N() == 0)
318  h_eff_phi->Sumw2();
319 
320  h_eff_phi->Divide(h_numerator_phi, h_denominator_phi, 1., 1., "B");
321  }
322 
324  string numpath_detIso_pt = inputdir + "/passProbes_detIsoID_pt";
325  string denpath_detIso_pt = inputdir + "/allProbes_ID_pt";
326 
327  MonitorElement *Numerator_detIso_pt = igetter.get(numpath_detIso_pt);
328  MonitorElement *Denominator_detIso_pt = igetter.get(denpath_detIso_pt);
329 
330  if (Numerator_detIso_pt && Denominator_detIso_pt) {
331  TH1F *h_numerator_detIso_pt = Numerator_detIso_pt->getTH1F();
332  TH1F *h_denominator_detIso_pt = Denominator_detIso_pt->getTH1F();
333 
334  TH1F *h_eff_detIso_pt = h_eff_pt_detIsoID->getTH1F();
335 
336  if (h_eff_detIso_pt->GetSumw2N() == 0)
337  h_eff_detIso_pt->Sumw2();
338  h_eff_detIso_pt->Divide(h_numerator_detIso_pt, h_denominator_detIso_pt, 1., 1., "B");
339  }
340 
342  string numpath_detIso_EB_pt = inputdir + "/passProbes_EB_detIsoID_pt";
343  string denpath_detIso_EB_pt = inputdir + "/allProbes_EB_ID_pt";
344 
345  MonitorElement *Numerator_detIso_EB_pt = igetter.get(numpath_detIso_EB_pt);
346  MonitorElement *Denominator_detIso_EB_pt = igetter.get(denpath_detIso_EB_pt);
347 
348  if (Numerator_detIso_EB_pt && Denominator_detIso_EB_pt) {
349  TH1F *h_numerator_detIso_EB_pt = Numerator_detIso_EB_pt->getTH1F();
350  TH1F *h_denominator_detIso_EB_pt = Denominator_detIso_EB_pt->getTH1F();
351 
352  TH1F *h_eff_detIso_EB_pt = h_eff_pt_EB_detIsoID->getTH1F();
353 
354  if (h_eff_detIso_EB_pt->GetSumw2N() == 0)
355  h_eff_detIso_EB_pt->Sumw2();
356 
357  h_eff_detIso_EB_pt->Divide(h_numerator_detIso_EB_pt, h_denominator_detIso_EB_pt, 1., 1., "B");
358  }
359 
361  string numpath_detIso_EE_pt = inputdir + "/passProbes_EE_detIsoID_pt";
362  string denpath_detIso_EE_pt = inputdir + "/allProbes_EE_ID_pt";
363 
364  MonitorElement *Numerator_detIso_EE_pt = igetter.get(numpath_detIso_EE_pt);
365  MonitorElement *Denominator_detIso_EE_pt = igetter.get(denpath_detIso_EE_pt);
366 
367  if (Numerator_detIso_EE_pt && Denominator_detIso_EE_pt) {
368  TH1F *h_numerator_detIso_EE_pt = Numerator_detIso_EE_pt->getTH1F();
369  TH1F *h_denominator_detIso_EE_pt = Denominator_detIso_EE_pt->getTH1F();
370 
371  TH1F *h_eff_detIso_EE_pt = h_eff_pt_EE_detIsoID->getTH1F();
372 
373  if (h_eff_detIso_EE_pt->GetSumw2N() == 0)
374  h_eff_detIso_EE_pt->Sumw2();
375 
376  h_eff_detIso_EE_pt->Divide(h_numerator_detIso_EE_pt, h_denominator_detIso_EE_pt, 1., 1., "B");
377  }
378 
380  string numpath_pfIso_pt = inputdir + "/passProbes_pfIsoID_pt";
381  string denpath_pfIso_pt = inputdir + "/allProbes_ID_pt";
382 
383  MonitorElement *Numerator_pfIso_pt = igetter.get(numpath_pfIso_pt);
384  MonitorElement *Denominator_pfIso_pt = igetter.get(denpath_pfIso_pt);
385 
386  if (Numerator_pfIso_pt && Denominator_pfIso_pt) {
387  TH1F *h_numerator_pfIso_pt = Numerator_pfIso_pt->getTH1F();
388  TH1F *h_denominator_pfIso_pt = Denominator_pfIso_pt->getTH1F();
389 
390  TH1F *h_eff_pfIso_pt = h_eff_pt_pfIsoID->getTH1F();
391 
392  if (h_eff_pfIso_pt->GetSumw2N() == 0)
393  h_eff_pfIso_pt->Sumw2();
394 
395  h_eff_pfIso_pt->Divide(h_numerator_pfIso_pt, h_denominator_pfIso_pt, 1., 1., "B");
396  }
397 
399 
400  string numpath_pfIso_EB_pt = inputdir + "/passProbes_EB_pfIsoID_pt";
401  string denpath_pfIso_EB_pt = inputdir + "/allProbes_EB_ID_pt";
402 
403  MonitorElement *Numerator_pfIso_EB_pt = igetter.get(numpath_pfIso_EB_pt);
404  MonitorElement *Denominator_pfIso_EB_pt = igetter.get(denpath_pfIso_EB_pt);
405 
406  if (Numerator_pfIso_EB_pt && Denominator_pfIso_EB_pt) {
407  TH1F *h_numerator_pfIso_EB_pt = Numerator_pfIso_EB_pt->getTH1F();
408  TH1F *h_denominator_pfIso_EB_pt = Denominator_pfIso_EB_pt->getTH1F();
409 
410  TH1F *h_eff_pfIso_EB_pt = h_eff_pt_EB_pfIsoID->getTH1F();
411 
412  if (h_eff_pfIso_EB_pt->GetSumw2N() == 0)
413  h_eff_pfIso_EB_pt->Sumw2();
414  h_eff_pfIso_EB_pt->Divide(h_numerator_pfIso_EB_pt, h_denominator_pfIso_EB_pt, 1., 1., "B");
415  }
416 
418  string numpath_pfIso_EE_pt = inputdir + "/passProbes_EE_pfIsoID_pt";
419  string denpath_pfIso_EE_pt = inputdir + "/allProbes_EE_ID_pt";
420 
421  MonitorElement *Numerator_pfIso_EE_pt = igetter.get(numpath_pfIso_EE_pt);
422  MonitorElement *Denominator_pfIso_EE_pt = igetter.get(denpath_pfIso_EE_pt);
423 
424  if (Numerator_pfIso_EE_pt && Denominator_pfIso_EE_pt) {
425  TH1F *h_numerator_pfIso_EE_pt = Numerator_pfIso_EE_pt->getTH1F();
426  TH1F *h_denominator_pfIso_EE_pt = Denominator_pfIso_EE_pt->getTH1F();
427 
428  TH1F *h_eff_pfIso_EE_pt = h_eff_pt_EE_pfIsoID->getTH1F();
429 
430  if (h_eff_pfIso_EE_pt->GetSumw2N() == 0)
431  h_eff_pfIso_EE_pt->Sumw2();
432 
433  h_eff_pfIso_EE_pt->Divide(h_numerator_pfIso_EE_pt, h_denominator_pfIso_EE_pt, 1., 1., "B");
434  }
435 
437 
438  string numpath_pfIso_nvtx = inputdir + "/passProbes_pfIsoID_nVtx";
439  string denpath_pfIso_nvtx = inputdir + "/allProbes_ID_nVtx";
440 
441  MonitorElement *Numerator_pfIso_nvtx = igetter.get(numpath_pfIso_nvtx);
442  MonitorElement *Denominator_pfIso_nvtx = igetter.get(denpath_pfIso_nvtx);
443 
444  if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx) {
445  TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->getTH1F();
446  TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->getTH1F();
447 
448  TH1F *h_eff_pfIso_nvtx = h_eff_vtx_pfIsoID->getTH1F();
449 
450  if (h_eff_pfIso_nvtx->GetSumw2N() == 0)
451  h_eff_pfIso_nvtx->Sumw2();
452 
453  h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1., "B");
454  }
455 
457  string numpath_detIso_nvtx = inputdir + "/passProbes_detIsoID_nVtx";
458  string denpath_detIso_nvtx = inputdir + "/allProbes_ID_nVtx";
459 
460  MonitorElement *Numerator_detIso_nvtx = igetter.get(numpath_detIso_nvtx);
461  MonitorElement *Denominator_detIso_nvtx = igetter.get(denpath_detIso_nvtx);
462 
463  if (Numerator_detIso_nvtx && Denominator_detIso_nvtx) {
464  TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->getTH1F();
465  TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->getTH1F();
466 
467  TH1F *h_eff_detIso_nvtx = h_eff_vtx_detIsoID->getTH1F();
468 
469  if (h_eff_detIso_nvtx->GetSumw2N() == 0)
470  h_eff_detIso_nvtx->Sumw2();
471 
472  h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1., "B");
473  }
474 
476  numpath_detIso_nvtx = inputdir + "/passProbes_EB_detIsoID_nVtx";
477  denpath_detIso_nvtx = inputdir + "/allProbes_EB_ID_nVtx";
478 
479  Numerator_detIso_nvtx = igetter.get(numpath_detIso_nvtx);
480  Denominator_detIso_nvtx = igetter.get(denpath_detIso_nvtx);
481 
482  if (Numerator_detIso_nvtx && Denominator_detIso_nvtx) {
483  TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->getTH1F();
484  TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->getTH1F();
485 
486  TH1F *h_eff_detIso_nvtx = h_eff_vtx_EB_detIsoID->getTH1F();
487 
488  if (h_eff_detIso_nvtx->GetSumw2N() == 0)
489  h_eff_detIso_nvtx->Sumw2();
490  h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1., "B");
491  }
492 
494  numpath_detIso_nvtx = inputdir + "/passProbes_EE_detIsoID_nVtx";
495  denpath_detIso_nvtx = inputdir + "/allProbes_EE_ID_nVtx";
496 
497  Numerator_detIso_nvtx = igetter.get(numpath_detIso_nvtx);
498  Denominator_detIso_nvtx = igetter.get(denpath_detIso_nvtx);
499 
500  if (Numerator_detIso_nvtx && Denominator_detIso_nvtx) {
501  TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->getTH1F();
502  TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->getTH1F();
503 
504  TH1F *h_eff_detIso_nvtx = h_eff_vtx_EE_detIsoID->getTH1F();
505 
506  if (h_eff_detIso_nvtx->GetSumw2N() == 0)
507  h_eff_detIso_nvtx->Sumw2();
508 
509  h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1., "B");
510  }
511 
513  numpath_pfIso_nvtx = inputdir + "/passProbes_EB_pfIsoID_nVtx";
514  denpath_pfIso_nvtx = inputdir + "/allProbes_EB_ID_nVtx";
515 
516  Numerator_pfIso_nvtx = igetter.get(numpath_pfIso_nvtx);
517  Denominator_pfIso_nvtx = igetter.get(denpath_pfIso_nvtx);
518 
519  if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx) {
520  TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->getTH1F();
521  TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->getTH1F();
522 
523  TH1F *h_eff_pfIso_nvtx = h_eff_vtx_EB_pfIsoID->getTH1F();
524 
525  if (h_eff_pfIso_nvtx->GetSumw2N() == 0)
526  h_eff_pfIso_nvtx->Sumw2();
527 
528  h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1., "B");
529  }
530 
532  numpath_pfIso_nvtx = inputdir + "/passProbes_EE_pfIsoID_nVtx";
533  denpath_pfIso_nvtx = inputdir + "/allProbes_EE_ID_nVtx";
534 
535  Numerator_pfIso_nvtx = igetter.get(numpath_pfIso_nvtx);
536  Denominator_pfIso_nvtx = igetter.get(denpath_pfIso_nvtx);
537 
538  if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx) {
539  TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->getTH1F();
540  TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->getTH1F();
541 
542  TH1F *h_eff_pfIso_nvtx = h_eff_vtx_EE_pfIsoID->getTH1F();
543 
544  if (h_eff_pfIso_nvtx->GetSumw2N() == 0)
545  h_eff_pfIso_nvtx->Sumw2();
546 
547  h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1., "B");
548  }
549 
551  string numpath_pfIsodB_pt = inputdir + "/passProbes_pfIsodBID_pt";
552  string denpath_pfIsodB_pt = inputdir + "/allProbes_ID_pt";
553 
554  MonitorElement *Numerator_pfIsodB_pt = igetter.get(numpath_pfIsodB_pt);
555  MonitorElement *Denominator_pfIsodB_pt = igetter.get(denpath_pfIsodB_pt);
556 
557  if (Numerator_pfIsodB_pt && Denominator_pfIsodB_pt) {
558  TH1F *h_numerator_pfIsodB_pt = Numerator_pfIsodB_pt->getTH1F();
559  TH1F *h_denominator_pfIsodB_pt = Denominator_pfIsodB_pt->getTH1F();
560 
561  TH1F *h_eff_pfIsodB_pt = h_eff_pt_pfIsodBID->getTH1F();
562 
563  if (h_eff_pfIsodB_pt->GetSumw2N() == 0)
564  h_eff_pfIsodB_pt->Sumw2();
565 
566  h_eff_pfIsodB_pt->Divide(h_numerator_pfIsodB_pt, h_denominator_pfIsodB_pt, 1., 1., "B");
567  }
568 
570 
571  string numpath_pfIsodB_EB_pt = inputdir + "/passProbes_EB_pfIsodBID_pt";
572  string denpath_pfIsodB_EB_pt = inputdir + "/allProbes_EB_ID_pt";
573 
574  MonitorElement *Numerator_pfIsodB_EB_pt = igetter.get(numpath_pfIsodB_EB_pt);
575  MonitorElement *Denominator_pfIsodB_EB_pt = igetter.get(denpath_pfIsodB_EB_pt);
576 
577  if (Numerator_pfIsodB_EB_pt && Denominator_pfIsodB_EB_pt) {
578  TH1F *h_numerator_pfIsodB_EB_pt = Numerator_pfIsodB_EB_pt->getTH1F();
579  TH1F *h_denominator_pfIsodB_EB_pt = Denominator_pfIsodB_EB_pt->getTH1F();
580 
581  TH1F *h_eff_pfIsodB_EB_pt = h_eff_pt_EB_pfIsodBID->getTH1F();
582 
583  if (h_eff_pfIsodB_EB_pt->GetSumw2N() == 0)
584  h_eff_pfIsodB_EB_pt->Sumw2();
585  h_eff_pfIsodB_EB_pt->Divide(h_numerator_pfIsodB_EB_pt, h_denominator_pfIsodB_EB_pt, 1., 1., "B");
586  }
587 
589  string numpath_pfIsodB_EE_pt = inputdir + "/passProbes_EE_pfIsodBID_pt";
590  string denpath_pfIsodB_EE_pt = inputdir + "/allProbes_EE_ID_pt";
591 
592  MonitorElement *Numerator_pfIsodB_EE_pt = igetter.get(numpath_pfIsodB_EE_pt);
593  MonitorElement *Denominator_pfIsodB_EE_pt = igetter.get(denpath_pfIsodB_EE_pt);
594 
595  if (Numerator_pfIsodB_EE_pt && Denominator_pfIsodB_EE_pt) {
596  TH1F *h_numerator_pfIsodB_EE_pt = Numerator_pfIsodB_EE_pt->getTH1F();
597  TH1F *h_denominator_pfIsodB_EE_pt = Denominator_pfIsodB_EE_pt->getTH1F();
598 
599  TH1F *h_eff_pfIsodB_EE_pt = h_eff_pt_EE_pfIsodBID->getTH1F();
600 
601  if (h_eff_pfIsodB_EE_pt->GetSumw2N() == 0)
602  h_eff_pfIsodB_EE_pt->Sumw2();
603 
604  h_eff_pfIsodB_EE_pt->Divide(h_numerator_pfIsodB_EE_pt, h_denominator_pfIsodB_EE_pt, 1., 1., "B");
605  }
606 
608 
609  string numpath_pfIsodB_nvtx = inputdir + "/passProbes_pfIsodBID_nVtx";
610  string denpath_pfIsodB_nvtx = inputdir + "/allProbes_ID_nVtx";
611 
612  MonitorElement *Numerator_pfIsodB_nvtx = igetter.get(numpath_pfIsodB_nvtx);
613  MonitorElement *Denominator_pfIsodB_nvtx = igetter.get(denpath_pfIsodB_nvtx);
614 
615  if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx) {
616  TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->getTH1F();
617  TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->getTH1F();
618 
619  TH1F *h_eff_pfIsodB_nvtx = h_eff_vtx_pfIsodBID->getTH1F();
620 
621  if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0)
622  h_eff_pfIsodB_nvtx->Sumw2();
623 
624  h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1., "B");
625  }
626 
628  numpath_pfIsodB_nvtx = inputdir + "/passProbes_EB_pfIsodBID_nVtx";
629  denpath_pfIsodB_nvtx = inputdir + "/allProbes_EB_ID_nVtx";
630 
631  Numerator_pfIsodB_nvtx = igetter.get(numpath_pfIsodB_nvtx);
632  Denominator_pfIsodB_nvtx = igetter.get(denpath_pfIsodB_nvtx);
633 
634  if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx) {
635  TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->getTH1F();
636  TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->getTH1F();
637 
638  TH1F *h_eff_pfIsodB_nvtx = h_eff_vtx_EB_pfIsodBID->getTH1F();
639 
640  if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0)
641  h_eff_pfIsodB_nvtx->Sumw2();
642 
643  h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1., "B");
644  }
645 
647  numpath_pfIsodB_nvtx = inputdir + "/passProbes_EE_pfIsodBID_nVtx";
648  denpath_pfIsodB_nvtx = inputdir + "/allProbes_EE_ID_nVtx";
649 
650  Numerator_pfIsodB_nvtx = igetter.get(numpath_pfIsodB_nvtx);
651  Denominator_pfIsodB_nvtx = igetter.get(denpath_pfIsodB_nvtx);
652 
653  if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx) {
654  TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->getTH1F();
655  TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->getTH1F();
656 
657  TH1F *h_eff_pfIsodB_nvtx = h_eff_vtx_EE_pfIsodBID->getTH1F();
658 
659  if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0)
660  h_eff_pfIsodB_nvtx->Sumw2();
661 
662  h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1., "B");
663  }
664 }

References dqm::implementation::IBooker::book1D(), etaBin(), ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, dqm::implementation::IGetter::get(), dqm::legacy::MonitorElement::getTH1F(), BeamMonitor_cff::phiBin, AlignmentTrackSelector_cfi::phiMax, AlignmentTrackSelector_cfi::phiMin, muonRecoAnalyzer_cfi::ptBin, AlignmentTrackSelector_cfi::ptMax, ptMin, dqm::implementation::NavigatorBase::setCurrentFolder(), dqm::impl::MonitorElement::setEfficiencyFlag(), EfficencyPlotter_cfi::vtxBin, EfficencyPlotter_cfi::vtxMax, and EfficencyPlotter_cfi::vtxMin.

Member Data Documentation

◆ etaBin

int EfficiencyPlotter::etaBin
private

Definition at line 48 of file EfficiencyPlotter.h.

◆ etaMax

double EfficiencyPlotter::etaMax
private

Definition at line 50 of file EfficiencyPlotter.h.

◆ etaMin

double EfficiencyPlotter::etaMin
private

Definition at line 49 of file EfficiencyPlotter.h.

◆ h_eff_eta_ID

MonitorElement* EfficiencyPlotter::h_eff_eta_ID
private

Definition at line 73 of file EfficiencyPlotter.h.

◆ h_eff_hp_eta_ID

MonitorElement* EfficiencyPlotter::h_eff_hp_eta_ID
private

Definition at line 74 of file EfficiencyPlotter.h.

◆ h_eff_inner_eta_ID

MonitorElement* EfficiencyPlotter::h_eff_inner_eta_ID
private

Definition at line 69 of file EfficiencyPlotter.h.

◆ h_eff_inner_phi_ID

MonitorElement* EfficiencyPlotter::h_eff_inner_phi_ID
private

Definition at line 70 of file EfficiencyPlotter.h.

◆ h_eff_inner_pt_ID

MonitorElement* EfficiencyPlotter::h_eff_inner_pt_ID
private

Definition at line 68 of file EfficiencyPlotter.h.

◆ h_eff_phi_ID

MonitorElement* EfficiencyPlotter::h_eff_phi_ID
private

Definition at line 75 of file EfficiencyPlotter.h.

◆ h_eff_pt_detIsoID

MonitorElement* EfficiencyPlotter::h_eff_pt_detIsoID
private

Definition at line 76 of file EfficiencyPlotter.h.

◆ h_eff_pt_EB_detIsoID

MonitorElement* EfficiencyPlotter::h_eff_pt_EB_detIsoID
private

Definition at line 77 of file EfficiencyPlotter.h.

◆ h_eff_pt_EB_ID

MonitorElement* EfficiencyPlotter::h_eff_pt_EB_ID
private

Definition at line 71 of file EfficiencyPlotter.h.

◆ h_eff_pt_EB_pfIsodBID

MonitorElement* EfficiencyPlotter::h_eff_pt_EB_pfIsodBID
private

Definition at line 91 of file EfficiencyPlotter.h.

◆ h_eff_pt_EB_pfIsoID

MonitorElement* EfficiencyPlotter::h_eff_pt_EB_pfIsoID
private

Definition at line 80 of file EfficiencyPlotter.h.

◆ h_eff_pt_EE_detIsoID

MonitorElement* EfficiencyPlotter::h_eff_pt_EE_detIsoID
private

Definition at line 78 of file EfficiencyPlotter.h.

◆ h_eff_pt_EE_ID

MonitorElement* EfficiencyPlotter::h_eff_pt_EE_ID
private

Definition at line 72 of file EfficiencyPlotter.h.

◆ h_eff_pt_EE_pfIsodBID

MonitorElement* EfficiencyPlotter::h_eff_pt_EE_pfIsodBID
private

Definition at line 92 of file EfficiencyPlotter.h.

◆ h_eff_pt_EE_pfIsoID

MonitorElement* EfficiencyPlotter::h_eff_pt_EE_pfIsoID
private

Definition at line 81 of file EfficiencyPlotter.h.

◆ h_eff_pt_ID

MonitorElement* EfficiencyPlotter::h_eff_pt_ID
private

Definition at line 67 of file EfficiencyPlotter.h.

◆ h_eff_pt_pfIsodBID

MonitorElement* EfficiencyPlotter::h_eff_pt_pfIsodBID
private

Definition at line 90 of file EfficiencyPlotter.h.

◆ h_eff_pt_pfIsoID

MonitorElement* EfficiencyPlotter::h_eff_pt_pfIsoID
private

Definition at line 79 of file EfficiencyPlotter.h.

◆ h_eff_vtx_detIsoID

MonitorElement* EfficiencyPlotter::h_eff_vtx_detIsoID
private

Definition at line 83 of file EfficiencyPlotter.h.

◆ h_eff_vtx_EB_detIsoID

MonitorElement* EfficiencyPlotter::h_eff_vtx_EB_detIsoID
private

Definition at line 85 of file EfficiencyPlotter.h.

◆ h_eff_vtx_EB_pfIsodBID

MonitorElement* EfficiencyPlotter::h_eff_vtx_EB_pfIsodBID
private

Definition at line 95 of file EfficiencyPlotter.h.

◆ h_eff_vtx_EB_pfIsoID

MonitorElement* EfficiencyPlotter::h_eff_vtx_EB_pfIsoID
private

Definition at line 87 of file EfficiencyPlotter.h.

◆ h_eff_vtx_EE_detIsoID

MonitorElement* EfficiencyPlotter::h_eff_vtx_EE_detIsoID
private

Definition at line 86 of file EfficiencyPlotter.h.

◆ h_eff_vtx_EE_pfIsodBID

MonitorElement* EfficiencyPlotter::h_eff_vtx_EE_pfIsodBID
private

Definition at line 96 of file EfficiencyPlotter.h.

◆ h_eff_vtx_EE_pfIsoID

MonitorElement* EfficiencyPlotter::h_eff_vtx_EE_pfIsoID
private

Definition at line 88 of file EfficiencyPlotter.h.

◆ h_eff_vtx_pfIsodBID

MonitorElement* EfficiencyPlotter::h_eff_vtx_pfIsodBID
private

Definition at line 94 of file EfficiencyPlotter.h.

◆ h_eff_vtx_pfIsoID

MonitorElement* EfficiencyPlotter::h_eff_vtx_pfIsoID
private

Definition at line 84 of file EfficiencyPlotter.h.

◆ ID_

std::string EfficiencyPlotter::ID_
private

Definition at line 64 of file EfficiencyPlotter.h.

◆ metname

std::string EfficiencyPlotter::metname
private

Definition at line 43 of file EfficiencyPlotter.h.

◆ nevents

int EfficiencyPlotter::nevents
private

Definition at line 38 of file EfficiencyPlotter.h.

◆ nLumiSegs

unsigned int EfficiencyPlotter::nLumiSegs
private

Definition at line 39 of file EfficiencyPlotter.h.

◆ parameters

edm::ParameterSet EfficiencyPlotter::parameters
private

Definition at line 45 of file EfficiencyPlotter.h.

◆ phiBin

int EfficiencyPlotter::phiBin
private

Definition at line 56 of file EfficiencyPlotter.h.

◆ phiMax

double EfficiencyPlotter::phiMax
private

Definition at line 58 of file EfficiencyPlotter.h.

◆ phiMin

double EfficiencyPlotter::phiMin
private

Definition at line 57 of file EfficiencyPlotter.h.

◆ prescaleFactor

int EfficiencyPlotter::prescaleFactor
private

Definition at line 40 of file EfficiencyPlotter.h.

◆ ptBin

int EfficiencyPlotter::ptBin
private

Definition at line 52 of file EfficiencyPlotter.h.

◆ ptMax

double EfficiencyPlotter::ptMax
private

Definition at line 54 of file EfficiencyPlotter.h.

◆ ptMin

double EfficiencyPlotter::ptMin
private

Definition at line 53 of file EfficiencyPlotter.h.

◆ run

int EfficiencyPlotter::run
private

◆ theFolder

std::string EfficiencyPlotter::theFolder
private

Definition at line 97 of file EfficiencyPlotter.h.

◆ vtxBin

int EfficiencyPlotter::vtxBin
private

Definition at line 60 of file EfficiencyPlotter.h.

◆ vtxMax

double EfficiencyPlotter::vtxMax
private

Definition at line 62 of file EfficiencyPlotter.h.

◆ vtxMin

double EfficiencyPlotter::vtxMin
private

Definition at line 61 of file EfficiencyPlotter.h.

EfficiencyPlotter::h_eff_pt_EE_ID
MonitorElement * h_eff_pt_EE_ID
Definition: EfficiencyPlotter.h:72
EfficiencyPlotter::h_eff_hp_eta_ID
MonitorElement * h_eff_hp_eta_ID
Definition: EfficiencyPlotter.h:74
EfficiencyPlotter::h_eff_pt_EB_ID
MonitorElement * h_eff_pt_EB_ID
Definition: EfficiencyPlotter.h:71
EfficiencyPlotter::h_eff_vtx_pfIsoID
MonitorElement * h_eff_vtx_pfIsoID
Definition: EfficiencyPlotter.h:84
EfficiencyPlotter::etaMax
double etaMax
Definition: EfficiencyPlotter.h:50
gather_cfg.cout
cout
Definition: gather_cfg.py:144
EfficiencyPlotter::h_eff_pt_EE_pfIsoID
MonitorElement * h_eff_pt_EE_pfIsoID
Definition: EfficiencyPlotter.h:81
dqm::impl::MonitorElement::setEfficiencyFlag
void setEfficiencyFlag()
Definition: MonitorElement.h:267
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
EfficiencyPlotter::phiMax
double phiMax
Definition: EfficiencyPlotter.h:58
EfficiencyPlotter::h_eff_vtx_EE_detIsoID
MonitorElement * h_eff_vtx_EE_detIsoID
Definition: EfficiencyPlotter.h:86
dqm::legacy::MonitorElement
Definition: MonitorElement.h:462
dqm::legacy::MonitorElement::getTH1F
virtual TH1F * getTH1F() const
Definition: MonitorElement.h:479
EfficiencyPlotter::h_eff_pt_EB_pfIsodBID
MonitorElement * h_eff_pt_EB_pfIsodBID
Definition: EfficiencyPlotter.h:91
EfficiencyPlotter::etaBin
int etaBin
Definition: EfficiencyPlotter.h:48
EfficiencyPlotter::h_eff_inner_eta_ID
MonitorElement * h_eff_inner_eta_ID
Definition: EfficiencyPlotter.h:69
EfficiencyPlotter::theFolder
std::string theFolder
Definition: EfficiencyPlotter.h:97
EfficiencyPlotter::h_eff_vtx_EB_pfIsodBID
MonitorElement * h_eff_vtx_EB_pfIsodBID
Definition: EfficiencyPlotter.h:95
EfficiencyPlotter::ptMax
double ptMax
Definition: EfficiencyPlotter.h:54
EfficiencyPlotter::h_eff_pt_pfIsoID
MonitorElement * h_eff_pt_pfIsoID
Definition: EfficiencyPlotter.h:79
EfficiencyPlotter::h_eff_pt_EB_pfIsoID
MonitorElement * h_eff_pt_EB_pfIsoID
Definition: EfficiencyPlotter.h:80
EfficiencyPlotter::h_eff_inner_phi_ID
MonitorElement * h_eff_inner_phi_ID
Definition: EfficiencyPlotter.h:70
EfficiencyPlotter::parameters
edm::ParameterSet parameters
Definition: EfficiencyPlotter.h:45
EfficiencyPlotter::h_eff_phi_ID
MonitorElement * h_eff_phi_ID
Definition: EfficiencyPlotter.h:75
EfficiencyPlotter::h_eff_inner_pt_ID
MonitorElement * h_eff_inner_pt_ID
Definition: EfficiencyPlotter.h:68
EfficiencyPlotter::h_eff_pt_ID
MonitorElement * h_eff_pt_ID
Definition: EfficiencyPlotter.h:67
EfficiencyPlotter::h_eff_vtx_EE_pfIsoID
MonitorElement * h_eff_vtx_EE_pfIsoID
Definition: EfficiencyPlotter.h:88
EfficiencyPlotter::h_eff_pt_EB_detIsoID
MonitorElement * h_eff_pt_EB_detIsoID
Definition: EfficiencyPlotter.h:77
EfficiencyPlotter::h_eff_eta_ID
MonitorElement * h_eff_eta_ID
Definition: EfficiencyPlotter.h:73
EfficiencyPlotter::h_eff_pt_detIsoID
MonitorElement * h_eff_pt_detIsoID
Definition: EfficiencyPlotter.h:76
EfficiencyPlotter::etaMin
double etaMin
Definition: EfficiencyPlotter.h:49
EfficiencyPlotter::ID_
std::string ID_
Definition: EfficiencyPlotter.h:64
EfficiencyPlotter::h_eff_vtx_detIsoID
MonitorElement * h_eff_vtx_detIsoID
Definition: EfficiencyPlotter.h:83
EfficiencyPlotter::vtxBin
int vtxBin
Definition: EfficiencyPlotter.h:60
EfficiencyPlotter::vtxMin
double vtxMin
Definition: EfficiencyPlotter.h:61
EfficiencyPlotter::h_eff_vtx_pfIsodBID
MonitorElement * h_eff_vtx_pfIsodBID
Definition: EfficiencyPlotter.h:94
EfficiencyPlotter::phiMin
double phiMin
Definition: EfficiencyPlotter.h:57
EfficiencyPlotter::h_eff_vtx_EB_pfIsoID
MonitorElement * h_eff_vtx_EB_pfIsoID
Definition: EfficiencyPlotter.h:87
EfficiencyPlotter::ptBin
int ptBin
Definition: EfficiencyPlotter.h:52
dqm::implementation::IGetter::get
virtual MonitorElement * get(std::string const &fullpath) const
Definition: DQMStore.cc:673
EfficiencyPlotter::h_eff_pt_EE_detIsoID
MonitorElement * h_eff_pt_EE_detIsoID
Definition: EfficiencyPlotter.h:78
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
EfficiencyPlotter::h_eff_pt_EE_pfIsodBID
MonitorElement * h_eff_pt_EE_pfIsodBID
Definition: EfficiencyPlotter.h:92
EfficiencyPlotter::h_eff_vtx_EE_pfIsodBID
MonitorElement * h_eff_vtx_EE_pfIsodBID
Definition: EfficiencyPlotter.h:96
dqm::impl::MonitorElement::setAxisTitle
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:800
EfficiencyPlotter::h_eff_pt_pfIsodBID
MonitorElement * h_eff_pt_pfIsodBID
Definition: EfficiencyPlotter.h:90
EfficiencyPlotter::h_eff_vtx_EB_detIsoID
MonitorElement * h_eff_vtx_EB_detIsoID
Definition: EfficiencyPlotter.h:85
EfficiencyPlotter::vtxMax
double vtxMax
Definition: EfficiencyPlotter.h:62
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
EfficiencyPlotter::phiBin
int phiBin
Definition: EfficiencyPlotter.h:56
EfficiencyPlotter::ptMin
double ptMin
Definition: EfficiencyPlotter.h:53