CMS 3D CMS Logo
List of all members | Public Member Functions | Protected Attributes
MatchCandidateBenchmark Class Reference

To plot Candidate quantities. More...

#include <MatchCandidateBenchmark.h>

Inheritance diagram for MatchCandidateBenchmark:
Benchmark

Public Member Functions

void fillOne (const reco::Candidate &candidate, const reco::Candidate &matchedCandidate)
 fill histograms with a given particle More...
 
void fillOne (const reco::Candidate &candidate, const reco::Candidate &matchedCandidate, const edm::ParameterSet &parameterSet)
 
 MatchCandidateBenchmark (Mode mode)
 
void setup (DQMStore::IBooker &b)
 book histograms More...
 
void setup (DQMStore::IBooker &b, const edm::ParameterSet &parameterSet)
 
 ~MatchCandidateBenchmark () override
 
- Public Member Functions inherited from Benchmark
 Benchmark (Mode mode=DEFAULT)
 
bool isInRange (float pt, float eta, float phi) const
 
virtual void setDirectory (TDirectory *dir)
 
void setParameters (Mode mode)
 
void setRange (float ptMin, float ptMax, float etaMin, float etaMax, float phiMin, float phiMax)
 
void write ()
 write to the TFile, in plain ROOT mode. No need to call this function in DQM mode More...
 
virtual ~Benchmark ()(false)
 

Protected Attributes

TH2F * BRdelta_et_Over_et_VS_et_
 
std::vector< TH1F * > BRpTRes_
 
TH2F * delta_et_Over_et_VS_et_
 
TH2F * delta_et_VS_et_
 
TH2F * delta_eta_VS_et_
 
TH2F * delta_phi_VS_et_
 
TH2F * ERdelta_et_Over_et_VS_et_
 
std::vector< TH1F * > ERpTRes_
 
bool histogramBooked_
 
std::vector< TH1F * > pTRes_
 
- Protected Attributes inherited from Benchmark
TDirectory * dir_
 
float etaMax_
 
float etaMin_
 
Mode mode_
 
float phiMax_
 
float phiMin_
 
float ptMax_
 
float ptMin_
 

Additional Inherited Members

- Public Types inherited from Benchmark
enum  Mode { DEFAULT, DQMOFFLINE, VALIDATION }
 
- Protected Member Functions inherited from Benchmark
TH1F * book1D (DQMStore::IBooker &b, const char *histname, const char *title, int nbins, float xmin, float xmax)
 book a 1D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not. More...
 
TH2F * book2D (DQMStore::IBooker &b, const char *histname, const char *title, int nbinsx, float xmin, float xmax, int nbinsy, float ymin, float ymax)
 book a 2D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not. More...
 
TH2F * book2D (DQMStore::IBooker &b, const char *histname, const char *title, int nbinsx, float *xbins, int nbinsy, float ymin, float ymax)
 book a 2D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not. More...
 
TProfile * bookProfile (DQMStore::IBooker &b, const char *histname, const char *title, int nbinsx, float xmin, float xmax, float ymin, float ymax, const char *option)
 book a TProfile histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not. More...
 
TProfile * bookProfile (DQMStore::IBooker &b, const char *histname, const char *title, int nbinsx, float *xbins, float ymin, float ymax, const char *option)
 book a TProfile histogram, either with DQM or plain root depending if DQM_ has been initialized in a child analyzer or not. More...
 

Detailed Description

To plot Candidate quantities.

Definition at line 13 of file MatchCandidateBenchmark.h.

Constructor & Destructor Documentation

MatchCandidateBenchmark::MatchCandidateBenchmark ( Mode  mode)

Definition at line 17 of file MatchCandidateBenchmark.cc.

References BRdelta_et_Over_et_VS_et_, delta_et_Over_et_VS_et_, delta_et_VS_et_, delta_eta_VS_et_, delta_phi_VS_et_, ERdelta_et_Over_et_VS_et_, and histogramBooked_.

