CMSSW_4_4_3_patch1/src/DQMServices/Core/src/MonitorElement.cc

Go to the documentation of this file.

#define __STDC_FORMAT_MACROS 1
#define DQM_ROOT_METHODS 1
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQMServices/Core/interface/QTest.h"
#include "DQMServices/Core/src/DQMError.h"
#include "TClass.h"
#include "TMath.h"
#include "TList.h"
#include <iostream>
#include <cassert>
#include <cfloat>
#include <inttypes.h>

static TH1 *
checkRootObject(const std::string &name, TObject *tobj, const char *func, int reqdim)
{
  if (! tobj)
    raiseDQMError("MonitorElement", "Method '%s' cannot be invoked on monitor"
                  " element '%s' because it is not a ROOT object.",
                  func, name.c_str());

  TH1 *h = static_cast<TH1 *>(tobj);
  int ndim = h->GetDimension();
  if (reqdim < 0 || reqdim > ndim)
    raiseDQMError("MonitorElement", "Method '%s' cannot be invoked on monitor"
                  " element '%s' because it requires %d dimensions; this"
                  " object of type '%s' has %d dimensions",
                  func, name.c_str(), reqdim, typeid(*h).name(), ndim);

  return h;
}

MonitorElement *
MonitorElement::initialise(Kind kind)
{
  switch (kind)
  {
  case DQM_KIND_INT:
  case DQM_KIND_REAL:
  case DQM_KIND_STRING:
  case DQM_KIND_TH1F:
  case DQM_KIND_TH1S:
  case DQM_KIND_TH1D:
  case DQM_KIND_TH2F:
  case DQM_KIND_TH2S:
  case DQM_KIND_TH2D:
  case DQM_KIND_TH3F:
  case DQM_KIND_TPROFILE:
  case DQM_KIND_TPROFILE2D:
    data_.flags &= ~DQMNet::DQM_PROP_TYPE_MASK;
    data_.flags |= kind;
    break;

  default:
    raiseDQMError("MonitorElement", "cannot initialise monitor element"
                  " to invalid type %d", (int) kind);
  }

  return this;
}

MonitorElement *
MonitorElement::initialise(Kind kind, TH1 *rootobj)
{
  initialise(kind);
  switch (kind)
  {
  case DQM_KIND_TH1F:
    assert(dynamic_cast<TH1F *>(rootobj));
    assert(! reference_ || dynamic_cast<TH1F *>(reference_));
    object_ = rootobj;
    break;

  case DQM_KIND_TH1S:
    assert(dynamic_cast<TH1S *>(rootobj));
    assert(! reference_ || dynamic_cast<TH1S *>(reference_));
    object_ = rootobj;
    break;

  case DQM_KIND_TH1D:
    assert(dynamic_cast<TH1D *>(rootobj));
    assert(! reference_ || dynamic_cast<TH1D *>(reference_));
    object_ = rootobj;
    break;

  case DQM_KIND_TH2F:
    assert(dynamic_cast<TH2F *>(rootobj));
    assert(! reference_ || dynamic_cast<TH2F *>(reference_));
    object_ = rootobj;
    break;

  case DQM_KIND_TH2S:
    assert(dynamic_cast<TH2S *>(rootobj));
    assert(! reference_ || dynamic_cast<TH2S *>(reference_));
    object_ = rootobj;
    break;

  case DQM_KIND_TH2D:
    assert(dynamic_cast<TH2D *>(rootobj));
    assert(! reference_ || dynamic_cast<TH1D *>(reference_));
    object_ = rootobj;
    break;

  case DQM_KIND_TH3F:
    assert(dynamic_cast<TH3F *>(rootobj));
    assert(! reference_ || dynamic_cast<TH3F *>(reference_));
    object_ = rootobj;
    break;

  case DQM_KIND_TPROFILE:
    assert(dynamic_cast<TProfile *>(rootobj));
    assert(! reference_ || dynamic_cast<TProfile *>(reference_));
    object_ = rootobj;
    break;

  case DQM_KIND_TPROFILE2D:
    assert(dynamic_cast<TProfile2D *>(rootobj));
    assert(! reference_ || dynamic_cast<TProfile2D *>(reference_));
    object_ = rootobj;
    break;

  default:
    raiseDQMError("MonitorElement", "cannot initialise monitor element"
                  " as a root object with type %d", (int) kind);
  }

  if (reference_)
    data_.flags |= DQMNet::DQM_PROP_HAS_REFERENCE;

  return this;
}

MonitorElement *
MonitorElement::initialise(Kind kind, const std::string &value)
{
  initialise(kind);
  if (kind == DQM_KIND_STRING)
    scalar_.str = value;
  else
    raiseDQMError("MonitorElement", "cannot initialise monitor element"
                  " as a string with type %d", (int) kind);

  return this;
}

MonitorElement::MonitorElement(void)
  : object_(0),
    reference_(0),
    refvalue_(0)
{
  data_.version = 0;
  data_.dirname = 0;
  data_.tag = 0;
  data_.flags = DQM_KIND_INVALID | DQMNet::DQM_PROP_NEW;
  scalar_.num = 0;
  scalar_.real = 0;
}

MonitorElement::MonitorElement(const std::string *path, const std::string &name)
  : object_(0),
    reference_(0),
    refvalue_(0)
{
  data_.version = 0;
  data_.dirname = path;
  data_.objname = name;
  data_.tag = 0;
  data_.flags = DQM_KIND_INVALID | DQMNet::DQM_PROP_NEW;
  scalar_.num = 0;
  scalar_.real = 0;
}

MonitorElement::MonitorElement(const MonitorElement &x)
  : data_(x.data_),
    scalar_(x.scalar_),
    object_(x.object_),
    reference_(x.reference_),
    refvalue_(x.refvalue_),
    qreports_(x.qreports_)
{
  if (object_)
    object_ = static_cast<TH1 *>(object_->Clone());

  if (refvalue_)
    refvalue_ = static_cast<TH1 *>(refvalue_->Clone());
}

