#include "TFile.h"
#include "TCanvas.h"
#include "TGraphAsymmErrors.h"
#include "ComputeTotalCrossSection.C"
#include <iostream>
#include <sstream>
#include <iomanip>
#include "AtlasUtils.h"
#include "AtlasStyle.C"
#include "AtlasLabels.C"

Double_t show_integrate(std::string gName, TFile* iFile, TFile* jFile)
{
  
  TGraphAsymmErrors* g = (TGraphAsymmErrors*)iFile->Get(gName.c_str());  
  TGraphAsymmErrors* gj = (TGraphAsymmErrors*)jFile->Get(gName.c_str());  
  TCanvas* c = new TCanvas(gName.c_str(),gName.c_str());
  g->Draw("ACP"); 
  gj->SetLineColor(2);
  gj->SetMarkerColor(2);
  gj->Draw("CP");
  std::cout<<"****************************** "<<gName<<" **************************"<<std::endl;
  Double_t xs = ComputeTotalCrossSection(g);
  Double_t xsj = ComputeTotalCrossSection(gj);
  std::cout<<"****** Diff = "<<xs-xsj<<std::endl;
  return xs-xsj;
}

void CompareFiles(std::string aa="pippo")
{
  //  TFile* iFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/rootFilesUCL/cnvHepMC/Zbb_all.v3.Rivet.root");
  //TFile* iFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/rootFilesUCL/cnvHepMC/Zbb_check.v3.Rivet.root");
  //TFile* iFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/events_eej_HERWIGPP_0.Rivet.v3.root");
  //TFile* jFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/rootFilesUCL/cnvHepMC/Zbb_checkv2.v2.Rivet.root");
    //  TFile* iFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSet4fZbb.v3.Rivet.root");
  //  TFile* jFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSet5fZb.v3.Rivet.root");
  //  TFile* kFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSet5fZj.v3.Rivet.root");
  //return;
// // // //     TFile* iFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSet4fZbb_1.v3.Rivet.root");
// // // //     TFile* jFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSet4fZbb_noFR.v3.Rivet.root");
// // // //     TFile* kFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSet4fZbb_MPI.v3.Rivet.root");

//    TFile* iFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSetNewCuts5fZj.v3.Rivet.root");
//    TFile* jFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSet5fZj.v3.Rivet.root");
    TFile* kFile = new TFile("../aMCatNLO/MPI_v3_fromGavin_Aug22.root");
//    TFile* jFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSet5fZj.v3.Rivet.root");

//// to evaluate the change from the scale for the Zbb (WEB vs MCFM)
//    TFile* iFile = new TFile("../aMCatNLO/lecceProd/histo_oldSetNewCuts_MCFMscales_4fZbb_1_1samples.v3.Rivet.root");
//    TFile* jFile = new TFile("../aMCatNLO/lecceProd/histo_oldSetNewCuts_WEBscales_4fZbb_1_1samples.v3.Rivet.root");

//    TFile* iFile = new TFile("../aMCatNLO/lecceProd/histo_oldSetNewCuts_MCFMscales_5fZbMPI_1_6samples.v3.Rivet.root");
//    TFile* jFile = new TFile("../aMCatNLO/lecceProd/histo_oldSetNewCuts_MCFMscales_5fZbNOMPI_1_6samples.v3.Rivet.root");


    // TFile* iFile = new TFile("../aMCatNLO/lecceProd/histo_oldSetNewCuts_MCFMscales_5fZbMPI_1_6samples.v3.Rivet.root");
    // TFile* jFile = new TFile("../aMCatNLO/lecceProd/histo_oldSetNewCuts_MCFMscales_5fZj_1_4samples.v3.Rivet.root");
    
    // TFile* iFile = new TFile("../aMCatNLO/lecceProd/histo_oldSetNewCuts_MCFMscales_5fZbMPI_1_6samples.v3.Rivet.root");
    // TFile* jFile = new TFile("../aMCatNLO/LouvainProd/histo_Louvain_5fZj_MCFMscales_1to5.v3.Rivet.root");

     TFile* iFile = new TFile("../aMCatNLO/Zbb_all.v3.Rivet.root");
     TFile* jFile = new TFile("../aMCatNLO/LouvainProd/histo_Louvain_4fZbb_MCFMscales.v3.Rivet.root");
    
//    TFile* jFile = new TFile("/afs/le.infn.it/project/atlasewk/Generators/aMC_at_NLO/lecceProd/oldsetup13Aug/histo_oldSet5fZj.v3.Rivet.root");
    

  std::string gname; 
  std::string gavinn;
  Double_t xs = 0.; 
  Double_t xs1 = 0.;
  TH1D* gg;

  std::cout<<"First  histo file = "<<iFile->GetName()<<std::endl;
  std::cout<<"Second histo file = "<<jFile->GetName()<<std::endl;
  
  gname = "j_bjet_Pt_binned";
  gavinn = "zb_jet_pT_MPI";
  xs = show_integrate(gname.c_str(), iFile, jFile);
  //xs = show_integrate(gname.c_str(), jFile);
  //xs = show_integrate(gname.c_str(), kFile);
  //xs = ComputeTotalCrossSection(ggg);

  gname = "j_bjet_Y_binned";
  // gavinn = "zb_jet_y_MPI";
  // gg = (TH1D*)kFile->Get(gavinn.c_str()); 
  // gg->Draw();
  // TGraphAsymmErrors* ggg = new TGraphAsymmErrors(gg);
  xs = show_integrate(gname.c_str(), iFile, jFile);

  gname = "j_bjet_Yboost_binned";
  // gavinn = "zb_yboost_MPI";
  // gg = (TH1D*)kFile->Get(gavinn.c_str()); 
  // TGraphAsymmErrors* ggg = new TGraphAsymmErrors(gg);
  xs = show_integrate(gname.c_str(), iFile, jFile);

  gname = "j_bjet_ZPt_binned";
  // gavinn = "zb_Z_pT_MPI";
  // gg = (TH1D*)kFile->Get(gavinn.c_str()); 
  // TGraphAsymmErrors* ggg = new TGraphAsymmErrors(gg);
  xs = show_integrate(gname.c_str(), iFile, jFile);


  gname = "j_bjet_ZY_binned";
  // gavinn = "zb_Z_y_MPI";
  // gg = (TH1D*)kFile->Get(gavinn.c_str()); 
  // ggg = new TGraphAsymmErrors(gg);
  xs = show_integrate(gname.c_str(), iFile, jFile);

  gname = "j_bjet_ZdPhiF20_binned";
  xs = show_integrate(gname.c_str(), iFile, jFile);

  gname = "j_bjet_ZdR620_binned";
  // gavinn = "zb_dR_MPI";
  // gg = (TH1D*)kFile->Get(gavinn.c_str()); 
  // ggg = new TGraphAsymmErrors(gg);
  xs = show_integrate(gname.c_str(), iFile, jFile);

  gname = "j_bjet_DY20_binned";
  // gavinn = "zb_dy_MPI";
  // gg = (TH1D*)kFile->Get(gavinn.c_str()); 
  // ggg = new TGraphAsymmErrors(gg);
  xs = show_integrate(gname.c_str(), iFile, jFile);

  // gname = "j_2bjet_ZPt_binned";
  // xs = show_integrate(gname.c_str(), iFile, jFile);

  gname = "j_2bjet_ZPtv3_binned";
  xs = show_integrate(gname.c_str(), iFile, jFile);

  gname = "j_2bjet_ZY_binned";
  gavinn = "zbb_Z_y_MPI";
  gg = (TH1D*)kFile->Get(gavinn.c_str()); 
  ggg = new TGraphAsymmErrors(gg);
  xs = show_integrate(gname.c_str(), iFile, jFile);

  gname = "j_2bjet_dR_binned";
  gavinn = "zbb_dR_MPI";
  gg = (TH1D*)kFile->Get(gavinn.c_str()); 
  ggg = new TGraphAsymmErrors(gg);
  xs = show_integrate(gname.c_str(), iFile, jFile);

  gname = "j_2bjet_Mbb_binned";
  xs = show_integrate(gname.c_str(), iFile, jFile);



  return;  
}