HistoCompare.cc

Go to the documentation of this file.

#include "Validation/RecoB/interface/HistoCompare.h"

#include <iostream>

HistoCompare::HistoCompare() {
  result_ = 0;
  do_nothing_ = false;
}

HistoCompare::HistoCompare(const TString &refFilename) {
  result_ = 0;
  do_nothing_ = false;
  // open reference file
  refFile_ = new TFile(refFilename);
  if (refFile_->IsZombie()) {
    std::cout << " Error openning file " << refFilename << std::endl;
    std::cout << " we will not compare histograms. " << std::endl;
    do_nothing_ = true;
  } else {
    std::cout << " open file " << refFilename << std::endl;
  }

  setChi2Test_ = true;
}

HistoCompare::~HistoCompare() {}

TH1 *HistoCompare::Compare(TH1 *h, const TString &hname) {
  if (do_nothing_)
    return nullptr;

  // std::cout << "histo: " << h->GetName() << " entries: "<< h->GetEntries() <<
  // std::endl; std::cout << "hname: " << hname << std::endl;

  // create a residual histogram
  // if (resHisto_) delete resHisto_;
  resHisto_ = (TH1 *)h->Clone(TString(h->GetName()) + "_residuals");
  resHisto_->Reset();

  // clone input histogram
  TH1 *htemp = (TH1 *)h->Clone(TString(h->GetName()));

  refFile_->cd();
  // std::cout << "get reference" << std::endl;
  refHisto_ = (TH1 *)gDirectory->Get(hname);
  // std::cout << "historef: " << refHisto_->GetName() << " entries: "<<
  // refHisto_->GetEntries() << std::endl; std::cout << "name: " <<
  // refHisto_->GetName() << std::endl;

  // normalize histograms
  htemp->Scale(1. / htemp->Integral());
  refHisto_->Scale(1. / refHisto_->Integral());

  if (setChi2Test_) {
    result_ = refHisto_->Chi2Test(htemp, "UFOF");
  }
  if (setKGTest_) {
    result_ = refHisto_->KolmogorovTest(htemp, "UO");
  }

  resHisto_->Add(htemp, refHisto_, 1., -1.);

  // std::cout << " residual done." << std::endl;

  return resHisto_;
}