MonitorElement &
MonitorElement::operator=(const MonitorElement &x)
{
  if (this != &x)
  {
    delete object_;
    delete refvalue_;

    data_ = x.data_;
    scalar_ = x.scalar_;
    object_ = x.object_;
    reference_ = x.reference_;
    refvalue_ = x.refvalue_;
    qreports_ = x.qreports_;

    if (object_)
      object_ = static_cast<TH1 *>(object_->Clone());

    if (refvalue_)
      refvalue_ = static_cast<TH1 *>(refvalue_->Clone());
  }

  return *this;
}

MonitorElement::~MonitorElement(void)
{
  delete object_;
  delete refvalue_;
}

void
MonitorElement::Fill(std::string &value)
{
  update();
  if (kind() == DQM_KIND_STRING)
    scalar_.str = value;
  else
    incompatible(__PRETTY_FUNCTION__);
}

void
MonitorElement::Fill(double x)
{
  update();
  if (kind() == DQM_KIND_INT)
    scalar_.num = static_cast<int64_t>(x);
  else if (kind() == DQM_KIND_REAL)
    scalar_.real = x;
  else if (kind() == DQM_KIND_TH1F)
    accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Fill(x, 1);
  else if (kind() == DQM_KIND_TH1S)
    accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Fill(x, 1);
  else if (kind() == DQM_KIND_TH1D)
    accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Fill(x, 1);
  else
    incompatible(__PRETTY_FUNCTION__);
}

void
MonitorElement::doFill(int64_t x)
{
  update();
  if (kind() == DQM_KIND_INT)
    scalar_.num = static_cast<int64_t>(x);
  else if (kind() == DQM_KIND_REAL)
    scalar_.real = static_cast<double>(x);
  else if (kind() == DQM_KIND_TH1F)
    accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Fill(static_cast<double>(x), 1);
  else if (kind() == DQM_KIND_TH1S)
    accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Fill(static_cast<double>(x), 1);
  else if (kind() == DQM_KIND_TH1D)
    accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Fill(static_cast<double>(x), 1);
  else
    incompatible(__PRETTY_FUNCTION__);
}

void
MonitorElement::Fill(double x, double yw)
{
  update();
  if (kind() == DQM_KIND_TH1F)
    accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Fill(x, yw);
  else if (kind() == DQM_KIND_TH1S)
    accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Fill(x, yw);
  else if (kind() == DQM_KIND_TH1D)
    accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Fill(x, yw);
  else if (kind() == DQM_KIND_TH2F)
    static_cast<TH2F *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, yw, 1);
  else if (kind() == DQM_KIND_TH2S)
    static_cast<TH2S *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, yw, 1);
  else if (kind() == DQM_KIND_TH2D)
    static_cast<TH2D *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, yw, 1);
  else if (kind() == DQM_KIND_TPROFILE)
    static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
      ->Fill(x, yw, 1);
  else
    incompatible(__PRETTY_FUNCTION__);
}

void
MonitorElement::ShiftFillLast(double y, double ye, int xscale)
{
  update();
  if (kind() == DQM_KIND_TH1F 
      || kind() == DQM_KIND_TH1S 
      || kind() == DQM_KIND_TH1D) 
  {
    int nbins = getNbinsX();
    int entries = (int)getEntries();
    // first fill bins from left to right
    int index = entries + 1 ;
    int xlow = 2 ; int xup = nbins ;
    // if more entries than bins then start shifting
    if ( entries >= nbins ) 
    {
      index = nbins;
      xlow = entries - nbins + 3 ; xup = entries+1 ;
      // average first bin
      double y1 = getBinContent(1);
      double y2 = getBinContent(2);
      double y1err = getBinError(1);
      double y2err = getBinError(2);
      double N = entries - nbins + 1.;
      if ( ye == 0. || y1err == 0. || y2err == 0.) 
      {
        // for errors zero calculate unweighted mean and its error
        double sum = N*y1 + y2;
        y1 = sum/(N+1.) ;
        // FIXME check if correct
        double s=(N+1.)*(N*y1*y1 + y2*y2) - sum*sum;
        if (s>=0.) 
          y1err = sqrt(s)/(N+1.);  
        else
          y1err = 0.;
      }
      else 
      {
        // for errors non-zero calculate weighted mean and its error
        double denom = (1./y1err + 1./y2err);
        double mean = (y1/y1err + y2/y2err)/denom;
        // FIXME check if correct
        y1err = sqrt(((y1-mean)*(y1-mean)/y1err +
                      (y2-mean)*(y2-mean)/y2err)/denom/2.);
        y1 = mean; // set y1 to mean for filling below
      }
      setBinContent(1,y1);
      setBinError(1,y1err);
      // shift remaining bins to the left
      for ( int i = 3; i <= nbins ; i++) 
      {
        setBinContent(i-1,getBinContent(i));
        setBinError(i-1,getBinError(i));
      }
    }
    // fill last bin with new values
    setBinContent(index,y);
    setBinError(index,ye); 
    // set entries
    setEntries(entries+1);
    // set axis labels and reset drawing option
    char buffer [10];
    sprintf (buffer, "%d", xlow*xscale); 
    std::string a(buffer); setBinLabel(2,a);
    sprintf (buffer, "%d", xup*xscale); 
    std::string b(buffer); setBinLabel(nbins,b);
    setBinLabel(1,"av.");
  }
  else
    incompatible(__PRETTY_FUNCTION__);
}
void
MonitorElement::Fill(double x, double y, double zw)
{
  update();
  if (kind() == DQM_KIND_TH2F)
    static_cast<TH2F *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, y, zw);
  else if (kind() == DQM_KIND_TH2S)
    static_cast<TH2S *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, y, zw);
  else if (kind() == DQM_KIND_TH2D)
    static_cast<TH2D *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, y, zw);
  else if (kind() == DQM_KIND_TH3F)
    static_cast<TH3F *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, y, zw, 1);
  else if (kind() == DQM_KIND_TPROFILE)
    static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, y, zw);
  else if (kind() == DQM_KIND_TPROFILE2D)
    static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, y, zw, 1);
  else
    incompatible(__PRETTY_FUNCTION__);
}

void
MonitorElement::Fill(double x, double y, double z, double w)
{
  update();
  if (kind() == DQM_KIND_TH3F)
    static_cast<TH3F *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, y, z, w);
  else if (kind() == DQM_KIND_TPROFILE2D)
    static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 2))
      ->Fill(x, y, z, w);
  else
    incompatible(__PRETTY_FUNCTION__);
}

