|
|
|
#include <EfficiencyPlotter.h>
Definition at line 26 of file EfficiencyPlotter.h.
| EfficiencyPlotter::EfficiencyPlotter | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 29 of file EfficiencyPlotter.cc.
References gather_cfg::cout, cppFunctionSkipper::operator, and Parameters::parameters.
{
#ifdef DEBUG
cout << "EfficiencyPlotter(): Constructor " << endl;
#endif
parameters = ps;
theDbe = edm::Service<DQMStore>().operator->();
}
| EfficiencyPlotter::~EfficiencyPlotter | ( | ) | [virtual] |
| void EfficiencyPlotter::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
| void EfficiencyPlotter::beginJob | ( | void | ) | [protected, virtual] |
BeginJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 38 of file EfficiencyPlotter.cc.
References gather_cfg::cout, jptDQMConfig_cff::etaMax, jptDQMConfig_cff::etaMin, LogTrace, metname, Parameters::parameters, jptDQMConfig_cff::phiMax, jptDQMConfig_cff::phiMin, jptDQMConfig_cff::ptMax, and PtMinSelector_cfg::ptMin.
{
#ifdef DEBUG
cout << "EfficiencyPlotter::beginJob " << endl;
#endif
theDbe->setCurrentFolder("Muons/EfficiencyAnalyzer");
metname = "EfficiencyAnalyzer";
LogTrace(metname)<<"[EfficiencyPlotter] beginJob: Parameters initialization";
// efficiency plot
etaBin = parameters.getParameter<int>("etaBin");
etaMin = parameters.getParameter<double>("etaMin");
etaMax = parameters.getParameter<double>("etaMax");
phiBin = parameters.getParameter<int>("phiBin");
phiMin = parameters.getParameter<double>("phiMin");
phiMax = parameters.getParameter<double>("phiMax");
ptBin = parameters.getParameter<int>("ptBin");
ptMin = parameters.getParameter<double>("ptMin");
ptMax = parameters.getParameter<double>("ptMax");
vtxBin = parameters.getParameter<int>("vtxBin");
vtxMin = parameters.getParameter<double>("vtxMin");
vtxMax = parameters.getParameter<double>("vtxMax");
h_eff_eta_TightMu = theDbe->book1D("Eff_eta_TightMu", "TightMu Eff. vs #eta", etaBin, etaMin, etaMax);
h_eff_hp_eta_TightMu = theDbe->book1D("Eff_hp_eta_TightMu", "High Pt TightMu Eff. vs #eta", etaBin, etaMin, etaMax);
h_eff_phi_TightMu = theDbe->book1D("Eff_phi_TightMu", "TightMu Eff. vs #phi", phiBin, phiMin, phiMax);
h_eff_pt_TightMu = theDbe->book1D("Eff_pt_TightMu", "TightMu Eff. vs Pt", ptBin, ptMin, ptMax);
h_eff_pt_EB_TightMu = theDbe->book1D("Eff_pt_EB_TightMu", "Barrel: TightMu Eff. vs Pt", ptBin, ptMin, ptMax);
h_eff_pt_EE_TightMu = theDbe->book1D("Eff_pt_EE_TightMu", "Endcap: TightMu Eff. vs Pt", ptBin, ptMin, ptMax);
h_eff_pt_detIsoTightMu = theDbe->book1D("Eff_pt_detIsoTightMu", "detIsoTightMu Efficiency vs Pt", ptBin, ptMin, ptMax);
h_eff_pt_EB_detIsoTightMu = theDbe->book1D("Eff_pt_EB_detIsoTightMu", "Barrel: detIsoTightMu Eff. vs Pt", ptBin, ptMin, ptMax);
h_eff_pt_EE_detIsoTightMu = theDbe->book1D("Eff_pt_EE_detIsoTightMu", "Endcap: detIsoTightMu Eff. vs Pt", ptBin, ptMin, ptMax);
h_eff_pt_pfIsoTightMu = theDbe->book1D("Eff_pt_pfIsoTightMu", "pfIsoTightMu Eff. vs Pt", ptBin, ptMin, ptMax);
h_eff_pt_EB_pfIsoTightMu = theDbe->book1D("Eff_pt_EB_pfIsoTightMu", "Barrel: pfIsoTightMu Eff. vs Pt", ptBin, ptMin, ptMax);
h_eff_pt_EE_pfIsoTightMu = theDbe->book1D("Eff_pt_EE_pfIsoTightMu", "Endcap: pfIsoTightMu Eff. vs Pt", ptBin, ptMin, ptMax);
h_eff_vtx_detIsoTightMu = theDbe->book1D("Eff_vtx_detIsoTightMu", "detIsoTightMu Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
h_eff_vtx_pfIsoTightMu = theDbe->book1D("Eff_vtx_pfIsoTightMu", "pfIsoTightMu Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
h_eff_vtx_EB_detIsoTightMu = theDbe->book1D("Eff_vtx_EB_detIsoTightMu", "Barrel: detIsoTightMu Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
h_eff_vtx_EB_pfIsoTightMu = theDbe->book1D("Eff_vtx_EB_pfIsoTightMu", "Barrel: pfIsoTightMu Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
h_eff_vtx_EE_detIsoTightMu = theDbe->book1D("Eff_vtx_EE_detIsoTightMu", "Endcap: detIsoTightMu Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
h_eff_vtx_EE_pfIsoTightMu = theDbe->book1D("Eff_vtx_EE_pfIsoTightMu", "Endcap: pfIsoTightMu Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
// This prevents this ME to be normalized when drawn into the GUI
h_eff_eta_TightMu ->setEfficiencyFlag();
h_eff_hp_eta_TightMu ->setEfficiencyFlag();
h_eff_phi_TightMu ->setEfficiencyFlag();
h_eff_pt_TightMu ->setEfficiencyFlag();
h_eff_pt_EB_TightMu ->setEfficiencyFlag();
h_eff_pt_EE_TightMu ->setEfficiencyFlag();
h_eff_pt_detIsoTightMu ->setEfficiencyFlag();
h_eff_pt_EB_detIsoTightMu ->setEfficiencyFlag();
h_eff_pt_EE_detIsoTightMu ->setEfficiencyFlag();
h_eff_pt_pfIsoTightMu ->setEfficiencyFlag();
h_eff_pt_EB_pfIsoTightMu ->setEfficiencyFlag();
h_eff_pt_EE_pfIsoTightMu ->setEfficiencyFlag();
h_eff_vtx_detIsoTightMu ->setEfficiencyFlag();
h_eff_vtx_pfIsoTightMu ->setEfficiencyFlag();
h_eff_vtx_EB_detIsoTightMu->setEfficiencyFlag();
h_eff_vtx_EB_pfIsoTightMu ->setEfficiencyFlag();
h_eff_vtx_EE_detIsoTightMu->setEfficiencyFlag();
h_eff_vtx_EE_pfIsoTightMu ->setEfficiencyFlag();
// AXIS TITLES....
h_eff_hp_eta_TightMu ->setAxisTitle("#eta", 1);
h_eff_eta_TightMu ->setAxisTitle("#eta", 1);
h_eff_phi_TightMu ->setAxisTitle("#phi", 1);
h_eff_pt_TightMu ->setAxisTitle("p_{T} (GeV)", 1);
h_eff_pt_EB_TightMu ->setAxisTitle("p_{T} (GeV)", 1);
h_eff_pt_EE_TightMu ->setAxisTitle("p_{T} (GeV)", 1);
h_eff_pt_detIsoTightMu ->setAxisTitle("p_{T} (GeV)", 1);
h_eff_pt_EB_detIsoTightMu ->setAxisTitle("p_{T} (GeV)", 1);
h_eff_pt_EE_detIsoTightMu ->setAxisTitle("p_{T} (GeV)", 1);
h_eff_pt_pfIsoTightMu ->setAxisTitle("p_{T} (GeV)", 1);
h_eff_pt_EB_pfIsoTightMu ->setAxisTitle("p_{T} (GeV)", 1);
h_eff_pt_EE_pfIsoTightMu ->setAxisTitle("p_{T} (GeV)", 1);
h_eff_vtx_detIsoTightMu ->setAxisTitle("Number of PV", 1);
h_eff_vtx_pfIsoTightMu ->setAxisTitle("Number of PV", 1);
h_eff_vtx_EB_detIsoTightMu->setAxisTitle("Number of PV", 1);
h_eff_vtx_EB_pfIsoTightMu ->setAxisTitle("Number of PV", 1);
h_eff_vtx_EE_detIsoTightMu->setAxisTitle("Number of PV", 1);
h_eff_vtx_EE_pfIsoTightMu ->setAxisTitle("Number of PV", 1);
// h_eff_eta_TightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_hp_eta_TightMu ->setAxisTitle("High p_{T} Mu Eff.", 2);
// h_eff_phi_TightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_pt_TightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_pt_EB_TightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_pt_EE_TightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_pt_detIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_pt_EB_detIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_pt_EE_detIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_pt_pfIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_pt_EB_pfIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_pt_EE_pfIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_vtx_detIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_vtx_pfIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_vtx_EB_detIsoTightMu->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_vtx_EB_pfIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_vtx_EE_detIsoTightMu->setAxisTitle("Tight Mu Eff.", 2);
// h_eff_vtx_EE_pfIsoTightMu ->setAxisTitle("Tight Mu Eff.", 2);
}
| void EfficiencyPlotter::beginLuminosityBlock | ( | edm::LuminosityBlock const & | lumiSeg, |
| edm::EventSetup const & | context | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 148 of file EfficiencyPlotter.cc.
References LogTrace, metname, edm::LuminosityBlockBase::run(), and DTTTrigCorrFirst::run.
| void EfficiencyPlotter::beginRun | ( | edm::Run const & | run, |
| edm::EventSetup const & | eSetup | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 145 of file EfficiencyPlotter.cc.
| void EfficiencyPlotter::endJob | ( | void | ) | [protected, virtual] |
| void EfficiencyPlotter::endLuminosityBlock | ( | edm::LuminosityBlock const & | lumiSeg, |
| edm::EventSetup const & | c | ||
| ) | [protected, virtual] |
DQM Client Diagnostic.
Reimplemented from edm::EDAnalyzer.
Definition at line 157 of file EfficiencyPlotter.cc.
References edm::LuminosityBlockBase::id(), and edm::LuminosityBlockID::luminosityBlock().
{
// LogTrace(metname)<<"[EfficiencyPlotter]: endLuminosityBlock, performing the DQM LS client operation";
// counts number of lumiSegs
nLumiSegs = lumiSeg.id().luminosityBlock();
}
| void EfficiencyPlotter::endRun | ( | edm::Run const & | run, |
| edm::EventSetup const & | eSetup | ||
| ) | [protected, virtual] |
--- Tight Muon efficiency vs muon Pt
--- 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]
Reimplemented from edm::EDAnalyzer.
Definition at line 162 of file EfficiencyPlotter.cc.
References MonitorElement::getTH1F(), LogTrace, and metname.
{
LogTrace(metname)<<"[EfficiencyPlotter]: endRun, performing the DQM end of run client operation";
string numpath_pt = "Muons/EfficiencyAnalyzer/passProbes_TightMu_pt";
string denpath_pt = "Muons/EfficiencyAnalyzer/allProbes_pt";
MonitorElement *Numerator_pt = theDbe->get(numpath_pt);
MonitorElement *Denominator_pt = theDbe->get(denpath_pt);
if (Numerator_pt && Denominator_pt){
TH1F *h_numerator_pt = Numerator_pt->getTH1F();
TH1F *h_denominator_pt = Denominator_pt->getTH1F();
TH1F *h_eff_pt = h_eff_pt_TightMu->getTH1F();
if (h_eff_pt->GetSumw2N() == 0) h_eff_pt->Sumw2();
h_eff_pt->Divide(h_numerator_pt, h_denominator_pt, 1., 1., "B");
}
string numpath_EB_pt = "Muons/EfficiencyAnalyzer/passProbes_TightMu_EB_pt";
string denpath_EB_pt = "Muons/EfficiencyAnalyzer/allProbes_EB_pt";
MonitorElement *Numerator_EB_pt = theDbe->get(numpath_EB_pt);
MonitorElement *Denominator_EB_pt = theDbe->get(denpath_EB_pt);
if (Numerator_EB_pt && Denominator_EB_pt){
TH1F *h_numerator_EB_pt = Numerator_EB_pt->getTH1F();
TH1F *h_denominator_EB_pt = Denominator_EB_pt->getTH1F();
TH1F *h_eff_EB_pt = h_eff_pt_EB_TightMu->getTH1F();
if (h_eff_EB_pt->GetSumw2N() == 0) h_eff_EB_pt->Sumw2();
h_eff_EB_pt->Divide(h_numerator_EB_pt, h_denominator_EB_pt, 1., 1., "B");
}
string numpath_EE_pt = "Muons/EfficiencyAnalyzer/passProbes_TightMu_EE_pt";
string denpath_EE_pt = "Muons/EfficiencyAnalyzer/allProbes_EE_pt";
MonitorElement *Numerator_EE_pt = theDbe->get(numpath_EE_pt);
MonitorElement *Denominator_EE_pt = theDbe->get(denpath_EE_pt);
if (Numerator_EE_pt && Denominator_EE_pt){
TH1F *h_numerator_EE_pt = Numerator_EE_pt->getTH1F();
TH1F *h_denominator_EE_pt = Denominator_EE_pt->getTH1F();
TH1F *h_eff_EE_pt = h_eff_pt_EE_TightMu->getTH1F();
if (h_eff_EE_pt->GetSumw2N() == 0) h_eff_EE_pt->Sumw2();
h_eff_EE_pt->Divide(h_numerator_EE_pt, h_denominator_EE_pt, 1., 1., "B");
}
string numpath_eta = "Muons/EfficiencyAnalyzer/passProbes_TightMu_eta";
string denpath_eta = "Muons/EfficiencyAnalyzer/allProbes_eta";
MonitorElement *Numerator_eta = theDbe->get(numpath_eta);
MonitorElement *Denominator_eta = theDbe->get(denpath_eta);
if (Numerator_eta && Denominator_eta){
TH1F *h_numerator_eta = Numerator_eta->getTH1F();
TH1F *h_denominator_eta = Denominator_eta->getTH1F();
TH1F *h_eff_eta = h_eff_eta_TightMu->getTH1F();
if (h_eff_eta->GetSumw2N() == 0) h_eff_eta->Sumw2();
h_eff_eta->Divide(h_numerator_eta, h_denominator_eta, 1., 1., "B");
}
string numpath_hp_eta = "Muons/EfficiencyAnalyzer/passProbes_TightMu_hp_eta";
string denpath_hp_eta = "Muons/EfficiencyAnalyzer/allProbes_hp_eta";
MonitorElement *Numerator_hp_eta = theDbe->get(numpath_hp_eta);
MonitorElement *Denominator_hp_eta = theDbe->get(denpath_hp_eta);
if (Numerator_hp_eta && Denominator_hp_eta){
TH1F *h_numerator_hp_eta = Numerator_hp_eta->getTH1F();
TH1F *h_denominator_hp_eta = Denominator_hp_eta->getTH1F();
TH1F *h_eff_hp_eta = h_eff_hp_eta_TightMu->getTH1F();
if (h_eff_hp_eta->GetSumw2N() == 0) h_eff_hp_eta->Sumw2();
h_eff_hp_eta->Divide(h_numerator_hp_eta, h_denominator_hp_eta, 1., 1., "B");
}
string numpath_phi = "Muons/EfficiencyAnalyzer/passProbes_TightMu_phi";
string denpath_phi = "Muons/EfficiencyAnalyzer/allProbes_phi";
MonitorElement *Numerator_phi = theDbe->get(numpath_phi);
MonitorElement *Denominator_phi = theDbe->get(denpath_phi);
if (Numerator_phi && Denominator_phi){
TH1F *h_numerator_phi = Numerator_phi->getTH1F();
TH1F *h_denominator_phi = Denominator_phi->getTH1F();
TH1F *h_eff_phi = h_eff_phi_TightMu->getTH1F();
if (h_eff_phi->GetSumw2N() == 0) h_eff_phi->Sumw2();
h_eff_phi->Divide(h_numerator_phi, h_denominator_phi, 1., 1., "B");
}
string numpath_detIso_pt = "Muons/EfficiencyAnalyzer/passProbes_detIsoTightMu_pt";
string denpath_detIso_pt = "Muons/EfficiencyAnalyzer/allProbes_TightMu_pt";
MonitorElement *Numerator_detIso_pt = theDbe->get(numpath_detIso_pt);
MonitorElement *Denominator_detIso_pt = theDbe->get(denpath_detIso_pt);
if (Numerator_detIso_pt && Denominator_detIso_pt){
TH1F *h_numerator_detIso_pt = Numerator_detIso_pt->getTH1F();
TH1F *h_denominator_detIso_pt = Denominator_detIso_pt->getTH1F();
TH1F *h_eff_detIso_pt = h_eff_pt_detIsoTightMu->getTH1F();
if (h_eff_detIso_pt->GetSumw2N() == 0) h_eff_detIso_pt->Sumw2();
h_eff_detIso_pt->Divide(h_numerator_detIso_pt, h_denominator_detIso_pt, 1., 1., "B");
}
string numpath_detIso_EB_pt = "Muons/EfficiencyAnalyzer/passProbes_EB_detIsoTightMu_pt";
string denpath_detIso_EB_pt = "Muons/EfficiencyAnalyzer/allProbes_EB_TightMu_pt";
MonitorElement *Numerator_detIso_EB_pt = theDbe->get(numpath_detIso_EB_pt);
MonitorElement *Denominator_detIso_EB_pt = theDbe->get(denpath_detIso_EB_pt);
if (Numerator_detIso_EB_pt && Denominator_detIso_EB_pt){
TH1F *h_numerator_detIso_EB_pt = Numerator_detIso_EB_pt->getTH1F();
TH1F *h_denominator_detIso_EB_pt = Denominator_detIso_EB_pt->getTH1F();
TH1F *h_eff_detIso_EB_pt = h_eff_pt_EB_detIsoTightMu->getTH1F();
if (h_eff_detIso_EB_pt->GetSumw2N() == 0) h_eff_detIso_EB_pt->Sumw2();
h_eff_detIso_EB_pt->Divide(h_numerator_detIso_EB_pt, h_denominator_detIso_EB_pt, 1., 1., "B");
}
string numpath_detIso_EE_pt = "Muons/EfficiencyAnalyzer/passProbes_EE_detIsoTightMu_pt";
string denpath_detIso_EE_pt = "Muons/EfficiencyAnalyzer/allProbes_EE_TightMu_pt";
MonitorElement *Numerator_detIso_EE_pt = theDbe->get(numpath_detIso_EE_pt);
MonitorElement *Denominator_detIso_EE_pt = theDbe->get(denpath_detIso_EE_pt);
if (Numerator_detIso_EE_pt && Denominator_detIso_EE_pt){
TH1F *h_numerator_detIso_EE_pt = Numerator_detIso_EE_pt->getTH1F();
TH1F *h_denominator_detIso_EE_pt = Denominator_detIso_EE_pt->getTH1F();
TH1F *h_eff_detIso_EE_pt = h_eff_pt_EE_detIsoTightMu->getTH1F();
if (h_eff_detIso_EE_pt->GetSumw2N() == 0) h_eff_detIso_EE_pt->Sumw2();
h_eff_detIso_EE_pt->Divide(h_numerator_detIso_EE_pt, h_denominator_detIso_EE_pt, 1., 1., "B");
}
string numpath_pfIso_pt = "Muons/EfficiencyAnalyzer/passProbes_pfIsoTightMu_pt";
string denpath_pfIso_pt = "Muons/EfficiencyAnalyzer/allProbes_TightMu_pt";
MonitorElement *Numerator_pfIso_pt = theDbe->get(numpath_pfIso_pt);
MonitorElement *Denominator_pfIso_pt = theDbe->get(denpath_pfIso_pt);
if (Numerator_pfIso_pt && Denominator_pfIso_pt){
TH1F *h_numerator_pfIso_pt = Numerator_pfIso_pt->getTH1F();
TH1F *h_denominator_pfIso_pt = Denominator_pfIso_pt->getTH1F();
TH1F *h_eff_pfIso_pt = h_eff_pt_pfIsoTightMu->getTH1F();
if (h_eff_pfIso_pt->GetSumw2N() == 0) h_eff_pfIso_pt->Sumw2();
h_eff_pfIso_pt->Divide(h_numerator_pfIso_pt, h_denominator_pfIso_pt, 1., 1., "B");
}
string numpath_pfIso_EB_pt = "Muons/EfficiencyAnalyzer/passProbes_EB_pfIsoTightMu_pt";
string denpath_pfIso_EB_pt = "Muons/EfficiencyAnalyzer/allProbes_EB_TightMu_pt";
MonitorElement *Numerator_pfIso_EB_pt = theDbe->get(numpath_pfIso_EB_pt);
MonitorElement *Denominator_pfIso_EB_pt = theDbe->get(denpath_pfIso_EB_pt);
if (Numerator_pfIso_EB_pt && Denominator_pfIso_EB_pt){
TH1F *h_numerator_pfIso_EB_pt = Numerator_pfIso_EB_pt->getTH1F();
TH1F *h_denominator_pfIso_EB_pt = Denominator_pfIso_EB_pt->getTH1F();
TH1F *h_eff_pfIso_EB_pt = h_eff_pt_EB_pfIsoTightMu->getTH1F();
if (h_eff_pfIso_EB_pt->GetSumw2N() == 0) h_eff_pfIso_EB_pt->Sumw2();
h_eff_pfIso_EB_pt->Divide(h_numerator_pfIso_EB_pt, h_denominator_pfIso_EB_pt, 1., 1., "B");
}
string numpath_pfIso_EE_pt = "Muons/EfficiencyAnalyzer/passProbes_EE_pfIsoTightMu_pt";
string denpath_pfIso_EE_pt = "Muons/EfficiencyAnalyzer/allProbes_EE_TightMu_pt";
MonitorElement *Numerator_pfIso_EE_pt = theDbe->get(numpath_pfIso_EE_pt);
MonitorElement *Denominator_pfIso_EE_pt = theDbe->get(denpath_pfIso_EE_pt);
if (Numerator_pfIso_EE_pt && Denominator_pfIso_EE_pt){
TH1F *h_numerator_pfIso_EE_pt = Numerator_pfIso_EE_pt->getTH1F();
TH1F *h_denominator_pfIso_EE_pt = Denominator_pfIso_EE_pt->getTH1F();
TH1F *h_eff_pfIso_EE_pt = h_eff_pt_EE_pfIsoTightMu->getTH1F();
if (h_eff_pfIso_EE_pt->GetSumw2N() == 0) h_eff_pfIso_EE_pt->Sumw2();
h_eff_pfIso_EE_pt->Divide(h_numerator_pfIso_EE_pt, h_denominator_pfIso_EE_pt, 1., 1., "B");
}
string numpath_pfIso_nvtx = "Muons/EfficiencyAnalyzer/passProbes_pfIsoTightMu_nVtx";
string denpath_pfIso_nvtx = "Muons/EfficiencyAnalyzer/allProbes_TightMu_nVtx";
MonitorElement *Numerator_pfIso_nvtx = theDbe->get(numpath_pfIso_nvtx);
MonitorElement *Denominator_pfIso_nvtx = theDbe->get(denpath_pfIso_nvtx);
if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->getTH1F();
TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->getTH1F();
TH1F *h_eff_pfIso_nvtx = h_eff_vtx_pfIsoTightMu->getTH1F();
if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1., "B");
}
string numpath_detIso_nvtx = "Muons/EfficiencyAnalyzer/passProbes_detIsoTightMu_nVtx";
string denpath_detIso_nvtx = "Muons/EfficiencyAnalyzer/allProbes_TightMu_nVtx";
MonitorElement *Numerator_detIso_nvtx = theDbe->get(numpath_detIso_nvtx);
MonitorElement *Denominator_detIso_nvtx = theDbe->get(denpath_detIso_nvtx);
if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->getTH1F();
TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->getTH1F();
TH1F *h_eff_detIso_nvtx = h_eff_vtx_detIsoTightMu->getTH1F();
if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1., "B");
}
numpath_detIso_nvtx = "Muons/EfficiencyAnalyzer/passProbes_EB_detIsoTightMu_nVtx";
denpath_detIso_nvtx = "Muons/EfficiencyAnalyzer/allProbes_EB_TightMu_nVtx";
Numerator_detIso_nvtx = theDbe->get(numpath_detIso_nvtx);
Denominator_detIso_nvtx = theDbe->get(denpath_detIso_nvtx);
if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->getTH1F();
TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->getTH1F();
TH1F *h_eff_detIso_nvtx = h_eff_vtx_EB_detIsoTightMu->getTH1F();
if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1., "B");
}
numpath_detIso_nvtx = "Muons/EfficiencyAnalyzer/passProbes_EE_detIsoTightMu_nVtx";
denpath_detIso_nvtx = "Muons/EfficiencyAnalyzer/allProbes_EE_TightMu_nVtx";
Numerator_detIso_nvtx = theDbe->get(numpath_detIso_nvtx);
Denominator_detIso_nvtx = theDbe->get(denpath_detIso_nvtx);
if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->getTH1F();
TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->getTH1F();
TH1F *h_eff_detIso_nvtx = h_eff_vtx_EE_detIsoTightMu->getTH1F();
if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1., "B");
}
numpath_pfIso_nvtx = "Muons/EfficiencyAnalyzer/passProbes_EB_pfIsoTightMu_nVtx";
denpath_pfIso_nvtx = "Muons/EfficiencyAnalyzer/allProbes_EB_TightMu_nVtx";
Numerator_pfIso_nvtx = theDbe->get(numpath_pfIso_nvtx);
Denominator_pfIso_nvtx = theDbe->get(denpath_pfIso_nvtx);
if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->getTH1F();
TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->getTH1F();
TH1F *h_eff_pfIso_nvtx = h_eff_vtx_EB_pfIsoTightMu->getTH1F();
if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1., "B");
}
numpath_pfIso_nvtx = "Muons/EfficiencyAnalyzer/passProbes_EE_pfIsoTightMu_nVtx";
denpath_pfIso_nvtx = "Muons/EfficiencyAnalyzer/allProbes_EE_TightMu_nVtx";
Numerator_pfIso_nvtx = theDbe->get(numpath_pfIso_nvtx);
Denominator_pfIso_nvtx = theDbe->get(denpath_pfIso_nvtx);
if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->getTH1F();
TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->getTH1F();
TH1F *h_eff_pfIso_nvtx = h_eff_vtx_EE_pfIsoTightMu->getTH1F();
if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1., "B");
}
}
int EfficiencyPlotter::etaBin [private] |
Definition at line 69 of file EfficiencyPlotter.h.
double EfficiencyPlotter::etaMax [private] |
Definition at line 71 of file EfficiencyPlotter.h.
double EfficiencyPlotter::etaMin [private] |
Definition at line 70 of file EfficiencyPlotter.h.
Definition at line 89 of file EfficiencyPlotter.h.
Definition at line 90 of file EfficiencyPlotter.h.
Definition at line 91 of file EfficiencyPlotter.h.
Definition at line 92 of file EfficiencyPlotter.h.
Definition at line 93 of file EfficiencyPlotter.h.
Definition at line 96 of file EfficiencyPlotter.h.
Definition at line 87 of file EfficiencyPlotter.h.
Definition at line 94 of file EfficiencyPlotter.h.
Definition at line 97 of file EfficiencyPlotter.h.
Definition at line 88 of file EfficiencyPlotter.h.
Definition at line 95 of file EfficiencyPlotter.h.
Definition at line 86 of file EfficiencyPlotter.h.
Definition at line 99 of file EfficiencyPlotter.h.
Definition at line 101 of file EfficiencyPlotter.h.
Definition at line 103 of file EfficiencyPlotter.h.
Definition at line 102 of file EfficiencyPlotter.h.
Definition at line 104 of file EfficiencyPlotter.h.
Definition at line 100 of file EfficiencyPlotter.h.
std::string EfficiencyPlotter::metname [private] |
Definition at line 63 of file EfficiencyPlotter.h.
int EfficiencyPlotter::nevents [private] |
Definition at line 58 of file EfficiencyPlotter.h.
unsigned int EfficiencyPlotter::nLumiSegs [private] |
Definition at line 59 of file EfficiencyPlotter.h.
Definition at line 66 of file EfficiencyPlotter.h.
int EfficiencyPlotter::phiBin [private] |
Definition at line 77 of file EfficiencyPlotter.h.
double EfficiencyPlotter::phiMax [private] |
Definition at line 79 of file EfficiencyPlotter.h.
double EfficiencyPlotter::phiMin [private] |
Definition at line 78 of file EfficiencyPlotter.h.
int EfficiencyPlotter::prescaleFactor [private] |
Definition at line 60 of file EfficiencyPlotter.h.
int EfficiencyPlotter::ptBin [private] |
Definition at line 73 of file EfficiencyPlotter.h.
double EfficiencyPlotter::ptMax [private] |
Definition at line 75 of file EfficiencyPlotter.h.
double EfficiencyPlotter::ptMin [private] |
Definition at line 74 of file EfficiencyPlotter.h.
int EfficiencyPlotter::run [private] |
Definition at line 61 of file EfficiencyPlotter.h.
DQMStore* EfficiencyPlotter::theDbe [private] |
Definition at line 65 of file EfficiencyPlotter.h.
int EfficiencyPlotter::vtxBin [private] |
Definition at line 81 of file EfficiencyPlotter.h.
double EfficiencyPlotter::vtxMax [private] |
Definition at line 83 of file EfficiencyPlotter.h.
double EfficiencyPlotter::vtxMin [private] |
Definition at line 82 of file EfficiencyPlotter.h.
1.7.3