CMS 3D CMS Logo
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
JetMonitor.cc
Go to the documentation of this file.
1 #include "DQMOffline/Trigger/plugins/JetMonitor.h"
2 
3 #include <utility>
4 
5 #include "FWCore/MessageLogger/interface/MessageLogger.h"
6 
7 #include "DQM/TrackingMonitor/interface/GetLumi.h"
8 
9 #include "CommonTools/TriggerUtils/interface/GenericTriggerEventFlag.h"
10 
11 // -----------------------------
12 // constructors and destructor
13 // -----------------------------
14 
15 JetMonitor::JetMonitor(const edm::ParameterSet& iConfig)
16  : num_genTriggerEventFlag_(new GenericTriggerEventFlag(
17  iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"), consumesCollector(), *this)),
18  den_genTriggerEventFlag_(new GenericTriggerEventFlag(
19  iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet"), consumesCollector(), *this)) {
20  folderName_ = iConfig.getParameter<std::string>("FolderName");
21  jetSrc_ = mayConsume<edm::View<reco::Jet> >(iConfig.getParameter<edm::InputTag>("jetSrc")); //jet
22  jetpT_variable_binning_ =
23  iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<std::vector<double> >("jetptBinning");
24  jetpT_binning =
25  getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("jetPSet"));
26  jetptThr_binning_ = getHistoPSet(
27  iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("jetPtThrPSet"));
28  ls_binning_ =
29  getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("lsPSet"));
30 
31  ptcut_ = iConfig.getParameter<double>("ptcut"); // for HLT Jet
32  isPFJetTrig = iConfig.getParameter<bool>("ispfjettrg");
33  isCaloJetTrig = iConfig.getParameter<bool>("iscalojettrg");
34 }
35 
36 JetMonitor::~JetMonitor() = default;
37 
38 JetMonitor::MEbinning JetMonitor::getHistoPSet(const edm::ParameterSet& pset) {
39  return JetMonitor::MEbinning{
40  pset.getParameter<unsigned int>("nbins"),
41  pset.getParameter<double>("xmin"),
42  pset.getParameter<double>("xmax"),
43  };
44 }
45 
46 JetMonitor::MEbinning JetMonitor::getHistoLSPSet(const edm::ParameterSet& pset) {
47  return JetMonitor::MEbinning{
48  pset.getParameter<unsigned int>("nbins"), 0., double(pset.getParameter<unsigned int>("nbins"))};
49 }
50 
51 void JetMonitor::setMETitle(JetME& me, const std::string& titleX, const std::string& titleY) {
52  me.numerator->setAxisTitle(titleX, 1);
53  me.numerator->setAxisTitle(titleY, 2);
54  me.denominator->setAxisTitle(titleX, 1);
55  me.denominator->setAxisTitle(titleY, 2);
56 }
57 void JetMonitor::bookME(DQMStore::IBooker& ibooker,
58  JetME& me,
59  std::string& histname,
60  std::string& histtitle,
61  unsigned int nbins,
62  double min,
63  double max) {
64  me.numerator = ibooker.book1D(histname + "_numerator", histtitle + " (numerator)", nbins, min, max);
65  me.denominator = ibooker.book1D(histname + "_denominator", histtitle + " (denominator)", nbins, min, max);
66 }
67 void JetMonitor::bookME(
68  DQMStore::IBooker& ibooker, JetME& me, std::string& histname, std::string& histtitle, std::vector<double> binning) {
69  int nbins = binning.size() - 1;
70  std::vector<float> fbinning(binning.begin(), binning.end());
71  float* arr = &fbinning[0];
72  me.numerator = ibooker.book1D(histname + "_numerator", histtitle + " (numerator)", nbins, arr);
73  me.denominator = ibooker.book1D(histname + "_denominator", histtitle + " (denominator)", nbins, arr);
74 }
75 void JetMonitor::bookME(DQMStore::IBooker& ibooker,
76  JetME& me,
77  std::string& histname,
78  std::string& histtitle,
79  int nbinsX,
80  double xmin,
81  double xmax,
82  double ymin,
83  double ymax) {
84  me.numerator =
85  ibooker.bookProfile(histname + "_numerator", histtitle + " (numerator)", nbinsX, xmin, xmax, ymin, ymax);
86  me.denominator =
87  ibooker.bookProfile(histname + "_denominator", histtitle + " (denominator)", nbinsX, xmin, xmax, ymin, ymax);
88 }
89 void JetMonitor::bookME(DQMStore::IBooker& ibooker,
90  JetME& me,
91  std::string& histname,
92  std::string& histtitle,
93  int nbinsX,
94  double xmin,
95  double xmax,
96  int nbinsY,
97  double ymin,
98  double ymax) {
99  me.numerator =
100  ibooker.book2D(histname + "_numerator", histtitle + " (numerator)", nbinsX, xmin, xmax, nbinsY, ymin, ymax);
101  me.denominator =
102  ibooker.book2D(histname + "_denominator", histtitle + " (denominator)", nbinsX, xmin, xmax, nbinsY, ymin, ymax);
103 }
104 void JetMonitor::bookME(DQMStore::IBooker& ibooker,
105  JetME& me,
106  std::string& histname,
107  std::string& histtitle,
108  std::vector<double> binningX,
109  std::vector<double> binningY) {
110  int nbinsX = binningX.size() - 1;
111  std::vector<float> fbinningX(binningX.begin(), binningX.end());
112  float* arrX = &fbinningX[0];
113  int nbinsY = binningY.size() - 1;
114  std::vector<float> fbinningY(binningY.begin(), binningY.end());
115  float* arrY = &fbinningY[0];
116 
117  me.numerator = ibooker.book2D(histname + "_numerator", histtitle + " (numerator)", nbinsX, arrX, nbinsY, arrY);
118  me.denominator = ibooker.book2D(histname + "_denominator", histtitle + " (denominator)", nbinsX, arrX, nbinsY, arrY);
119 }
120 
121 void JetMonitor::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const& iRun, edm::EventSetup const& iSetup) {
122  std::string histname, histtitle;
123  std::string hist_obtag = "";
124  std::string histtitle_obtag = "";
125  std::string currentFolder = folderName_;
126  ibooker.setCurrentFolder(currentFolder);
127 
128  if (isPFJetTrig) {
129  hist_obtag = "pfjet";
130  histtitle_obtag = "PFJet";
131  } else if (isCaloJetTrig) {
132  hist_obtag = "calojet";
133  histtitle_obtag = "CaloJet";
134  } else {
135  hist_obtag = "pfjet";
136  histtitle_obtag = "PFJet";
137  } //default is pfjet
138 
139  bookMESub(ibooker, a_ME, sizeof(a_ME) / sizeof(a_ME[0]), hist_obtag, histtitle_obtag, "", "");
140  bookMESub(ibooker,
141  a_ME_HB,
142  sizeof(a_ME_HB) / sizeof(a_ME_HB[0]),
143  hist_obtag,
144  histtitle_obtag,
145  "HB",
146  "(HB)",
147  true,
148  true,
149  true,
150  false);
151  bookMESub(ibooker,
152  a_ME_HE,
153  sizeof(a_ME_HE) / sizeof(a_ME_HE[0]),
154  hist_obtag,
155  histtitle_obtag,
156  "HE",
157  "(HE)",
158  true,
159  true,
160  true,
161  false);
162  bookMESub(ibooker,
163  a_ME_HF,
164  sizeof(a_ME_HF) / sizeof(a_ME_HF[0]),
165  hist_obtag,
166  histtitle_obtag,
167  "HF",
168  "(HF)",
169  true,
170  true,
171  true,
172  false);
173  bookMESub(ibooker,
174  a_ME_HE_p,
175  sizeof(a_ME_HE_p) / sizeof(a_ME_HE_p[0]),
176  hist_obtag,
177  histtitle_obtag,
178  "HE_p",
179  "(HE+)",
180  true,
181  true,
182  true,
183  false);
184  bookMESub(ibooker,
185  a_ME_HE_m,
186  sizeof(a_ME_HE_m) / sizeof(a_ME_HE_m[0]),
187  hist_obtag,
188  histtitle_obtag,
189  "HE_m",
190  "(HE-)",
191  true,
192  true,
193  true,
194  false);
195  bookMESub(ibooker,
196  a_ME_HEP17,
197  sizeof(a_ME_HEP17) / sizeof(a_ME_HEP17[0]),
198  hist_obtag,
199  histtitle_obtag,
200  "HEP17",
201  "(HEP17)",
202  true,
203  false,
204  false,
205  false);
206  bookMESub(ibooker,
207  a_ME_HEM17,
208  sizeof(a_ME_HEM17) / sizeof(a_ME_HEM17[0]),
209  hist_obtag,
210  histtitle_obtag,
211  "HEM17",
212  "(HEM17)",
213  true,
214  false,
215  false,
216  false);
217  bookMESub(ibooker,
218  a_ME_HEP18,
219  sizeof(a_ME_HEP18) / sizeof(a_ME_HEP18[0]),
220  hist_obtag,
221  histtitle_obtag,
222  "HEP18",
223  "(HEP18)",
224  false,
225  false,
226  false,
227  false);
228 
229  /*
230  WE WOULD NEED TURNON CURVES TO BE COMPARED NOT JUST THE ZOOM OF A 2D MAP !!!
231 
232  histname = hist_obtag +"AbsEtaVsPhi_HEP17"; histtitle = histtitle_obtag + " |eta| Vs phi (HEP17) ";
233  bookME(ibooker,jetHEP17_AbsEtaVsPhi_,histname,histtitle, eta_binning_hep17_.nbins, eta_binning_hep17_.xmin, eta_binning_hep17_.xmax, phi_binning_hep17_.nbins,phi_binning_hep17_.xmin,phi_binning_hep17_.xmax);
234  setMETitle(jetHEP17_AbsEtaVsPhi_,histtitle_obtag + " |#eta|","#phi");
235 
236  histname = hist_obtag +"AbsEtaVsPhi_HEM17"; histtitle = histtitle_obtag + " |eta| Vs phi (HEM17) ";
237  bookME(ibooker,jetHEM17_AbsEtaVsPhi_,histname,histtitle, eta_binning_hep17_.nbins, eta_binning_hep17_.xmin, eta_binning_hep17_.xmax, phi_binning_hep17_.nbins,phi_binning_hep17_.xmin,phi_binning_hep17_.xmax);
238  setMETitle(jetHEM17_AbsEtaVsPhi_,histtitle_obtag + " |#eta|","#phi");
239  */
240 
241  histname = hist_obtag + "abseta_HEP17";
242  histtitle = histtitle_obtag + " |#eta| (HEP17) ";
243  bookME(ibooker,
244  jetHEP17_AbsEta_,
245  histname,
246  histtitle,
247  eta_binning_hep17_.nbins,
248  eta_binning_hep17_.xmin,
249  eta_binning_hep17_.xmax);
250  setMETitle(jetHEP17_AbsEta_, histtitle_obtag + " |#eta|", "events / |#eta|");
251 
252  histname = hist_obtag + "abseta_HEM17";
253  histtitle = histtitle_obtag + " |eta| (HEM17) ";
254  bookME(ibooker,
255  jetHEM17_AbsEta_,
256  histname,
257  histtitle,
258  eta_binning_hep17_.nbins,
259  eta_binning_hep17_.xmin,
260  eta_binning_hep17_.xmax);
261  setMETitle(jetHEM17_AbsEta_, histtitle_obtag + " |#eta|", "events / |#eta|");
262 
263  // Initialize the GenericTriggerEventFlag
264  if (num_genTriggerEventFlag_ && num_genTriggerEventFlag_->on())
265  num_genTriggerEventFlag_->initRun(iRun, iSetup);
266  if (den_genTriggerEventFlag_ && den_genTriggerEventFlag_->on())
267  den_genTriggerEventFlag_->initRun(iRun, iSetup);
268 }
269 
270 #include "FWCore/Framework/interface/ESHandle.h"
271 #include "FWCore/Framework/interface/EventSetup.h"
272 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
273 #include "Geometry/Records/interface/TrackerTopologyRcd.h"
274 #include "DataFormats/Math/interface/deltaR.h" // For Delta R
275 void JetMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) {
276  // Filter out events if Trigger Filtering is requested
277  if (den_genTriggerEventFlag_->on() && !den_genTriggerEventFlag_->accept(iEvent, iSetup))
278  return;
279  // edm::Handle<reco::PFJetCollection> pfjetHandle;
280  // iEvent.getByToken( pfjetToken_, pfjetHandle );
281 
282  // edm::Handle<reco::CaloJetCollection> calojetHandle;
283  // iEvent.getByToken( calojetToken_, calojetHandle );
284 
285  int ls = iEvent.id().luminosityBlock();
286  v_jetpt.clear();
287  v_jeteta.clear();
288  v_jetphi.clear();
289 
290  edm::Handle<edm::View<reco::Jet> > offjets;
291  iEvent.getByToken(jetSrc_, offjets);
292  if (!offjets.isValid()) {
293  edm::LogWarning("JetMonitor") << "Jet handle not valid \n";
294  return;
295  }
296  for (edm::View<reco::Jet>::const_iterator ibegin = offjets->begin(), iend = offjets->end(), ijet = ibegin;
297  ijet != iend;
298  ++ijet) {
299  //if (ijet->pt()< 20) {continue;}
300  if (ijet->pt() < ptcut_) {
301  continue;
302  }
303  v_jetpt.push_back(ijet->pt());
304  v_jeteta.push_back(ijet->eta());
305  v_jetphi.push_back(ijet->phi());
306  // cout << "jetpt (view ) : " << ijet->pt() << endl;
307  }
308 
309  if (v_jetpt.empty())
310  return;
311  double jetpt_ = v_jetpt[0];
312  double jeteta_ = v_jeteta[0];
313  double jetphi_ = v_jetphi[0];
314 
315  FillME(a_ME, jetpt_, jetphi_, jeteta_, ls, "denominator");
316  if (isBarrel(jeteta_)) {
317  FillME(a_ME_HB, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false);
318  } else if (isEndCapP(jeteta_)) {
319  FillME(a_ME_HE, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false);
320  FillME(a_ME_HE_p, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false);
321  } else if (isEndCapM(jeteta_)) {
322  FillME(a_ME_HE, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false);
323  FillME(a_ME_HE_m, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false);
324  } else if (isForward(jeteta_)) {
325  FillME(a_ME_HF, jetpt_, jetphi_, jeteta_, ls, "denominator", true, true, true, false);
326  }
327 
328  if (isHEP17(jeteta_, jetphi_)) {
329  FillME(a_ME_HEP17,
330  jetpt_,
331  jetphi_,
332  jeteta_,
333  ls,
334  "denominator",
335  true,
336  false,
337  false,
338  false); // doPhi, doEta, doEtaPhi, doVsLS
339  jetHEP17_AbsEta_.denominator->Fill(abs(jeteta_));
340  } else if (isHEM17(jeteta_, jetphi_)) {
341  FillME(
342  a_ME_HEM17, jetpt_, jetphi_, jeteta_, ls, "denominator", true, false, false, false); // doPhi, doEta, doEtaPhi
343  jetHEM17_AbsEta_.denominator->Fill(abs(jeteta_));
344  } else if (isHEP18(jeteta_, jetphi_)) {
345  FillME(
346  a_ME_HEP18, jetpt_, jetphi_, jeteta_, ls, "denominator", false, false, false, false); // doPhi, doEta, doEtaPhi
347  }
348 
349  if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup))
350  return; // Require Numerator //
351 
352  FillME(a_ME, jetpt_, jetphi_, jeteta_, ls, "numerator");
353  if (isBarrel(jeteta_)) {
354  FillME(a_ME_HB, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false);
355  } else if (isEndCapP(jeteta_)) {
356  FillME(a_ME_HE, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false);
357  FillME(a_ME_HE_p, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false);
358  } else if (isEndCapM(jeteta_)) {
359  FillME(a_ME_HE, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false);
360  FillME(a_ME_HE_m, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false);
361  } else if (isForward(jeteta_)) {
362  FillME(a_ME_HF, jetpt_, jetphi_, jeteta_, ls, "numerator", true, true, true, false);
363  }
364 
365  if (isHEP17(jeteta_, jetphi_)) {
366  FillME(a_ME_HEP17,
367  jetpt_,
368  jetphi_,
369  jeteta_,
370  ls,
371  "numerator",
372  true,
373  false,
374  false,
375  false); // doPhi, doEta, doEtaPhi, doVsLS
376  jetHEP17_AbsEta_.numerator->Fill(abs(jeteta_));
377  } else if (isHEM17(jeteta_, jetphi_)) {
378  FillME(a_ME_HEM17,
379  jetpt_,
380  jetphi_,
381  jeteta_,
382  ls,
383  "numerator",
384  true,
385  false,
386  false,
387  false); // doPhi, doEta, doEtaPhi, doVsLS
388  jetHEM17_AbsEta_.numerator->Fill(abs(jeteta_));
389  } else if (isHEP18(jeteta_, jetphi_)) {
390  FillME(a_ME_HEP18,
391  jetpt_,
392  jetphi_,
393  jeteta_,
394  ls,
395  "numerator",
396  false,
397  false,
398  false,
399  false); // doPhi, doEta, doEtaPhi, doVsLS
400  }
401 }
402 
403 void JetMonitor::fillHistoPSetDescription(edm::ParameterSetDescription& pset) {
404  pset.add<unsigned int>("nbins");
405  pset.add<double>("xmin");
406  pset.add<double>("xmax");
407 }
408 
409 void JetMonitor::fillHistoLSPSetDescription(edm::ParameterSetDescription& pset) {
410  pset.add<unsigned int>("nbins", 2500);
411  pset.add<double>("xmin", 0.);
412  pset.add<double>("xmax", 2500.);
413 }
414 
415 void JetMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
416  edm::ParameterSetDescription desc;
417  desc.add<std::string>("FolderName", "HLT/Jet");
418 
419  desc.add<edm::InputTag>("met", edm::InputTag("pfMet"));
420  // desc.add<edm::InputTag>( "pfjets", edm::InputTag("ak4PFJetsCHS") );
421  // desc.add<edm::InputTag>( "calojets", edm::InputTag("ak4CaloJets") );
422  desc.add<edm::InputTag>("jetSrc", edm::InputTag("ak4PFJetsCHS"));
423  desc.add<edm::InputTag>("electrons", edm::InputTag("gedGsfElectrons"));
424  desc.add<edm::InputTag>("muons", edm::InputTag("muons"));
425  desc.add<int>("njets", 0);
426  desc.add<int>("nelectrons", 0);
427  desc.add<double>("ptcut", 20);
428  desc.add<bool>("ispfjettrg", true);
429  desc.add<bool>("iscalojettrg", false);
430 
431  edm::ParameterSetDescription genericTriggerEventPSet;
432  genericTriggerEventPSet.add<bool>("andOr");
433  genericTriggerEventPSet.add<edm::InputTag>("dcsInputTag", edm::InputTag("scalersRawToDigi"));
434  genericTriggerEventPSet.add<std::vector<int> >("dcsPartitions", {});
435  genericTriggerEventPSet.add<bool>("andOrDcs", false);
436  genericTriggerEventPSet.add<bool>("errorReplyDcs", true);
437  genericTriggerEventPSet.add<std::string>("dbLabel", "");
438  genericTriggerEventPSet.add<bool>("andOrHlt", true);
439  genericTriggerEventPSet.add<edm::InputTag>("hltInputTag", edm::InputTag("TriggerResults::HLT"));
440  genericTriggerEventPSet.add<std::vector<std::string> >("hltPaths", {});
441  // genericTriggerEventPSet.add<std::string>("hltDBKey","");
442  genericTriggerEventPSet.add<bool>("errorReplyHlt", false);
443  genericTriggerEventPSet.add<unsigned int>("verbosityLevel", 1);
444 
445  desc.add<edm::ParameterSetDescription>("numGenericTriggerEventPSet", genericTriggerEventPSet);
446  desc.add<edm::ParameterSetDescription>("denGenericTriggerEventPSet", genericTriggerEventPSet);
447 
448  edm::ParameterSetDescription histoPSet;
449  edm::ParameterSetDescription jetPSet;
450  edm::ParameterSetDescription jetPtThrPSet;
451  fillHistoPSetDescription(jetPSet);
452  fillHistoPSetDescription(jetPtThrPSet);
453  histoPSet.add<edm::ParameterSetDescription>("jetPSet", jetPSet);
454  histoPSet.add<edm::ParameterSetDescription>("jetPtThrPSet", jetPtThrPSet);
455  std::vector<double> bins = {
456  0., 20., 40., 60., 80., 90., 100., 110., 120., 130., 140., 150., 160.,
457  170., 180., 190., 200., 220., 240., 260., 280., 300., 350., 400., 450., 1000.}; // Jet pT Binning
458  histoPSet.add<std::vector<double> >("jetptBinning", bins);
459 
460  edm::ParameterSetDescription lsPSet;
461  fillHistoLSPSetDescription(lsPSet);
462  histoPSet.add<edm::ParameterSetDescription>("lsPSet", lsPSet);
463 
464  desc.add<edm::ParameterSetDescription>("histoPSet", histoPSet);
465 
466  descriptions.add("jetMonitoring", desc);
467 }
468 
469 bool JetMonitor::isBarrel(double eta) {
470  bool output = false;
471  if (fabs(eta) <= 1.3)
472  output = true;
473  return output;
474 }
475 
476 //------------------------------------------------------------------------//
477 bool JetMonitor::isEndCapM(double eta) {
478  bool output = false;
479  if (fabs(eta) <= 3.0 && fabs(eta) > 1.3 && (eta < 0))
480  output = true; // (mia) this magic number should come from some file in CMSSW !!!
481  return output;
482 }
484 bool JetMonitor::isEndCapP(double eta) {
485  bool output = false;
486  //if ( eta<=3.0 && eta >1.3) output=true;
487  if (fabs(eta) <= 3.0 && fabs(eta) > 1.3 && (eta > 0))
488  output = true; // (mia) this magic number should come from some file in CMSSW !!!
489  return output;
490 }
492 bool JetMonitor::isForward(double eta) {
493  bool output = false;
494  if (fabs(eta) > 3.0)
495  output = true;
496  return output;
497 }
499 bool JetMonitor::isHEP17(double eta, double phi) {
500  bool output = false;
501  // phi -0.87 to -0.52
502  if (fabs(eta) <= 3.0 && fabs(eta) > 1.3 && (eta > 0) && phi > -0.87 && phi <= -0.52) {
503  output = true;
504  } // (mia) this magic number should come from some file in CMSSW !!!
505  return output;
506 }
508 bool JetMonitor::isHEM17(double eta, double phi) {
509  bool output = false;
510  if (fabs(eta) <= 3.0 && fabs(eta) > 1.3 && (eta < 0) && phi > -0.87 && phi <= -0.52) {
511  output = true;
512  } // (mia) this magic number should come from some file in CMSSW !!!
513  return output;
514 }
516 bool JetMonitor::isHEP18(double eta, double phi) {
517  bool output = false;
518  // phi -0.87 to -0.52
519  if (fabs(eta) <= 3.0 && fabs(eta) > 1.3 && (eta > 0) && phi > -0.52 && phi <= -0.17) {
520  output = true;
521  } // (mia) this magic number should come from some file in CMSSW !!!
522  return output;
523 }
524 /*void JetMonitor::AutoNullPtr(JetME* a_me,const int len_){
525  for (int i =0; i < len_; ++i)
526  {
527  a_me[i].denominator = nullptr;
528  a_me[i].numerator = nullptr;
529  }
530 }*/
531 void JetMonitor::FillME(JetME* a_me,
532  double pt_,
533  double phi_,
534  double eta_,
535  int ls_,
536  const std::string& denu,
537  bool doPhi,
538  bool doEta,
539  bool doEtaPhi,
540  bool doVsLS) {
541  std::string isDeno = "";
542  isDeno = denu;
543  std::string DenoOrNume = "";
544  DenoOrNume = denu;
545 
546  if (DenoOrNume == "denominator") {
547  // index 0 = pt, 1 = ptThreshold , 2 = pt vs ls , 3 = phi, 4 = eta,
548  // 5 = eta vs phi, 6 = eta vs pt , 7 = abs(eta) , 8 = abs(eta) vs phi
549  a_me[0].denominator->Fill(pt_); // pt
550  a_me[1].denominator->Fill(pt_); // jetpT Threshold binning for pt
551  if (doVsLS)
552  a_me[2].denominator->Fill(ls_, pt_); // pt vs ls
553  if (doPhi)
554  a_me[3].denominator->Fill(phi_); // phi
555  if (doEta)
556  a_me[4].denominator->Fill(eta_); // eta
557  if (doEtaPhi)
558  a_me[5].denominator->Fill(eta_, phi_); // eta vs phi
559  if (doEta)
560  a_me[6].denominator->Fill(eta_, pt_); // eta vs pT
561  } else if (DenoOrNume == "numerator") {
562  a_me[0].numerator->Fill(pt_); // pt
563  a_me[1].numerator->Fill(pt_); // jetpT Threshold binning for pt
564  if (doVsLS)
565  a_me[2].numerator->Fill(ls_, pt_); // pt vs ls
566  if (doPhi)
567  a_me[3].numerator->Fill(phi_); // phi
568  if (doEta)
569  a_me[4].numerator->Fill(eta_); // eat
570  if (doEtaPhi)
571  a_me[5].numerator->Fill(eta_, phi_); // eta vs phi
572  if (doEta)
573  a_me[6].numerator->Fill(eta_, pt_); // eta vs pT
574  } else {
575  edm::LogWarning("JetMonitor") << "CHECK OUT denu option in FillME !!! DenoOrNume ? : " << DenoOrNume << std::endl;
576  }
577 }
578 void JetMonitor::bookMESub(DQMStore::IBooker& Ibooker,
579  JetME* a_me,
580  const int len_,
581  const std::string& h_Name,
582  const std::string& h_Title,
583  const std::string& h_subOptName,
584  std::string h_suOptTitle,
585  bool doPhi,
586  bool doEta,
587  bool doEtaPhi,
588  bool doVsLS) {
589  std::string hName = h_Name;
590  std::string hTitle = h_Title;
591  std::string hSubN = "";
592  std::string hSubT = "";
593  hSubT = std::move(h_suOptTitle);
594 
595  int nbin_phi = jet_phi_binning_.nbins;
596  double maxbin_phi = jet_phi_binning_.xmax;
597  double minbin_phi = jet_phi_binning_.xmin;
598 
599  int nbin_eta = jet_eta_binning_.nbins;
600  double maxbin_eta = jet_eta_binning_.xmax;
601  double minbin_eta = jet_eta_binning_.xmin;
602 
603  if (!h_subOptName.empty()) {
604  hSubN = "_" + h_subOptName;
605  }
606 
607  if (h_subOptName == "HEP17") {
608  nbin_phi = phi_binning_hep17_.nbins;
609  maxbin_phi = phi_binning_hep17_.xmax;
610  minbin_phi = phi_binning_hep17_.xmin;
611 
612  nbin_eta = eta_binning_hep17_.nbins;
613  maxbin_eta = eta_binning_hep17_.xmax;
614  minbin_eta = eta_binning_hep17_.xmin;
615  }
616  if (h_subOptName == "HEM17") {
617  nbin_phi = phi_binning_hep17_.nbins;
618  maxbin_phi = phi_binning_hep17_.xmax;
619  minbin_phi = phi_binning_hep17_.xmin;
620 
621  nbin_eta = eta_binning_hem17_.nbins;
622  maxbin_eta = eta_binning_hem17_.xmax;
623  minbin_eta = eta_binning_hem17_.xmin;
624  }
625  if (h_subOptName == "HEP18") {
626  nbin_phi = phi_binning_hep18_.nbins;
627  maxbin_phi = phi_binning_hep18_.xmax;
628  minbin_phi = phi_binning_hep18_.xmin;
629 
630  nbin_eta = eta_binning_hep17_.nbins;
631  maxbin_eta = eta_binning_hep17_.xmax;
632  minbin_eta = eta_binning_hep17_.xmin;
633  }
634  hName = h_Name + "pT" + hSubN;
635  hTitle = h_Title + " pT " + hSubT;
636  bookME(Ibooker, a_me[0], hName, hTitle, jetpT_binning.nbins, jetpT_binning.xmin, jetpT_binning.xmax);
637  setMETitle(a_me[0], h_Title + " pT [GeV]", "events / [GeV]");
638 
639  hName = h_Name + "pT_pTThresh" + hSubN;
640  hTitle = h_Title + " pT " + hSubT;
641  bookME(Ibooker, a_me[1], hName, hTitle, jetptThr_binning_.nbins, jetptThr_binning_.xmin, jetptThr_binning_.xmax);
642  setMETitle(a_me[1], h_Title + "pT [GeV]", "events / [GeV]");
643 
644  if (doVsLS) {
645  hName = h_Name + "pTVsLS" + hSubN;
646  hTitle = h_Title + " vs LS " + hSubT;
647  bookME(Ibooker,
648  a_me[2],
649  hName,
650  hTitle,
651  ls_binning_.nbins,
652  ls_binning_.xmin,
653  ls_binning_.xmax,
654  jetpT_binning.xmin,
655  jetpT_binning.xmax);
656  setMETitle(a_me[2], "LS", h_Title + "pT [GeV]");
657  }
658 
659  if (doPhi) {
660  hName = h_Name + "phi" + hSubN;
661  hTitle = h_Title + " phi " + hSubT;
662  bookME(Ibooker, a_me[3], hName, hTitle, nbin_phi, minbin_phi, maxbin_phi);
663  setMETitle(a_me[3], h_Title + " #phi", "events / 0.1 rad");
664  }
665 
666  if (doEta) {
667  hName = h_Name + "eta" + hSubN;
668  hTitle = h_Title + " eta " + hSubT;
669  bookME(Ibooker, a_me[4], hName, hTitle, nbin_eta, minbin_eta, maxbin_eta);
670  setMETitle(a_me[4], h_Title + " #eta", "events / #eta");
671  }
672 
673  if (doEtaPhi) {
674  hName = h_Name + "EtaVsPhi" + hSubN;
675  hTitle = h_Title + " eta Vs phi " + hSubT;
676  bookME(Ibooker, a_me[5], hName, hTitle, nbin_eta, minbin_eta, maxbin_eta, nbin_phi, minbin_phi, maxbin_phi);
677  setMETitle(a_me[5], h_Title + " #eta", "#phi");
678  }
679 
680  if (doEta) {
681  hName = h_Name + "EtaVspT" + hSubN;
682  hTitle = h_Title + " eta Vs pT " + hSubT;
683  bookME(Ibooker,
684  a_me[6],
685  hName,
686  hTitle,
687  nbin_eta,
688  minbin_eta,
689  maxbin_eta,
690  jetpT_binning.nbins,
691  jetpT_binning.xmin,
692  jetpT_binning.xmax);
693  setMETitle(a_me[6], h_Title + " #eta", "Leading Jet pT [GeV]");
694  }
695 }
696 // Define this as a plug-in
697 #include "FWCore/Framework/interface/MakerMacros.h"
698 DEFINE_FWK_MODULE(JetMonitor);
dqm::dqmstoreimpl::DQMStore::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX)
Definition: DQMStore.cc:239
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
JetMonitor::isBarrel
bool isBarrel(double eta)
Definition: JetMonitor.cc:469
JetMonitor::a_ME_HB
JetME a_ME_HB[7]
Definition: JetMonitor.h:174
HLT_2018_cff.InputTag
InputTag
Definition: HLT_2018_cff.py:78913
JetMonitor::isHEP17
bool isHEP17(double eta, double phi)
For Hcal HEP17 Area.
Definition: JetMonitor.cc:499
JetMonitor::a_ME_HE_m
JetME a_ME_HE_m[7]
Definition: JetMonitor.h:178
hlt_dqm_clientPB-live_cfg.nbinsY
nbinsY
Definition: hlt_dqm_clientPB-live_cfg.py:55
JetMonitor::jetHEP17_AbsEta_
JetME jetHEP17_AbsEta_
Definition: JetMonitor.h:186
JetMonitor::eta_binning_hem17_
MEbinning eta_binning_hem17_
Definition: JetMonitor.h:230
JetMonitor::a_ME_HEM17
JetME a_ME_HEM17[7]
Definition: JetMonitor.h:179
JetMonitor::phi_binning_hep17_
MEbinning phi_binning_hep17_
Definition: JetMonitor.h:217
LumiMonitor_cff.histoPSet
histoPSet
Definition: LumiMonitor_cff.py:12
JetMonitor::ls_binning_
MEbinning ls_binning_
Definition: JetMonitor.h:171
JetMonitor::JetMonitor
JetMonitor(const edm::ParameterSet &)
Definition: JetMonitor.cc:15
JetMonitor::getHistoPSet
static MEbinning getHistoPSet(const edm::ParameterSet &pset)
Definition: JetMonitor.cc:38
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
L1TOccupancyClient_cfi.ymax
ymax
Definition: L1TOccupancyClient_cfi.py:43
SiStripSourceConfigTier0_cff.genericTriggerEventPSet
genericTriggerEventPSet
Definition: SiStripSourceConfigTier0_cff.py:96
JetMonitor::fillHistoLSPSetDescription
static void fillHistoLSPSetDescription(edm::ParameterSetDescription &pset)
Definition: JetMonitor.cc:409
TrackerOfflineValidation_Dqm_cff.xmin
xmin
Definition: TrackerOfflineValidation_Dqm_cff.py:10
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
JetMonitor::folderName_
std::string folderName_
Definition: JetMonitor.h:157
dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:418
JetMonitor::a_ME_HE_p
JetME a_ME_HE_p[7]
Definition: JetMonitor.h:177
LumiMonitor_cff.lsPSet
lsPSet
Definition: LumiMonitor_cff.py:33
LaserClient_cfi.nbins
nbins
Definition: LaserClient_cfi.py:51
JetMonitor::FillME
void FillME(JetME *a_me, double pt_, double phi_, double eta_, int ls_, const std::string &denu, bool doPhi=true, bool doEta=true, bool doEtaPhi=true, bool doVsLS=true)
Definition: JetMonitor.cc:531
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
convertSQLitetoXML_cfg.output
output
Definition: convertSQLitetoXML_cfg.py:32
GenericTriggerEventFlag
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
Definition: GenericTriggerEventFlag.h:38
JetMonitor::bookME
void bookME(DQMStore::IBooker &, JetME &me, std::string &histname, std::string &histtitle, unsigned int nbins, double xmin, double xmax)
Definition: JetMonitor.cc:57
JetMonitor::a_ME_HEP18
JetME a_ME_HEP18[7]
Definition: JetMonitor.h:181
JetMonitor::isCaloJetTrig
bool isCaloJetTrig
Definition: JetMonitor.h:195
edm::EventID::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: EventID.h:39
JetMonitor::bookMESub
void bookMESub(DQMStore::IBooker &, JetME *a_me, const int len_, const std::string &h_Name, const std::string &h_Title, const std::string &h_subOptName, std::string h_subOptTitle, bool doPhi=true, bool doEta=true, bool doEtaPhi=true, bool doVsLS=true)
Definition: JetMonitor.cc:578
JetMonitor::jetpT_binning
MEbinning jetpT_binning
Definition: JetMonitor.h:169
trigObjTnPSource_cfi.bins
bins
Definition: trigObjTnPSource_cfi.py:20
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:242
JetMonitor::eta_binning_hep17_
MEbinning eta_binning_hep17_
Definition: JetMonitor.h:228
JetMonitor::phi_binning_hep18_
MEbinning phi_binning_hep18_
Definition: JetMonitor.h:223
iEvent
int iEvent
Definition: GenABIO.cc:224
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
JetMonitor::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: JetMonitor.cc:415
JetMonitor::jet_phi_binning_
MEbinning jet_phi_binning_
Definition: JetMonitor.h:205
JetMonitor::den_genTriggerEventFlag_
std::unique_ptr< GenericTriggerEventFlag > den_genTriggerEventFlag_
Definition: JetMonitor.h:190
JetMonitor_cfi.jetPSet
jetPSet
Definition: JetMonitor_cfi.py:11
JetMonitor::isEndCapM
bool isEndCapM(double eta)
Definition: JetMonitor.cc:477
hlt_dqm_clientPB-live_cfg.me
me
Definition: hlt_dqm_clientPB-live_cfg.py:47
JetMonitor::analyze
void analyze(edm::Event const &iEvent, edm::EventSetup const &iSetup) override
Definition: JetMonitor.cc:275
dqm::dqmstoreimpl::DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, char const *option="s")
Definition: DQMStore.cc:333
JetMonitor::jetpT_variable_binning_
std::vector< double > jetpT_variable_binning_
Definition: JetMonitor.h:168
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
hlt_dqm_clientPB-live_cfg.nbinsX
nbinsX
Definition: hlt_dqm_clientPB-live_cfg.py:51
JetMonitor::a_ME_HF
JetME a_ME_HF[7]
Definition: JetMonitor.h:176
JetMonitor::a_ME_HEP17
JetME a_ME_HEP17[7]
Definition: JetMonitor.h:180
min
T min(T a, T b)
Definition: MathUtil.h:58
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
JetMonitor::isHEM17
bool isHEM17(double eta, double phi)
For Hcal HEM17 Area.
Definition: JetMonitor.cc:508
JetMonitor::jetptThr_binning_
MEbinning jetptThr_binning_
Definition: JetMonitor.h:170
JetMonitor::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: JetMonitor.cc:121
JetMonitor::jetHEM17_AbsEta_
JetME jetHEM17_AbsEta_
Definition: JetMonitor.h:187
JetMonitor::jetSrc_
edm::EDGetTokenT< edm::View< reco::Jet > > jetSrc_
Definition: JetMonitor.h:166
eostools.ls
def ls(path, rec=False)
Definition: eostools.py:349
TrackerOfflineValidation_Dqm_cff.xmax
xmax
Definition: TrackerOfflineValidation_Dqm_cff.py:11
JetMonitor::a_ME_HE
JetME a_ME_HE[7]
Definition: JetMonitor.h:175
JetMonitor::v_jetpt
std::vector< double > v_jetpt
Definition: JetMonitor.h:197
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
JetMonitor::num_genTriggerEventFlag_
std::unique_ptr< GenericTriggerEventFlag > num_genTriggerEventFlag_
Definition: JetMonitor.h:189
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:59
edm
HLT enums.
Definition: AlignableModifier.h:19
edm::View::const_iterator
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
dqm::dqmstoreimpl::DQMStore::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Definition: DQMStore.cc:266
JetMonitor::isPFJetTrig
bool isPFJetTrig
Definition: JetMonitor.h:194
JetMonitor::a_ME
JetME a_ME[7]
Definition: JetMonitor.h:173
JetMonitor::~JetMonitor
~JetMonitor() override
JetMonitor::fillHistoPSetDescription
static void fillHistoPSetDescription(edm::ParameterSetDescription &pset)
Definition: JetMonitor.cc:403
JetMonitor::v_jeteta
std::vector< double > v_jeteta
Definition: JetMonitor.h:198
JetMonitor::JetME::denominator
MonitorElement * denominator
Definition: JetMonitor.h:57
JetMonitor::isForward
bool isForward(double eta)
For Hcal Forward Plus Area.
Definition: JetMonitor.cc:492
JetMonitor::v_jetphi
std::vector< double > v_jetphi
Definition: JetMonitor.h:199
eostools.move
def move(src, dest)
Definition: eostools.py:511
JetMonitor::jet_eta_binning_
MEbinning jet_eta_binning_
Definition: JetMonitor.h:211
JetMonitor::getHistoLSPSet
static MEbinning getHistoLSPSet(const edm::ParameterSet &pset)
Definition: JetMonitor.cc:46
edm::Run
Definition: Run.h:45
JetMonitor::setMETitle
void setMETitle(JetME &me, const std::string &titleX, const std::string &titleY)
Definition: JetMonitor.cc:51
JetMonitor::JetME::numerator
MonitorElement * numerator
Definition: JetMonitor.h:56
JetMonitor::isHEP18
bool isHEP18(double eta, double phi)
For Hcal HEP18 Area.
Definition: JetMonitor.cc:516
JetMonitor_cff.jetPtThrPSet
jetPtThrPSet
Definition: JetMonitor_cff.py:9
L1TOccupancyClient_cfi.ymin
ymin
Definition: L1TOccupancyClient_cfi.py:43
JetMonitor::isEndCapP
bool isEndCapP(double eta)
For Hcal Endcap Plus Area.
Definition: JetMonitor.cc:484
JetMonitor::MEbinning::nbins
unsigned int nbins
Definition: JetMonitor.h:50
JetMonitor::ptcut_
double ptcut_
Definition: JetMonitor.h:193
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:813
Generated for CMSSW Reference Manual by doxygen 1.8.11