void
MonitorElement::Reset(void)
{
  update();
  if (kind() == DQM_KIND_INT)
    scalar_.num = 0;
  else if (kind() == DQM_KIND_REAL)
    scalar_.real = 0;
  else if (kind() == DQM_KIND_STRING)
    scalar_.str.clear();
  else
    return accessRootObject(__PRETTY_FUNCTION__, 1)
      ->Reset();
}

void
MonitorElement::packScalarData(std::string &into, const char *prefix) const
{
  char buf[64];
  if (kind() == DQM_KIND_INT)
  {
    snprintf(buf, sizeof(buf), "%s%" PRId64, prefix, scalar_.num);
    into = buf;
  }
  else if (kind() == DQM_KIND_REAL)
  {
    snprintf(buf, sizeof(buf), "%s%.*g", prefix, DBL_DIG+2, scalar_.real);
    into = buf;
  }
  else if (kind() == DQM_KIND_STRING)
  {
    into.reserve(strlen(prefix) + scalar_.str.size());
    into += prefix;
    into += scalar_.str;
  }
  else
    incompatible(__PRETTY_FUNCTION__);
}

void
MonitorElement::packQualityData(std::string &into) const
{
  DQMNet::packQualityData(into, data_.qreports);
}

std::string
MonitorElement::valueString(void) const
{
  std::string result;
  if (kind() == DQM_KIND_INT)
    packScalarData(result, "i=");
  else if (kind() == DQM_KIND_REAL)
    packScalarData(result, "f=");
  else if (kind() == DQM_KIND_STRING)
    packScalarData(result, "s=");
  else
    incompatible(__PRETTY_FUNCTION__);

  return result;
}

std::string
MonitorElement::tagString(void) const
{
  std::string result;
  std::string val(valueString());
  result.reserve(6 + 2*data_.objname.size() + val.size());
  result += '<'; result += data_.objname; result += '>';
  result += val;
  result += '<'; result += '/'; result += data_.objname; result += '>';
  return result;
}

std::string
MonitorElement::tagLabelString(void) const
{
  char buf[32];
  std::string result;
  size_t len = sprintf(buf, "t=%" PRIu32, data_.tag);

  result.reserve(6 + 2*data_.objname.size() + len);
  result += '<'; result += data_.objname; result += '>';
  result += buf;
  result += '<'; result += '/'; result += data_.objname; result += '>';
  return result;
}

std::string
MonitorElement::qualityTagString(const DQMNet::QValue &qv) const
{
  char buf[64];
  std::string result;
  size_t titlelen = data_.objname.size() + qv.qtname.size() + 1;
  size_t buflen = sprintf(buf, "qr=st:%d:%.*g:", qv.code, DBL_DIG+2, qv.qtresult);

  result.reserve(7 + 2*titlelen + buflen + qv.algorithm.size() + qv.message.size());
  result += '<'; result += data_.objname; result += '.'; result += qv.qtname; result += '>';
  result += buf; result += qv.algorithm; result += ':'; result += qv.message;
  result += '<'; result += '/'; result += data_.objname; result += '.'; result += qv.qtname; result += '>';
  return result;
}

const QReport *
MonitorElement::getQReport(const std::string &qtname) const
{
  QReport *qr;
  DQMNet::QValue *qv;
  const_cast<MonitorElement *>(this)->getQReport(false, qtname, qr, qv);
  return qr;
}

std::vector<QReport *>
MonitorElement::getQReports(void) const
{
  std::vector<QReport *> result;
  result.reserve(qreports_.size());
  for (size_t i = 0, e = qreports_.size(); i != e; ++i)
  {
    const_cast<MonitorElement *>(this)->qreports_[i].qvalue_
      = const_cast<DQMNet::QValue *>(&data_.qreports[i]);
    result.push_back(const_cast<QReport *>(&qreports_[i]));
  }
  return result;
}

std::vector<QReport *>
MonitorElement::getQWarnings(void) const
{
  std::vector<QReport *> result;
  result.reserve(qreports_.size());
  for (size_t i = 0, e = qreports_.size(); i != e; ++i)
    if (data_.qreports[i].code == dqm::qstatus::WARNING)
    {
      const_cast<MonitorElement *>(this)->qreports_[i].qvalue_
        = const_cast<DQMNet::QValue *>(&data_.qreports[i]);
      result.push_back(const_cast<QReport *>(&qreports_[i]));
    }
  return result;
}

std::vector<QReport *>
MonitorElement::getQErrors(void) const
{
  std::vector<QReport *> result;
  result.reserve(qreports_.size());
  for (size_t i = 0, e = qreports_.size(); i != e; ++i)
    if (data_.qreports[i].code == dqm::qstatus::ERROR)
    {
      const_cast<MonitorElement *>(this)->qreports_[i].qvalue_
        = const_cast<DQMNet::QValue *>(&data_.qreports[i]);
      result.push_back(const_cast<QReport *>(&qreports_[i]));
    }
  return result;
}

std::vector<QReport *>
MonitorElement::getQOthers(void) const
{
  std::vector<QReport *> result;
  result.reserve(qreports_.size());
  for (size_t i = 0, e = qreports_.size(); i != e; ++i)
    if (data_.qreports[i].code != dqm::qstatus::STATUS_OK
        && data_.qreports[i].code != dqm::qstatus::WARNING
        && data_.qreports[i].code != dqm::qstatus::ERROR)
    {
      const_cast<MonitorElement *>(this)->qreports_[i].qvalue_
        = const_cast<DQMNet::QValue *>(&data_.qreports[i]);
      result.push_back(const_cast<QReport *>(&qreports_[i]));
    }
  return result;
}

void
MonitorElement::runQTests(void)
{
  assert(qreports_.size() == data_.qreports.size());

  // Rerun quality tests where the ME or the quality algorithm was modified.
  bool dirty = wasUpdated();
  for (size_t i = 0, e = data_.qreports.size(); i < e; ++i)
  {
    DQMNet::QValue &qv = data_.qreports[i];
    QReport &qr = qreports_[i];
    QCriterion *qc = qr.qcriterion_;
    qr.qvalue_ = &qv;

    // if (qc && (dirty || qc->wasModified()))  // removed for new QTest (abm-090503)
    if (qc && dirty)
    {
      assert(qc->getName() == qv.qtname);
      std::string oldMessage = qv.message;
      int oldStatus = qv.code;

      qc->runTest(this, qr, qv);

      if (oldStatus != qv.code || oldMessage != qv.message)
        update();
    }
  }

  // Update QReport statistics.
  updateQReportStats();
}