17  : Benchmark(mode) {
18 
19  delta_et_Over_et_VS_et_ = nullptr;
20  delta_et_VS_et_ = nullptr;
21  delta_eta_VS_et_ = nullptr;
22  delta_phi_VS_et_ = nullptr;
23 
24  BRdelta_et_Over_et_VS_et_ = nullptr;
25  ERdelta_et_Over_et_VS_et_ = nullptr;
26  // pTRes are initialzied in the setup since ptBinsPS.size() is needed
27 
28  histogramBooked_ = false;
29 }
ALCARECOPromptCalibProdSiPixelAli0T_cff.mode
mode
Definition: ALCARECOPromptCalibProdSiPixelAli0T_cff.py:96
Benchmark::Benchmark
Benchmark(Mode mode=DEFAULT)
Definition: Benchmark.h:41
MatchCandidateBenchmark::~MatchCandidateBenchmark ( )
override

Definition at line 32 of file MatchCandidateBenchmark.cc.

32 {}

Member Function Documentation

void MatchCandidateBenchmark::fillOne ( const reco::Candidate &  candidate,
const reco::Candidate &  matchedCandidate 
)

fill histograms with a given particle

Definition at line 268 of file MatchCandidateBenchmark.cc.

References BRdelta_et_Over_et_VS_et_, delta_et_Over_et_VS_et_, delta_et_VS_et_, delta_eta_VS_et_, delta_phi_VS_et_, ERdelta_et_Over_et_VS_et_, reco::Candidate::eta(), histogramBooked_, Benchmark::isInRange(), reco::Candidate::phi(), and reco::Candidate::pt().

Referenced by PFCandidateMonitor::fill(), PFJetMonitor::fill(), PFCandidateManager::fill(), and PFMETMonitor::fillOne().

269  {
270 
271  if( !isInRange(cand.pt(), cand.eta(), cand.phi() ) ) return;
272 
273  if (histogramBooked_) {
274  if (delta_et_Over_et_VS_et_) delta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() );
275  if ( fabs(cand.eta()) <= 1.4 )
276  if (BRdelta_et_Over_et_VS_et_) BRdelta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() );
277  if ( fabs(cand.eta()) >= 1.6 && fabs(cand.eta()) <= 2.4 )
278  if (ERdelta_et_Over_et_VS_et_) ERdelta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() );
279  if (delta_et_VS_et_) delta_et_VS_et_->Fill( matchedCand.pt(), cand.pt() - matchedCand.pt() );
280  if (delta_eta_VS_et_) delta_eta_VS_et_->Fill( matchedCand.pt(), cand.eta() - matchedCand.eta() );
281  if (delta_phi_VS_et_) delta_phi_VS_et_->Fill( matchedCand.pt(), cand.phi() - matchedCand.phi() );
282  /*
283  // TProfile
284  if (profile_delta_et_Over_et_VS_et_) {
285  profile_delta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() );
286  profileRMS_delta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() ); }
287  if (profile_delta_et_VS_et_) {
288  profile_delta_et_VS_et_->Fill( matchedCand.pt(), cand.pt() - matchedCand.pt() );
289  profileRMS_delta_et_VS_et_->Fill( matchedCand.pt(), cand.pt() - matchedCand.pt() ); }
290  if (profile_delta_eta_VS_et_) {
291  profile_delta_eta_VS_et_->Fill( matchedCand.pt(), cand.eta() - matchedCand.eta() );
292  profileRMS_delta_eta_VS_et_->Fill( matchedCand.pt(), cand.eta() - matchedCand.eta() ); }
293  if (profile_delta_phi_VS_et_) {
294  profile_delta_phi_VS_et_->Fill( matchedCand.pt(), cand.phi() - matchedCand.phi() );
295  profileRMS_delta_phi_VS_et_->Fill( matchedCand.pt(), cand.phi() - matchedCand.phi() ); }
296  */
297  }
298 }
Benchmark::isInRange
bool isInRange(float pt, float eta, float phi) const
Definition: Benchmark.h:58
void MatchCandidateBenchmark::fillOne ( const reco::Candidate &  candidate,
const reco::Candidate &  matchedCandidate,
const edm::ParameterSet &  parameterSet 
)

Definition at line 300 of file MatchCandidateBenchmark.cc.

