CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
MonitorElement.cc
Go to the documentation of this file.
1 #define __STDC_FORMAT_MACROS 1
2 #define DQM_ROOT_METHODS 1
3 #include "DQMServices/Core/interface/MonitorElement.h"
4 #include "DQMServices/Core/interface/QTest.h"
5 #include "DQMServices/Core/src/DQMError.h"
6 #include "TClass.h"
7 #include "TMath.h"
8 #include "TList.h"
9 #include "THashList.h"
10 #include <iostream>
11 #include <cassert>
12 #include <cfloat>
13 #include <inttypes.h>
14 
15 static TH1 *
16 checkRootObject(const std::string &name, TObject *tobj, const char *func, int reqdim)
17 {
18  if (! tobj)
19  raiseDQMError("MonitorElement", "Method '%s' cannot be invoked on monitor"
20  " element '%s' because it is not a ROOT object.",
21  func, name.c_str());
22 
23  TH1 *h = static_cast<TH1 *>(tobj);
24  int ndim = h->GetDimension();
25  if (reqdim < 0 || reqdim > ndim)
26  raiseDQMError("MonitorElement", "Method '%s' cannot be invoked on monitor"
27  " element '%s' because it requires %d dimensions; this"
28  " object of type '%s' has %d dimensions",
29  func, name.c_str(), reqdim, typeid(*h).name(), ndim);
30 
31  return h;
32 }
33 
34 MonitorElement *
35 MonitorElement::initialise(Kind kind)
36 {
37  switch (kind)
38  {
39  case DQM_KIND_INT:
40  case DQM_KIND_REAL:
41  case DQM_KIND_STRING:
42  case DQM_KIND_TH1F:
43  case DQM_KIND_TH1S:
44  case DQM_KIND_TH1D:
45  case DQM_KIND_TH2F:
46  case DQM_KIND_TH2S:
47  case DQM_KIND_TH2D:
48  case DQM_KIND_TH3F:
49  case DQM_KIND_TPROFILE:
50  case DQM_KIND_TPROFILE2D:
51  data_.flags &= ~DQMNet::DQM_PROP_TYPE_MASK;
52  data_.flags |= kind;
53  break;
54 
55  default:
56  raiseDQMError("MonitorElement", "cannot initialise monitor element"
57  " to invalid type %d", (int) kind);
58  }
59 
60  return this;
61 }
62 
63 MonitorElement *
64 MonitorElement::initialise(Kind kind, TH1 *rootobj)
65 {
66  initialise(kind);
67  switch (kind)
68  {
69  case DQM_KIND_TH1F:
70  assert(dynamic_cast<TH1F *>(rootobj));
71  assert(! reference_ || dynamic_cast<TH1F *>(reference_));
72  object_ = rootobj;
73  break;
74 
75  case DQM_KIND_TH1S:
76  assert(dynamic_cast<TH1S *>(rootobj));
77  assert(! reference_ || dynamic_cast<TH1S *>(reference_));
78  object_ = rootobj;
79  break;
80 
81  case DQM_KIND_TH1D:
82  assert(dynamic_cast<TH1D *>(rootobj));
83  assert(! reference_ || dynamic_cast<TH1D *>(reference_));
84  object_ = rootobj;
85  break;
86 
87  case DQM_KIND_TH2F:
88  assert(dynamic_cast<TH2F *>(rootobj));
89  assert(! reference_ || dynamic_cast<TH2F *>(reference_));
90  object_ = rootobj;
91  break;
92 
93  case DQM_KIND_TH2S:
94  assert(dynamic_cast<TH2S *>(rootobj));
95  assert(! reference_ || dynamic_cast<TH2S *>(reference_));
96  object_ = rootobj;
97  break;
98 
99  case DQM_KIND_TH2D:
100  assert(dynamic_cast<TH2D *>(rootobj));
101  assert(! reference_ || dynamic_cast<TH1D *>(reference_));
102  object_ = rootobj;
103  break;
104 
105  case DQM_KIND_TH3F:
106  assert(dynamic_cast<TH3F *>(rootobj));
107  assert(! reference_ || dynamic_cast<TH3F *>(reference_));
108  object_ = rootobj;
109  break;
110 
111  case DQM_KIND_TPROFILE:
112  assert(dynamic_cast<TProfile *>(rootobj));
113  assert(! reference_ || dynamic_cast<TProfile *>(reference_));
114  object_ = rootobj;
115  break;
116 
117  case DQM_KIND_TPROFILE2D:
118  assert(dynamic_cast<TProfile2D *>(rootobj));
119  assert(! reference_ || dynamic_cast<TProfile2D *>(reference_));
120  object_ = rootobj;
121  break;
122 
123  default:
124  raiseDQMError("MonitorElement", "cannot initialise monitor element"
125  " as a root object with type %d", (int) kind);
126  }
127 
128  if (reference_)
129  data_.flags |= DQMNet::DQM_PROP_HAS_REFERENCE;
130 
131  return this;
132 }
133 
134 MonitorElement *
135 MonitorElement::initialise(Kind kind, const std::string &value)
136 {
137  initialise(kind);
138  if (kind == DQM_KIND_STRING)
139  scalar_.str = value;
140  else
141  raiseDQMError("MonitorElement", "cannot initialise monitor element"
142  " as a string with type %d", (int) kind);
143 
144  return this;
145 }
146 
147 MonitorElement::MonitorElement(void)
148  : object_(0),
149  reference_(0),
150  refvalue_(0)
151 {
152  data_.version = 0;
153  data_.dirname = 0;
154  data_.run = 0;
155  data_.lumi = 0;
156  data_.streamId = 0;
157  data_.moduleId = 0;
158  data_.tag = 0;
159  data_.flags = DQM_KIND_INVALID | DQMNet::DQM_PROP_NEW;
160  scalar_.num = 0;
161  scalar_.real = 0;
162 }
163 
164 MonitorElement::MonitorElement(const std::string *path,
165  const std::string &name,
166  uint32_t run /* = 0 */,
167  uint32_t streamId /* = 0 */,
168  uint32_t moduleId /* = 0 */)
169  : object_(0),
170  reference_(0),
171  refvalue_(0)
172 {
173  data_.version = 0;
174  data_.run = run;
175  data_.lumi = 0;
176  data_.streamId = streamId;
177  data_.moduleId = moduleId;
178  data_.dirname = path;
179  data_.objname = name;
180  data_.tag = 0;
181  data_.flags = DQM_KIND_INVALID | DQMNet::DQM_PROP_NEW;
182  scalar_.num = 0;
183  scalar_.real = 0;
184 }
185 
186 MonitorElement::MonitorElement(const MonitorElement &x)
187  : data_(x.data_),
188  scalar_(x.scalar_),
189  object_(x.object_),
190  reference_(x.reference_),
191  refvalue_(x.refvalue_),
192  qreports_(x.qreports_)
193 {
194  if (object_)
195  object_ = static_cast<TH1 *>(object_->Clone());
196 
197  if (refvalue_)
198  refvalue_ = static_cast<TH1 *>(refvalue_->Clone());
199 }
200 
201 MonitorElement &
202 MonitorElement::operator=(const MonitorElement &x)
203 {
204  if (this != &x)
205  {
206  delete object_;
207  delete refvalue_;
208 
209  data_ = x.data_;
210  scalar_ = x.scalar_;
211  object_ = x.object_;
212  reference_ = x.reference_;
213  refvalue_ = x.refvalue_;
214  qreports_ = x.qreports_;
215 
216  if (object_)
217  object_ = static_cast<TH1 *>(object_->Clone());
218 
219  if (refvalue_)
220  refvalue_ = static_cast<TH1 *>(refvalue_->Clone());
221  }
222 
223  return *this;
224 }
225 
226 MonitorElement::~MonitorElement(void)
227 {
228  delete object_;
229  delete refvalue_;
230 }
231 
232 //utility function to check the consistency of the axis labels
233 //taken from TH1::CheckBinLabels which is not public
234 bool
235 MonitorElement::CheckBinLabels(const TAxis* a1, const TAxis * a2)
236 {
237  // check that axis have same labels
238  THashList *l1 = (const_cast<TAxis*>(a1))->GetLabels();
239  THashList *l2 = (const_cast<TAxis*>(a2))->GetLabels();
240 
241  if (!l1 && !l2 )
242  return true;
243  if (!l1 || !l2 ) {
244  return false;
245  }
246  // check now labels sizes are the same
247  if (l1->GetSize() != l2->GetSize() ) {
248  return false;
249  }
250  for (int i = 1; i <= a1->GetNbins(); ++i) {
251  TString label1 = a1->GetBinLabel(i);
252  TString label2 = a2->GetBinLabel(i);
253  if (label1 != label2) {
254  return false;
255  }
256  }
257  return true;
258 }
259 
261 void
262 MonitorElement::Fill(std::string &value)
263 {
264  update();
265  if (kind() == DQM_KIND_STRING)
266  scalar_.str = value;
267  else
268  incompatible(__PRETTY_FUNCTION__);
269 }
270 
272 void
273 MonitorElement::Fill(double x)
274 {
275  update();
276  if (kind() == DQM_KIND_INT)
277  scalar_.num = static_cast<int64_t>(x);
278  else if (kind() == DQM_KIND_REAL)
279  scalar_.real = x;
280  else if (kind() == DQM_KIND_TH1F)
281  accessRootObject(__PRETTY_FUNCTION__, 1)
282  ->Fill(x, 1);
283  else if (kind() == DQM_KIND_TH1S)
284  accessRootObject(__PRETTY_FUNCTION__, 1)
285  ->Fill(x, 1);
286  else if (kind() == DQM_KIND_TH1D)
287  accessRootObject(__PRETTY_FUNCTION__, 1)
288  ->Fill(x, 1);
289  else
290  incompatible(__PRETTY_FUNCTION__);
291 }
292 
294 void
295 MonitorElement::doFill(int64_t x)
296 {
297  update();
298  if (kind() == DQM_KIND_INT)
299  scalar_.num = static_cast<int64_t>(x);
300  else if (kind() == DQM_KIND_REAL)
301  scalar_.real = static_cast<double>(x);
302  else if (kind() == DQM_KIND_TH1F)
303  accessRootObject(__PRETTY_FUNCTION__, 1)
304  ->Fill(static_cast<double>(x), 1);
305  else if (kind() == DQM_KIND_TH1S)
306  accessRootObject(__PRETTY_FUNCTION__, 1)
307  ->Fill(static_cast<double>(x), 1);
308  else if (kind() == DQM_KIND_TH1D)
309  accessRootObject(__PRETTY_FUNCTION__, 1)
310  ->Fill(static_cast<double>(x), 1);
311  else
312  incompatible(__PRETTY_FUNCTION__);
313 }
314 
316 void
317 MonitorElement::Fill(double x, double yw)
318 {
319  update();
320  if (kind() == DQM_KIND_TH1F)
321  accessRootObject(__PRETTY_FUNCTION__, 1)
322  ->Fill(x, yw);
323  else if (kind() == DQM_KIND_TH1S)
324  accessRootObject(__PRETTY_FUNCTION__, 1)
325  ->Fill(x, yw);
326  else if (kind() == DQM_KIND_TH1D)
327  accessRootObject(__PRETTY_FUNCTION__, 1)
328  ->Fill(x, yw);
329  else if (kind() == DQM_KIND_TH2F)
330  static_cast<TH2F *>(accessRootObject(__PRETTY_FUNCTION__, 2))
331  ->Fill(x, yw, 1);
332  else if (kind() == DQM_KIND_TH2S)
333  static_cast<TH2S *>(accessRootObject(__PRETTY_FUNCTION__, 2))
334  ->Fill(x, yw, 1);
335  else if (kind() == DQM_KIND_TH2D)
336  static_cast<TH2D *>(accessRootObject(__PRETTY_FUNCTION__, 2))
337  ->Fill(x, yw, 1);
338  else if (kind() == DQM_KIND_TPROFILE)
339  static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
340  ->Fill(x, yw, 1);
341  else
342  incompatible(__PRETTY_FUNCTION__);
343 }
344 
348 void
349 MonitorElement::ShiftFillLast(double y, double ye, int xscale)
350 {
351  update();
352  if (kind() == DQM_KIND_TH1F
353  || kind() == DQM_KIND_TH1S
354  || kind() == DQM_KIND_TH1D)
355  {
356  int nbins = getNbinsX();
357  int entries = (int)getEntries();
358  // first fill bins from left to right
359  int index = entries + 1 ;
360  int xlow = 2 ; int xup = nbins ;
361  // if more entries than bins then start shifting
362  if ( entries >= nbins )
363  {
364  index = nbins;
365  xlow = entries - nbins + 3 ; xup = entries+1 ;
366  // average first bin
367  double y1 = getBinContent(1);
368  double y2 = getBinContent(2);
369  double y1err = getBinError(1);
370  double y2err = getBinError(2);
371  double N = entries - nbins + 1.;
372  if ( ye == 0. || y1err == 0. || y2err == 0.)
373  {
374  // for errors zero calculate unweighted mean and its error
375  double sum = N*y1 + y2;
376  y1 = sum/(N+1.) ;
377  // FIXME check if correct
378  double s=(N+1.)*(N*y1*y1 + y2*y2) - sum*sum;
379  if (s>=0.)
380  y1err = sqrt(s)/(N+1.);
381  else
382  y1err = 0.;
383  }
384  else
385  {
386  // for errors non-zero calculate weighted mean and its error
387  double denom = (1./y1err + 1./y2err);
388  double mean = (y1/y1err + y2/y2err)/denom;
389  // FIXME check if correct
390  y1err = sqrt(((y1-mean)*(y1-mean)/y1err +
391  (y2-mean)*(y2-mean)/y2err)/denom/2.);
392  y1 = mean; // set y1 to mean for filling below
393  }
394  setBinContent(1,y1);
395  setBinError(1,y1err);
396  // shift remaining bins to the left
397  for ( int i = 3; i <= nbins ; i++)
398  {
399  setBinContent(i-1,getBinContent(i));
400  setBinError(i-1,getBinError(i));
401  }
402  }
403  // fill last bin with new values
404  setBinContent(index,y);
405  setBinError(index,ye);
406  // set entries
407  setEntries(entries+1);
408  // set axis labels and reset drawing option
409  char buffer [10];
410  sprintf (buffer, "%d", xlow*xscale);
411  std::string a(buffer); setBinLabel(2,a);
412  sprintf (buffer, "%d", xup*xscale);
413  std::string b(buffer); setBinLabel(nbins,b);
414  setBinLabel(1,"av.");
415  }
416  else
417  incompatible(__PRETTY_FUNCTION__);
418 }
420 void
421 MonitorElement::Fill(double x, double y, double zw)
422 {
423  update();
424  if (kind() == DQM_KIND_TH2F)
425  static_cast<TH2F *>(accessRootObject(__PRETTY_FUNCTION__, 2))
426  ->Fill(x, y, zw);
427  else if (kind() == DQM_KIND_TH2S)
428  static_cast<TH2S *>(accessRootObject(__PRETTY_FUNCTION__, 2))
429  ->Fill(x, y, zw);
430  else if (kind() == DQM_KIND_TH2D)
431  static_cast<TH2D *>(accessRootObject(__PRETTY_FUNCTION__, 2))
432  ->Fill(x, y, zw);
433  else if (kind() == DQM_KIND_TH3F)
434  static_cast<TH3F *>(accessRootObject(__PRETTY_FUNCTION__, 2))
435  ->Fill(x, y, zw, 1);
436  else if (kind() == DQM_KIND_TPROFILE)
437  static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 2))
438  ->Fill(x, y, zw);
439  else if (kind() == DQM_KIND_TPROFILE2D)
440  static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 2))
441  ->Fill(x, y, zw, 1);
442  else
443  incompatible(__PRETTY_FUNCTION__);
444 }
445 
447 void
448 MonitorElement::Fill(double x, double y, double z, double w)
449 {
450  update();
451  if (kind() == DQM_KIND_TH3F)
452  static_cast<TH3F *>(accessRootObject(__PRETTY_FUNCTION__, 2))
453  ->Fill(x, y, z, w);
454  else if (kind() == DQM_KIND_TPROFILE2D)
455  static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 2))
456  ->Fill(x, y, z, w);
457  else
458  incompatible(__PRETTY_FUNCTION__);
459 }
460 
462 void
463 MonitorElement::Reset(void)
464 {
465  update();
466  if (kind() == DQM_KIND_INT)
467  scalar_.num = 0;
468  else if (kind() == DQM_KIND_REAL)
469  scalar_.real = 0;
470  else if (kind() == DQM_KIND_STRING)
471  scalar_.str.clear();
472  else
473  return accessRootObject(__PRETTY_FUNCTION__, 1)
474  ->Reset();
475 }
476 
478 void
479 MonitorElement::packScalarData(std::string &into, const char *prefix) const
480 {
481  char buf[64];
482  if (kind() == DQM_KIND_INT)
483  {
484  snprintf(buf, sizeof(buf), "%s%" PRId64, prefix, scalar_.num);
485  into = buf;
486  }
487  else if (kind() == DQM_KIND_REAL)
488  {
489  snprintf(buf, sizeof(buf), "%s%.*g", prefix, DBL_DIG+2, scalar_.real);
490  into = buf;
491  }
492  else if (kind() == DQM_KIND_STRING)
493  {
494  into.reserve(strlen(prefix) + scalar_.str.size());
495  into += prefix;
496  into += scalar_.str;
497  }
498  else
499  incompatible(__PRETTY_FUNCTION__);
500 }
501 
503 void
504 MonitorElement::packQualityData(std::string &into) const
505 {
506  DQMNet::packQualityData(into, data_.qreports);
507 }
508 
511 std::string
512 MonitorElement::valueString(void) const
513 {
514  std::string result;
515  if (kind() == DQM_KIND_INT)
516  packScalarData(result, "i=");
517  else if (kind() == DQM_KIND_REAL)
518  packScalarData(result, "f=");
519  else if (kind() == DQM_KIND_STRING)
520  packScalarData(result, "s=");
521  else
522  incompatible(__PRETTY_FUNCTION__);
523 
524  return result;
525 }
526 
530 std::string
531 MonitorElement::tagString(void) const
532 {
533  std::string result;
534  std::string val(valueString());
535  result.reserve(6 + 2*data_.objname.size() + val.size());
536  result += '<'; result += data_.objname; result += '>';
537  result += val;
538  result += '<'; result += '/'; result += data_.objname; result += '>';
539  return result;
540 }
541 
543 std::string
544 MonitorElement::tagLabelString(void) const
545 {
546  char buf[32];
547  std::string result;
548  size_t len = sprintf(buf, "t=%" PRIu32, data_.tag);
549 
550  result.reserve(6 + 2*data_.objname.size() + len);
551  result += '<'; result += data_.objname; result += '>';
552  result += buf;
553  result += '<'; result += '/'; result += data_.objname; result += '>';
554  return result;
555 }
556 
558 std::string
559 MonitorElement::effLabelString(void) const
560 {
561  std::string result;
562 
563  result.reserve(6 + 2*data_.objname.size() + 3);
564  result += '<'; result += data_.objname; result += '>';
565  result += "e=1";
566  result += '<'; result += '/'; result += data_.objname; result += '>';
567  return result;
568 }
569 
570 std::string
571 MonitorElement::qualityTagString(const DQMNet::QValue &qv) const
572 {
573  char buf[64];
574  std::string result;
575  size_t titlelen = data_.objname.size() + qv.qtname.size() + 1;
576  size_t buflen = sprintf(buf, "qr=st:%d:%.*g:", qv.code, DBL_DIG+2, qv.qtresult);
577 
578  result.reserve(7 + 2*titlelen + buflen + qv.algorithm.size() + qv.message.size());
579  result += '<'; result += data_.objname; result += '.'; result += qv.qtname; result += '>';
580  result += buf; result += qv.algorithm; result += ':'; result += qv.message;
581  result += '<'; result += '/'; result += data_.objname; result += '.'; result += qv.qtname; result += '>';
582  return result;
583 }
584 
585 const QReport *
586 MonitorElement::getQReport(const std::string &qtname) const
587 {
588  QReport *qr;
589  DQMNet::QValue *qv;
590  const_cast<MonitorElement *>(this)->getQReport(false, qtname, qr, qv);
591  return qr;
592 }
593 
594 std::vector<QReport *>
595 MonitorElement::getQReports(void) const
596 {
597  std::vector<QReport *> result;
598  result.reserve(qreports_.size());
599  for (size_t i = 0, e = qreports_.size(); i != e; ++i)
600  {
601  const_cast<MonitorElement *>(this)->qreports_[i].qvalue_
602  = const_cast<DQMNet::QValue *>(&data_.qreports[i]);
603  result.push_back(const_cast<QReport *>(&qreports_[i]));
604  }
605  return result;
606 }
607 
608 std::vector<QReport *>
609 MonitorElement::getQWarnings(void) const
610 {
611  std::vector<QReport *> result;
612  result.reserve(qreports_.size());
613  for (size_t i = 0, e = qreports_.size(); i != e; ++i)
614  if (data_.qreports[i].code == dqm::qstatus::WARNING)
615  {
616  const_cast<MonitorElement *>(this)->qreports_[i].qvalue_
617  = const_cast<DQMNet::QValue *>(&data_.qreports[i]);
618  result.push_back(const_cast<QReport *>(&qreports_[i]));
619  }
620  return result;
621 }
622 
623 std::vector<QReport *>
624 MonitorElement::getQErrors(void) const
625 {
626  std::vector<QReport *> result;
627  result.reserve(qreports_.size());
628  for (size_t i = 0, e = qreports_.size(); i != e; ++i)
629  if (data_.qreports[i].code == dqm::qstatus::ERROR)
630  {
631  const_cast<MonitorElement *>(this)->qreports_[i].qvalue_
632  = const_cast<DQMNet::QValue *>(&data_.qreports[i]);
633  result.push_back(const_cast<QReport *>(&qreports_[i]));
634  }
635  return result;
636 }
637 
638 std::vector<QReport *>
639 MonitorElement::getQOthers(void) const
640 {
641  std::vector<QReport *> result;
642  result.reserve(qreports_.size());
643  for (size_t i = 0, e = qreports_.size(); i != e; ++i)
644  if (data_.qreports[i].code != dqm::qstatus::STATUS_OK
645  && data_.qreports[i].code != dqm::qstatus::WARNING
646  && data_.qreports[i].code != dqm::qstatus::ERROR)
647  {
648  const_cast<MonitorElement *>(this)->qreports_[i].qvalue_
649  = const_cast<DQMNet::QValue *>(&data_.qreports[i]);
650  result.push_back(const_cast<QReport *>(&qreports_[i]));
651  }
652  return result;
653 }
654 
656 void
657 MonitorElement::runQTests(void)
658 {
659  assert(qreports_.size() == data_.qreports.size());
660 
661  // Rerun quality tests where the ME or the quality algorithm was modified.
662  bool dirty = wasUpdated();
663  for (size_t i = 0, e = data_.qreports.size(); i < e; ++i)
664  {
665  DQMNet::QValue &qv = data_.qreports[i];
666  QReport &qr = qreports_[i];
667  QCriterion *qc = qr.qcriterion_;
668  qr.qvalue_ = &qv;
669 
670  // if (qc && (dirty || qc->wasModified())) // removed for new QTest (abm-090503)
671  if (qc && dirty)
672  {
673  assert(qc->getName() == qv.qtname);
674  std::string oldMessage = qv.message;
675  int oldStatus = qv.code;
676 
677  qc->runTest(this, qr, qv);
678 
679  if (oldStatus != qv.code || oldMessage != qv.message)
680  update();
681  }
682  }
683 
684  // Update QReport statistics.
685  updateQReportStats();
686 }
687 
688 void
689 MonitorElement::incompatible(const char *func) const
690 {
691  raiseDQMError("MonitorElement", "Method '%s' cannot be invoked on monitor"
692  " element '%s'", func, data_.objname.c_str());
693 }
694 
695 TH1 *
696 MonitorElement::accessRootObject(const char *func, int reqdim) const
697 {
698  if (kind() < DQM_KIND_TH1F)
699  raiseDQMError("MonitorElement", "Method '%s' cannot be invoked on monitor"
700  " element '%s' because it is not a root object",
701  func, data_.objname.c_str());
702 
703  return checkRootObject(data_.objname, object_, func, reqdim);
704 }
705 
706 /*** getter methods (wrapper around ROOT methods) ****/
707 //
709 double
710 MonitorElement::getMean(int axis /* = 1 */) const
711 { return accessRootObject(__PRETTY_FUNCTION__, axis-1)
712  ->GetMean(axis); }
713 
716 double
717 MonitorElement::getMeanError(int axis /* = 1 */) const
718 { return accessRootObject(__PRETTY_FUNCTION__, axis-1)
719  ->GetMeanError(axis); }
720 
722 double
723 MonitorElement::getRMS(int axis /* = 1 */) const
724 { return accessRootObject(__PRETTY_FUNCTION__, axis-1)
725  ->GetRMS(axis); }
726 
728 double
729 MonitorElement::getRMSError(int axis /* = 1 */) const
730 { return accessRootObject(__PRETTY_FUNCTION__, axis-1)
731  ->GetRMSError(axis); }
732 
734 int
735 MonitorElement::getNbinsX(void) const
736 { return accessRootObject(__PRETTY_FUNCTION__, 1)
737  ->GetNbinsX(); }
738 
740 int
741 MonitorElement::getNbinsY(void) const
742 { return accessRootObject(__PRETTY_FUNCTION__, 2)
743  ->GetNbinsY(); }
744 
746 int
747 MonitorElement::getNbinsZ(void) const
748 { return accessRootObject(__PRETTY_FUNCTION__, 3)
749  ->GetNbinsZ(); }
750 
752 double
753 MonitorElement::getBinContent(int binx) const
754 { return accessRootObject(__PRETTY_FUNCTION__, 1)
755  ->GetBinContent(binx); }
756 
758 double
759 MonitorElement::getBinContent(int binx, int biny) const
760 { return accessRootObject(__PRETTY_FUNCTION__, 2)
761  ->GetBinContent(binx, biny); }
762 
764 double
765 MonitorElement::getBinContent(int binx, int biny, int binz) const
766 { return accessRootObject(__PRETTY_FUNCTION__, 3)
767  ->GetBinContent(binx, biny, binz); }
768 
770 double
771 MonitorElement::getBinError(int binx) const
772 { return accessRootObject(__PRETTY_FUNCTION__, 1)
773  ->GetBinError(binx); }
774 
776 double
777 MonitorElement::getBinError(int binx, int biny) const
778 { return accessRootObject(__PRETTY_FUNCTION__, 2)
779  ->GetBinError(binx, biny); }
780 
782 double
783 MonitorElement::getBinError(int binx, int biny, int binz) const
784 { return accessRootObject(__PRETTY_FUNCTION__, 3)
785  ->GetBinError(binx, biny, binz); }
786 
788 double
789 MonitorElement::getEntries(void) const
790 { return accessRootObject(__PRETTY_FUNCTION__, 1)
791  ->GetEntries(); }
792 
794 double
795 MonitorElement::getBinEntries(int bin) const
796 {
797  if (kind() == DQM_KIND_TPROFILE)
798  return static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
799  ->GetBinEntries(bin);
800  else if (kind() == DQM_KIND_TPROFILE2D)
801  return static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 1))
802  ->GetBinEntries(bin);
803  else
804  {
805  incompatible(__PRETTY_FUNCTION__);
806  return 0;
807  }
808 }
809 
811 double
812 MonitorElement::getYmin(void) const
813 {
814  if (kind() == DQM_KIND_TPROFILE)
815  return static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
816  ->GetYmin();
817  else
818  {
819  incompatible(__PRETTY_FUNCTION__);
820  return 0;
821  }
822 }
823 
825 double
826 MonitorElement::getYmax(void) const
827 {
828  if (kind() == DQM_KIND_TPROFILE)
829  return static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
830  ->GetYmax();
831  else
832  {
833  incompatible(__PRETTY_FUNCTION__);
834  return 0;
835  }
836 }
837 
839 std::string
840 MonitorElement::getAxisTitle(int axis /* = 1 */) const
841 { return getAxis(__PRETTY_FUNCTION__, axis)
842  ->GetTitle(); }
843 
845 std::string
846 MonitorElement::getTitle(void) const
847 { return accessRootObject(__PRETTY_FUNCTION__, 1)
848  ->GetTitle(); }
849 
850 /*** setter methods (wrapper around ROOT methods) ****/
851 //
853 void
854 MonitorElement::setBinContent(int binx, double content)
855 {
856  update();
857  accessRootObject(__PRETTY_FUNCTION__, 1)
858  ->SetBinContent(binx, content);
859 }
860 
862 void
863 MonitorElement::setBinContent(int binx, int biny, double content)
864 {
865  update();
866  accessRootObject(__PRETTY_FUNCTION__, 2)
867  ->SetBinContent(binx, biny, content); }
868 
870 void
871 MonitorElement::setBinContent(int binx, int biny, int binz, double content)
872 {
873  update();
874  accessRootObject(__PRETTY_FUNCTION__, 3)
875  ->SetBinContent(binx, biny, binz, content); }
876 
878 void
879 MonitorElement::setBinError(int binx, double error)
880 {
881  update();
882  accessRootObject(__PRETTY_FUNCTION__, 1)
883  ->SetBinError(binx, error);
884 }
885 
887 void
888 MonitorElement::setBinError(int binx, int biny, double error)
889 {
890  update();
891  accessRootObject(__PRETTY_FUNCTION__, 2)
892  ->SetBinError(binx, biny, error);
893 }
894 
896 void
897 MonitorElement::setBinError(int binx, int biny, int binz, double error)
898 {
899  update();
900  accessRootObject(__PRETTY_FUNCTION__, 3)
901  ->SetBinError(binx, biny, binz, error);
902 }
903 
905 void
906 MonitorElement::setBinEntries(int bin, double nentries)
907 {
908  update();
909  if (kind() == DQM_KIND_TPROFILE)
910  static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
911  ->SetBinEntries(bin, nentries);
912  else if (kind() == DQM_KIND_TPROFILE2D)
913  static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 1))
914  ->SetBinEntries(bin, nentries);
915  else
916  incompatible(__PRETTY_FUNCTION__);
917 }
918 
920 void
921 MonitorElement::setEntries(double nentries)
922 {
923  update();
924  accessRootObject(__PRETTY_FUNCTION__, 1)
925  ->SetEntries(nentries);
926 }
927 
929 void
930 MonitorElement::setBinLabel(int bin, const std::string &label, int axis /* = 1 */)
931 {
932  update();
933  if ( getAxis(__PRETTY_FUNCTION__, axis)->GetNbins() >= bin )
934  {
935  getAxis(__PRETTY_FUNCTION__, axis)
936  ->SetBinLabel(bin, label.c_str());
937  }
938  else
939  {
940  // edm::LogWarning ("MonitorElement")
941  std::cout << "*** MonitorElement: WARNING:"
942  <<"setBinLabel: attempting to set label of non-existent bin number for ME: "<< getFullname() << " \n";
943  }
944 }
945 
947 void
948 MonitorElement::setAxisRange(double xmin, double xmax, int axis /* = 1 */)
949 {
950  update();
951  getAxis(__PRETTY_FUNCTION__, axis)
952  ->SetRangeUser(xmin, xmax);
953 }
954 
956 void
957 MonitorElement::setAxisTitle(const std::string &title, int axis /* = 1 */)
958 {
959  update();
960  getAxis(__PRETTY_FUNCTION__, axis)
961  ->SetTitle(title.c_str());
962 }
963 
965 void
966 MonitorElement::setAxisTimeDisplay(int value, int axis /* = 1 */)
967 {
968  update();
969  getAxis(__PRETTY_FUNCTION__, axis)
970  ->SetTimeDisplay(value);
971 }
972 
974 void
975 MonitorElement::setAxisTimeFormat(const char *format /* = "" */, int axis /* = 1 */)
976 {
977  update();
978  getAxis(__PRETTY_FUNCTION__, axis)
979  ->SetTimeFormat(format);
980 }
981 
983 void
984 MonitorElement::setAxisTimeOffset(double toffset, const char *option /* ="local" */, int axis /* = 1 */)
985 {
986  update();
987  getAxis(__PRETTY_FUNCTION__, axis)
988  ->SetTimeOffset(toffset, option);
989 }
990 
992 void
993 MonitorElement::setTitle(const std::string &title)
994 {
995  update();
996  accessRootObject(__PRETTY_FUNCTION__, 1)
997  ->SetTitle(title.c_str());
998 }
999 
1000 TAxis *
1001 MonitorElement::getAxis(const char *func, int axis) const
1002 {
1003  TH1 *h = accessRootObject(func, axis-1);
1004  TAxis *a = 0;
1005  if (axis == 1)
1006  a = h->GetXaxis();
1007  else if (axis == 2)
1008  a = h->GetYaxis();
1009  else if (axis == 3)
1010  a = h->GetZaxis();
1011 
1012  if (! a)
1013  raiseDQMError("MonitorElement", "No such axis %d in monitor element"
1014  " '%s' of type '%s'", axis, data_.objname.c_str(),
1015  typeid(*h).name());
1016 
1017  return a;
1018 }
1019 
1020 // ------------ Operations for MEs that are normally never reset ---------
1021 
1024 void
1025 MonitorElement::softReset(void)
1026 {
1027  update();
1028 
1029  // Create the reference object the first time this is called.
1030  // On subsequent calls accumulate the current value to the
1031  // reference, and then reset the current value. This way the
1032  // future contents will have the reference "subtracted".
1033  if (kind() == DQM_KIND_TH1F)
1034  {
1035  TH1F *orig = static_cast<TH1F *>(object_);
1036  TH1F *r = static_cast<TH1F *>(refvalue_);
1037  if (! r)
1038  {
1039  refvalue_ = r = new TH1F((std::string(orig->GetName()) + "_ref").c_str(),
1040  orig->GetTitle(),
1041  orig->GetNbinsX(),
1042  orig->GetXaxis()->GetXmin(),
1043  orig->GetXaxis()->GetXmax());
1044  r->SetDirectory(0);
1045  r->Reset();
1046  }
1047 
1048  r->Add(orig);
1049  orig->Reset();
1050  }
1051  else if (kind() == DQM_KIND_TH1S)
1052  {
1053  TH1S *orig = static_cast<TH1S *>(object_);
1054  TH1S *r = static_cast<TH1S *>(refvalue_);
1055  if (! r)
1056  {
1057  refvalue_ = r = new TH1S((std::string(orig->GetName()) + "_ref").c_str(),
1058  orig->GetTitle(),
1059  orig->GetNbinsX(),
1060  orig->GetXaxis()->GetXmin(),
1061  orig->GetXaxis()->GetXmax());
1062  r->SetDirectory(0);
1063  r->Reset();
1064  }
1065 
1066  r->Add(orig);
1067  orig->Reset();
1068  }
1069  else if (kind() == DQM_KIND_TH1D)
1070  {
1071  TH1D *orig = static_cast<TH1D *>(object_);
1072  TH1D *r = static_cast<TH1D *>(refvalue_);
1073  if (! r)
1074  {
1075  refvalue_ = r = new TH1D((std::string(orig->GetName()) + "_ref").c_str(),
1076  orig->GetTitle(),
1077  orig->GetNbinsX(),
1078  orig->GetXaxis()->GetXmin(),
1079  orig->GetXaxis()->GetXmax());
1080  r->SetDirectory(0);
1081  r->Reset();
1082  }
1083 
1084  r->Add(orig);
1085  orig->Reset();
1086  }
1087  else if (kind() == DQM_KIND_TH2F)
1088  {
1089  TH2F *orig = static_cast<TH2F *>(object_);
1090  TH2F *r = static_cast<TH2F *>(refvalue_);
1091  if (! r)
1092  {
1093  refvalue_ = r = new TH2F((std::string(orig->GetName()) + "_ref").c_str(),
1094  orig->GetTitle(),
1095  orig->GetNbinsX(),
1096  orig->GetXaxis()->GetXmin(),
1097  orig->GetXaxis()->GetXmax(),
1098  orig->GetNbinsY(),
1099  orig->GetYaxis()->GetXmin(),
1100  orig->GetYaxis()->GetXmax());
1101  r->SetDirectory(0);
1102  r->Reset();
1103  }
1104 
1105  r->Add(orig);
1106  orig->Reset();
1107  }
1108  else if (kind() == DQM_KIND_TH2S)
1109  {
1110  TH2S *orig = static_cast<TH2S *>(object_);
1111  TH2S *r = static_cast<TH2S *>(refvalue_);
1112  if (! r)
1113  {
1114  refvalue_ = r = new TH2S((std::string(orig->GetName()) + "_ref").c_str(),
1115  orig->GetTitle(),
1116  orig->GetNbinsX(),
1117  orig->GetXaxis()->GetXmin(),
1118  orig->GetXaxis()->GetXmax(),
1119  orig->GetNbinsY(),
1120  orig->GetYaxis()->GetXmin(),
1121  orig->GetYaxis()->GetXmax());
1122  r->SetDirectory(0);
1123  r->Reset();
1124  }
1125 
1126  r->Add(orig);
1127  orig->Reset();
1128  }
1129  else if (kind() == DQM_KIND_TH2D)
1130  {
1131  TH2D *orig = static_cast<TH2D *>(object_);
1132  TH2D *r = static_cast<TH2D *>(refvalue_);
1133  if (! r)
1134  {
1135  refvalue_ = r = new TH2D((std::string(orig->GetName()) + "_ref").c_str(),
1136  orig->GetTitle(),
1137  orig->GetNbinsX(),
1138  orig->GetXaxis()->GetXmin(),
1139  orig->GetXaxis()->GetXmax(),
1140  orig->GetNbinsY(),
1141  orig->GetYaxis()->GetXmin(),
1142  orig->GetYaxis()->GetXmax());
1143  r->SetDirectory(0);
1144  r->Reset();
1145  }
1146 
1147  r->Add(orig);
1148  orig->Reset();
1149  }
1150  else if (kind() == DQM_KIND_TH3F)
1151  {
1152  TH3F *orig = static_cast<TH3F *>(object_);
1153  TH3F *r = static_cast<TH3F *>(refvalue_);
1154  if (! r)
1155  {
1156  refvalue_ = r = new TH3F((std::string(orig->GetName()) + "_ref").c_str(),
1157  orig->GetTitle(),
1158  orig->GetNbinsX(),
1159  orig->GetXaxis()->GetXmin(),
1160  orig->GetXaxis()->GetXmax(),
1161  orig->GetNbinsY(),
1162  orig->GetYaxis()->GetXmin(),
1163  orig->GetYaxis()->GetXmax(),
1164  orig->GetNbinsZ(),
1165  orig->GetZaxis()->GetXmin(),
1166  orig->GetZaxis()->GetXmax());
1167  r->SetDirectory(0);
1168  r->Reset();
1169  }
1170 
1171  r->Add(orig);
1172  orig->Reset();
1173  }
1174  else if (kind() == DQM_KIND_TPROFILE)
1175  {
1176  TProfile *orig = static_cast<TProfile *>(object_);
1177  TProfile *r = static_cast<TProfile *>(refvalue_);
1178  if (! r)
1179  {
1180  refvalue_ = r = new TProfile((std::string(orig->GetName()) + "_ref").c_str(),
1181  orig->GetTitle(),
1182  orig->GetNbinsX(),
1183  orig->GetXaxis()->GetXmin(),
1184  orig->GetXaxis()->GetXmax(),
1185  orig->GetYaxis()->GetXmin(),
1186  orig->GetYaxis()->GetXmax(),
1187  orig->GetErrorOption());
1188  r->SetDirectory(0);
1189  r->Reset();
1190  }
1191 
1192  addProfiles(r, orig, r, 1, 1);
1193  orig->Reset();
1194  }
1195  else if (kind() == DQM_KIND_TPROFILE2D)
1196  {
1197  TProfile2D *orig = static_cast<TProfile2D *>(object_);
1198  TProfile2D *r = static_cast<TProfile2D *>(refvalue_);
1199  if (! r)
1200  {
1201  refvalue_ = r = new TProfile2D((std::string(orig->GetName()) + "_ref").c_str(),
1202  orig->GetTitle(),
1203  orig->GetNbinsX(),
1204  orig->GetXaxis()->GetXmin(),
1205  orig->GetXaxis()->GetXmax(),
1206  orig->GetNbinsY(),
1207  orig->GetYaxis()->GetXmin(),
1208  orig->GetYaxis()->GetXmax(),
1209  orig->GetZaxis()->GetXmin(),
1210  orig->GetZaxis()->GetXmax(),
1211  orig->GetErrorOption());
1212  r->SetDirectory(0);
1213  r->Reset();
1214  }
1215 
1216  addProfiles(r, orig, r, 1, 1);
1217  orig->Reset();
1218  }
1219  else
1220  incompatible(__PRETTY_FUNCTION__);
1221 }
1222 
1224 void
1225 MonitorElement::disableSoftReset(void)
1226 {
1227  if (refvalue_)
1228  {
1229  if (kind() == DQM_KIND_TH1F
1230  || kind() == DQM_KIND_TH1S
1231  || kind() == DQM_KIND_TH1D
1232  || kind() == DQM_KIND_TH2F
1233  || kind() == DQM_KIND_TH2S
1234  || kind() == DQM_KIND_TH2D
1235  || kind() == DQM_KIND_TH3F)
1236  {
1237  TH1 *orig = static_cast<TH1 *>(object_);
1238  orig->Add(refvalue_);
1239  }
1240  else if (kind() == DQM_KIND_TPROFILE)
1241  {
1242  TProfile *orig = static_cast<TProfile *>(object_);
1243  TProfile *r = static_cast<TProfile *>(refvalue_);
1244  addProfiles(orig, r, orig, 1, 1);
1245  }
1246  else if (kind() == DQM_KIND_TPROFILE2D)
1247  {
1248  TProfile2D *orig = static_cast<TProfile2D *>(object_);
1249  TProfile2D *r = static_cast<TProfile2D *>(refvalue_);
1250  addProfiles(orig, r, orig, 1, 1);
1251  }
1252  else
1253  incompatible(__PRETTY_FUNCTION__);
1254 
1255  delete refvalue_;
1256  refvalue_ = 0;
1257  }
1258 }
1259 
1260 // implementation: Giuseppe.Della-Ricca@ts.infn.it
1261 // Can be called with sum = h1 or sum = h2
1262 void
1263 MonitorElement::addProfiles(TProfile *h1, TProfile *h2, TProfile *sum, float c1, float c2)
1264 {
1265  assert(h1);
1266  assert(h2);
1267  assert(sum);
1268 
1269  static const Int_t NUM_STAT = 6;
1270  Double_t stats1[NUM_STAT];
1271  Double_t stats2[NUM_STAT];
1272  Double_t stats3[NUM_STAT];
1273 
1274  bool isRebinOn = sum->TestBit(TH1::kCanRebin);
1275  sum->ResetBit(TH1::kCanRebin);
1276 
1277  for (Int_t i = 0; i < NUM_STAT; ++i)
1278  stats1[i] = stats2[i] = stats3[i] = 0;
1279 
1280  h1->GetStats(stats1);
1281  h2->GetStats(stats2);
1282 
1283  for (Int_t i = 0; i < NUM_STAT; ++i)
1284  stats3[i] = c1*stats1[i] + c2*stats2[i];
1285 
1286  stats3[1] = c1*TMath::Abs(c1)*stats1[1]
1287  + c2*TMath::Abs(c2)*stats2[1];
1288 
1289  Double_t entries = c1*h1->GetEntries() + c2* h2->GetEntries();
1290  TArrayD* h1sumw2 = h1->GetSumw2();
1291  TArrayD* h2sumw2 = h2->GetSumw2();
1292  for (Int_t bin = 0, nbin = sum->GetNbinsX()+1; bin <= nbin; ++bin)
1293  {
1294  Double_t entries = c1*h1->GetBinEntries(bin)
1295  + c2*h2->GetBinEntries(bin);
1296  Double_t content = c1*h1->GetBinEntries(bin)*h1->GetBinContent(bin)
1297  + c2*h2->GetBinEntries(bin)*h2->GetBinContent(bin);
1298  Double_t error = TMath::Sqrt(c1*TMath::Abs(c1)*h1sumw2->fArray[bin]
1299  + c2*TMath::Abs(c2)*h2sumw2->fArray[bin]);
1300  sum->SetBinContent(bin, content);
1301  sum->SetBinError(bin, error);
1302  sum->SetBinEntries(bin, entries);
1303  }
1304 
1305  sum->SetEntries(entries);
1306  sum->PutStats(stats3);
1307  if (isRebinOn) sum->SetBit(TH1::kCanRebin);
1308 }
1309 
1310 // implementation: Giuseppe.Della-Ricca@ts.infn.it
1311 // Can be called with sum = h1 or sum = h2
1312 void
1313 MonitorElement::addProfiles(TProfile2D *h1, TProfile2D *h2, TProfile2D *sum, float c1, float c2)
1314 {
1315  assert(h1);
1316  assert(h2);
1317  assert(sum);
1318 
1319  static const Int_t NUM_STAT = 9;
1320  Double_t stats1[NUM_STAT];
1321  Double_t stats2[NUM_STAT];
1322  Double_t stats3[NUM_STAT];
1323 
1324  bool isRebinOn = sum->TestBit(TH1::kCanRebin);
1325  sum->ResetBit(TH1::kCanRebin);
1326 
1327  for (Int_t i = 0; i < NUM_STAT; ++i)
1328  stats1[i] = stats2[i] = stats3[i] = 0;
1329 
1330  h1->GetStats(stats1);
1331  h2->GetStats(stats2);
1332 
1333  for (Int_t i = 0; i < NUM_STAT; i++)
1334  stats3[i] = c1*stats1[i] + c2*stats2[i];
1335 
1336  stats3[1] = c1*TMath::Abs(c1)*stats1[1]
1337  + c2*TMath::Abs(c2)*stats2[1];
1338 
1339  Double_t entries = c1*h1->GetEntries() + c2*h2->GetEntries();
1340  TArrayD *h1sumw2 = h1->GetSumw2();
1341  TArrayD *h2sumw2 = h2->GetSumw2();
1342  for (Int_t xbin = 0, nxbin = sum->GetNbinsX()+1; xbin <= nxbin; ++xbin)
1343  for (Int_t ybin = 0, nybin = sum->GetNbinsY()+1; ybin <= nybin; ++ybin)
1344  {
1345  Int_t bin = sum->GetBin(xbin, ybin);
1346  Double_t entries = c1*h1->GetBinEntries(bin)
1347  + c2*h2->GetBinEntries(bin);
1348  Double_t content = c1*h1->GetBinEntries(bin)*h1->GetBinContent(bin)
1349  + c2*h2->GetBinEntries(bin)*h2->GetBinContent(bin);
1350  Double_t error = TMath::Sqrt(c1*TMath::Abs(c1)*h1sumw2->fArray[bin]
1351  + c2*TMath::Abs(c2)*h2sumw2->fArray[bin]);
1352 
1353  sum->SetBinContent(bin, content);
1354  sum->SetBinError(bin, error);
1355  sum->SetBinEntries(bin, entries);
1356  }
1357  sum->SetEntries(entries);
1358  sum->PutStats(stats3);
1359  if (isRebinOn) sum->SetBit(TH1::kCanRebin);
1360 }
1361 
1362 void
1363 MonitorElement::copyFunctions(TH1 *from, TH1 *to)
1364 {
1365  // will copy functions only if local-copy and original-object are equal
1366  // (ie. no soft-resetting or accumulating is enabled)
1367  if (isSoftResetEnabled() || isAccumulateEnabled())
1368  return;
1369 
1370  update();
1371  TList *fromf = from->GetListOfFunctions();
1372  TList *tof = to->GetListOfFunctions();
1373  for (int i = 0, nfuncs = fromf ? fromf->GetSize() : 0; i < nfuncs; ++i)
1374  {
1375  TObject *obj = fromf->At(i);
1376  // not interested in statistics
1377  if (!strcmp(obj->IsA()->GetName(), "TPaveStats"))
1378  continue;
1379 
1380  if (TF1 *fn = dynamic_cast<TF1 *>(obj))
1381  tof->Add(new TF1(*fn));
1382  //else if (dynamic_cast<TPaveStats *>(obj))
1383  // ; // FIXME? tof->Add(new TPaveStats(*stats));
1384  else
1385  raiseDQMError("MonitorElement", "Cannot extract function '%s' of type"
1386  " '%s' from monitor element '%s' for a copy",
1387  obj->GetName(), obj->IsA()->GetName(), data_.objname.c_str());
1388  }
1389 }
1390 
1391 void
1392 MonitorElement::copyFrom(TH1 *from)
1393 {
1394  TH1 *orig = accessRootObject(__PRETTY_FUNCTION__, 1);
1395  if (orig->GetTitle() != from->GetTitle())
1396  orig->SetTitle(from->GetTitle());
1397 
1398  if (!isAccumulateEnabled())
1399  orig->Reset();
1400 
1401  if (isSoftResetEnabled())
1402  {
1403  if (kind() == DQM_KIND_TH1F
1404  || kind() == DQM_KIND_TH1S
1405  || kind() == DQM_KIND_TH1D
1406  || kind() == DQM_KIND_TH2F
1407  || kind() == DQM_KIND_TH2S
1408  || kind() == DQM_KIND_TH2D
1409  || kind() == DQM_KIND_TH3F)
1410  // subtract "reference"
1411  orig->Add(from, refvalue_, 1, -1);
1412  else if (kind() == DQM_KIND_TPROFILE)
1413  // subtract "reference"
1414  addProfiles(static_cast<TProfile *>(from),
1415  static_cast<TProfile *>(refvalue_),
1416  static_cast<TProfile *>(orig),
1417  1, -1);
1418  else if (kind() == DQM_KIND_TPROFILE2D)
1419  // subtract "reference"
1420  addProfiles(static_cast<TProfile2D *>(from),
1421  static_cast<TProfile2D *>(refvalue_),
1422  static_cast<TProfile2D *>(orig),
1423  1, -1);
1424  else
1425  incompatible(__PRETTY_FUNCTION__);
1426  }
1427  else
1428  orig->Add(from);
1429 
1430  copyFunctions(from, orig);
1431 }
1432 
1433 // --- Operations on MEs that are normally reset at end of monitoring cycle ---
1434 void
1435 MonitorElement::getQReport(bool create, const std::string &qtname, QReport *&qr, DQMNet::QValue *&qv)
1436 {
1437  assert(qreports_.size() == data_.qreports.size());
1438 
1439  qr = 0;
1440  qv = 0;
1441 
1442  size_t pos = 0, end = qreports_.size();
1443  while (pos < end && data_.qreports[pos].qtname != qtname)
1444  ++pos;
1445 
1446  if (pos == end && ! create)
1447  return;
1448  else if (pos == end)
1449  {
1450  data_.qreports.push_back(DQMNet::QValue());
1451  qreports_.push_back(QReport(0, 0));
1452 
1453  DQMNet::QValue &q = data_.qreports.back();
1454  q.code = dqm::qstatus::DID_NOT_RUN;
1455  q.qtresult = 0;
1456  q.qtname = qtname;
1457  q.message = "NO_MESSAGE_ASSIGNED";
1458  q.algorithm = "UNKNOWN_ALGORITHM";
1459  }
1460 
1461  qr = &qreports_[pos];
1462  qv = &data_.qreports[pos];
1463 }
1464 
1466 void
1467 MonitorElement::addQReport(const DQMNet::QValue &desc, QCriterion *qc)
1468 {
1469  QReport *qr;
1470  DQMNet::QValue *qv;
1471  getQReport(true, desc.qtname, qr, qv);
1472  qr->qcriterion_ = qc;
1473  *qv = desc;
1474  update();
1475 }
1476 
1477 void
1478 MonitorElement::addQReport(QCriterion *qc)
1479 {
1480  QReport *qr;
1481  DQMNet::QValue *qv;
1482  getQReport(true, qc->getName(), qr, qv);
1483  qv->code = dqm::qstatus::DID_NOT_RUN;
1484  qv->message = "NO_MESSAGE_ASSIGNED";
1485  qr->qcriterion_ = qc;
1486  update();
1487 }
1488 
1490 void
1491 MonitorElement::updateQReportStats(void)
1492 {
1493  data_.flags &= ~DQMNet::DQM_PROP_REPORT_ALARM;
1494  for (size_t i = 0, e = data_.qreports.size(); i < e; ++i)
1495  switch (data_.qreports[i].code)
1496  {
1497  case dqm::qstatus::STATUS_OK:
1498  break;
1499  case dqm::qstatus::WARNING:
1500  data_.flags |= DQMNet::DQM_PROP_REPORT_WARN;
1501  break;
1502  case dqm::qstatus::ERROR:
1503  data_.flags |= DQMNet::DQM_PROP_REPORT_ERROR;
1504  break;
1505  default:
1506  data_.flags |= DQMNet::DQM_PROP_REPORT_OTHER;
1507  break;
1508  }
1509 }
1510 
1511 // -------------------------------------------------------------------
1512 TObject *
1513 MonitorElement::getRootObject(void) const
1514 {
1515  const_cast<MonitorElement *>(this)->update();
1516  return object_;
1517 }
1518 
1519 TH1 *
1520 MonitorElement::getTH1(void) const
1521 {
1522  const_cast<MonitorElement *>(this)->update();
1523  return accessRootObject(__PRETTY_FUNCTION__, 0);
1524 }
1525 
1526 TH1F *
1527 MonitorElement::getTH1F(void) const
1528 {
1529  assert(kind() == DQM_KIND_TH1F);
1530  const_cast<MonitorElement *>(this)->update();
1531  return static_cast<TH1F *>(accessRootObject(__PRETTY_FUNCTION__, 1));
1532 }
1533 
1534 TH1S *
1535 MonitorElement::getTH1S(void) const
1536 {
1537  assert(kind() == DQM_KIND_TH1S);
1538  const_cast<MonitorElement *>(this)->update();
1539  return static_cast<TH1S *>(accessRootObject(__PRETTY_FUNCTION__, 1));
1540 }
1541 
1542 TH1D *
1543 MonitorElement::getTH1D(void) const
1544 {
1545  assert(kind() == DQM_KIND_TH1D);
1546  const_cast<MonitorElement *>(this)->update();
1547  return static_cast<TH1D *>(accessRootObject(__PRETTY_FUNCTION__, 1));
1548 }
1549 
1550 TH2F *
1551 MonitorElement::getTH2F(void) const
1552 {
1553  assert(kind() == DQM_KIND_TH2F);
1554  const_cast<MonitorElement *>(this)->update();
1555  return static_cast<TH2F *>(accessRootObject(__PRETTY_FUNCTION__, 2));
1556 }
1557 
1558 TH2S *
1559 MonitorElement::getTH2S(void) const
1560 {
1561  assert(kind() == DQM_KIND_TH2S);
1562  const_cast<MonitorElement *>(this)->update();
1563  return static_cast<TH2S *>(accessRootObject(__PRETTY_FUNCTION__, 2));
1564 }
1565 
1566 TH2D *
1567 MonitorElement::getTH2D(void) const
1568 {
1569  assert(kind() == DQM_KIND_TH2D);
1570  const_cast<MonitorElement *>(this)->update();
1571  return static_cast<TH2D *>(accessRootObject(__PRETTY_FUNCTION__, 2));
1572 }
1573 
1574 TH3F *
1575 MonitorElement::getTH3F(void) const
1576 {
1577  assert(kind() == DQM_KIND_TH3F);
1578  const_cast<MonitorElement *>(this)->update();
1579  return static_cast<TH3F *>(accessRootObject(__PRETTY_FUNCTION__, 3));
1580 }
1581 
1582 TProfile *
1583 MonitorElement::getTProfile(void) const
1584 {
1585  assert(kind() == DQM_KIND_TPROFILE);
1586  const_cast<MonitorElement *>(this)->update();
1587  return static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1));
1588 }
1589 
1590 TProfile2D *
1591 MonitorElement::getTProfile2D(void) const
1592 {
1593  assert(kind() == DQM_KIND_TPROFILE2D);
1594  const_cast<MonitorElement *>(this)->update();
1595  return static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 2));
1596 }
1597 
1598 // -------------------------------------------------------------------
1599 TObject *
1600 MonitorElement::getRefRootObject(void) const
1601 {
1602  const_cast<MonitorElement *>(this)->update();
1603  return reference_;
1604 }
1605 
1606 TH1 *
1607 MonitorElement::getRefTH1(void) const
1608 {
1609  const_cast<MonitorElement *>(this)->update();
1610  return checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 0);
1611 }
1612 
1613 TH1F *
1614 MonitorElement::getRefTH1F(void) const
1615 {
1616  assert(kind() == DQM_KIND_TH1F);
1617  const_cast<MonitorElement *>(this)->update();
1618  return static_cast<TH1F *>
1619  (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 1));
1620 }
1621 
1622 TH1S *
1623 MonitorElement::getRefTH1S(void) const
1624 {
1625  assert(kind() == DQM_KIND_TH1S);
1626  const_cast<MonitorElement *>(this)->update();
1627  return static_cast<TH1S *>
1628  (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 1));
1629 }
1630 
1631 TH1D *
1632 MonitorElement::getRefTH1D(void) const
1633 {
1634  assert(kind() == DQM_KIND_TH1D);
1635  const_cast<MonitorElement *>(this)->update();
1636  return static_cast<TH1D *>
1637  (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 1));
1638 }
1639 
1640 TH2F *
1641 MonitorElement::getRefTH2F(void) const
1642 {
1643  assert(kind() == DQM_KIND_TH2F);
1644  const_cast<MonitorElement *>(this)->update();
1645  return static_cast<TH2F *>
1646  (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 2));
1647 }
1648 
1649 TH2S *
1650 MonitorElement::getRefTH2S(void) const
1651 {
1652  assert(kind() == DQM_KIND_TH2S);
1653  const_cast<MonitorElement *>(this)->update();
1654  return static_cast<TH2S *>
1655  (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 2));
1656 }
1657 
1658 TH2D *
1659 MonitorElement::getRefTH2D(void) const
1660 {
1661  assert(kind() == DQM_KIND_TH2D);
1662  const_cast<MonitorElement *>(this)->update();
1663  return static_cast<TH2D *>
1664  (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 2));
1665 }
1666 
1667 TH3F *
1668 MonitorElement::getRefTH3F(void) const
1669 {
1670  assert(kind() == DQM_KIND_TH3F);
1671  const_cast<MonitorElement *>(this)->update();
1672  return static_cast<TH3F *>
1673  (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 3));
1674 }
1675 
1676 TProfile *
1677 MonitorElement::getRefTProfile(void) const
1678 {
1679  assert(kind() == DQM_KIND_TPROFILE);
1680  const_cast<MonitorElement *>(this)->update();
1681  return static_cast<TProfile *>
1682  (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 1));
1683 }
1684 
1685 TProfile2D *
1686 MonitorElement::getRefTProfile2D(void) const
1687 {
1688  assert(kind() == DQM_KIND_TPROFILE2D);
1689  const_cast<MonitorElement *>(this)->update();
1690  return static_cast<TProfile2D *>
1691  (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 2));
1692 }
MonitorElement::getRefTH1F
TH1F * getRefTH1F(void) const
Definition: MonitorElement.cc:1614
MonitorElement::getTH2S
TH2S * getTH2S(void) const
Definition: MonitorElement.cc:1559
MonitorElement::getTH1S
TH1S * getTH1S(void) const
Definition: MonitorElement.cc:1535
QReport::qcriterion_
QCriterion * qcriterion_
Definition: QReport.h:51
DQMNet::DQM_PROP_REPORT_WARN
static const uint32_t DQM_PROP_REPORT_WARN
Definition: DQMNet.h:47
i
int i
Definition: DBlmapReader.cc:9
MonitorElement::incompatible
void incompatible(const char *func) const
Definition: MonitorElement.cc:689
DQMNet::CoreObject::moduleId
uint32_t moduleId
Definition: DQMNet.h:104
MonitorElement::getQReport
const QReport * getQReport(const std::string &qtname) const
get QReport corresponding to &lt;qtname&gt; (null pointer if QReport does not exist)
Definition: MonitorElement.cc:586
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:854
DQMNet::CoreObject::qreports
QReports qreports
Definition: DQMNet.h:107
MonitorElement::accessRootObject
TH1 * accessRootObject(const char *func, int reqdim) const
Definition: MonitorElement.cc:696
QCriterion::getName
std::string getName(void) const
get name of quality test
Definition: QTest.h:57
MonitorElement::copyFrom
void copyFrom(TH1 *from)
Definition: MonitorElement.cc:1392
MonitorElement::initialise
MonitorElement * initialise(Kind kind)
Definition: MonitorElement.cc:35
DQMNet::QValue::algorithm
std::string algorithm
Definition: DQMNet.h:93
w
const double w
Definition: UKUtility.cc:23
MonitorElement::operator=
MonitorElement & operator=(const MonitorElement &)
Definition: MonitorElement.cc:202
MonitorElement::getTProfile2D
TProfile2D * getTProfile2D(void) const
Definition: MonitorElement.cc:1591
MonitorElement::getRefTH1
TH1 * getRefTH1(void) const
Definition: MonitorElement.cc:1607
DQMNet::CoreObject::version
uint64_t version
Definition: DQMNet.h:100
MonitorElement::setAxisRange
void setAxisRange(double xmin, double xmax, int axis=1)
set x-, y- or z-axis range (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:948
MonitorElement::setAxisTimeFormat
void setAxisTimeFormat(const char *format="", int axis=1)
set the format of the time values that are displayed on an axis
Definition: MonitorElement.cc:975
MonitorElement::updateQReportStats
void updateQReportStats(void)
Refresh QReport stats, usually after MEs were read in from a file.
Definition: MonitorElement.cc:1491
MonitorElement::qreports_
std::vector< QReport > qreports_
Definition: MonitorElement.h:71
MonitorElement::tagLabelString
std::string tagLabelString(void) const
return label string for the monitor element tag (eg. &lt;name&gt;t=12345&lt;/name&gt;)
Definition: MonitorElement.cc:544
DQMNet
Definition: DQMNet.h:22
MonitorElement::qualityTagString
std::string qualityTagString(const DQMNet::QValue &qv) const
Definition: MonitorElement.cc:571
MonitorElement::getAxisTitle
std::string getAxisTitle(int axis=1) const
get x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:840
MonitorElement::getRefTProfile2D
TProfile2D * getRefTProfile2D(void) const
Definition: MonitorElement.cc:1686
zw
auto zw(V v) -> Vec2< typenamestd::remove_reference< decltype(v[0])>::type >
Definition: ExtVec.h:36
DQMNet::CoreObject::flags
uint32_t flags
Definition: DQMNet.h:98
MonitorElement::disableSoftReset
void disableSoftReset(void)
reverts action of softReset
Definition: MonitorElement.cc:1225
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
dqm::qstatus::WARNING
static const int WARNING
Definition: DQMDefinitions.h:18
MonitorElement::getTH3F
TH3F * getTH3F(void) const
Definition: MonitorElement.cc:1575
MonitorElement::getTH1D
TH1D * getTH1D(void) const
Definition: MonitorElement.cc:1543
MonitorElement::runQTests
void runQTests(void)
run all quality tests
Definition: MonitorElement.cc:657
MonitorElement::softReset
void softReset(void)
Definition: MonitorElement.cc:1025
MonitorElement::getTH2D
TH2D * getTH2D(void) const
Definition: MonitorElement.cc:1567
DQMNet::CoreObject::tag
uint32_t tag
Definition: DQMNet.h:99
DQMNet::CoreObject::dirname
const std::string * dirname
Definition: DQMNet.h:105
MonitorElement::update
void update(void)
Mark the object updated.
Definition: MonitorElement.h:135
MonitorElement::getEntries
double getEntries(void) const
get # of entries
Definition: MonitorElement.cc:789
MonitorElement::getMean
double getMean(int axis=1) const
get mean value of histogram along x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:710
DQMNet::CoreObject::run
uint32_t run
Definition: DQMNet.h:101
detailsBasic3DVector::z
float float float z
Definition: extBasic3DVector.h:15
SiStripMonitorClusterAlca_cfi.xmin
tuple xmin
Definition: SiStripMonitorClusterAlca_cfi.py:34
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
MonitorElement::getMeanError
double getMeanError(int axis=1) const
Definition: MonitorElement.cc:717
MonitorElement::packScalarData
void packScalarData(std::string &into, const char *prefix) const
convert scalar data into a string.
Definition: MonitorElement.cc:479
MonitorElement::getNbinsY
int getNbinsY(void) const
get # of bins in Y-axis
Definition: MonitorElement.cc:741
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:164
MonitorElement::getQErrors
std::vector< QReport * > getQErrors(void) const
get errors from last set of quality tests
Definition: MonitorElement.cc:624
MonitorElement::wasUpdated
bool wasUpdated(void) const
true if ME was updated in last monitoring cycle
Definition: MonitorElement.h:131
MonitorElement::getRefTH3F
TH3F * getRefTH3F(void) const
Definition: MonitorElement.cc:1668
counter.c2
tuple c2
Definition: counter.py:145
MonitorElement::ShiftFillLast
void ShiftFillLast(double y, double ye=0., int32_t xscale=1)
Definition: MonitorElement.cc:349
checkRootObject
static TH1 * checkRootObject(const std::string &name, TObject *tobj, const char *func, int reqdim)
Definition: MonitorElement.cc:16
ws_sso_content_reader.content
tuple content
Definition: ws_sso_content_reader.py:117
MonitorElement::packQualityData
void packQualityData(std::string &into) const
serialise quality report information into a string.
Definition: MonitorElement.cc:504
DQMNet::DQM_PROP_HAS_REFERENCE
static const uint32_t DQM_PROP_HAS_REFERENCE
Definition: DQMNet.h:53
cmsHarvester.path
tuple path
else: Piece not in the list, fine.
Definition: cmsHarvester.py:1591
MonitorElement::getRefTH2F
TH2F * getRefTH2F(void) const
Definition: MonitorElement.cc:1641
MonitorElement::setAxisTimeDisplay
void setAxisTimeDisplay(int value, int axis=1)
set x-, y-, or z-axis to display time values
Definition: MonitorElement.cc:966
MonitorElement::doFill
void doFill(int64_t x)
&quot;Fill&quot; ME method for int64_t
Definition: MonitorElement.cc:295
MonitorElement::getNbinsZ
int getNbinsZ(void) const
get # of bins in Z-axis
Definition: MonitorElement.cc:747
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
DQMNet::CoreObject::lumi
uint32_t lumi
Definition: DQMNet.h:102
query.result
tuple result
Definition: query.py:137
MonitorElement::addProfiles
void addProfiles(TProfile *h1, TProfile *h2, TProfile *sum, float c1, float c2)
Definition: MonitorElement.cc:1263
MonitorElement::getTitle
std::string getTitle(void) const
get MonitorElement title
Definition: MonitorElement.cc:846
Abs
T Abs(T a)
Definition: MathUtil.h:49
dqm::qstatus::DID_NOT_RUN
static const int DID_NOT_RUN
Definition: DQMDefinitions.h:16
MonitorElement::getTH1
TH1 * getTH1(void) const
Definition: MonitorElement.cc:1520
DQMNet::DQM_PROP_REPORT_ERROR
static const uint32_t DQM_PROP_REPORT_ERROR
Definition: DQMNet.h:46
h
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
Definition: Activities.doc:4
MonitorElement::setBinError
void setBinError(int binx, double error)
set uncertainty on content of bin (1-D)
Definition: MonitorElement.cc:879
end
#define end
Definition: vmac.h:37
DQMNet::DQM_PROP_REPORT_OTHER
static const uint32_t DQM_PROP_REPORT_OTHER
Definition: DQMNet.h:48
MonitorElement::kind
Kind kind(void) const
Get the type of the monitor element.
Definition: MonitorElement.h:103
MonitorElement::setEntries
void setEntries(double nentries)
set # of entries
Definition: MonitorElement.cc:921
MonitorElement::getFullname
const std::string getFullname(void) const
get full name of ME including Pathname
Definition: MonitorElement.h:119
MonitorElement::effLabelString
std::string effLabelString(void) const
return label string for the monitor element tag (eg. &lt;name&gt;t=12345&lt;/name&gt;)
Definition: MonitorElement.cc:559
MonitorElement::isAccumulateEnabled
bool isAccumulateEnabled(void) const
whether ME contents should be accumulated over multiple monitoring periods; default: false ...
Definition: MonitorElement.h:286
MonitorElement::streamId
const uint32_t streamId(void) const
Definition: MonitorElement.h:385
DQMNet::CoreObject::objname
std::string objname
Definition: DQMNet.h:106
MonitorElement::valueString
std::string valueString(void) const
Definition: MonitorElement.cc:512
MonitorElement::getRefTH2D
TH2D * getRefTH2D(void) const
Definition: MonitorElement.cc:1659
DQMNet::QValue::qtname
std::string qtname
Definition: DQMNet.h:92
MonitorElement::data_
DQMNet::CoreObject data_
Definition: MonitorElement.h:66
MonitorElement::getAxis
TAxis * getAxis(const char *func, int axis) const
Definition: MonitorElement.cc:1001
MonitorElement::setTitle
void setTitle(const std::string &title)
set (ie. change) histogram/profile title
Definition: MonitorElement.cc:993
MonitorElement::getRMSError
double getRMSError(int axis=1) const
get RMS uncertainty of histogram along x, y or z axis(axis=1,2,3 respectively)
Definition: MonitorElement.cc:729
MonitorElement::getRootObject
TObject * getRootObject(void) const
Definition: MonitorElement.cc:1513
N
#define N
Definition: blowfish.cc:9
MonitorElement::getQReports
std::vector< QReport * > getQReports(void) const
get map of QReports
Definition: MonitorElement.cc:595
MonitorElement::setAxisTimeOffset
void setAxisTimeOffset(double toffset, const char *option="local", int axis=1)
set the time offset, if option = &quot;gmt&quot; then the offset is treated as a GMT time
Definition: MonitorElement.cc:984
MonitorElement::getQOthers
std::vector< QReport * > getQOthers(void) const
Definition: MonitorElement.cc:639
MonitorElement::getBinError
double getBinError(int binx) const
get uncertainty on content of bin (1-D) - See TH1::GetBinError for details
Definition: MonitorElement.cc:771
QCriterion::runTest
virtual float runTest(const MonitorElement *me)
Definition: QTest.cc:27
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1527
MonitorElement::tagString
std::string tagString(void) const
Definition: MonitorElement.cc:531
b
double b
Definition: hdecay.h:120
DQMNet::packQualityData
static void packQualityData(std::string &into, const QReports &qr)
Definition: DQMNet.cc:177
SiStripMonitorClusterAlca_cfi.xmax
tuple xmax
Definition: SiStripMonitorClusterAlca_cfi.py:35
MonitorElement::getRefTProfile
TProfile * getRefTProfile(void) const
Definition: MonitorElement.cc:1677
MonitorElement::copyFunctions
void copyFunctions(TH1 *from, TH1 *to)
Definition: MonitorElement.cc:1363
MonitorElement::moduleId
const uint32_t moduleId(void) const
Definition: MonitorElement.h:386
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:753
MonitorElement::getRMS
double getRMS(int axis=1) const
get RMS of histogram along x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:723
MonitorElement::getYmax
double getYmax(void) const
get max Y value (for profiles)
Definition: MonitorElement.cc:826
a
double a
Definition: hdecay.h:121
QReport::qvalue_
DQMNet::QValue * qvalue_
Definition: QReport.h:50
MonitorElement::getTProfile
TProfile * getTProfile(void) const
Definition: MonitorElement.cc:1583
MonitorElement::getRefTH2S
TH2S * getRefTH2S(void) const
Definition: MonitorElement.cc:1650
MonitorElement::isSoftResetEnabled
bool isSoftResetEnabled(void) const
whether soft-reset is enabled; default is false
Definition: MonitorElement.h:282
MonitorElement::getBinEntries
double getBinEntries(int bin) const
get # of bin entries (for profiles)
Definition: MonitorElement.cc:795
MonitorElement::getQWarnings
std::vector< QReport * > getQWarnings(void) const
get warnings from last set of quality tests
Definition: MonitorElement.cc:609
MonitorElement::getNbinsX
int getNbinsX(void) const
get # of bins in X-axis
Definition: MonitorElement.cc:735
DQMNet::QValue::message
std::string message
Definition: DQMNet.h:91
dqm::qstatus::STATUS_OK
static const int STATUS_OK
Definition: DQMDefinitions.h:17
MonitorElement::run
const uint32_t run(void) const
Definition: MonitorElement.h:383
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
DQMNet::DQM_PROP_NEW
static const uint32_t DQM_PROP_NEW
Definition: DQMNet.h:58
x
Definition: DDAxes.h:10
MonitorElement::getTH2F
TH2F * getTH2F(void) const
Definition: MonitorElement.cc:1551
if
if(conf.exists("allCellsPositionCalc"))
Definition: Basic2DGenericPFlowClusterizer.cc:58
DQMNet::CoreObject::streamId
uint32_t streamId
Definition: DQMNet.h:103
DQMNet::QValue::qtresult
float qtresult
Definition: DQMNet.h:90
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:957
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:463
MonitorElement::getRefRootObject
TObject * getRefRootObject(void) const
Definition: MonitorElement.cc:1600
MonitorElement::setBinEntries
void setBinEntries(int bin, double nentries)
set # of bin entries (to be used for profiles)
Definition: MonitorElement.cc:906
MonitorElement::getRefTH1S
TH1S * getRefTH1S(void) const
Definition: MonitorElement.cc:1623
SurfaceDeformationFactory::create
SurfaceDeformation * create(int type, const std::vector< double > &params)
Definition: SurfaceDeformationFactory.cc:27
MonitorElement::addQReport
void addQReport(const DQMNet::QValue &desc, QCriterion *qc)
Add quality report, from DQMStore.
Definition: MonitorElement.cc:1467
dqm::qstatus::ERROR
static const int ERROR
Definition: DQMDefinitions.h:19
MonitorElement::getYmin
double getYmin(void) const
get min Y value (for profiles)
Definition: MonitorElement.cc:812
MonitorElement::getRefTH1D
TH1D * getRefTH1D(void) const
Definition: MonitorElement.cc:1632
MonitorElement::CheckBinLabels
static bool CheckBinLabels(const TAxis *a1, const TAxis *a2)
Check the consistency of the axis labels.
Definition: MonitorElement.cc:235
raiseDQMError
void raiseDQMError(const char *context, const char *fmt,...)
Definition: DQMError.cc:11
Generated for CMSSW Reference Manual by doxygen 1.8.5