void
MonitorElement::incompatible(const char *func) const
{
  raiseDQMError("MonitorElement", "Method '%s' cannot be invoked on monitor"
                " element '%s'", func, data_.objname.c_str());
}

TH1 *
MonitorElement::accessRootObject(const char *func, int reqdim) const
{
  if (kind() < DQM_KIND_TH1F)
    raiseDQMError("MonitorElement", "Method '%s' cannot be invoked on monitor"
                  " element '%s' because it is not a root object",
                  func, data_.objname.c_str());

  return checkRootObject(data_.objname, object_, func, reqdim);
}

/*** getter methods (wrapper around ROOT methods) ****/
// 
double
MonitorElement::getMean(int axis /* = 1 */) const
{ return accessRootObject(__PRETTY_FUNCTION__, axis-1)
    ->GetMean(axis); }

double
MonitorElement::getMeanError(int axis /* = 1 */) const
{ return accessRootObject(__PRETTY_FUNCTION__, axis-1)
    ->GetMeanError(axis); }

double
MonitorElement::getRMS(int axis /* = 1 */) const
{ return accessRootObject(__PRETTY_FUNCTION__, axis-1)
    ->GetRMS(axis); }

double
MonitorElement::getRMSError(int axis /* = 1 */) const
{ return accessRootObject(__PRETTY_FUNCTION__, axis-1)
    ->GetRMSError(axis); }

int
MonitorElement::getNbinsX(void) const
{ return accessRootObject(__PRETTY_FUNCTION__, 1)
    ->GetNbinsX(); }

int
MonitorElement::getNbinsY(void) const
{ return accessRootObject(__PRETTY_FUNCTION__, 2)
    ->GetNbinsY(); }

int
MonitorElement::getNbinsZ(void) const
{ return accessRootObject(__PRETTY_FUNCTION__, 3)
    ->GetNbinsZ(); }

double
MonitorElement::getBinContent(int binx) const
{ return accessRootObject(__PRETTY_FUNCTION__, 1)
    ->GetBinContent(binx); }

double
MonitorElement::getBinContent(int binx, int biny) const
{ return accessRootObject(__PRETTY_FUNCTION__, 2)
    ->GetBinContent(binx, biny); }

double
MonitorElement::getBinContent(int binx, int biny, int binz) const
{ return accessRootObject(__PRETTY_FUNCTION__, 3)
    ->GetBinContent(binx, biny, binz); }

double
MonitorElement::getBinError(int binx) const
{ return accessRootObject(__PRETTY_FUNCTION__, 1)
    ->GetBinError(binx); }

double
MonitorElement::getBinError(int binx, int biny) const
{ return accessRootObject(__PRETTY_FUNCTION__, 2)
    ->GetBinError(binx, biny); }

double
MonitorElement::getBinError(int binx, int biny, int binz) const
{ return accessRootObject(__PRETTY_FUNCTION__, 3)
    ->GetBinError(binx, biny, binz); }

double
MonitorElement::getEntries(void) const
{ return accessRootObject(__PRETTY_FUNCTION__, 1)
    ->GetEntries(); }

double
MonitorElement::getBinEntries(int bin) const
{
  if (kind() == DQM_KIND_TPROFILE)
    return static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
      ->GetBinEntries(bin);
  else if (kind() == DQM_KIND_TPROFILE2D)
    return static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 1))
      ->GetBinEntries(bin);
  else
  {
    incompatible(__PRETTY_FUNCTION__);
    return 0;
  }
}

double
MonitorElement::getYmin(void) const
{
  if (kind() == DQM_KIND_TPROFILE)
    return static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
      ->GetYmin();
  else
  {
    incompatible(__PRETTY_FUNCTION__);
    return 0;
  }
}

double
MonitorElement::getYmax(void) const
{
  if (kind() == DQM_KIND_TPROFILE)
    return static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
      ->GetYmax();
  else
  {
    incompatible(__PRETTY_FUNCTION__);
    return 0;
  }
}

std::string
MonitorElement::getAxisTitle(int axis /* = 1 */) const
{ return getAxis(__PRETTY_FUNCTION__, axis)
    ->GetTitle(); }

std::string
MonitorElement::getTitle(void) const
{ return accessRootObject(__PRETTY_FUNCTION__, 1)
    ->GetTitle(); }

/*** setter methods (wrapper around ROOT methods) ****/
// 
void
MonitorElement::setBinContent(int binx, double content)
{
  update();
  accessRootObject(__PRETTY_FUNCTION__, 1)
    ->SetBinContent(binx, content);
}

void
MonitorElement::setBinContent(int binx, int biny, double content)
{
  update();
  accessRootObject(__PRETTY_FUNCTION__, 2)
    ->SetBinContent(binx, biny, content); }

void
MonitorElement::setBinContent(int binx, int biny, int binz, double content)
{
  update();
  accessRootObject(__PRETTY_FUNCTION__, 3)
    ->SetBinContent(binx, biny, binz, content); }

void
MonitorElement::setBinError(int binx, double error)
{
  update();
  accessRootObject(__PRETTY_FUNCTION__, 1)
    ->SetBinError(binx, error);
}

void
MonitorElement::setBinError(int binx, int biny, double error)
{
  update();
  accessRootObject(__PRETTY_FUNCTION__, 2)
    ->SetBinError(binx, biny, error);
}

void
MonitorElement::setBinError(int binx, int biny, int binz, double error)
{
  update();
  accessRootObject(__PRETTY_FUNCTION__, 3)
    ->SetBinError(binx, biny, binz, error);
}

void
MonitorElement::setBinEntries(int bin, double nentries)
{
  update();
  if (kind() == DQM_KIND_TPROFILE)
    static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1))
      ->SetBinEntries(bin, nentries);
  else if (kind() == DQM_KIND_TPROFILE2D)
    static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 1))
      ->SetBinEntries(bin, nentries);
  else
    incompatible(__PRETTY_FUNCTION__);
}