References BRdelta_et_Over_et_VS_et_, BRpTRes_, delta_et_Over_et_VS_et_, delta_et_VS_et_, delta_eta_VS_et_, delta_phi_VS_et_, ERdelta_et_Over_et_VS_et_, ERpTRes_, reco::Candidate::eta(), edm::ParameterSet::getParameter(), histogramBooked_, mps_fire::i, Benchmark::isInRange(), reco::Candidate::phi(), reco::Candidate::pt(), and pTRes_.

302  {
303  if( !isInRange(cand.pt(), cand.eta(), cand.phi() ) ) return;
304 
305  if (histogramBooked_) {
306 
307  std::vector<double> ptBinsPS = parameterSet.getParameter< std::vector<double> >( "VariablePtBins" );
308  edm::ParameterSet ptPS = parameterSet.getParameter<edm::ParameterSet>("PtHistoParameter");
309  std::vector<float> ptBins;
310  if (ptBinsPS.size() > 1) {
311  ptBins.reserve(ptBinsPS.size());
312  for (size_t i = 0; i < ptBinsPS.size(); i++) {
313  ptBins.push_back(ptBinsPS[i]);
314  }
315  } else {
316  Int_t nFixedBins = ptPS.getParameter<int32_t>("nBin");
317  ptBins.reserve(nFixedBins + 1);
318  for (Int_t i = 0; i <= nFixedBins; i++) {
319  ptBins.push_back( ptPS.getParameter<double>("xMin") + i*((ptPS.getParameter<double>("xMax") - ptPS.getParameter<double>("xMin")) / nFixedBins) );
320  }
321  ptBinsPS.resize(nFixedBins);
322  }
323 
324  edm::ParameterSet dptOvptPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPtOvPtHistoParameter");
325  if (matchedCand.pt() > ptBins.at(0)) { // underflow problem
326  if (delta_et_Over_et_VS_et_) delta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt()) / matchedCand.pt() );
327  if ( fabs(cand.eta()) >= dptOvptPS.getParameter<double>("BREtaMin") && fabs(cand.eta()) <= dptOvptPS.getParameter<double>("BREtaMax"))
328  if (BRdelta_et_Over_et_VS_et_) BRdelta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() );
329  if ( fabs(cand.eta()) >= dptOvptPS.getParameter<double>("EREtaMin") && fabs(cand.eta()) <= dptOvptPS.getParameter<double>("EREtaMax"))
330  if (ERdelta_et_Over_et_VS_et_) ERdelta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() );
331  if (delta_et_VS_et_) delta_et_VS_et_->Fill( matchedCand.pt(), cand.pt() - matchedCand.pt() );
332  if (delta_eta_VS_et_) delta_eta_VS_et_->Fill( matchedCand.pt(), cand.eta() - matchedCand.eta() );
333  if (delta_phi_VS_et_) delta_phi_VS_et_->Fill( matchedCand.pt(), cand.phi() - matchedCand.phi() );
334  }
335  /*
336  // TProfile
337  if (profile_delta_et_Over_et_VS_et_) {
338  profile_delta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() );
339  profileRMS_delta_et_Over_et_VS_et_->Fill( matchedCand.pt(), (cand.pt() - matchedCand.pt())/matchedCand.pt() ); }
340  if (profile_delta_et_VS_et_) {
341  profile_delta_et_VS_et_->Fill( matchedCand.pt(), cand.pt() - matchedCand.pt() );
342  profileRMS_delta_et_VS_et_->Fill( matchedCand.pt(), cand.pt() - matchedCand.pt() ); }
343  if (profile_delta_eta_VS_et_) {
344  profile_delta_eta_VS_et_->Fill( matchedCand.pt(), cand.eta() - matchedCand.eta() );
345  profileRMS_delta_eta_VS_et_->Fill( matchedCand.pt(), cand.eta() - matchedCand.eta() ); }
346  if (profile_delta_phi_VS_et_) {
347  profile_delta_phi_VS_et_->Fill( matchedCand.pt(), cand.phi() - matchedCand.phi() );
348  profileRMS_delta_phi_VS_et_->Fill( matchedCand.pt(), cand.phi() - matchedCand.phi() ); }
349  */
350 
351  for (size_t i = 0; i < pTRes_.size(); i++) {
352  if (matchedCand.pt() >= ptBins.at(i) && matchedCand.pt() < ptBins.at(i+1)) {
353  if (pTRes_[i]) pTRes_[i]->Fill( (cand.pt() - matchedCand.pt()) / matchedCand.pt() ) ;
354  if ( fabs(cand.eta()) >= dptOvptPS.getParameter<double>("BREtaMin") && fabs(cand.eta()) <= dptOvptPS.getParameter<double>("BREtaMax"))
355  if (BRpTRes_[i]) BRpTRes_[i]->Fill( (cand.pt() - matchedCand.pt()) / matchedCand.pt() ) ; // Fill Barrel
356  if ( fabs(cand.eta()) >= dptOvptPS.getParameter<double>("EREtaMin") && fabs(cand.eta()) <= dptOvptPS.getParameter<double>("EREtaMax"))
357  if (ERpTRes_[i]) ERpTRes_[i]->Fill( (cand.pt() - matchedCand.pt()) / matchedCand.pt() ) ; // Fill Endcap
358  }
359  }
360  }
361 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
mps_fire.i
i
Definition: mps_fire.py:269
MatchCandidateBenchmark::pTRes_
std::vector< TH1F * > pTRes_
Definition: MatchCandidateBenchmark.h:42
MatchCandidateBenchmark::BRpTRes_
std::vector< TH1F * > BRpTRes_
Definition: MatchCandidateBenchmark.h:43
MatchCandidateBenchmark::ERpTRes_
std::vector< TH1F * > ERpTRes_
Definition: MatchCandidateBenchmark.h:44
Benchmark::isInRange
bool isInRange(float pt, float eta, float phi) const
Definition: Benchmark.h:58
void MatchCandidateBenchmark::setup ( DQMStore::IBooker &  b)

