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