void
MonitorElement::setEntries(double nentries)
{
  update();
  accessRootObject(__PRETTY_FUNCTION__, 1)
    ->SetEntries(nentries);
}

void
MonitorElement::setBinLabel(int bin, const std::string &label, int axis /* = 1 */)
{
  update();
  if ( getAxis(__PRETTY_FUNCTION__, axis)->GetNbins() >= bin ) 
  {
    getAxis(__PRETTY_FUNCTION__, axis)
      ->SetBinLabel(bin, label.c_str());
  }
  else
  {
    //  edm::LogWarning ("MonitorElement") 
    std::cout << "*** MonitorElement: WARNING:"
              <<"setBinLabel: attempting to set label of non-existent bin number \n";
  }
}

void
MonitorElement::setAxisRange(double xmin, double xmax, int axis /* = 1 */)
{
  update();
  getAxis(__PRETTY_FUNCTION__, axis)
    ->SetRangeUser(xmin, xmax);
}

void
MonitorElement::setAxisTitle(const std::string &title, int axis /* = 1 */)
{
  update();
  getAxis(__PRETTY_FUNCTION__, axis)
    ->SetTitle(title.c_str());
}

void
MonitorElement::setAxisTimeDisplay(int value, int axis /* = 1 */)
{
  update();
  getAxis(__PRETTY_FUNCTION__, axis)
    ->SetTimeDisplay(value);
}

void
MonitorElement::setAxisTimeFormat(const char *format /* = "" */, int axis /* = 1 */)
{
  update();
  getAxis(__PRETTY_FUNCTION__, axis)
    ->SetTimeFormat(format);
}

void
MonitorElement::setAxisTimeOffset(double toffset, const char *option /* ="local" */, int axis /* = 1 */)
{
  update();
  getAxis(__PRETTY_FUNCTION__, axis)
    ->SetTimeOffset(toffset, option);
}

void
MonitorElement::setTitle(const std::string &title)
{
  update();
  accessRootObject(__PRETTY_FUNCTION__, 1)
    ->SetTitle(title.c_str());
}

TAxis *
MonitorElement::getAxis(const char *func, int axis) const
{
  TH1 *h = accessRootObject(func, axis-1);
  TAxis *a = 0;
  if (axis == 1)
    a = h->GetXaxis();
  else if (axis == 2)
    a = h->GetYaxis();
  else if (axis == 3)
    a = h->GetZaxis();

  if (! a)
    raiseDQMError("MonitorElement", "No such axis %d in monitor element"
                  " '%s' of type '%s'", axis, data_.objname.c_str(),
                  typeid(*h).name());

  return a;
}

// ------------ Operations for MEs that are normally never reset ---------