book histograms

Definition at line 35 of file MatchCandidateBenchmark.cc.

References Benchmark::book2D(), BRdelta_et_Over_et_VS_et_, BRpTRes_, delta_et_Over_et_VS_et_, delta_et_VS_et_, delta_eta_VS_et_, delta_phi_VS_et_, Benchmark::DQMOFFLINE, ERdelta_et_Over_et_VS_et_, ERpTRes_, histogramBooked_, mps_fire::i, Benchmark::PhaseSpace::m, Benchmark::PhaseSpace::M, Benchmark::mode_, Benchmark::PhaseSpace::n, pTRes_, findQualityFiles::size, and Benchmark::VALIDATION.

Referenced by PFMETMonitor::setup(), PFCandidateMonitor::setup(), and PFJetMonitor::setup().

35  {
36  if (!histogramBooked_) {
37  PhaseSpace ptPS;
38  PhaseSpace dptOvptPS;
39  PhaseSpace dptPS;
40  PhaseSpace detaPS;
41  PhaseSpace dphiPS;
42  switch(mode_) {
43  case VALIDATION:
44  ptPS = PhaseSpace(100,0,1000);
45  dptOvptPS = PhaseSpace( 200, -1, 1);
46  dphiPS = PhaseSpace( 200, -1, 1);
47  detaPS = PhaseSpace( 200, -1, 1);
48  dptPS = PhaseSpace( 100, -100, 100);
49  break;
50  case DQMOFFLINE:
51  default:
52  ptPS = PhaseSpace(50,0,100);
53  dptOvptPS = PhaseSpace( 50, -1, 1);
54  dphiPS = PhaseSpace( 50, -1, 1);
55  detaPS = PhaseSpace( 50, -1, 1);
56  dptPS = PhaseSpace( 50, -50, 50);
57  break;
58  }
59  float ptBins[11] = {0, 1, 2, 5, 10, 20, 50, 100, 200, 400, 1000};
60  int size = sizeof(ptBins)/sizeof(*ptBins);
61 
62  delta_et_Over_et_VS_et_ = book2D(b, "delta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
63  size, ptBins,
64  dptOvptPS.n, dptOvptPS.m, dptOvptPS.M );
65 
66  BRdelta_et_Over_et_VS_et_ = book2D(b, "BRdelta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
67  size, ptBins,
68  dptOvptPS.n, dptOvptPS.m, dptOvptPS.M );
69  ERdelta_et_Over_et_VS_et_ = book2D(b, "ERdelta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
70  size, ptBins,
71  dptOvptPS.n, dptOvptPS.m, dptOvptPS.M );
72 
73  delta_et_VS_et_ = book2D(b, "delta_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}",
74  size, ptBins,
75  dptPS.n, dptPS.m, dptPS.M );
76 
77  delta_eta_VS_et_ = book2D(b, "delta_eta_VS_et_", ";#E_{T, true} (GeV);#Delta#eta",
78  size, ptBins,
79  detaPS.n, detaPS.m, detaPS.M );
80 
81  delta_phi_VS_et_ = book2D(b, "delta_phi_VS_et_", ";E_{T, true} (GeV);#Delta#phi",
82  size, ptBins,
83  dphiPS.n, dphiPS.m, dphiPS.M );
84  /*
85  // TProfile
86  profile_delta_et_Over_et_VS_et_ = bookProfile("profile_delta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
87  size, ptBins,
88  dptOvptPS.m, dptOvptPS.M, "" );
89 
90  profile_delta_et_VS_et_ = bookProfile("profile_delta_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}",
91  size, ptBins,
92  dptPS.m, dptPS.M, "" );
93 
94  profile_delta_eta_VS_et_ = bookProfile("profile_delta_eta_VS_et_", ";#E_{T, true} (GeV);#Delta#eta",
95  size, ptBins,
96  detaPS.m, detaPS.M, "" );
97 
98  profile_delta_phi_VS_et_ = bookProfile("profile_delta_phi_VS_et_", ";E_{T, true} (GeV);#Delta#phi",
99  size, ptBins,
100  dphiPS.m, dphiPS.M, "" );
101  // TProfile RMS
102  profileRMS_delta_et_Over_et_VS_et_ = bookProfile("profileRMS_delta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
103  size, ptBins,
104  dptOvptPS.m, dptOvptPS.M, "s" );
105 
106  profileRMS_delta_et_VS_et_ = bookProfile("profileRMS_delta_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}",
107  size, ptBins,
108  dptPS.m, dptPS.M, "s" );
109 
110  profileRMS_delta_eta_VS_et_ = bookProfile("profileRMS_delta_eta_VS_et_", ";#E_{T, true} (GeV);#Delta#eta",
111  size, ptBins,
112  detaPS.m, detaPS.M, "s" );
113 
114  profileRMS_delta_phi_VS_et_ = bookProfile("profileRMS_delta_phi_VS_et_", ";E_{T, true} (GeV);#Delta#phi",
115  size, ptBins,
116  dphiPS.m, dphiPS.M, "s" );
117  */
118  pTRes_.resize(size); BRpTRes_.resize(size); ERpTRes_.resize(size);
119  for (size_t i = 0; i < pTRes_.size(); i++) {
120  pTRes_[i] = nullptr; BRpTRes_[i] = nullptr; ERpTRes_[i] = nullptr; }
121 
122  histogramBooked_ = true;
123  }
124 }
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:442
mps_fire.i
i
Definition: mps_fire.py:269
MatchCandidateBenchmark::pTRes_
std::vector< TH1F * > pTRes_
Definition: MatchCandidateBenchmark.h:42
MatchCandidateBenchmark::BRpTRes_
std::vector< TH1F * > BRpTRes_
Definition: MatchCandidateBenchmark.h:43
Benchmark::book2D
TH2F * book2D(DQMStore::IBooker &b, const char *histname, const char *title, int nbinsx, float xmin, float xmax, int nbinsy, float ymin, float ymax)
book a 2D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child ...
Definition: Benchmark.cc:29
MatchCandidateBenchmark::ERpTRes_
std::vector< TH1F * > ERpTRes_
Definition: MatchCandidateBenchmark.h:44
Benchmark::mode_
Mode mode_
Definition: Benchmark.h:105
void MatchCandidateBenchmark::setup ( DQMStore::IBooker &  b,
const edm::ParameterSet &  parameterSet 
)

