/*
	Compare root files
*/
{
  gROOT->Reset();

  gROOT->ProcessLine(".L CompareAdvanced.C");

  TCanvas *c1 = new TCanvas("c1236","Muon",200,10,800,700);

  const string Cut = "(abs(jgenEta[0])<5.0)&&(abs(jgenEta[1])<5.0)&&(jgenPt[0]>25)&&(jgenPt[1]>25)&&(abs(jgenEta[0]-jgenEta[1])>4)";

//  double deltaPhi(double phi1, double phi2)
//  {
//          double deltaphi = fabs(phi1-phi2);
//          if (deltaphi > TMath::PI ) deltaphi = 2 * TMath::PI - deltaphi;
//          if (deltaphi > 2 * TMath::PI ) cout<< "Delta Phi = "<<deltaphi<<endl;
//          return deltaphi;
//  };

c1->Divide(1,2);
c1->cd(1);
compareQuantities("genMuonPhi","",10,-4,4,"",6,"pp_WWJJ_phantom.root","phantom","pp_WWJJ_notop_LL_8TeV_GENSIM.root","LL","pp_WWJJ_notop_LT_8TeV_GENSIM.root","LT","pp_WWJJ_notop_TT_8TeV_GENSIM.root","TT","pp_WWJJ_notop_TOT_8TeV_GENSIM.root","TOT","pp_WWJJ_phantom_qcd.root","phan_bkg");
c1->cd(2);
compareQuantities("genMuonPhi","",10,-4,4,Cut,6,"pp_WWJJ_phantom.root","phantom","pp_WWJJ_notop_LL_8TeV_GENSIM.root","LL","pp_WWJJ_notop_LT_8TeV_GENSIM.root","LT","pp_WWJJ_notop_TT_8TeV_GENSIM.root","TT","pp_WWJJ_notop_TOT_8TeV_GENSIM.root","TOT","pp_WWJJ_phantom_qcd.root","phan_bkg");
c1->SaveAs("Phi.pdf");
c1->Clear();
}