/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/DQM/EcalCommon/src/MESetNonObject.cc

Go to the documentation of this file.

#include "DQM/EcalCommon/interface/MESetNonObject.h"

#include "FWCore/Utilities/interface/Exception.h"

namespace ecaldqm
{
  MESetNonObject::MESetNonObject(std::string const& _fullpath, MEData const& _data, bool _readOnly/* = false*/) :
    MESet(_fullpath, _data, _readOnly)
  {
  }

  void
  MESetNonObject::book()
  {
    dqmStore_->setCurrentFolder(dir_);

    mes_.clear();

    MonitorElement* me(0);

    BinService::AxisSpecs* xaxis(data_->xaxis);
    BinService::AxisSpecs* yaxis(data_->yaxis);
    BinService::AxisSpecs* zaxis(data_->zaxis);

    switch(data_->kind) {
    case MonitorElement::DQM_KIND_REAL :
      me = dqmStore_->bookFloat(name_);
      break;

    case MonitorElement::DQM_KIND_TH1F :
      {
        if(!xaxis)
          throw cms::Exception("InvalidCall") << "No xaxis found for MESetNonObject" << std::endl;

        if(!xaxis->edges)
          me = dqmStore_->book1D(name_, name_, xaxis->nbins, xaxis->low, xaxis->high);
        else{
          float* edges(new float[xaxis->nbins + 1]);
          for(int i(0); i < xaxis->nbins + 1; i++)
            edges[i] = xaxis->edges[i];
          me = dqmStore_->book1D(name_, name_, xaxis->nbins, edges);
          delete [] edges;
        }
      }
      break;

    case MonitorElement::DQM_KIND_TPROFILE :
      {
        if(!xaxis)
          throw cms::Exception("InvalidCall") << "No xaxis found for MESetNonObject" << std::endl;

        double ylow, yhigh;
        if(!yaxis){
          ylow = -std::numeric_limits<double>::max();
          yhigh = std::numeric_limits<double>::max();
        }
        else{
          ylow = yaxis->low;
          yhigh = yaxis->high;
        }
        if(xaxis->edges)
          me = dqmStore_->bookProfile(name_, name_, xaxis->nbins, xaxis->edges, ylow, yhigh, "");
        else
          me = dqmStore_->bookProfile(name_, name_, xaxis->nbins, xaxis->low, xaxis->high, ylow, yhigh, "");

      }
      break;

    case MonitorElement::DQM_KIND_TH2F :
      {
        if(!xaxis || !yaxis)
          throw cms::Exception("InvalidCall") << "No x/yaxis found for MESetNonObject" << std::endl;

        if(!xaxis->edges || !yaxis->edges)
          me = dqmStore_->book2D(name_, name_, xaxis->nbins, xaxis->low, xaxis->high, yaxis->nbins, yaxis->low, yaxis->high);
        else{
          float* xedges(new float[xaxis->nbins + 1]);
          for(int i(0); i < xaxis->nbins + 1; i++)
            xedges[i] = xaxis->edges[i];
          float* yedges(new float[yaxis->nbins + 1]);
          for(int i(0); i < yaxis->nbins + 1; i++)
            yedges[i] = yaxis->edges[i];
          me = dqmStore_->book2D(name_, name_, xaxis->nbins, xedges, yaxis->nbins, yedges);
          delete [] xedges;
          delete [] yedges;
        }
      }
      break;

    case MonitorElement::DQM_KIND_TPROFILE2D :
      {
        if(!xaxis || !yaxis)
          throw cms::Exception("InvalidCall") << "No x/yaxis found for MESetNonObject" << std::endl;
        double high(0.), low(0.);
        if(zaxis){
          low = zaxis->low;
          high = zaxis->high;
        }
        else{
          low = -std::numeric_limits<double>::max();
          high = std::numeric_limits<double>::max();
        }

        me = dqmStore_->bookProfile2D(name_, name_, xaxis->nbins, xaxis->low, xaxis->high, yaxis->nbins, yaxis->low, yaxis->high, low, high, "");
      }
      break;

    default :
      throw cms::Exception("InvalidCall") << "MESetNonObject of type " << data_->kind << " not implemented" << std::endl;
    }

    mes_.push_back(me);

    active_ = true;
  }

  bool
  MESetNonObject::retrieve() const
  {
    mes_.clear();

    MonitorElement* me(dqmStore_->get(dir_ + "/" + name_));
    if(!me) return false;

    mes_.push_back(me);

    active_ = true;
    return true;
  }

  MESetNonObject::~MESetNonObject()
  {
  }

  void
  MESetNonObject::fill(double _x, double _wy/* = 1.*/, double _w/* = 1.*/)
  {
    if(mes_.size() == 0 || !mes_[0]) return;

    switch(data_->kind) {
    case MonitorElement::DQM_KIND_REAL :
      mes_[0]->Fill(_x);
      break;
    case MonitorElement::DQM_KIND_TH1F :
    case MonitorElement::DQM_KIND_TPROFILE :
      mes_[0]->Fill(_x, _wy);
      break;
    case MonitorElement::DQM_KIND_TH2F :
    case MonitorElement::DQM_KIND_TPROFILE2D :
      mes_[0]->Fill(_x, _wy, _w);
      break;
    default :
      break;
    }
  }

}