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

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

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

const string Cut = "(abs(j0genEta)<5.0)&&(abs(j1genEta)<5.0)&&(j0genPt>25)&&(j1genPt>25)&&(abs(j0genEta-j1genEta)>4)&&(cos(j0genTheta-j1genTheta)<0)&&(((j0genEta+0.7<genMuonEta)&&(j1genEta-0.7>genMuonEta))||((j0genEta-0.7>genMuonEta)&&(j1genEta+0.7<genMuonEta)))&&((genMuonPt>20)||(genElePt>20))&&(mgenPt>20)&&((abs(genMuonEta)<2)||(abs(genEleEta)<2))&&(sqrt((genMuonEta-j0genEta)**2+(genMuonPhi-j0genPhi)**2)>0.7)&&(sqrt((genMuonEta-j1genEta)**2+(genMuonPhi-j1genPhi)**2)>0.7)&&(sqrt((genMuonEta-j2genEta)**2+(genMuonPhi-j2genPhi)**2)>0.7)&&(sqrt((genMuonEta-j3genEta)**2+(genMuonPhi-j3genPhi)**2)>0.7)&&(bQuark==-10)&&(tQuark==-10)";

///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"ParticleIntermediate","",56,-28,28,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"ParticleIntermediate","",56,-28,28,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("ParticleIntermediate.pdf");
	c1->Clear();	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"ParticleIn","",14,-7,7,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"ParticleIn","",14,-7,7,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("ParticleIn.pdf");
	c1->Clear();	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"mWW","mWW",20,100,800.0,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"mWW","mWW",20,100,800.0,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("Masses.pdf");
	c1->Clear();
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"mWLep","",20,50,110,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"mWLep","",20,50,110,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("Masses1.pdf");
	c1->Clear();
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"mWHad","",20,50,210,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"mWHad","",20,50,210,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("Masses2.pdf");
	c1->Clear();
//	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"costheta1","",20,-1,1,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"costheta1","",20,-1,1,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("costheta1.pdf");
	c1->Clear();
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"costheta2","",20,-1,1,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"costheta2","",20,-1,1,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("costheta2.pdf");
	c1->Clear();	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"costhetastar","",20,-1,1,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"costhetastar","",20,-1,1,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("costhetastar.pdf");
	c1->Clear();	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"phi","",20,-3.5,3.5,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"phi","",20,-3.5,3.5,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("phi.pdf");
	c1->Clear();	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"phi1","",20,-3.5,3.5,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"phi1","",20,-3.5,3.5,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("phi1.pdf");
	c1->Clear();	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"dEtajj","",20,0,9,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"dEtajj","",20,0,9,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("dEtajj.pdf");
	c1->Clear();	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"dPhijj","",20,0,3.5,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"dPhijj","",20,0,3.5,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("dPhijj.pdf");
	c1->Clear();
	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"j0genEta-((j0genEta+j1genEta)/2)","",15,-3.5,3.5,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"j0genEta-((j0genEta+j1genEta)/2)","",15,-3.5,3.5,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("zeppenfeld1.pdf");
	c1->Clear();	
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"j1genEta-((j0genEta+j1genEta)/2)","",15,-3.5,3.5,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"j1genEta-((j0genEta+j1genEta)/2)","",15,-3.5,3.5,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("zeppenfeld2.pdf");
	c1->Clear();
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"j2genEta-((j0genEta+j1genEta)/2)","",15,-3.5,3.5,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"j2genEta-((j0genEta+j1genEta)/2)","",15,-3.5,3.5,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("zeppenfeld3.pdf");
	c1->Clear();
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////	
	c1->Divide(2,1);
	c1->cd(1);
	compareQuantities(1,1,"j3genEta-((j0genEta+j1genEta)/2)","",15,-3.5,3.5,"",6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->cd(2);
	compareQuantities(0,1,"j3genEta-((j0genEta+j1genEta)/2)","",15,-3.5,3.5,Cut,6,"lhe_LL.root","LL","lhe_LT.root","LT","lhe_TT.root","TT","lhe_TOT.root","TOT","lhe_Phantom.root","phantom","lhe_Phantom_qcd.root","phantom_qcd");	
	c1->SaveAs("zeppenfeld4.pdf");
	c1->Clear();

	/*
	*/
}