Definition at line 127 of file MatchCandidateBenchmark.cc.

References Benchmark::book1D(), Benchmark::book2D(), BRdelta_et_Over_et_VS_et_, BRpTRes_, delta_et_Over_et_VS_et_, delta_et_VS_et_, delta_eta_VS_et_, delta_phi_VS_et_, ERdelta_et_Over_et_VS_et_, ERpTRes_, edm::ParameterSet::getParameter(), histogramBooked_, mps_fire::i, and pTRes_.

127  {
128 
129  std::vector<double> ptBinsPS = parameterSet.getParameter< std::vector<double> >( "VariablePtBins" );
130  pTRes_.resize(ptBinsPS.size()-1); BRpTRes_.resize(ptBinsPS.size()-1); ERpTRes_.resize(ptBinsPS.size()-1);
131  if (!pTRes_.empty())
132  for (size_t i = 0; i < pTRes_.size(); i++) {
133  pTRes_[i] = nullptr; BRpTRes_[i] = nullptr; ERpTRes_[i] = nullptr; }
134 
135  if (!histogramBooked_) {
136 
137  edm::ParameterSet ptPS = parameterSet.getParameter<edm::ParameterSet>("PtHistoParameter");
138  edm::ParameterSet dptPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPtHistoParameter");
139  edm::ParameterSet dptOvptPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPtOvPtHistoParameter");
140  edm::ParameterSet detaPS = parameterSet.getParameter<edm::ParameterSet>("DeltaEtaHistoParameter");
141  edm::ParameterSet dphiPS = parameterSet.getParameter<edm::ParameterSet>("DeltaPhiHistoParameter");
142 
143  std::vector<float> ptBins;
144  if (ptBinsPS.size() > 1) {
145  ptBins.reserve(ptBinsPS.size());
146  for (size_t i = 0; i < ptBinsPS.size(); i++)
147  ptBins.push_back(ptBinsPS[i]);
148  } else {
149  Int_t nFixedBins = ptPS.getParameter<int32_t>("nBin");
150  ptBins.reserve(nFixedBins+1);
151  for (Int_t i = 0; i <= nFixedBins; i++)
152  ptBins.push_back(ptPS.getParameter<double>("xMin") + i*((ptPS.getParameter<double>("xMax") - ptPS.getParameter<double>("xMin")) / nFixedBins)) ;
153  ptBinsPS.resize(nFixedBins);
154  }
155 
156  if (dptOvptPS.getParameter<bool>("switchOn")) {
157  delta_et_Over_et_VS_et_ = book2D(b, "delta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
158  ptBinsPS.size()-1, &(ptBins.front()),
159  dptOvptPS.getParameter<int32_t>("nBin"),
160  dptOvptPS.getParameter<double>("xMin"),
161  dptOvptPS.getParameter<double>("xMax"));
162  }
163  if (dptOvptPS.getParameter<bool>("slicingOn")) {
164  for (size_t i = 0; i < pTRes_.size(); i++) {
165  pTRes_[i] = book1D( b, TString::Format("Pt%d_%d", (int)ptBins[i], (int)ptBins[i+1]), ";#Deltap_{T}/p_{T};Entries",
166  dptOvptPS.getParameter<int32_t>("nBin"),
167  dptOvptPS.getParameter<double>("xMin"),
168  dptOvptPS.getParameter<double>("xMax"));
169  BRpTRes_[i] = book1D( b, TString::Format("BRPt%d_%d", (int)ptBins[i], (int)ptBins[i+1]), ";#Deltap_{T}/p_{T};Entries",
170  dptOvptPS.getParameter<int32_t>("nBin"),
171  dptOvptPS.getParameter<double>("xMin"),
172  dptOvptPS.getParameter<double>("xMax"));
173  ERpTRes_[i] = book1D( b, TString::Format("ERPt%d_%d", (int)ptBins[i], (int)ptBins[i+1]), ";#Deltap_{T}/p_{T};Entries",
174  dptOvptPS.getParameter<int32_t>("nBin"),
175  dptOvptPS.getParameter<double>("xMin"),
176  dptOvptPS.getParameter<double>("xMax"));
177  }
178  }
179  if (dptOvptPS.getParameter<bool>("BROn")) {
180  BRdelta_et_Over_et_VS_et_ = book2D(b, "BRdelta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
181  ptBinsPS.size()-1, &(ptBins.front()),
182  dptOvptPS.getParameter<int32_t>("nBin"),
183  dptOvptPS.getParameter<double>("xMin"),
184  dptOvptPS.getParameter<double>("xMax"));
185  }
186  if (dptOvptPS.getParameter<bool>("EROn")) {
187  ERdelta_et_Over_et_VS_et_ = book2D(b, "ERdelta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
188  ptBinsPS.size()-1, &(ptBins.front()),
189  dptOvptPS.getParameter<int32_t>("nBin"),
190  dptOvptPS.getParameter<double>("xMin"),
191  dptOvptPS.getParameter<double>("xMax"));
192  }
193 
194  if (dptPS.getParameter<bool>("switchOn")) {
195  delta_et_VS_et_ = book2D(b, "delta_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}",
196  ptBinsPS.size()-1, &(ptBins.front()),
197  dptPS.getParameter<int32_t>("nBin"),
198  dptPS.getParameter<double>("xMin"),
199  dptPS.getParameter<double>("xMax"));
200  }
201 
202  if (detaPS.getParameter<bool>("switchOn")) {
203  delta_eta_VS_et_ = book2D(b, "delta_eta_VS_et_", ";E_{T, true} (GeV);#Delta#eta",
204  ptBinsPS.size()-1, &(ptBins.front()),
205  detaPS.getParameter<int32_t>("nBin"),
206  detaPS.getParameter<double>("xMin"),
207  detaPS.getParameter<double>("xMax"));
208  }
209 
210  if (dphiPS.getParameter<bool>("switchOn")) {
211  delta_phi_VS_et_ = book2D(b, "delta_phi_VS_et_", ";E_{T, true} (GeV);#Delta#phi",
212  ptBinsPS.size()-1, &(ptBins.front()),
213  dphiPS.getParameter<int32_t>("nBin"),
214  dphiPS.getParameter<double>("xMin"),
215  dphiPS.getParameter<double>("xMax"));
216  }
217  /*
218  // TProfile
219  if (dptOvptPS.getParameter<bool>("switchOn")) {
220  profile_delta_et_Over_et_VS_et_ = bookProfile("profile_delta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
221  ptBinsPS.size()-1, &(ptBins.front()),
222  dptOvptPS.getParameter<double>("xMin"),
223  dptOvptPS.getParameter<double>("xMax"), "" );
224  profileRMS_delta_et_Over_et_VS_et_ = bookProfile("profileRMS_delta_et_Over_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}/E_{T}",
225  ptBinsPS.size()-1, &(ptBins.front()),
226  dptOvptPS.getParameter<double>("xMin"),
227  dptOvptPS.getParameter<double>("xMax"), "s" );
228  }
229 
230  if (dptPS.getParameter<bool>("switchOn")) {
231  profile_delta_et_VS_et_ = bookProfile("profile_delta_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}",
232  ptBinsPS.size()-1, &(ptBins.front()),
233  dptPS.getParameter<double>("xMin"),
234  dptPS.getParameter<double>("xMax"), "" );
235  profileRMS_delta_et_VS_et_ = bookProfile("profileRMS_delta_et_VS_et_", ";E_{T, true} (GeV);#DeltaE_{T}",
236  ptBinsPS.size()-1, &(ptBins.front()),
237  dptPS.getParameter<double>("xMin"),
238  dptPS.getParameter<double>("xMax"), "s" );
239  }
240 
241  if (detaPS.getParameter<bool>("switchOn")) {
242  profile_delta_eta_VS_et_ = bookProfile("profile_delta_eta_VS_et_", ";E_{T, true} (GeV);#Delta#eta",
243  ptBinsPS.size()-1, &(ptBins.front()),
244  detaPS.getParameter<double>("xMin"),
245  detaPS.getParameter<double>("xMax"), "" );
246  profileRMS_delta_eta_VS_et_ = bookProfile("profileRMS_delta_eta_VS_et_", ";E_{T, true} (GeV);#Delta#eta",
247  ptBinsPS.size()-1, &(ptBins.front()),
248  detaPS.getParameter<double>("xMin"),
249  detaPS.getParameter<double>("xMax"), "s" );
250  }
251 
252  if (dphiPS.getParameter<bool>("switchOn")) {
253  profile_delta_phi_VS_et_ = bookProfile("profile_delta_phi_VS_et_", ";E_{T, true} (GeV);#Delta#phi",
254  ptBinsPS.size()-1, &(ptBins.front()),
255  dphiPS.getParameter<double>("xMin"),
256  dphiPS.getParameter<double>("xMax"), "" );
257  profileRMS_delta_phi_VS_et_ = bookProfile("profileRMS_delta_phi_VS_et_", ";E_{T, true} (GeV);#Delta#phi",
258  ptBinsPS.size()-1, &(ptBins.front()),
259  dphiPS.getParameter<double>("xMin"),
260  dphiPS.getParameter<double>("xMax"), "s" );
261  }
262  */
263 
264  histogramBooked_ = true;
265  }
266 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Benchmark::book1D
TH1F * book1D(DQMStore::IBooker &b, const char *histname, const char *title, int nbins, float xmin, float xmax)
book a 1D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child ...
Definition: Benchmark.cc:23
mps_fire.i
i
Definition: mps_fire.py:269
MatchCandidateBenchmark::pTRes_
std::vector< TH1F * > pTRes_
Definition: MatchCandidateBenchmark.h:42
MatchCandidateBenchmark::BRpTRes_
std::vector< TH1F * > BRpTRes_
Definition: MatchCandidateBenchmark.h:43
Benchmark::book2D
TH2F * book2D(DQMStore::IBooker &b, const char *histname, const char *title, int nbinsx, float xmin, float xmax, int nbinsy, float ymin, float ymax)
book a 2D histogram, either with DQM or plain root depending if DQM_ has been initialized in a child ...
Definition: Benchmark.cc:29
MatchCandidateBenchmark::ERpTRes_
std::vector< TH1F * > ERpTRes_
Definition: MatchCandidateBenchmark.h:44

Member Data Documentation

TH2F* MatchCandidateBenchmark::BRdelta_et_Over_et_VS_et_
protected

Definition at line 40 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), MatchCandidateBenchmark(), and setup().

