HistsCompare.cpp

Go to the documentation of this file.

// Compare all histograms in two root files with same structure.
 
#include <iostream>
#include <cstdlib> // For getenv
#include "AutoHistoCompare.h"
#include <TH1F.h>
#include <TCanvas.h>
 
int main(int argc, char** argv) {
  
  const char* reffilename = std::getenv("REFFILE");
  const char* curfilename = std::getenv("CURFILE");
 
  if ( argc > 1 ) {
    reffilename = argv[1];
  }
  if ( argc > 2 ) {
    curfilename = argv[2];
  }
 
  if ( reffilename == nullptr ) {
    std::cerr << "Error: reffilename should provide." << std::endl;
    return 1;
  }
  if ( curfilename == nullptr ) {
    std::cerr << "Error: curfilename should provide." << std::endl;
    return 1;
  }
 
  TH1F *PvHisto = new TH1F("PV", "PV Distribution", 200, -0.1, 1.1);
 
  AutoHistoCompare *myPV = new AutoHistoCompare(PvHisto, reffilename, curfilename, 2, 0);
 
  TCanvas* canvas = new TCanvas("canvas", "Histogram Canvas", 800, 600);
  myPV->DoCompare();
  PvHisto->Draw();
  canvas->SaveAs("PVdistribution.eps");
 
  delete PvHisto;
  delete myPV;
  delete canvas;
 
  return 0;
}