#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)
{
  
  TGraphAsymmErrors* g = (TGraphAsymmErrors*)iFile->Get(gName.c_str());  
  TCanvas* c = new TCanvas(gName.c_str(),gName.c_str());
  g->Draw("ACP");  
  std::cout<<"****************************** "<<gName<<" **************************"<<std::endl;
  Double_t xs = ComputeTotalCrossSection(g);
  return xs;
}

void getCrossSections(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* 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;
  std::cout<<"Third  histo file = "<<kFile->GetName()<<std::endl;
  
  gname = "j_bjet_Pt_binned";
  gavinn = "zb_jet_pT_MPI";
  xs = show_integrate(gname.c_str(), iFile);
  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);
  xs = show_integrate(gname.c_str(), jFile);
  //xs = show_integrate(gname.c_str(), kFile);
  xs = ComputeTotalCrossSection(ggg);

  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);
  xs1 = show_integrate(gname.c_str(), jFile);
  std::cout<<"MPI on top of Zb ME = "<<xs-xs1<<std::endl; 
  //xs = show_integrate(gname.c_str(), kFile);
  xs = ComputeTotalCrossSection(ggg);

  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);
  xs1 = show_integrate(gname.c_str(), jFile);
  std::cout<<"MPI on top of Zb ME = "<<xs-xs1<<std::endl; 
  //xs = show_integrate(gname.c_str(), kFile);
  xs = ComputeTotalCrossSection(ggg);

  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);
  xs = show_integrate(gname.c_str(), jFile);
  //xs = show_integrate(gname.c_str(), kFile);
  xs = ComputeTotalCrossSection(ggg);

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

  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);
  xs1 = show_integrate(gname.c_str(), jFile);
  std::cout<<"MPI on top of Zb ME = "<<xs-xs1<<std::endl; 
  //xs = show_integrate(gname.c_str(), kFile);
  xs = ComputeTotalCrossSection(ggg);

  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);
  xs1 = show_integrate(gname.c_str(), jFile);
  std::cout<<"MPI on top of Zb ME = "<<xs-xs1<<std::endl; 
  //xs = show_integrate(gname.c_str(), kFile);
  xs = ComputeTotalCrossSection(ggg);

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

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

  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);
  xs = show_integrate(gname.c_str(), jFile);
  //xs = show_integrate(gname.c_str(), kFile);
  xs = ComputeTotalCrossSection(ggg);

  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);
  xs1 = show_integrate(gname.c_str(), jFile);
  std::cout<<"MPI on top of Zb ME = "<<xs-xs1<<std::endl; 
  //xs = show_integrate(gname.c_str(), kFile);
  xs = ComputeTotalCrossSection(ggg);

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



  return;  
}