MTVHistoProducerAlgo.cc

Go to the documentation of this file.

#include "Validation/RecoTrack/interface/MTVHistoProducerAlgo.h"

void MTVHistoProducerAlgo::doProfileX(TH2 * th2, MonitorElement* me){
  if (th2->GetNbinsX()==me->getNbinsX()){
    TProfile * p1 = (TProfile*) th2->ProfileX();
    p1->Copy(*me->getTProfile());
    delete p1;
  } else {
    throw cms::Exception("MultiTrackValidator") << "Different number of bins!";
  }
}

void MTVHistoProducerAlgo::fillPlotFromVector(MonitorElement* h, std::vector<int>& vec) {
  for (unsigned int j=0; j<vec.size(); j++){
    h->setBinContent(j+1, vec[j]);
  }
}

void MTVHistoProducerAlgo::fillPlotFromVectors(MonitorElement* h, 
                           std::vector<int>& numerator, 
                           std::vector<int>& denominator,
                           std::string type){
  double value,err;
  for (unsigned int j=0; j<numerator.size(); j++){
    if (denominator[j]!=0){
      if (type=="effic"){
    value = ((double) numerator[j])*1./((double) denominator[j]);
        err = sqrt( value*(1-value)/(double) denominator[j] );
      } else if (type=="fakerate"){
    value = 1-((double) numerator[j])*1./((double) denominator[j]);
        err = sqrt( value*(1-value)/(double) denominator[j] );
      } else if (type=="pileup"){
    value = ((double) numerator[j])*1./((double) denominator[j]);
        err = sqrt( value*(1+value)/(double) denominator[j] );
      } else return;
      h->setBinContent(j+1, value);
      h->setBinError(j+1, err);
    }
    else {
      h->setBinContent(j+1, 0.);
      h->setBinError(j+1, 0.);
    }
  }
}

void MTVHistoProducerAlgo::fillPlotFromPlots(MonitorElement* h, 
                         TH1* numerator, 
                         TH1* denominator,
                         std::string type){
  assert(h->getNbinsX()==numerator->GetNbinsX());
  assert(h->getNbinsX()==denominator->GetNbinsX());
  double value,err;
  for (int bin=1;bin<=h->getNbinsX();++bin) {
    if (denominator->GetBinContent(bin)!=0) {
      if (type=="effic"){
    value = ((double) numerator->GetBinContent(bin))/((double) denominator->GetBinContent(bin));
    err = sqrt( value*(1.-value)/((double) denominator->GetBinContent(bin)) );
      } else if (type=="fakerate"){
    value = 1.-((double) numerator->GetBinContent(bin))/((double) denominator->GetBinContent(bin));
        err = sqrt( value*(1.-value)/(double)denominator->GetBinContent(bin) );
      } else if (type=="pileup"){
    value = ((double) numerator->GetBinContent(bin))/((double) denominator->GetBinContent(bin));
        err = sqrt( value*(1.+value)/(double) denominator->GetBinContent(bin) );
      } else return;
      h->setBinContent(bin,value);
      h->setBinError(bin,err);
    }
  }
}

void MTVHistoProducerAlgo::BinLogX(TH1*h){  
  TAxis *axis = h->GetXaxis();
  int bins = axis->GetNbins();
  
  float from = axis->GetXmin();
  float to = axis->GetXmax();
  float width = (to - from) / bins;
  float *new_bins = new float[bins + 1];
  
  for (int i = 0; i <= bins; i++) {
    new_bins[i] = TMath::Power(10, from + i * width);
    
  }
  axis->Set(bins, new_bins);
  delete[] new_bins;
}


//void MTVHistoProducerAlgo::
//void MTVHistoProducerAlgo::