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