CMS 3D CMS Logo
MESetDet1D.cc
Go to the documentation of this file.
1 #include "DQM/EcalCommon/interface/MESetDet1D.h"
2 
3 #include "DQM/EcalCommon/interface/EcalDQMCommonUtils.h"
4 
5 namespace ecaldqm {
6 
7  MESetDet1D::MESetDet1D(std::string const &_fullPath,
8  binning::ObjectType _otype,
9  binning::BinningType _btype,
10  MonitorElement::Kind _kind,
11  binning::AxisSpecs const *_yaxis /* = 0*/)
12  : MESetEcal(_fullPath, _otype, _btype, _kind, 1, nullptr, _yaxis, nullptr) {
13  switch (kind_) {
14  case MonitorElement::DQM_KIND_TH1F:
15  case MonitorElement::DQM_KIND_TPROFILE:
16  case MonitorElement::DQM_KIND_TH2F:
17  case MonitorElement::DQM_KIND_TPROFILE2D:
18  break;
19  default:
20  throw_("Unsupported MonitorElement kind");
21  }
22  }
23 
24  MESetDet1D::MESetDet1D(MESetDet1D const &_orig) : MESetEcal(_orig) {}
25 
26  MESetDet1D::~MESetDet1D() {}
27 
28  MESet *MESetDet1D::clone(std::string const &_path /* = ""*/) const {
29  std::string path(path_);
30  if (!_path.empty())
31  path_ = _path;
32  MESet *copy(new MESetDet1D(*this));
33  path_ = path;
34  return copy;
35  }
36 
37  void MESetDet1D::book(DQMStore::IBooker &_ibooker) {
38  MESetEcal::book(_ibooker);
39 
40  if (btype_ == binning::kDCC) {
41  for (unsigned iME(0); iME < mes_.size(); iME++) {
42  MonitorElement *me(mes_[iME]);
43 
44  binning::ObjectType actualObject(binning::getObject(otype_, iME));
45  if (actualObject == binning::kEB) {
46  for (int iBin(1); iBin <= me->getNbinsX(); iBin++)
47  me->setBinLabel(iBin, binning::channelName(iBin + kEBmLow));
48  } else if (actualObject == binning::kEE) {
49  for (int iBin(1); iBin <= me->getNbinsX() / 2; iBin++) {
50  me->setBinLabel(iBin, binning::channelName(iBin));
51  me->setBinLabel(iBin + me->getNbinsX() / 2, binning::channelName(iBin + 45));
52  }
53  } else if (actualObject == binning::kEEm) {
54  for (int iBin(1); iBin <= me->getNbinsX(); iBin++)
55  me->setBinLabel(iBin, binning::channelName(iBin));
56  } else if (actualObject == binning::kEEp) {
57  for (int iBin(1); iBin <= me->getNbinsX(); iBin++)
58  me->setBinLabel(iBin, binning::channelName(iBin + 45));
59  }
60  }
61  } else if (btype_ == binning::kTriggerTower) {
62  for (unsigned iME(0); iME < mes_.size(); iME++) {
63  MonitorElement *me(mes_[iME]);
64 
65  binning::ObjectType actualObject(binning::getObject(otype_, iME));
66  unsigned dccid(0);
67  if (actualObject == binning::kSM && (iME <= kEEmHigh || iME >= kEEpLow))
68  dccid = iME + 1;
69  else if (actualObject == binning::kEESM)
70  dccid = iME <= kEEmHigh ? iME + 1 : iME + 37;
71 
72  if (dccid > 0) {
73  std::stringstream ss;
74  std::pair<unsigned, unsigned> inner(innerTCCs(iME + 1));
75  std::pair<unsigned, unsigned> outer(outerTCCs(iME + 1));
76  ss << "TCC" << inner.first << " TT1";
77  me->setBinLabel(1, ss.str());
78  ss.str("");
79  ss << "TCC" << inner.second << " TT1";
80  me->setBinLabel(1 + nTTInner, ss.str());
81  ss.str("");
82  ss << "TCC" << outer.first << " TT1";
83  me->setBinLabel(1 + 2 * nTTInner, ss.str());
84  ss.str("");
85  ss << "TCC" << outer.second << " TT1";
86  me->setBinLabel(1 + 2 * nTTInner + nTTOuter, ss.str());
87  // Bins are numbered:
88  // inner1:(1)-->(nTTInner)
89  // inner2:(1+nTTInner)-->(1+nTTInner + nTTInner-1 = 2*nTTInner)
90  // outer1:(1+2*nTTInner)-->(1+2*nTTInner+nTTOuter-1=2*nTTInner+nTTOuter)
91  // outer2:(1+2*nTTInner+nTTOuter)-->(1+2*nTTInner+nTTOuter + nTTOuter-1
92  // = 2*nTTInner+2*nTTOuter)
93  int offset(0);
94  for (int iBin(4); iBin <= (2 * nTTOuter + 2 * nTTInner); iBin += 4) {
95  if (iBin == 4 + nTTInner)
96  offset = nTTInner;
97  else if (iBin == 4 + 2 * nTTInner)
98  offset = 2 * nTTInner;
99  else if (iBin == 4 + 2 * nTTInner + nTTOuter)
100  offset = 2 * nTTInner + nTTOuter;
101  ss.str("");
102  ss << iBin - offset;
103  me->setBinLabel(iBin, ss.str());
104  }
105  }
106  }
107  }
108  }
109 
110  void MESetDet1D::fill(DetId const &_id, double _wy /* = 1.*/, double _w /* = 1.*/, double) {
111  if (!active_)
112  return;
113 
114  unsigned iME(binning::findPlotIndex(otype_, _id));
115  checkME_(iME);
116 
117  binning::ObjectType obj(binning::getObject(otype_, iME));
118  int xbin(binning::findBin1D(obj, btype_, _id));
119 
120  if (kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D)
121  fill_(iME, xbin, _wy, _w);
122  else
123  fill_(iME, xbin, _wy);
124  }
125 
126  void MESetDet1D::fill(EcalElectronicsId const &_id, double _wy /* = 1.*/, double _w /* = 1.*/, double) {
127  if (!active_)
128  return;
129 
130  unsigned iME(binning::findPlotIndex(otype_, _id));
131  checkME_(iME);
132 
133  binning::ObjectType obj(binning::getObject(otype_, iME));
134  int xbin(binning::findBin1D(obj, btype_, _id));
135 
136  if (kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D)
137  fill_(iME, xbin, _wy, _w);
138  else
139  fill_(iME, xbin, _wy);
140  }
141 
142  void MESetDet1D::fill(int _dcctccid, double _wy /* = 1.*/, double _w /* = 1.*/, double) {
143  if (!active_)
144  return;
145 
146  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
147  checkME_(iME);
148 
149  binning::ObjectType obj(binning::getObject(otype_, iME));
150  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
151 
152  if (kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D)
153  fill_(iME, xbin, _wy, _w);
154  else
155  fill_(iME, xbin, _wy);
156  }
157 
158  void MESetDet1D::setBinContent(DetId const &_id, double _content) {
159  if (!active_)
160  return;
161 
162  unsigned iME(binning::findPlotIndex(otype_, _id));
163  checkME_(iME);
164 
165  MonitorElement *me(mes_[iME]);
166 
167  binning::ObjectType obj(binning::getObject(otype_, iME));
168  int xbin(binning::findBin1D(obj, btype_, _id));
169 
170  if (kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D) {
171  int nbinsY(me->getTH1()->GetNbinsY());
172  for (int iY(1); iY <= nbinsY; iY++)
173  me->setBinContent(xbin, iY, _content);
174  } else
175  me->setBinContent(xbin, _content);
176  }
177 
178  void MESetDet1D::setBinContent(EcalElectronicsId const &_id, double _content) {
179  if (!active_)
180  return;
181 
182  unsigned iME(binning::findPlotIndex(otype_, _id));
183  checkME_(iME);
184 
185  MonitorElement *me(mes_[iME]);
186 
187  binning::ObjectType obj(binning::getObject(otype_, iME));
188  int xbin(binning::findBin1D(obj, btype_, _id));
189 
190  if (kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D) {
191  int nbinsY(me->getTH1()->GetNbinsY());
192  for (int iY(1); iY <= nbinsY; iY++)
193  me->setBinContent(xbin, iY, _content);
194  } else
195  me->setBinContent(xbin, _content);
196  }
197 
198  void MESetDet1D::setBinContent(int _dcctccid, double _content) {
199  if (!active_)
200  return;
201 
202  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
203  checkME_(iME);
204 
205  MonitorElement *me(mes_[iME]);
206 
207  binning::ObjectType obj(binning::getObject(otype_, iME));
208  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
209 
210  if (kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D) {
211  int nbinsY(me->getTH1()->GetNbinsY());
212  for (int iY(1); iY <= nbinsY; iY++)
213  me->setBinContent(xbin, iY, _content);
214  } else
215  me->setBinContent(xbin, _content);
216  }
217 
218  void MESetDet1D::setBinContent(DetId const &_id, int _ybin, double _content) {
219  if (!active_)
220  return;
221  if (kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
222  return;
223 
224  unsigned iME(binning::findPlotIndex(otype_, _id));
225  checkME_(iME);
226 
227  MonitorElement *me(mes_[iME]);
228 
229  binning::ObjectType obj(binning::getObject(otype_, iME));
230  int xbin(binning::findBin1D(obj, btype_, _id));
231  me->setBinContent(xbin, _ybin, _content);
232  }
233 
234  void MESetDet1D::setBinContent(EcalElectronicsId const &_id, int _ybin, double _content) {
235  if (!active_)
236  return;
237  if (kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
238  return;
239 
240  unsigned iME(binning::findPlotIndex(otype_, _id));
241  checkME_(iME);
242 
243  MonitorElement *me(mes_[iME]);
244 
245  binning::ObjectType obj(binning::getObject(otype_, iME));
246  int xbin(binning::findBin1D(obj, btype_, _id));
247  me->setBinContent(xbin, _ybin, _content);
248  }
249 
250  void MESetDet1D::setBinContent(int _dcctccid, int _ybin, double _content) {
251  if (!active_)
252  return;
253  if (kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
254  return;
255 
256  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
257  checkME_(iME);
258 
259  MonitorElement *me(mes_[iME]);
260 
261  binning::ObjectType obj(binning::getObject(otype_, iME));
262  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
263  me->setBinContent(xbin, _ybin, _content);
264  }
265 
266  void MESetDet1D::setBinError(DetId const &_id, double _error) {
267  if (!active_)
268  return;
269 
270  unsigned iME(binning::findPlotIndex(otype_, _id));
271  checkME_(iME);
272 
273  MonitorElement *me(mes_[iME]);
274 
275  binning::ObjectType obj(binning::getObject(otype_, iME));
276  int xbin(binning::findBin1D(obj, btype_, _id));
277 
278  if (kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D) {
279  int nbinsY(me->getTH1()->GetNbinsY());
280  for (int iY(1); iY <= nbinsY; iY++)
281  me->setBinError(xbin, iY, _error);
282  } else
283  me->setBinError(xbin, _error);
284  }
285 
286  void MESetDet1D::setBinError(EcalElectronicsId const &_id, double _error) {
287  if (!active_)
288  return;
289 
290  unsigned iME(binning::findPlotIndex(otype_, _id));
291  checkME_(iME);
292 
293  MonitorElement *me(mes_[iME]);
294 
295  binning::ObjectType obj(binning::getObject(otype_, iME));
296  int xbin(binning::findBin1D(obj, btype_, _id));
297 
298  if (kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D) {
299  int nbinsY(me->getTH1()->GetNbinsY());
300  for (int iY(1); iY <= nbinsY; iY++)
301  me->setBinError(xbin, iY, _error);
302  } else
303  me->setBinError(xbin, _error);
304  }
305 
306  void MESetDet1D::setBinError(int _dcctccid, double _error) {
307  if (!active_)
308  return;
309 
310  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
311  checkME_(iME);
312 
313  MonitorElement *me(mes_[iME]);
314 
315  binning::ObjectType obj(binning::getObject(otype_, iME));
316  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
317 
318  if (kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D) {
319  int nbinsY(me->getTH1()->GetNbinsY());
320  for (int iY(1); iY <= nbinsY; iY++)
321  me->setBinError(xbin, iY, _error);
322  } else
323  me->setBinError(xbin, _error);
324  }
325 
326  void MESetDet1D::setBinError(DetId const &_id, int _ybin, double _error) {
327  if (!active_)
328  return;
329  if (kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
330  return;
331 
332  unsigned iME(binning::findPlotIndex(otype_, _id));
333  checkME_(iME);
334 
335  MonitorElement *me(mes_[iME]);
336 
337  binning::ObjectType obj(binning::getObject(otype_, iME));
338  int xbin(binning::findBin1D(obj, btype_, _id));
339  me->setBinError(xbin, _ybin, _error);
340  }
341 
342  void MESetDet1D::setBinError(EcalElectronicsId const &_id, int _ybin, double _error) {
343  if (!active_)
344  return;
345  if (kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
346  return;
347 
348  unsigned iME(binning::findPlotIndex(otype_, _id));
349  checkME_(iME);
350 
351  MonitorElement *me(mes_[iME]);
352 
353  binning::ObjectType obj(binning::getObject(otype_, iME));
354  int xbin(binning::findBin1D(obj, btype_, _id));
355  me->setBinError(xbin, _ybin, _error);
356  }
357 
358  void MESetDet1D::setBinError(int _dcctccid, int _ybin, double _error) {
359  if (!active_)
360  return;
361  if (kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
362  return;
363 
364  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
365  checkME_(iME);
366 
367  MonitorElement *me(mes_[iME]);
368 
369  binning::ObjectType obj(binning::getObject(otype_, iME));
370  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
371  me->setBinError(xbin, _ybin, _error);
372  }
373 
374  void MESetDet1D::setBinEntries(DetId const &_id, double _entries) {
375  if (!active_)
376  return;
377  if (kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
378  return;
379 
380  unsigned iME(binning::findPlotIndex(otype_, _id));
381  checkME_(iME);
382 
383  MonitorElement *me(mes_[iME]);
384 
385  binning::ObjectType obj(binning::getObject(otype_, iME));
386  int xbin(binning::findBin1D(obj, btype_, _id));
387 
388  if (kind_ == MonitorElement::DQM_KIND_TPROFILE2D) {
389  int nbinsX(me->getTH1()->GetNbinsX());
390  int nbinsY(me->getTH1()->GetNbinsY());
391  for (int iY(1); iY <= nbinsY; iY++)
392  me->setBinEntries((nbinsX + 2) * iY + xbin, _entries);
393  } else
394  me->setBinEntries(xbin, _entries);
395  }
396 
397  void MESetDet1D::setBinEntries(EcalElectronicsId const &_id, double _entries) {
398  if (!active_)
399  return;
400  if (kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
401  return;
402 
403  unsigned iME(binning::findPlotIndex(otype_, _id));
404  checkME_(iME);
405 
406  MonitorElement *me(mes_[iME]);
407 
408  binning::ObjectType obj(binning::getObject(otype_, iME));
409  int xbin(binning::findBin1D(obj, btype_, _id));
410 
411  if (kind_ == MonitorElement::DQM_KIND_TPROFILE2D) {
412  int nbinsX(me->getTH1()->GetNbinsX());
413  int nbinsY(me->getTH1()->GetNbinsY());
414  for (int iY(1); iY <= nbinsY; iY++)
415  me->setBinEntries((nbinsX + 2) * iY + xbin, _entries);
416  } else
417  me->setBinEntries(xbin, _entries);
418  }
419 
420  void MESetDet1D::setBinEntries(int _dcctccid, double _entries) {
421  if (!active_)
422  return;
423  if (kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
424  return;
425 
426  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
427  checkME_(iME);
428 
429  MonitorElement *me(mes_[iME]);
430 
431  binning::ObjectType obj(binning::getObject(otype_, iME));
432  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
433 
434  if (kind_ == MonitorElement::DQM_KIND_TPROFILE2D) {
435  int nbinsX(me->getTH1()->GetNbinsX());
436  int nbinsY(me->getTH1()->GetNbinsY());
437  for (int iY(1); iY <= nbinsY; iY++)
438  me->setBinEntries((nbinsX + 2) * iY + xbin, _entries);
439  } else
440  me->setBinEntries(xbin, _entries);
441  }
442 
443  void MESetDet1D::setBinEntries(DetId const &_id, int _ybin, double _entries) {
444  if (!active_)
445  return;
446  if (kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
447  return;
448 
449  unsigned iME(binning::findPlotIndex(otype_, _id));
450  checkME_(iME);
451 
452  MonitorElement *me(mes_[iME]);
453 
454  binning::ObjectType obj(binning::getObject(otype_, iME));
455  int xbin(binning::findBin1D(obj, btype_, _id));
456  int nbinsX(me->getTH1()->GetNbinsX());
457  me->setBinEntries((nbinsX + 2) * _ybin + xbin, _entries);
458  }
459 
460  void MESetDet1D::setBinEntries(EcalElectronicsId const &_id, int _ybin, double _entries) {
461  if (!active_)
462  return;
463  if (kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
464  return;
465 
466  unsigned iME(binning::findPlotIndex(otype_, _id));
467  checkME_(iME);
468 
469  MonitorElement *me(mes_[iME]);
470 
471  binning::ObjectType obj(binning::getObject(otype_, iME));
472  int xbin(binning::findBin1D(obj, btype_, _id));
473  int nbinsX(me->getTH1()->GetNbinsX());
474  me->setBinEntries((nbinsX + 2) * _ybin + xbin, _entries);
475  }
476 
477  void MESetDet1D::setBinEntries(int _dcctccid, int _ybin, double _entries) {
478  if (!active_)
479  return;
480  if (kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
481  return;
482 
483  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
484  checkME_(iME);
485 
486  MonitorElement *me(mes_[iME]);
487 
488  binning::ObjectType obj(binning::getObject(otype_, iME));
489  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
490  int nbinsX(me->getTH1()->GetNbinsX());
491  me->setBinEntries((nbinsX + 2) * _ybin + xbin, _entries);
492  }
493 
494  double MESetDet1D::getBinContent(DetId const &_id, int _ybin /* = 0*/) const {
495  if (!active_)
496  return 0.;
497 
498  unsigned iME(binning::findPlotIndex(otype_, _id));
499  checkME_(iME);
500 
501  MonitorElement *me(mes_[iME]);
502 
503  binning::ObjectType obj(binning::getObject(otype_, iME));
504  int xbin(binning::findBin1D(obj, btype_, _id));
505  int nbinsX(me->getTH1()->GetNbinsX());
506 
507  return me->getBinContent((nbinsX + 2) * _ybin + xbin);
508  }
509 
510  double MESetDet1D::getBinContent(EcalElectronicsId const &_id, int _ybin /* = 0*/) const {
511  if (!active_)
512  return 0.;
513 
514  unsigned iME(binning::findPlotIndex(otype_, _id));
515  checkME_(iME);
516 
517  MonitorElement *me(mes_[iME]);
518 
519  binning::ObjectType obj(binning::getObject(otype_, iME));
520  int xbin(binning::findBin1D(obj, btype_, _id));
521  int nbinsX(me->getTH1()->GetNbinsX());
522 
523  return me->getBinContent((nbinsX + 2) * _ybin + xbin);
524  }
525 
526  double MESetDet1D::getBinContent(int _dcctccid, int _ybin /* = 0*/) const {
527  if (!active_)
528  return 0.;
529 
530  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
531  checkME_(iME);
532 
533  MonitorElement *me(mes_[iME]);
534 
535  binning::ObjectType obj(binning::getObject(otype_, iME));
536  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
537  int nbinsX(me->getTH1()->GetNbinsX());
538 
539  return me->getBinContent((nbinsX + 2) * _ybin + xbin);
540  }
541 
542  double MESetDet1D::getBinError(DetId const &_id, int _ybin /* = 0*/) const {
543  if (!active_)
544  return 0.;
545 
546  unsigned iME(binning::findPlotIndex(otype_, _id));
547  checkME_(iME);
548 
549  MonitorElement *me(mes_[iME]);
550 
551  binning::ObjectType obj(binning::getObject(otype_, iME));
552  int xbin(binning::findBin1D(obj, btype_, _id));
553  int nbinsX(me->getTH1()->GetNbinsX());
554 
555  return me->getBinError((nbinsX + 2) * _ybin + xbin);
556  }
557 
558  double MESetDet1D::getBinError(EcalElectronicsId const &_id, int _ybin /* = 0*/) const {
559  if (!active_)
560  return 0.;
561 
562  unsigned iME(binning::findPlotIndex(otype_, _id));
563  checkME_(iME);
564 
565  MonitorElement *me(mes_[iME]);
566 
567  binning::ObjectType obj(binning::getObject(otype_, iME));
568  int xbin(binning::findBin1D(obj, btype_, _id));
569  int nbinsX(me->getTH1()->GetNbinsX());
570 
571  return me->getBinError((nbinsX + 2) * _ybin + xbin);
572  }
573 
574  double MESetDet1D::getBinError(int _dcctccid, int _ybin /* = 0*/) const {
575  if (!active_)
576  return 0.;
577 
578  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
579  checkME_(iME);
580 
581  MonitorElement *me(mes_[iME]);
582 
583  binning::ObjectType obj(binning::getObject(otype_, iME));
584  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
585  int nbinsX(me->getTH1()->GetNbinsX());
586 
587  return me->getBinError((nbinsX + 2) * _ybin + xbin);
588  }
589 
590  double MESetDet1D::getBinEntries(DetId const &_id, int _ybin /* = 0*/) const {
591  if (!active_)
592  return 0.;
593  if (kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
594  return 0.;
595 
596  unsigned iME(binning::findPlotIndex(otype_, _id));
597  checkME_(iME);
598 
599  MonitorElement *me(mes_[iME]);
600 
601  binning::ObjectType obj(binning::getObject(otype_, iME));
602  int xbin(binning::findBin1D(obj, btype_, _id));
603  int nbinsX(me->getTH1()->GetNbinsX());
604 
605  return me->getBinEntries((nbinsX + 2) * _ybin + xbin);
606  }
607 
608  double MESetDet1D::getBinEntries(EcalElectronicsId const &_id, int _ybin /* = 0*/) const {
609  if (!active_)
610  return 0.;
611  if (kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
612  return 0.;
613 
614  unsigned iME(binning::findPlotIndex(otype_, _id));
615  checkME_(iME);
616 
617  MonitorElement *me(mes_[iME]);
618 
619  binning::ObjectType obj(binning::getObject(otype_, iME));
620  int xbin(binning::findBin1D(obj, btype_, _id));
621  int nbinsX(me->getTH1()->GetNbinsX());
622 
623  return me->getBinEntries((nbinsX + 2) * _ybin + xbin);
624  }
625 
626  double MESetDet1D::getBinEntries(int _dcctccid, int _ybin /* = 0*/) const {
627  if (!active_)
628  return 0.;
629  if (kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
630  return 0.;
631 
632  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
633  checkME_(iME);
634 
635  MonitorElement *me(mes_[iME]);
636 
637  binning::ObjectType obj(binning::getObject(otype_, iME));
638  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
639  int nbinsX(me->getTH1()->GetNbinsX());
640 
641  return me->getBinEntries((nbinsX + 2) * _ybin + xbin);
642  }
643 
644  int MESetDet1D::findBin(DetId const &_id) const {
645  if (!active_)
646  return -1;
647  if (kind_ == MonitorElement::DQM_KIND_TPROFILE || kind_ == MonitorElement::DQM_KIND_TPROFILE2D)
648  return -1;
649 
650  unsigned iME(binning::findPlotIndex(otype_, _id));
651  checkME_(iME);
652 
653  binning::ObjectType obj(binning::getObject(otype_, iME));
654  return binning::findBin1D(obj, btype_, _id);
655  }
656 
657  int MESetDet1D::findBin(EcalElectronicsId const &_id) const {
658  if (!active_)
659  return -1;
660  if (kind_ == MonitorElement::DQM_KIND_TPROFILE || kind_ == MonitorElement::DQM_KIND_TPROFILE2D)
661  return -1;
662 
663  unsigned iME(binning::findPlotIndex(otype_, _id));
664  checkME_(iME);
665 
666  binning::ObjectType obj(binning::getObject(otype_, iME));
667  return binning::findBin1D(obj, btype_, _id);
668  }
669 
670  int MESetDet1D::findBin(int _dcctccid) const {
671  if (!active_)
672  return -1;
673  if (kind_ == MonitorElement::DQM_KIND_TPROFILE || kind_ == MonitorElement::DQM_KIND_TPROFILE2D)
674  return -1;
675 
676  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
677  checkME_(iME);
678 
679  binning::ObjectType obj(binning::getObject(otype_, iME));
680  return binning::findBin1D(obj, btype_, _dcctccid);
681  }
682 
683  int MESetDet1D::findBin(DetId const &_id, double _y, double) const {
684  if (!active_)
685  return -1;
686  if (kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
687  return -1;
688 
689  unsigned iME(binning::findPlotIndex(otype_, _id));
690  checkME_(iME);
691 
692  MonitorElement *me(mes_[iME]);
693 
694  binning::ObjectType obj(binning::getObject(otype_, iME));
695  int xbin(binning::findBin1D(obj, btype_, _id));
696  int nbinsX(me->getTH1()->GetNbinsX());
697  return xbin + (nbinsX + 2) * me->getTH1()->GetYaxis()->FindBin(_y);
698  }
699 
700  int MESetDet1D::findBin(EcalElectronicsId const &_id, double _y, double) const {
701  if (!active_)
702  return -1;
703  if (kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
704  return -1;
705 
706  unsigned iME(binning::findPlotIndex(otype_, _id));
707  checkME_(iME);
708 
709  MonitorElement *me(mes_[iME]);
710 
711  binning::ObjectType obj(binning::getObject(otype_, iME));
712  int xbin(binning::findBin1D(obj, btype_, _id));
713  int nbinsX(me->getTH1()->GetNbinsX());
714  return xbin + (nbinsX + 2) * me->getTH1()->GetYaxis()->FindBin(_y);
715  }
716 
717  int MESetDet1D::findBin(int _dcctccid, double _y, double) const {
718  if (!active_)
719  return -1;
720  if (kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D)
721  return -1;
722 
723  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
724  checkME_(iME);
725 
726  MonitorElement *me(mes_[iME]);
727 
728  binning::ObjectType obj(binning::getObject(otype_, iME));
729  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
730  int nbinsX(me->getTH1()->GetNbinsX());
731  return xbin + (nbinsX + 2) * me->getTH1()->GetYaxis()->FindBin(_y);
732  }
733 
734  void MESetDet1D::reset(double _content /* = 0.*/, double _err /* = 0.*/, double _entries /* = 0.*/) {
735  unsigned nME(binning::getNObjects(otype_));
736 
737  bool isProfile(kind_ == MonitorElement::DQM_KIND_TPROFILE || kind_ == MonitorElement::DQM_KIND_TPROFILE2D);
738  bool is2D(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D);
739 
740  for (unsigned iME(0); iME < nME; iME++) {
741  MonitorElement *me(mes_[iME]);
742 
743  int nbinsX(me->getTH1()->GetNbinsX());
744  int nbinsY(me->getTH1()->GetNbinsY());
745  for (int ix(1); ix <= nbinsX; ix++) {
746  for (int iy(1); iy <= nbinsY; iy++) {
747  int bin(is2D ? (nbinsX + 2) * iy + ix : ix);
748  me->setBinContent(bin, _content);
749  me->setBinError(bin, _err);
750  if (isProfile)
751  me->setBinEntries(bin, _entries);
752  }
753  }
754  }
755  }
756 } // namespace ecaldqm
ecaldqm::MESet::kind_
MonitorElement::Kind kind_
Definition: MESet.h:129
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:866
ecaldqm::MESetDet1D::reset
void reset(double=0., double=0., double=0.) override
Definition: MESetDet1D.cc:734
popcon2dropbox.copy
def copy(args, dbName)
Definition: popcon2dropbox.py:122
ecaldqm::MESet::checkME_
virtual void checkME_(unsigned _iME) const
Definition: MESet.h:114
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
EcalElectronicsId
Ecal readout channel identification [32:20] Unused (so far) [19:13] DCC id [12:6] tower [5:3] strip [...
Definition: EcalElectronicsId.h:18
ecaldqm::MESetDet1D::getBinContent
double getBinContent(DetId const &, int=0) const override
Definition: MESetDet1D.cc:494
ecaldqm::binning::getNObjects
unsigned getNObjects(ObjectType)
Definition: MESetBinningUtils.cc:511
nullptr
#define nullptr
Definition: GCC11Compatibility.h:37
MonitorElement::getTH1
TH1 * getTH1() const
Definition: MonitorElement.cc:1476
MonitorElement::setBinLabel
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:942
ecaldqm::MESetDet1D::findBin
int findBin(DetId const &) const
Definition: MESetDet1D.cc:644
ecaldqm::outerTCCs
std::pair< unsigned, unsigned > outerTCCs(unsigned)
Definition: EcalDQMCommonUtils.cc:153
ecaldqm::binning::getObject
ObjectType getObject(ObjectType, unsigned)
Definition: MESetBinningUtils.cc:466
ecaldqm::MESetDet1D::clone
MESet * clone(std::string const &="") const override
Definition: MESetDet1D.cc:28
ecaldqm::MESetDet1D::setBinContent
void setBinContent(DetId const &, double) override
Definition: MESetDet1D.cc:158
ecaldqm::MESetDet1D::getBinError
double getBinError(DetId const &, int=0) const override
Definition: MESetDet1D.cc:542
ecaldqm::MESet::throw_
void throw_(std::string const &_message) const
Definition: MESet.h:122
fastTrackerRecHitType::is2D
bool is2D(HitType hitType)
Definition: FastTrackerRecHit.h:31
ecaldqm::MESetDet1D::setBinError
void setBinError(DetId const &, double) override
Definition: MESetDet1D.cc:266
photons_cff._id
_id
Definition: photons_cff.py:52
MonitorElement::setBinError
void setBinError(int binx, double error)
set uncertainty on content of bin (1-D)
Definition: MonitorElement.cc:891
ecaldqm::MESetDet1D::MESetDet1D
MESetDet1D(std::string const &, binning::ObjectType, binning::BinningType, MonitorElement::Kind, binning::AxisSpecs const *=0)
Definition: MESetDet1D.cc:7
ecaldqm::binning::findBin1D
int findBin1D(ObjectType, BinningType, DetId const &)
ecaldqm::binning::channelName
std::string channelName(uint32_t, BinningType _btype=kDCC)
Definition: MESetBinningUtils.cc:590
ecaldqm::MESet::fill_
virtual void fill_(unsigned, int, double)
Definition: MESet.cc:266
stringResolutionProvider_cfi.bin
bin
set the eta bin as selection string.
Definition: stringResolutionProvider_cfi.py:11
DetId
Definition: DetId.h:18
PFRecoTauDiscriminationByIsolation_cfi.offset
offset
Definition: PFRecoTauDiscriminationByIsolation_cfi.py:62
ecaldqm::binning::findPlotIndex
unsigned findPlotIndex(ObjectType, DetId const &)
MonitorElement::getBinError
double getBinError(int binx) const
get uncertainty on content of bin (1-D) - See TH1::GetBinError for details
Definition: MonitorElement.cc:783
ecaldqm::MESet::active_
bool active_
Definition: MESet.h:134
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:765
ecaldqm::MESetDet1D::book
void book(DQMStore::IBooker &) override
Definition: MESetDet1D.cc:37
MonitorElement::getBinEntries
double getBinEntries(int bin) const
get # of bin entries (for profiles)
Definition: MonitorElement.cc:807
ecaldqm::innerTCCs
std::pair< unsigned, unsigned > innerTCCs(unsigned)
Definition: EcalDQMCommonUtils.cc:133
ecaldqm::MESet::otype_
binning::ObjectType otype_
Definition: MESet.h:127
ecaldqm::MESetDet1D::~MESetDet1D
~MESetDet1D() override
Definition: MESetDet1D.cc:26
ecaldqm::MESetDet1D::fill
void fill(DetId const &, double=1., double=1., double=0.) override
Definition: MESetDet1D.cc:110
MonitorElement::getNbinsX
int getNbinsX() const
get # of bins in X-axis
Definition: MonitorElement.cc:747
ecaldqm::MESet::btype_
binning::BinningType btype_
Definition: MESet.h:128
ecaldqm::MESet::path_
std::string path_
Definition: MESet.h:126
MonitorElement::setBinEntries
void setBinEntries(int bin, double nentries)
set # of bin entries (to be used for profiles)
Definition: MonitorElement.cc:918
ecaldqm::MESetDet1D::getBinEntries
double getBinEntries(DetId const &, int=0) const override
Definition: MESetDet1D.cc:590
ecaldqm::MESetDet1D::setBinEntries
void setBinEntries(DetId const &, double) override
Definition: MESetDet1D.cc:374
ecaldqm::MESet::mes_
std::vector< MonitorElement * > mes_
Definition: MESet.h:124
ecaldqm::MESetEcal::book
void book(DQMStore::IBooker &) override
Definition: MESetEcal.cc:70
btvTracks_cfi._content
_content
Definition: btvTracks_cfi.py:5
Generated for CMSSW Reference Manual by doxygen 1.8.11