std::vector<TH1F*> MatchCandidateBenchmark::BRpTRes_
protected

Definition at line 43 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), and setup().

TH2F* MatchCandidateBenchmark::delta_et_Over_et_VS_et_
protected

Definition at line 35 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), MatchCandidateBenchmark(), and setup().

TH2F* MatchCandidateBenchmark::delta_et_VS_et_
protected

Definition at line 36 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), MatchCandidateBenchmark(), and setup().

TH2F* MatchCandidateBenchmark::delta_eta_VS_et_
protected

Definition at line 37 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), MatchCandidateBenchmark(), and setup().

TH2F* MatchCandidateBenchmark::delta_phi_VS_et_
protected

Definition at line 38 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), MatchCandidateBenchmark(), and setup().

TH2F* MatchCandidateBenchmark::ERdelta_et_Over_et_VS_et_
protected

Definition at line 41 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), MatchCandidateBenchmark(), and setup().

std::vector<TH1F*> MatchCandidateBenchmark::ERpTRes_
protected

Definition at line 44 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), and setup().

bool MatchCandidateBenchmark::histogramBooked_
protected

Definition at line 46 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), MatchCandidateBenchmark(), and setup().

std::vector<TH1F*> MatchCandidateBenchmark::pTRes_
protected

Definition at line 42 of file MatchCandidateBenchmark.h.

Referenced by fillOne(), and setup().

Generated for CMSSW Reference Manual by doxygen 1.8.11