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 
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, nullptr, _yaxis, nullptr)
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(1+nTTInner, ss.str());
89  ss.str("");
90  ss << "TCC" << outer.first << " TT1";
91  me->setBinLabel(1+2*nTTInner, ss.str());
92  ss.str("");
93  ss << "TCC" << outer.second << " TT1";
94  me->setBinLabel(1+2*nTTInner+nTTOuter, ss.str());
95  // Bins are numbered:
96  // inner1:(1)-->(nTTInner)
97  // inner2:(1+nTTInner)-->(1+nTTInner + nTTInner-1 = 2*nTTInner)
98  // outer1:(1+2*nTTInner)-->(1+2*nTTInner+nTTOuter-1=2*nTTInner+nTTOuter)
99  // outer2:(1+2*nTTInner+nTTOuter)-->(1+2*nTTInner+nTTOuter + nTTOuter-1 = 2*nTTInner+2*nTTOuter)
100  int offset(0);
101  for(int iBin(4); iBin <= (2*nTTOuter + 2*nTTInner); iBin += 4){
102  if(iBin == 4+nTTInner) offset = nTTInner;
103  else if(iBin == 4+2*nTTInner) offset = 2*nTTInner;
104  else if(iBin == 4+2*nTTInner+nTTOuter) offset = 2*nTTInner+nTTOuter;
105  ss.str("");
106  ss << iBin - offset;
107  me->setBinLabel(iBin, ss.str());
108  }
109  }
110  }
111  }
112  }
113 
114  void
115  MESetDet1D::fill(DetId const& _id, double _wy/* = 1.*/, double _w/* = 1.*/, double)
116  {
117  if(!active_) return;
118 
119  unsigned iME(binning::findPlotIndex(otype_, _id));
120  checkME_(iME);
121 
122  binning::ObjectType obj(binning::getObject(otype_, iME));
123  int xbin(binning::findBin1D(obj, btype_, _id));
124 
125  if(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D)
126  fill_(iME, xbin, _wy, _w);
127  else
128  fill_(iME, xbin, _wy);
129  }
130 
131  void
132  MESetDet1D::fill(EcalElectronicsId const& _id, double _wy/* = 1.*/, double _w/* = 1.*/, double)
133  {
134  if(!active_) return;
135 
136  unsigned iME(binning::findPlotIndex(otype_, _id));
137  checkME_(iME);
138 
139  binning::ObjectType obj(binning::getObject(otype_, iME));
140  int xbin(binning::findBin1D(obj, btype_, _id));
141 
142  if(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D)
143  fill_(iME, xbin, _wy, _w);
144  else
145  fill_(iME, xbin, _wy);
146  }
147 
148  void
149  MESetDet1D::fill(int _dcctccid, double _wy/* = 1.*/, double _w/* = 1.*/, double)
150  {
151  if(!active_) return;
152 
153  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
154  checkME_(iME);
155 
156  binning::ObjectType obj(binning::getObject(otype_, iME));
157  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
158 
159  if(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D)
160  fill_(iME, xbin, _wy, _w);
161  else
162  fill_(iME, xbin, _wy);
163  }
164 
165  void
166  MESetDet1D::setBinContent(DetId const& _id, double _content)
167  {
168  if(!active_) return;
169 
170  unsigned iME(binning::findPlotIndex(otype_, _id));
171  checkME_(iME);
172 
173  MonitorElement* me(mes_[iME]);
174 
175  binning::ObjectType obj(binning::getObject(otype_, iME));
176  int xbin(binning::findBin1D(obj, btype_, _id));
177 
178  if(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D){
179  int nbinsY(me->getTH1()->GetNbinsY());
180  for(int iY(1); iY <= nbinsY; iY++)
181  me->setBinContent(xbin, iY, _content);
182  }
183  else
184  me->setBinContent(xbin, _content);
185  }
186 
187  void
188  MESetDet1D::setBinContent(EcalElectronicsId const& _id, double _content)
189  {
190  if(!active_) 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 
200  if(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D){
201  int nbinsY(me->getTH1()->GetNbinsY());
202  for(int iY(1); iY <= nbinsY; iY++)
203  me->setBinContent(xbin, iY, _content);
204  }
205  else
206  me->setBinContent(xbin, _content);
207  }
208 
209  void
210  MESetDet1D::setBinContent(int _dcctccid, double _content)
211  {
212  if(!active_) return;
213 
214  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
215  checkME_(iME);
216 
217  MonitorElement* me(mes_[iME]);
218 
219  binning::ObjectType obj(binning::getObject(otype_, iME));
220  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
221 
222  if(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D){
223  int nbinsY(me->getTH1()->GetNbinsY());
224  for(int iY(1); iY <= nbinsY; iY++)
225  me->setBinContent(xbin, iY, _content);
226  }
227  else
228  me->setBinContent(xbin, _content);
229  }
230 
231  void
232  MESetDet1D::setBinContent(DetId const& _id, int _ybin, double _content)
233  {
234  if(!active_) return;
235  if(kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
236 
237  unsigned iME(binning::findPlotIndex(otype_, _id));
238  checkME_(iME);
239 
240  MonitorElement* me(mes_[iME]);
241 
242  binning::ObjectType obj(binning::getObject(otype_, iME));
243  int xbin(binning::findBin1D(obj, btype_, _id));
244  me->setBinContent(xbin, _ybin, _content);
245  }
246 
247  void
248  MESetDet1D::setBinContent(EcalElectronicsId const& _id, int _ybin, double _content)
249  {
250  if(!active_) return;
251  if(kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
252 
253  unsigned iME(binning::findPlotIndex(otype_, _id));
254  checkME_(iME);
255 
256  MonitorElement* me(mes_[iME]);
257 
258  binning::ObjectType obj(binning::getObject(otype_, iME));
259  int xbin(binning::findBin1D(obj, btype_, _id));
260  me->setBinContent(xbin, _ybin, _content);
261  }
262 
263  void
264  MESetDet1D::setBinContent(int _dcctccid, int _ybin, double _content)
265  {
266  if(!active_) return;
267  if(kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
268 
269  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
270  checkME_(iME);
271 
272  MonitorElement* me(mes_[iME]);
273 
274  binning::ObjectType obj(binning::getObject(otype_, iME));
275  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
276  me->setBinContent(xbin, _ybin, _content);
277  }
278 
279  void
280  MESetDet1D::setBinError(DetId const& _id, double _error)
281  {
282  if(!active_) return;
283 
284  unsigned iME(binning::findPlotIndex(otype_, _id));
285  checkME_(iME);
286 
287  MonitorElement* me(mes_[iME]);
288 
289  binning::ObjectType obj(binning::getObject(otype_, iME));
290  int xbin(binning::findBin1D(obj, btype_, _id));
291 
292  if(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D){
293  int nbinsY(me->getTH1()->GetNbinsY());
294  for(int iY(1); iY <= nbinsY; iY++)
295  me->setBinError(xbin, iY, _error);
296  }
297  else
298  me->setBinError(xbin, _error);
299  }
300 
301  void
302  MESetDet1D::setBinError(EcalElectronicsId const& _id, double _error)
303  {
304  if(!active_) 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 
314  if(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D){
315  int nbinsY(me->getTH1()->GetNbinsY());
316  for(int iY(1); iY <= nbinsY; iY++)
317  me->setBinError(xbin, iY, _error);
318  }
319  else
320  me->setBinError(xbin, _error);
321  }
322 
323  void
324  MESetDet1D::setBinError(int _dcctccid, double _error)
325  {
326  if(!active_) return;
327 
328  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
329  checkME_(iME);
330 
331  MonitorElement* me(mes_[iME]);
332 
333  binning::ObjectType obj(binning::getObject(otype_, iME));
334  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
335 
336  if(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D){
337  int nbinsY(me->getTH1()->GetNbinsY());
338  for(int iY(1); iY <= nbinsY; iY++)
339  me->setBinError(xbin, iY, _error);
340  }
341  else
342  me->setBinError(xbin, _error);
343  }
344 
345  void
346  MESetDet1D::setBinError(DetId const& _id, int _ybin, double _error)
347  {
348  if(!active_) return;
349  if(kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
350 
351  unsigned iME(binning::findPlotIndex(otype_, _id));
352  checkME_(iME);
353 
354  MonitorElement* me(mes_[iME]);
355 
356  binning::ObjectType obj(binning::getObject(otype_, iME));
357  int xbin(binning::findBin1D(obj, btype_, _id));
358  me->setBinError(xbin, _ybin, _error);
359  }
360 
361  void
362  MESetDet1D::setBinError(EcalElectronicsId const& _id, int _ybin, double _error)
363  {
364  if(!active_) return;
365  if(kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
366 
367  unsigned iME(binning::findPlotIndex(otype_, _id));
368  checkME_(iME);
369 
370  MonitorElement* me(mes_[iME]);
371 
372  binning::ObjectType obj(binning::getObject(otype_, iME));
373  int xbin(binning::findBin1D(obj, btype_, _id));
374  me->setBinError(xbin, _ybin, _error);
375  }
376 
377  void
378  MESetDet1D::setBinError(int _dcctccid, int _ybin, double _error)
379  {
380  if(!active_) return;
381  if(kind_ != MonitorElement::DQM_KIND_TH2F && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
382 
383  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
384  checkME_(iME);
385 
386  MonitorElement* me(mes_[iME]);
387 
388  binning::ObjectType obj(binning::getObject(otype_, iME));
389  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
390  me->setBinError(xbin, _ybin, _error);
391  }
392 
393  void
394  MESetDet1D::setBinEntries(DetId const& _id, double _entries)
395  {
396  if(!active_) return;
397  if(kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
398 
399  unsigned iME(binning::findPlotIndex(otype_, _id));
400  checkME_(iME);
401 
402  MonitorElement* me(mes_[iME]);
403 
404  binning::ObjectType obj(binning::getObject(otype_, iME));
405  int xbin(binning::findBin1D(obj, btype_, _id));
406 
407  if(kind_ == MonitorElement::DQM_KIND_TPROFILE2D){
408  int nbinsX(me->getTH1()->GetNbinsX());
409  int nbinsY(me->getTH1()->GetNbinsY());
410  for(int iY(1); iY <= nbinsY; iY++)
411  me->setBinEntries((nbinsX + 2) * iY + xbin, _entries);
412  }
413  else
414  me->setBinEntries(xbin, _entries);
415  }
416 
417  void
418  MESetDet1D::setBinEntries(EcalElectronicsId const& _id, double _entries)
419  {
420  if(!active_) return;
421  if(kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
422 
423  unsigned iME(binning::findPlotIndex(otype_, _id));
424  checkME_(iME);
425 
426  MonitorElement* me(mes_[iME]);
427 
428  binning::ObjectType obj(binning::getObject(otype_, iME));
429  int xbin(binning::findBin1D(obj, btype_, _id));
430 
431  if(kind_ == MonitorElement::DQM_KIND_TPROFILE2D){
432  int nbinsX(me->getTH1()->GetNbinsX());
433  int nbinsY(me->getTH1()->GetNbinsY());
434  for(int iY(1); iY <= nbinsY; iY++)
435  me->setBinEntries((nbinsX + 2) * iY + xbin, _entries);
436  }
437  else
438  me->setBinEntries(xbin, _entries);
439  }
440 
441  void
442  MESetDet1D::setBinEntries(int _dcctccid, double _entries)
443  {
444  if(!active_) return;
445  if(kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
446 
447  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
448  checkME_(iME);
449 
450  MonitorElement* me(mes_[iME]);
451 
452  binning::ObjectType obj(binning::getObject(otype_, iME));
453  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
454 
455  if(kind_ == MonitorElement::DQM_KIND_TPROFILE2D){
456  int nbinsX(me->getTH1()->GetNbinsX());
457  int nbinsY(me->getTH1()->GetNbinsY());
458  for(int iY(1); iY <= nbinsY; iY++)
459  me->setBinEntries((nbinsX + 2) * iY + xbin, _entries);
460  }
461  else
462  me->setBinEntries(xbin, _entries);
463  }
464 
465  void
466  MESetDet1D::setBinEntries(DetId const& _id, int _ybin, double _entries)
467  {
468  if(!active_) return;
469  if(kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
470 
471  unsigned iME(binning::findPlotIndex(otype_, _id));
472  checkME_(iME);
473 
474  MonitorElement* me(mes_[iME]);
475 
476  binning::ObjectType obj(binning::getObject(otype_, iME));
477  int xbin(binning::findBin1D(obj, btype_, _id));
478  int nbinsX(me->getTH1()->GetNbinsX());
479  me->setBinEntries((nbinsX + 2) * _ybin + xbin, _entries);
480  }
481 
482  void
483  MESetDet1D::setBinEntries(EcalElectronicsId const& _id, int _ybin, double _entries)
484  {
485  if(!active_) return;
486  if(kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
487 
488  unsigned iME(binning::findPlotIndex(otype_, _id));
489  checkME_(iME);
490 
491  MonitorElement* me(mes_[iME]);
492 
493  binning::ObjectType obj(binning::getObject(otype_, iME));
494  int xbin(binning::findBin1D(obj, btype_, _id));
495  int nbinsX(me->getTH1()->GetNbinsX());
496  me->setBinEntries((nbinsX + 2) * _ybin + xbin, _entries);
497  }
498 
499  void
500  MESetDet1D::setBinEntries(int _dcctccid, int _ybin, double _entries)
501  {
502  if(!active_) return;
503  if(kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return;
504 
505  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
506  checkME_(iME);
507 
508  MonitorElement* me(mes_[iME]);
509 
510  binning::ObjectType obj(binning::getObject(otype_, iME));
511  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
512  int nbinsX(me->getTH1()->GetNbinsX());
513  me->setBinEntries((nbinsX + 2) * _ybin + xbin, _entries);
514  }
515 
516  double
517  MESetDet1D::getBinContent(DetId const& _id, int _ybin/* = 0*/) const
518  {
519  if(!active_) return 0.;
520 
521  unsigned iME(binning::findPlotIndex(otype_, _id));
522  checkME_(iME);
523 
524  MonitorElement* me(mes_[iME]);
525 
526  binning::ObjectType obj(binning::getObject(otype_, iME));
527  int xbin(binning::findBin1D(obj, btype_, _id));
528  int nbinsX(me->getTH1()->GetNbinsX());
529 
530  return me->getBinContent((nbinsX + 2) * _ybin + xbin);
531  }
532 
533  double
534  MESetDet1D::getBinContent(EcalElectronicsId const& _id, int _ybin/* = 0*/) const
535  {
536  if(!active_) return 0.;
537 
538  unsigned iME(binning::findPlotIndex(otype_, _id));
539  checkME_(iME);
540 
541  MonitorElement* me(mes_[iME]);
542 
543  binning::ObjectType obj(binning::getObject(otype_, iME));
544  int xbin(binning::findBin1D(obj, btype_, _id));
545  int nbinsX(me->getTH1()->GetNbinsX());
546 
547  return me->getBinContent((nbinsX + 2) * _ybin + xbin);
548  }
549 
550  double
551  MESetDet1D::getBinContent(int _dcctccid, int _ybin/* = 0*/) const
552  {
553  if(!active_) return 0.;
554 
555  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
556  checkME_(iME);
557 
558  MonitorElement* me(mes_[iME]);
559 
560  binning::ObjectType obj(binning::getObject(otype_, iME));
561  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
562  int nbinsX(me->getTH1()->GetNbinsX());
563 
564  return me->getBinContent((nbinsX + 2) * _ybin + xbin);
565  }
566 
567  double
568  MESetDet1D::getBinError(DetId const& _id, int _ybin/* = 0*/) const
569  {
570  if(!active_) return 0.;
571 
572  unsigned iME(binning::findPlotIndex(otype_, _id));
573  checkME_(iME);
574 
575  MonitorElement* me(mes_[iME]);
576 
577  binning::ObjectType obj(binning::getObject(otype_, iME));
578  int xbin(binning::findBin1D(obj, btype_, _id));
579  int nbinsX(me->getTH1()->GetNbinsX());
580 
581  return me->getBinError((nbinsX + 2) * _ybin + xbin);
582  }
583 
584  double
585  MESetDet1D::getBinError(EcalElectronicsId const& _id, int _ybin/* = 0*/) const
586  {
587  if(!active_) return 0.;
588 
589  unsigned iME(binning::findPlotIndex(otype_, _id));
590  checkME_(iME);
591 
592  MonitorElement* me(mes_[iME]);
593 
594  binning::ObjectType obj(binning::getObject(otype_, iME));
595  int xbin(binning::findBin1D(obj, btype_, _id));
596  int nbinsX(me->getTH1()->GetNbinsX());
597 
598  return me->getBinError((nbinsX + 2) * _ybin + xbin);
599  }
600 
601  double
602  MESetDet1D::getBinError(int _dcctccid, int _ybin/* = 0*/) const
603  {
604  if(!active_) return 0.;
605 
606  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
607  checkME_(iME);
608 
609  MonitorElement* me(mes_[iME]);
610 
611  binning::ObjectType obj(binning::getObject(otype_, iME));
612  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
613  int nbinsX(me->getTH1()->GetNbinsX());
614 
615  return me->getBinError((nbinsX + 2) * _ybin + xbin);
616  }
617 
618  double
619  MESetDet1D::getBinEntries(DetId const& _id, int _ybin/* = 0*/) const
620  {
621  if(!active_) return 0.;
622  if(kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return 0.;
623 
624  unsigned iME(binning::findPlotIndex(otype_, _id));
625  checkME_(iME);
626 
627  MonitorElement* me(mes_[iME]);
628 
629  binning::ObjectType obj(binning::getObject(otype_, iME));
630  int xbin(binning::findBin1D(obj, btype_, _id));
631  int nbinsX(me->getTH1()->GetNbinsX());
632 
633  return me->getBinEntries((nbinsX + 2) * _ybin + xbin);
634  }
635 
636  double
637  MESetDet1D::getBinEntries(EcalElectronicsId const& _id, int _ybin/* = 0*/) const
638  {
639  if(!active_) return 0.;
640  if(kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return 0.;
641 
642  unsigned iME(binning::findPlotIndex(otype_, _id));
643  checkME_(iME);
644 
645  MonitorElement* me(mes_[iME]);
646 
647  binning::ObjectType obj(binning::getObject(otype_, iME));
648  int xbin(binning::findBin1D(obj, btype_, _id));
649  int nbinsX(me->getTH1()->GetNbinsX());
650 
651  return me->getBinEntries((nbinsX + 2) * _ybin + xbin);
652  }
653 
654  double
655  MESetDet1D::getBinEntries(int _dcctccid, int _ybin/* = 0*/) const
656  {
657  if(!active_) return 0.;
658  if(kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return 0.;
659 
660  unsigned iME(binning::findPlotIndex(otype_, _dcctccid, btype_));
661  checkME_(iME);
662 
663  MonitorElement* me(mes_[iME]);
664 
665  binning::ObjectType obj(binning::getObject(otype_, iME));
666  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
667  int nbinsX(me->getTH1()->GetNbinsX());
668 
669  return me->getBinEntries((nbinsX + 2) * _ybin + xbin);
670  }
671 
672  int
673  MESetDet1D::findBin(DetId const& _id) const
674  {
675  if(!active_) return -1;
676  if(kind_ == MonitorElement::DQM_KIND_TPROFILE || kind_ == MonitorElement::DQM_KIND_TPROFILE2D) return -1;
677 
678  unsigned iME(binning::findPlotIndex(otype_, _id));
679  checkME_(iME);
680 
681  binning::ObjectType obj(binning::getObject(otype_, iME));
682  return binning::findBin1D(obj, btype_, _id);
683  }
684 
685  int
686  MESetDet1D::findBin(EcalElectronicsId const& _id) const
687  {
688  if(!active_) return -1;
689  if(kind_ == MonitorElement::DQM_KIND_TPROFILE || kind_ == MonitorElement::DQM_KIND_TPROFILE2D) return -1;
690 
691  unsigned iME(binning::findPlotIndex(otype_, _id));
692  checkME_(iME);
693 
694  binning::ObjectType obj(binning::getObject(otype_, iME));
695  return binning::findBin1D(obj, btype_, _id);
696  }
697 
698  int
699  MESetDet1D::findBin(int _dcctccid) const
700  {
701  if(!active_) return -1;
702  if(kind_ == MonitorElement::DQM_KIND_TPROFILE || kind_ == MonitorElement::DQM_KIND_TPROFILE2D) return -1;
703 
704  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
705  checkME_(iME);
706 
707  binning::ObjectType obj(binning::getObject(otype_, iME));
708  return binning::findBin1D(obj, btype_, _dcctccid);
709  }
710 
711  int
712  MESetDet1D::findBin(DetId const& _id, double _y, double) const
713  {
714  if(!active_) return -1;
715  if(kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return -1;
716 
717  unsigned iME(binning::findPlotIndex(otype_, _id));
718  checkME_(iME);
719 
720  MonitorElement* me(mes_[iME]);
721 
722  binning::ObjectType obj(binning::getObject(otype_, iME));
723  int xbin(binning::findBin1D(obj, btype_, _id));
724  int nbinsX(me->getTH1()->GetNbinsX());
725  return xbin + (nbinsX + 2) * me->getTH1()->GetYaxis()->FindBin(_y);
726  }
727 
728  int
729  MESetDet1D::findBin(EcalElectronicsId const& _id, double _y, double) const
730  {
731  if(!active_) return -1;
732  if(kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return -1;
733 
734  unsigned iME(binning::findPlotIndex(otype_, _id));
735  checkME_(iME);
736 
737  MonitorElement* me(mes_[iME]);
738 
739  binning::ObjectType obj(binning::getObject(otype_, iME));
740  int xbin(binning::findBin1D(obj, btype_, _id));
741  int nbinsX(me->getTH1()->GetNbinsX());
742  return xbin + (nbinsX + 2) * me->getTH1()->GetYaxis()->FindBin(_y);
743  }
744 
745  int
746  MESetDet1D::findBin(int _dcctccid, double _y, double) const
747  {
748  if(!active_) return -1;
749  if(kind_ != MonitorElement::DQM_KIND_TPROFILE && kind_ != MonitorElement::DQM_KIND_TPROFILE2D) return -1;
750 
751  unsigned iME(binning::findPlotIndex(otype_, _dcctccid));
752  checkME_(iME);
753 
754  MonitorElement* me(mes_[iME]);
755 
756  binning::ObjectType obj(binning::getObject(otype_, iME));
757  int xbin(binning::findBin1D(obj, btype_, _dcctccid));
758  int nbinsX(me->getTH1()->GetNbinsX());
759  return xbin + (nbinsX + 2) * me->getTH1()->GetYaxis()->FindBin(_y);
760  }
761 
762  void
763  MESetDet1D::reset(double _content/* = 0.*/, double _err/* = 0.*/, double _entries/* = 0.*/)
764  {
765  unsigned nME(binning::getNObjects(otype_));
766 
767  bool isProfile(kind_ == MonitorElement::DQM_KIND_TPROFILE || kind_ == MonitorElement::DQM_KIND_TPROFILE2D);
768  bool is2D(kind_ == MonitorElement::DQM_KIND_TH2F || kind_ == MonitorElement::DQM_KIND_TPROFILE2D);
769 
770  for(unsigned iME(0); iME < nME; iME++) {
771  MonitorElement* me(mes_[iME]);
772 
773  int nbinsX(me->getTH1()->GetNbinsX());
774  int nbinsY(me->getTH1()->GetNbinsY());
775  for(int ix(1); ix <= nbinsX; ix++){
776  for(int iy(1); iy <= nbinsY; iy++){
777  int bin(is2D ? (nbinsX + 2) * iy + ix : ix);
778  me->setBinContent(bin, _content);
779  me->setBinError(bin, _err);
780  if(isProfile) me->setBinEntries(bin, _entries);
781  }
782  }
783  }
784  }
785 }
786 
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:866
ecaldqm::MESetDet1D::reset
void reset(double=0., double=0., double=0.) override
Definition: MESetDet1D.cc:763
popcon2dropbox.copy
def copy(args, dbName)
Definition: popcon2dropbox.py:122
ecaldqm::MESet::checkME_
virtual void checkME_(unsigned _iME) const
Definition: MESet.h:115
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:20
ecaldqm::MESetDet1D::getBinContent
double getBinContent(DetId const &, int=0) const override
Definition: MESetDet1D.cc:517
ecaldqm::binning::getNObjects
unsigned getNObjects(ObjectType)
Definition: MESetBinningUtils.cc:476
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:673
ecaldqm::outerTCCs
std::pair< unsigned, unsigned > outerTCCs(unsigned)
Definition: EcalDQMCommonUtils.cc:171
nullptr
#define nullptr
ecaldqm::binning::getObject
ObjectType getObject(ObjectType, unsigned)
Definition: MESetBinningUtils.cc:442
ecaldqm::MESetDet1D::clone
MESet * clone(std::string const &="") const override
Definition: MESetDet1D.cc:32
ecaldqm::MESetDet1D::setBinContent
void setBinContent(DetId const &, double) override
Definition: MESetDet1D.cc:166
ecaldqm::MESet::mes_
std::vector< MonitorElement * > mes_
Definition: MESet.h:129
ecaldqm::MESetDet1D::getBinError
double getBinError(DetId const &, int=0) const override
Definition: MESetDet1D.cc:568
ecaldqm::MESet::throw_
void throw_(std::string const &_message) const
Definition: MESet.h:124
fastTrackerRecHitType::is2D
bool is2D(HitType hitType)
Definition: FastTrackerRecHit.h:31
ecaldqm::MESetDet1D::setBinError
void setBinError(DetId const &, double) override
Definition: MESetDet1D.cc:280
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:8
ecaldqm::binning::findBin1D
int findBin1D(ObjectType, BinningType, DetId const &)
ecaldqm::binning::channelName
std::string channelName(uint32_t, BinningType _btype=kDCC)
Definition: MESetBinningUtils.cc:555
ecaldqm::MESet::fill_
virtual void fill_(unsigned, int, double)
Definition: MESet.cc:278
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:61
GetRecoTauVFromDQM_MC_cff.obj
obj
Definition: GetRecoTauVFromDQM_MC_cff.py:88
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:138
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:42
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:150
ecaldqm::MESet::otype_
binning::ObjectType otype_
Definition: MESet.h:132
ecaldqm::MESetDet1D::~MESetDet1D
~MESetDet1D() override
Definition: MESetDet1D.cc:27
ecaldqm::MESetDet1D::fill
void fill(DetId const &, double=1., double=1., double=0.) override
Definition: MESetDet1D.cc:115
MonitorElement::getNbinsX
int getNbinsX() const
get # of bins in X-axis
Definition: MonitorElement.cc:747
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:918
ecaldqm::MESetDet1D::getBinEntries
double getBinEntries(DetId const &, int=0) const override
Definition: MESetDet1D.cc:619
ecaldqm::MESetDet1D::setBinEntries
void setBinEntries(DetId const &, double) override
Definition: MESetDet1D.cc:394
ecaldqm::MESetEcal::book
void book(DQMStore::IBooker &) override
Definition: MESetEcal.cc:69
btvTracks_cfi._content
_content
Definition: btvTracks_cfi.py:5
Generated for CMSSW Reference Manual by doxygen 1.8.11