void
MonitorElement::softReset(void)
{
  update();

  // Create the reference object the first time this is called.
  // On subsequent calls accumulate the current value to the
  // reference, and then reset the current value.  This way the
  // future contents will have the reference "subtracted".
  if (kind() == DQM_KIND_TH1F)
  {
    TH1F *orig = static_cast<TH1F *>(object_);
    TH1F *r = static_cast<TH1F *>(refvalue_);
    if (! r)
    {
      refvalue_ = r = new TH1F((std::string(orig->GetName()) + "_ref").c_str(),
                               orig->GetTitle(),
                               orig->GetNbinsX(),
                               orig->GetXaxis()->GetXmin(),
                               orig->GetXaxis()->GetXmax());
      r->SetDirectory(0);
      r->Reset();
    }

    r->Add(orig);
    orig->Reset();
  }
  else if (kind() == DQM_KIND_TH1S)
  {
    TH1S *orig = static_cast<TH1S *>(object_);
    TH1S *r = static_cast<TH1S *>(refvalue_);
    if (! r)
    {
      refvalue_ = r = new TH1S((std::string(orig->GetName()) + "_ref").c_str(),
                               orig->GetTitle(),
                               orig->GetNbinsX(),
                               orig->GetXaxis()->GetXmin(),
                               orig->GetXaxis()->GetXmax());
      r->SetDirectory(0);
      r->Reset();
    }

    r->Add(orig);
    orig->Reset();
  }
  else if (kind() == DQM_KIND_TH1D)
  {
    TH1D *orig = static_cast<TH1D *>(object_);
    TH1D *r = static_cast<TH1D *>(refvalue_);
    if (! r)
    {
      refvalue_ = r = new TH1D((std::string(orig->GetName()) + "_ref").c_str(),
                               orig->GetTitle(),
                               orig->GetNbinsX(),
                               orig->GetXaxis()->GetXmin(),
                               orig->GetXaxis()->GetXmax());
      r->SetDirectory(0);
      r->Reset();
    }

    r->Add(orig);
    orig->Reset();
  }
  else if (kind() == DQM_KIND_TH2F)
  {
    TH2F *orig = static_cast<TH2F *>(object_);
    TH2F *r = static_cast<TH2F *>(refvalue_);
    if (! r)
    {
      refvalue_ = r = new TH2F((std::string(orig->GetName()) + "_ref").c_str(),
                               orig->GetTitle(),
                               orig->GetNbinsX(),
                               orig->GetXaxis()->GetXmin(),
                               orig->GetXaxis()->GetXmax(),
                               orig->GetNbinsY(),
                               orig->GetYaxis()->GetXmin(),
                               orig->GetYaxis()->GetXmax());
      r->SetDirectory(0);
      r->Reset();
    }

    r->Add(orig);
    orig->Reset();
  }
  else if (kind() == DQM_KIND_TH2S)
  {
    TH2S *orig = static_cast<TH2S *>(object_);
    TH2S *r = static_cast<TH2S *>(refvalue_);
    if (! r)
    {
      refvalue_ = r = new TH2S((std::string(orig->GetName()) + "_ref").c_str(),
                               orig->GetTitle(),
                               orig->GetNbinsX(),
                               orig->GetXaxis()->GetXmin(),
                               orig->GetXaxis()->GetXmax(),
                               orig->GetNbinsY(),
                               orig->GetYaxis()->GetXmin(),
                               orig->GetYaxis()->GetXmax());
      r->SetDirectory(0);
      r->Reset();
    }

    r->Add(orig);
    orig->Reset();
  }
  else if (kind() == DQM_KIND_TH2D)
  {
    TH2D *orig = static_cast<TH2D *>(object_);
    TH2D *r = static_cast<TH2D *>(refvalue_);
    if (! r)
    {
      refvalue_ = r = new TH2D((std::string(orig->GetName()) + "_ref").c_str(),
                               orig->GetTitle(),
                               orig->GetNbinsX(),
                               orig->GetXaxis()->GetXmin(),
                               orig->GetXaxis()->GetXmax(),
                               orig->GetNbinsY(),
                               orig->GetYaxis()->GetXmin(),
                               orig->GetYaxis()->GetXmax());
      r->SetDirectory(0);
      r->Reset();
    }

    r->Add(orig);
    orig->Reset();
  }
  else if (kind() == DQM_KIND_TH3F)
  {
    TH3F *orig = static_cast<TH3F *>(object_);
    TH3F *r = static_cast<TH3F *>(refvalue_);
    if (! r)
    {
      refvalue_ = r = new TH3F((std::string(orig->GetName()) + "_ref").c_str(),
                               orig->GetTitle(),
                               orig->GetNbinsX(),
                               orig->GetXaxis()->GetXmin(),
                               orig->GetXaxis()->GetXmax(),
                               orig->GetNbinsY(),
                               orig->GetYaxis()->GetXmin(),
                               orig->GetYaxis()->GetXmax(),
                               orig->GetNbinsZ(),
                               orig->GetZaxis()->GetXmin(),
                               orig->GetZaxis()->GetXmax());
      r->SetDirectory(0);
      r->Reset();
    }

    r->Add(orig);
    orig->Reset();
  }
  else if (kind() == DQM_KIND_TPROFILE)
  {
    TProfile *orig = static_cast<TProfile *>(object_);
    TProfile *r = static_cast<TProfile *>(refvalue_);
    if (! r)
    {
      refvalue_ = r = new TProfile((std::string(orig->GetName()) + "_ref").c_str(),
                                   orig->GetTitle(),
                                   orig->GetNbinsX(),
                                   orig->GetXaxis()->GetXmin(),
                                   orig->GetXaxis()->GetXmax(),
                                   orig->GetYaxis()->GetXmin(),
                                   orig->GetYaxis()->GetXmax(),
                                   orig->GetErrorOption());
      r->SetDirectory(0);
      r->Reset();
    }

    addProfiles(r, orig, r, 1, 1);
    orig->Reset();
  }
  else if (kind() == DQM_KIND_TPROFILE2D)
  {
    TProfile2D *orig = static_cast<TProfile2D *>(object_);
    TProfile2D *r = static_cast<TProfile2D *>(refvalue_);
    if (! r)
    {
      refvalue_ = r = new TProfile2D((std::string(orig->GetName()) + "_ref").c_str(),
                                     orig->GetTitle(),
                                     orig->GetNbinsX(),
                                     orig->GetXaxis()->GetXmin(),
                                     orig->GetXaxis()->GetXmax(),
                                     orig->GetNbinsY(),
                                     orig->GetYaxis()->GetXmin(),
                                     orig->GetYaxis()->GetXmax(),
                                     orig->GetZaxis()->GetXmin(),
                                     orig->GetZaxis()->GetXmax(),
                                     orig->GetErrorOption());
      r->SetDirectory(0);
      r->Reset();
    }

    addProfiles(r, orig, r, 1, 1);
    orig->Reset();
  }
  else
    incompatible(__PRETTY_FUNCTION__);
}

void
MonitorElement::disableSoftReset(void)
{
  if (refvalue_)
  {
    if (kind() == DQM_KIND_TH1F
        || kind() == DQM_KIND_TH1S
        || kind() == DQM_KIND_TH1D
        || kind() == DQM_KIND_TH2F
        || kind() == DQM_KIND_TH2S
        || kind() == DQM_KIND_TH2D
        || kind() == DQM_KIND_TH3F)
    {
      TH1 *orig = static_cast<TH1 *>(object_);
      orig->Add(refvalue_);
    }
    else if (kind() == DQM_KIND_TPROFILE)
    {
      TProfile *orig = static_cast<TProfile *>(object_);
      TProfile *r = static_cast<TProfile *>(refvalue_);
      addProfiles(orig, r, orig, 1, 1);
    }
    else if (kind() == DQM_KIND_TPROFILE2D)
    {
      TProfile2D *orig = static_cast<TProfile2D *>(object_);
      TProfile2D *r = static_cast<TProfile2D *>(refvalue_);
      addProfiles(orig, r, orig, 1, 1);
    }
    else
      incompatible(__PRETTY_FUNCTION__);

    delete refvalue_;
    refvalue_ = 0;
  }
}
  
// implementation: Giuseppe.Della-Ricca@ts.infn.it
// Can be called with sum = h1 or sum = h2
void
MonitorElement::addProfiles(TProfile *h1, TProfile *h2, TProfile *sum, float c1, float c2)
{
  assert(h1);
  assert(h2);
  assert(sum);

  static const Int_t NUM_STAT = 6;
  Double_t stats1[NUM_STAT];
  Double_t stats2[NUM_STAT];
  Double_t stats3[NUM_STAT];

  bool isRebinOn = sum->TestBit(TH1::kCanRebin);
  sum->ResetBit(TH1::kCanRebin);

  for (Int_t i = 0; i < NUM_STAT; ++i)
    stats1[i] = stats2[i] = stats3[i] = 0;

  h1->GetStats(stats1);
  h2->GetStats(stats2);

  for (Int_t i = 0; i < NUM_STAT; ++i)
    stats3[i] = c1*stats1[i] + c2*stats2[i];

  stats3[1] = c1*TMath::Abs(c1)*stats1[1]
              + c2*TMath::Abs(c2)*stats2[1];

  Double_t entries = c1*h1->GetEntries() + c2* h2->GetEntries();
  TArrayD* h1sumw2 = h1->GetSumw2();
  TArrayD* h2sumw2 = h2->GetSumw2();
  for (Int_t bin = 0, nbin = sum->GetNbinsX()+1; bin <= nbin; ++bin) 
  {
    Double_t entries = c1*h1->GetBinEntries(bin)
                       + c2*h2->GetBinEntries(bin);
    Double_t content = c1*h1->GetBinEntries(bin)*h1->GetBinContent(bin)
                       + c2*h2->GetBinEntries(bin)*h2->GetBinContent(bin);
    Double_t error = TMath::Sqrt(c1*TMath::Abs(c1)*h1sumw2->fArray[bin]
                                 + c2*TMath::Abs(c2)*h2sumw2->fArray[bin]);
    sum->SetBinContent(bin, content);
    sum->SetBinError(bin, error);
    sum->SetBinEntries(bin, entries);
  }
  
  sum->SetEntries(entries);
  sum->PutStats(stats3);
  if (isRebinOn) sum->SetBit(TH1::kCanRebin);
}

