RunHistogramManager.cc

Go to the documentation of this file.

#include "DPGAnalysis/SiStripTools/interface/RunHistogramManager.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CommonTools/UtilAlgos/interface/TFileService.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TProfile.h"
#include "TProfile2D.h"

BaseHistoParams::BaseHistoParams() {}

BaseHistoParams::~BaseHistoParams() {}

/*
void BaseHistoParams::beginRun(const edm::Run& iRun, TFileDirectory& subrun) {

  beginRun(iRun.run(),subrun);

}
*/

RunHistogramManager::RunHistogramManager(edm::ConsumesCollector& iC, const bool fillHistograms)
    : _fillHistograms(fillHistograms),
      _histograms(),
      _conditionsInRunToken(iC.consumes<edm::ConditionsInRunBlock, edm::InRun>(edm::InputTag("conditionsInEdm"))) {}

RunHistogramManager::RunHistogramManager(edm::ConsumesCollector&& iC, const bool fillHistograms)
    : _fillHistograms(fillHistograms),
      _histograms(),
      _conditionsInRunToken(iC.consumes<edm::ConditionsInRunBlock, edm::InRun>(edm::InputTag("conditionsInEdm"))) {}

TH1F** RunHistogramManager::makeTH1F(
    const char* name, const char* title, const unsigned int nbinx, const double xmin, const double xmax) {
  TH1F** pointer = new TH1F*(nullptr);

  BaseHistoParams* hp = new HistoParams<TH1F>(pointer, "TH1F", name, title, nbinx, xmin, xmax);
  _histograms.push_back(hp);

  LogDebug("TH1Fmade") << "Histogram " << name << " " << title << " pre-booked:" << _histograms.size();

  return pointer;
}

RunHistogramManager::~RunHistogramManager() {
  for (std::vector<BaseHistoParams*>::const_iterator hp = _histograms.begin(); hp != _histograms.end(); ++hp) {
    delete *hp;
  }
  LogDebug("Destructor") << "All BaseHistoParams destroyed ";
}

TProfile** RunHistogramManager::makeTProfile(
    const char* name, const char* title, const unsigned int nbinx, const double xmin, const double xmax) {
  TProfile** pointer = new TProfile*(nullptr);

  BaseHistoParams* hp = new HistoParams<TProfile>(pointer, "TProfile", name, title, nbinx, xmin, xmax);
  _histograms.push_back(hp);

  LogDebug("TProfilemade") << "Histogram " << name << " " << title << " pre-booked:" << _histograms.size();

  return pointer;
}

TH2F** RunHistogramManager::makeTH2F(const char* name,
                                     const char* title,
                                     const unsigned int nbinx,
                                     const double xmin,
                                     const double xmax,
                                     const unsigned int nbiny,
                                     const double ymin,
                                     const double ymax) {
  TH2F** pointer = new TH2F*(nullptr);

  BaseHistoParams* hp = new HistoParams<TH2F>(pointer, "TH2F", name, title, nbinx, xmin, xmax, nbiny, ymin, ymax);
  _histograms.push_back(hp);

  LogDebug("TH2Fmade") << "Histogram " << name << " " << title << " pre-booked :" << _histograms.size();

  return pointer;
}

TProfile2D** RunHistogramManager::makeTProfile2D(const char* name,
                                                 const char* title,
                                                 const unsigned int nbinx,
                                                 const double xmin,
                                                 const double xmax,
                                                 const unsigned int nbiny,
                                                 const double ymin,
                                                 const double ymax) {
  TProfile2D** pointer = new TProfile2D*(nullptr);

  BaseHistoParams* hp =
      new HistoParams<TProfile2D>(pointer, "TProfile2D", name, title, nbinx, xmin, xmax, nbiny, ymin, ymax);
  _histograms.push_back(hp);

  LogDebug("TProfile2Dmade") << "Histogram " << name << " " << title << " pre-booked :" << _histograms.size();

  return pointer;
}

void RunHistogramManager::beginRun(const edm::Run& iRun) {
  edm::Service<TFileService> tfserv;
  beginRun(iRun, tfserv->tFileDirectory());
}

void RunHistogramManager::beginRun(const edm::Run& iRun, TFileDirectory& subdir) {
  if (!_fillHistograms) {
    beginRun(iRun.run(), subdir);
  } else {
    unsigned int fillnum = 0;

    edm::Handle<edm::ConditionsInRunBlock> cirb;
    iRun.getByToken(_conditionsInRunToken, cirb);

    if (!cirb.failedToGet() && cirb.isValid())
      fillnum = cirb->lhcFillNumber;

    beginRun(fillnum, subdir);
  }
}

void RunHistogramManager::beginRun(const unsigned int irun) {
  edm::Service<TFileService> tfserv;
  beginRun(irun, tfserv->tFileDirectory());
}

void RunHistogramManager::beginRun(const unsigned int irun, TFileDirectory& subdir) {
  // create/go to the run subdirectory

  char fillrun[30];

  if (!_fillHistograms) {
    sprintf(fillrun, "%s", "run");
  } else {
    sprintf(fillrun, "%s", "fill");
  }

  char dirname[300];
  sprintf(dirname, "%s_%d", fillrun, irun);
  TFileDirectory subrun = subdir.mkdir(dirname);

  // loop on the histograms and update the pointer references

  for (unsigned int ih = 0; ih < _histograms.size(); ++ih) {
    _histograms[ih]->beginRun(irun, subrun, fillrun);
  }
}