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);

  }
}