// implementation: Giuseppe.Della-Ricca@ts.infn.it
// Can be called with sum = h1 or sum = h2
void
MonitorElement::addProfiles(TProfile2D *h1, TProfile2D *h2, TProfile2D *sum, float c1, float c2)
{
  assert(h1);
  assert(h2);
  assert(sum);

  static const Int_t NUM_STAT = 9;
  Double_t stats1[NUM_STAT];
  Double_t stats2[NUM_STAT];
  Double_t stats3[NUM_STAT];

  bool isRebinOn = sum->TestBit(TH1::kCanRebin);
  sum->ResetBit(TH1::kCanRebin);

  for (Int_t i = 0; i < NUM_STAT; ++i)
    stats1[i] = stats2[i] = stats3[i] = 0;

  h1->GetStats(stats1);
  h2->GetStats(stats2);

  for (Int_t i = 0; i < NUM_STAT; i++)
    stats3[i] = c1*stats1[i] + c2*stats2[i];

  stats3[1] = c1*TMath::Abs(c1)*stats1[1]
              + c2*TMath::Abs(c2)*stats2[1];

  Double_t entries = c1*h1->GetEntries() + c2*h2->GetEntries();
  TArrayD *h1sumw2 = h1->GetSumw2();
  TArrayD *h2sumw2 = h2->GetSumw2();
  for (Int_t xbin = 0, nxbin = sum->GetNbinsX()+1; xbin <= nxbin; ++xbin)
    for (Int_t ybin = 0, nybin = sum->GetNbinsY()+1; ybin <= nybin; ++ybin)
    {
      Int_t bin = sum->GetBin(xbin, ybin);
      Double_t entries = c1*h1->GetBinEntries(bin)
                         + c2*h2->GetBinEntries(bin);
      Double_t content = c1*h1->GetBinEntries(bin)*h1->GetBinContent(bin)
                         + c2*h2->GetBinEntries(bin)*h2->GetBinContent(bin);
      Double_t error = TMath::Sqrt(c1*TMath::Abs(c1)*h1sumw2->fArray[bin]
                                   + c2*TMath::Abs(c2)*h2sumw2->fArray[bin]);

      sum->SetBinContent(bin, content);
      sum->SetBinError(bin, error);
      sum->SetBinEntries(bin, entries);
    }
  sum->SetEntries(entries);
  sum->PutStats(stats3);
  if (isRebinOn) sum->SetBit(TH1::kCanRebin);
}

void
MonitorElement::copyFunctions(TH1 *from, TH1 *to)
{
  // will copy functions only if local-copy and original-object are equal
  // (ie. no soft-resetting or accumulating is enabled)
  if (isSoftResetEnabled() || isAccumulateEnabled())
    return;

  update();
  TList *fromf = from->GetListOfFunctions();
  TList *tof   = to->GetListOfFunctions();
  for (int i = 0, nfuncs = fromf ? fromf->GetSize() : 0; i < nfuncs; ++i)
  {
    TObject *obj = fromf->At(i);
    // not interested in statistics
    if (!strcmp(obj->IsA()->GetName(), "TPaveStats"))
      continue;

    if (TF1 *fn = dynamic_cast<TF1 *>(obj))
      tof->Add(new TF1(*fn));
    //else if (dynamic_cast<TPaveStats *>(obj))
    //  ; // FIXME? tof->Add(new TPaveStats(*stats));
    else
      raiseDQMError("MonitorElement", "Cannot extract function '%s' of type"
                    " '%s' from monitor element '%s' for a copy",
                    obj->GetName(), obj->IsA()->GetName(), data_.objname.c_str());
  }
}

void
MonitorElement::copyFrom(TH1 *from)
{
  TH1 *orig = accessRootObject(__PRETTY_FUNCTION__, 1);
  if (orig->GetTitle() != from->GetTitle())
    orig->SetTitle(from->GetTitle());

  if (!isAccumulateEnabled())
    orig->Reset();

  if (isSoftResetEnabled())
  {
    if (kind() == DQM_KIND_TH1F
        || kind() == DQM_KIND_TH1S
        || kind() == DQM_KIND_TH1D
        || kind() == DQM_KIND_TH2F
        || kind() == DQM_KIND_TH2S
        || kind() == DQM_KIND_TH2D
        || kind() == DQM_KIND_TH3F)
      // subtract "reference"
      orig->Add(from, refvalue_, 1, -1);
    else if (kind() == DQM_KIND_TPROFILE)
      // subtract "reference"
      addProfiles(static_cast<TProfile *>(from),
                  static_cast<TProfile *>(refvalue_),
                  static_cast<TProfile *>(orig),
                  1, -1);
    else if (kind() == DQM_KIND_TPROFILE2D)
      // subtract "reference"
      addProfiles(static_cast<TProfile2D *>(from),
                  static_cast<TProfile2D *>(refvalue_),
                  static_cast<TProfile2D *>(orig),
                  1, -1);
    else
      incompatible(__PRETTY_FUNCTION__);
  }
  else
    orig->Add(from);

  copyFunctions(from, orig);
}

// --- Operations on MEs that are normally reset at end of monitoring cycle ---
void
MonitorElement::getQReport(bool create, const std::string &qtname, QReport *&qr, DQMNet::QValue *&qv)
{
  assert(qreports_.size() == data_.qreports.size());

  qr = 0;
  qv = 0;

  size_t pos = 0, end = qreports_.size();
  while (pos < end && data_.qreports[pos].qtname != qtname)
    ++pos;

  if (pos == end && ! create)
    return;
  else if (pos == end)
  {
    data_.qreports.push_back(DQMNet::QValue());
    qreports_.push_back(QReport(0, 0));

    DQMNet::QValue &q = data_.qreports.back();
    q.code = dqm::qstatus::DID_NOT_RUN;
    q.qtresult = 0;
    q.qtname = qtname;
    q.message = "NO_MESSAGE_ASSIGNED";
    q.algorithm = "UNKNOWN_ALGORITHM";
  }

  qr = &qreports_[pos];
  qv = &data_.qreports[pos];
}
  
void
MonitorElement::addQReport(const DQMNet::QValue &desc, QCriterion *qc)
{
  QReport *qr;
  DQMNet::QValue *qv;
  getQReport(true, desc.qtname, qr, qv);
  qr->qcriterion_ = qc;
  *qv = desc;
  update();
}

void
MonitorElement::addQReport(QCriterion *qc)
{
  QReport *qr;
  DQMNet::QValue *qv;
  getQReport(true, qc->getName(), qr, qv);
  qv->code = dqm::qstatus::DID_NOT_RUN;
  qv->message = "NO_MESSAGE_ASSIGNED";
  qr->qcriterion_ = qc;
  update();
}

void
MonitorElement::updateQReportStats(void)
{
  data_.flags &= ~DQMNet::DQM_PROP_REPORT_ALARM;
  for (size_t i = 0, e = data_.qreports.size(); i < e; ++i)
    switch (data_.qreports[i].code)
    {
    case dqm::qstatus::STATUS_OK:
      break;
    case dqm::qstatus::WARNING:
      data_.flags |= DQMNet::DQM_PROP_REPORT_WARN;
      break;
    case dqm::qstatus::ERROR:
      data_.flags |= DQMNet::DQM_PROP_REPORT_ERROR;
      break;
    default:
      data_.flags |= DQMNet::DQM_PROP_REPORT_OTHER;
      break;
    }
}

// -------------------------------------------------------------------
TObject *
MonitorElement::getRootObject(void) const
{
  const_cast<MonitorElement *>(this)->update();
  return object_;
}

TH1 *
MonitorElement::getTH1(void) const
{
  const_cast<MonitorElement *>(this)->update();
  return accessRootObject(__PRETTY_FUNCTION__, 0);
}

TH1F *
MonitorElement::getTH1F(void) const
{
  assert(kind() == DQM_KIND_TH1F);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH1F *>(accessRootObject(__PRETTY_FUNCTION__, 1));
}

TH1S *
MonitorElement::getTH1S(void) const
{
  assert(kind() == DQM_KIND_TH1S);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH1S *>(accessRootObject(__PRETTY_FUNCTION__, 1));
}

TH1D *
MonitorElement::getTH1D(void) const
{
  assert(kind() == DQM_KIND_TH1D);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH1D *>(accessRootObject(__PRETTY_FUNCTION__, 1));
}

TH2F *
MonitorElement::getTH2F(void) const
{
  assert(kind() == DQM_KIND_TH2F);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH2F *>(accessRootObject(__PRETTY_FUNCTION__, 2));
}

TH2S *
MonitorElement::getTH2S(void) const
{
  assert(kind() == DQM_KIND_TH2S);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH2S *>(accessRootObject(__PRETTY_FUNCTION__, 2));
}

TH2D *
MonitorElement::getTH2D(void) const
{
  assert(kind() == DQM_KIND_TH2D);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH2D *>(accessRootObject(__PRETTY_FUNCTION__, 2));
}

TH3F *
MonitorElement::getTH3F(void) const
{
  assert(kind() == DQM_KIND_TH3F);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH3F *>(accessRootObject(__PRETTY_FUNCTION__, 3));
}

TProfile *
MonitorElement::getTProfile(void) const
{
  assert(kind() == DQM_KIND_TPROFILE);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TProfile *>(accessRootObject(__PRETTY_FUNCTION__, 1));
}

TProfile2D *
MonitorElement::getTProfile2D(void) const
{
  assert(kind() == DQM_KIND_TPROFILE2D);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TProfile2D *>(accessRootObject(__PRETTY_FUNCTION__, 2));
}

// -------------------------------------------------------------------
TObject *
MonitorElement::getRefRootObject(void) const
{
  const_cast<MonitorElement *>(this)->update();
  return reference_;
}

TH1 *
MonitorElement::getRefTH1(void) const
{
  const_cast<MonitorElement *>(this)->update();
  return checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 0);
}

TH1F *
MonitorElement::getRefTH1F(void) const
{
  assert(kind() == DQM_KIND_TH1F);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH1F *>
    (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 1));
}

TH1S *
MonitorElement::getRefTH1S(void) const
{
  assert(kind() == DQM_KIND_TH1S);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH1S *>
    (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 1));
}

TH1D *
MonitorElement::getRefTH1D(void) const
{
  assert(kind() == DQM_KIND_TH1D);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH1D *>
    (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 1));
}

TH2F *
MonitorElement::getRefTH2F(void) const
{
  assert(kind() == DQM_KIND_TH2F);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH2F *>
    (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 2));
}

TH2S *
MonitorElement::getRefTH2S(void) const
{
  assert(kind() == DQM_KIND_TH2S);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH2S *>
    (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 2));
}

TH2D *
MonitorElement::getRefTH2D(void) const
{
  assert(kind() == DQM_KIND_TH2D);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH2D *>
    (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 2));
}

TH3F *
MonitorElement::getRefTH3F(void) const
{
  assert(kind() == DQM_KIND_TH3F);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TH3F *>
    (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 3));
}

TProfile *
MonitorElement::getRefTProfile(void) const
{
  assert(kind() == DQM_KIND_TPROFILE);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TProfile *>
    (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 1));
}

TProfile2D *
MonitorElement::getRefTProfile2D(void) const
{
  assert(kind() == DQM_KIND_TPROFILE2D);
  const_cast<MonitorElement *>(this)->update();
  return static_cast<TProfile2D *>
    (checkRootObject(data_.objname, reference_, __PRETTY_FUNCTION__, 2));
}