The Z Analysis ROOTbookMany analyses selecting leptons suffer from Z + jets as a contributing background due to its large production cross section. It is therefore vital to check the correct modelling of this process by the Monte-Carlo simulated data. It is important to measure well known Standard Model particles, to confirm that we understand properly the detector and software. We are then ready to search for new physics. |
%jsroot on
TChain *InputSignal_Zee = new TChain("mini");
TChain *InputSignal_Zmumu = new TChain("mini");
TChain *InputSignal_Ztautau = new TChain("mini");
TChain *MonteCarlo = new TChain("mini");
TChain *Data = new TChain("mini");
// Data
Data->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/Data/DataEgamma.root");
Data->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/Data/DataMuons.root");
//Z
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_147770.Zee.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_147771.Zmumu.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_147772.Ztautau.root");
//Diboson
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_105985.WW.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_105986.ZZ.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_105987.WZ.root");
//Drell Yan
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_173041.DYeeM08to15.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_173042.DYeeM15to40.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_173043.DYmumuM08to15.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_173044.DYmumuM15to40.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_173045.DYtautauM08to15.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_173046.DYtautauM15to40.root");
//Stop
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_110090.stop_tchan_top.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_110091.stop_tchan_antitop.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_110119.stop_schan.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_110140.stop_wtchan.root");
//TTBar
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_117049.ttbar_had.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/MC/mc_117050.ttbar_lep.root");
//Data
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/Data/DataEgamma.root");
MonteCarlo->Add("http://atlas-opendata.web.cern.ch/atlas-opendata/release/samples/Data/DataMuons.root");
Float_t Luminosity = 1000;
//MonteCarlo
Int_t MonteCarlo_lep_n = -1,
MonteCarlo_lep_type[5],
MonteCarlo_runNumber=-1;
Float_t MonteCarlo_lep_pt[5],
MonteCarlo_lep_charge[5],
MonteCarlo_lep_eta[5],
MonteCarlo_lep_phi[5],
MonteCarlo_lep_E[6],
MonteCarlo_scaleFactor_PILEUP=-999,
MonteCarlo_scaleFactor_ELE=-999,
MonteCarlo_scaleFactor_MUON=-999,
MonteCarlo_scaleFactor_BTAG=-999,
MonteCarlo_scaleFactor_TRIGGER=-999,
MonteCarlo_scaleFactor_JVFSF=-999,
MonteCarlo_scaleFactor_ZVERTEX=-999,
MonteCarlo_mcWeight=-999;
MonteCarlo->SetBranchAddress("lep_pt", &MonteCarlo_lep_pt);
MonteCarlo->SetBranchAddress("lep_charge", &MonteCarlo_lep_charge);
MonteCarlo->SetBranchAddress("lep_n", &MonteCarlo_lep_n);
MonteCarlo->SetBranchAddress("lep_type", &MonteCarlo_lep_type);
MonteCarlo->SetBranchAddress("lep_eta", &MonteCarlo_lep_eta);
MonteCarlo->SetBranchAddress("lep_phi", &MonteCarlo_lep_phi);
MonteCarlo->SetBranchAddress("lep_E", &MonteCarlo_lep_E);
MonteCarlo->SetBranchAddress("channelNumber", &MonteCarlo_runNumber);
MonteCarlo->SetBranchAddress("scaleFactor_PILEUP", &MonteCarlo_scaleFactor_PILEUP);
MonteCarlo->SetBranchAddress("scaleFactor_ELE", &MonteCarlo_scaleFactor_ELE);
MonteCarlo->SetBranchAddress("scaleFactor_MUON", &MonteCarlo_scaleFactor_MUON);
MonteCarlo->SetBranchAddress("scaleFactor_BTAG", &MonteCarlo_scaleFactor_BTAG);
MonteCarlo->SetBranchAddress("scaleFactor_TRIGGER", &MonteCarlo_scaleFactor_TRIGGER);
MonteCarlo->SetBranchAddress("scaleFactor_JVFSF", &MonteCarlo_scaleFactor_JVFSF);
MonteCarlo->SetBranchAddress("scaleFactor_JVFSF", &MonteCarlo_scaleFactor_ZVERTEX);
MonteCarlo->SetBranchAddress("mcWeight", &MonteCarlo_mcWeight);
//Data
Int_t Data_lep_n = -1,
Data_lep_type[5],
Data_runNumber=-1;
Float_t Data_lep_pt[5],
Data_lep_charge[5],
Data_lep_eta[5],
Data_lep_phi[5],
Data_lep_E[6];
Data->SetBranchAddress("lep_pt", &Data_lep_pt);
Data->SetBranchAddress("lep_charge", &Data_lep_charge);
Data->SetBranchAddress("lep_n", &Data_lep_n);
Data->SetBranchAddress("lep_type", &Data_lep_type);
Data->SetBranchAddress("lep_eta", &Data_lep_eta);
Data->SetBranchAddress("lep_phi", &Data_lep_phi);
Data->SetBranchAddress("lep_E", &Data_lep_E);
Data->SetBranchAddress("channelNumber", &Data_runNumber);
//PT HISTOGRAMS
//Zee histrogram - last cut only
TH1F *h_lep_pt_Zee = new TH1F("h_lep_pt_Zee","h_lep_pt_Zee",20,0,200);
h_lep_pt_Zee->SetFillColor(kRed);
//Zmumu histrogram - last cut only
TH1F *h_lep_pt_Zmumu = new TH1F("h_lep_pt_Zmumu","h_lep_pt_Zmumu",20,0,200);
h_lep_pt_Zmumu->SetFillColor(kYellow);
//Ztautau histrogram - last cut only
TH1F *h_lep_pt_Ztautau = new TH1F("h_lep_pt_Ztautau","h_lep_pt_Ztautau",20,0,200);
h_lep_pt_Ztautau->SetFillColor(kBlue);
//Merged Zee,Zmumu,Ztautau
TH1F *h_lep_pt_mergedZ = new TH1F("h_lep_pt_mergedZ","h_lep_pt_mergedZ",20,0,200);
h_lep_pt_mergedZ->SetFillColor(kGreen);
//Diboson WW histrogram
TH1F *h_lep_pt_Dww = new TH1F("h_lep_pt_Dww","h_lep_pt_Dww",20,0,200);
h_lep_pt_Dww->SetFillColor(kRed);
//Diboson ZZ histrogram
TH1F *h_lep_pt_Dzz = new TH1F("h_lep_pt_Dzz","h_lep_pt_Dzz",20,0,200);
h_lep_pt_Dzz->SetFillColor(kYellow);
//Diboson WZ histrogram
TH1F *h_lep_pt_Dwz = new TH1F("h_lep_pt_Dwz","h_lep_pt_Dwz",20,0,200);
h_lep_pt_Dwz->SetFillColor(kBlue);
//Merged Diboson WW,ZZ,WZ histogram
TH1F *h_lep_pt_mergedD = new TH1F("h_lep_pt_mergedD","h_lep_pt_mergedD",20,0,200);
h_lep_pt_mergedD->SetFillColor(kGreen);
//DYee1 histrogram
TH1F *h_lep_pt_DYee1 = new TH1F("h_lep_pt_DYee1","h_lep_pt_DYee1",20,0,200);
h_lep_pt_DYee1->SetFillColor(kRed);
//DYee2 histrogram
TH1F *h_lep_pt_DYee2 = new TH1F("h_lep_pt_DYee2","h_lep_pt_DYee2",20,0,200);
h_lep_pt_DYee2->SetFillColor(kYellow);
//DYmumu1 histrogram
TH1F *h_lep_pt_DYmumu1 = new TH1F("h_lep_pt_DYmumu1 ","h_lep_pt_DYmumu1",20,0,200);
h_lep_pt_DYmumu1->SetFillColor(kBlue);
//DYmumu2 histrogram
TH1F *h_lep_pt_DYmumu2 = new TH1F("h_lep_pt_DYmumu2","h_lep_pt_DYmumu2",20,0,200);
h_lep_pt_DYmumu2->SetFillColor(kGreen);
//DYtautau1 histrogram
TH1F *h_lep_pt_DYtautau1 = new TH1F("h_lep_pt_DYtautau1","h_lep_pt_DYtautau1",20,0,200);
h_lep_pt_DYtautau1->SetFillColor(kOrange);
//DYtautau2 histrogram
TH1F *h_lep_pt_DYtautau2 = new TH1F("h_lep_pt_DYtautau2","h_lep_pt_DYtautau2",20,0,200);
h_lep_pt_DYtautau2->SetFillColor(kPink);
//Merged DY EE,MUMU,TAUTAU histogram
TH1F *h_lep_pt_mergedDY = new TH1F("h_lep_pt_mergedDY","h_lep_pt_mergedDY",20,0,200);
h_lep_pt_mergedDY->SetFillColor(kGray);
//Stop tchan_top histrogram
TH1F *h_lep_pt_Stchan_top = new TH1F("h_lep_pt_Stchan_top","h_lep_pt_Stchan_top",20,0,200);
h_lep_pt_Stchan_top->SetFillColor(kRed);
//Stop tchan_antitop histrogram
TH1F *h_lep_pt_Stchan_antitop = new TH1F("h_lep_pt_Stchan_antitop","h_lep_pt_Stchan_antitop",20,0,200);
h_lep_pt_Stchan_antitop->SetFillColor(kYellow);
//Stop schan histogram
TH1F *h_lep_pt_Sschan = new TH1F("h_lep_pt_Sschan ","h_lep_pt_Sschan",20,0,200);
h_lep_pt_Sschan->SetFillColor(kBlue);
//Stop wtchan histogram
TH1F *h_lep_pt_Swtchan = new TH1F("h_lep_pt_Swtchan","h_lep_pt_Swtchan",20,0,200);
h_lep_pt_Swtchan->SetFillColor(kGreen);
//Merged stop
TH1F *h_lep_pt_mergedS = new TH1F("h_lep_pt_mergedS","h_lep_pt_mergedS",20,0,200);
h_lep_pt_mergedS->SetFillColor(kGray);
//TTBar had histrogram
TH1F *h_lep_pt_TThad = new TH1F("h_lep_pt_TThad","h_lep_pt_TThad",20,0,200);
h_lep_pt_TThad->SetFillColor(kRed);
//TTBar lep histrogram
TH1F *h_lep_pt_TTlep = new TH1F("h_lep_pt_TTlep","h_lep_pt_TTlep",20,0,200);
h_lep_pt_TTlep->SetFillColor(kYellow);
//Merged TTBar
TH1F *h_lep_pt_mergedTT = new TH1F("h_lep_pt_mergedTT","h_lep_pt_mergedTT",20,0,200);
h_lep_pt_mergedTT->SetFillColor(kGray);
//ETA HISTOGRAMS
//Zee histrogram - last cut only
TH1F *h_lep_eta_Zee = new TH1F("h_lep_eta_Zee","h_lep_eta_Zee",33,-3,3);
h_lep_eta_Zee->SetFillColor(kRed);
//Zmumu histrogram - last cut only
TH1F *h_lep_eta_Zmumu = new TH1F("h_lep_eta_Zmumu","h_lep_eta_Zmumu",33,-3,3);
h_lep_eta_Zmumu->SetFillColor(kYellow);
//Ztautau histrogram - last cut only
TH1F *h_lep_eta_Ztautau = new TH1F("h_lep_eta_Ztautau","h_lep_eta_Ztautau",33,-3,3);
h_lep_eta_Ztautau->SetFillColor(kBlue);
//Merged Zee,Zmumu,Ztautau
TH1F *h_lep_eta_mergedZ = new TH1F("h_lep_eta_mergedZ","h_lep_eta_mergedZ",33,-3,3);
h_lep_eta_mergedZ->SetFillColor(kGreen);
//Diboson WW histrogram
TH1F *h_lep_eta_Dww = new TH1F("h_lep_eta_Dww","h_lep_eta_Dww",33,-3,3);
h_lep_eta_Dww->SetFillColor(kRed);
//Diboson ZZ histrogram
TH1F *h_lep_eta_Dzz = new TH1F("h_lep_eta_Dzz","h_lep_eta_Dzz",33,-3,3);
h_lep_eta_Dzz->SetFillColor(kYellow);
//Diboson WZ histrogram
TH1F *h_lep_eta_Dwz = new TH1F("h_lep_eta_Dwz","h_lep_eta_Dwz",33,-3,3);
h_lep_eta_Dwz->SetFillColor(kBlue);
//Merged Dibson WW,ZZ,WZ histogram
TH1F *h_lep_eta_mergedD = new TH1F("h_lep_eta_mergedD","h_lep_eta_mergedD",33,-3,3);
h_lep_eta_mergedD->SetFillColor(kGreen);
//DYee1 histrogram
TH1F *h_lep_eta_DYee1 = new TH1F("h_lep_eta_DYee1","h_lep_eta_DYee1",33,-3,3);
h_lep_eta_DYee1->SetFillColor(kRed);
//DYee2 histrogram
TH1F *h_lep_eta_DYee2 = new TH1F("h_lep_eta_DYee2","h_lep_eta_DYee2",33,-3,3);
h_lep_eta_DYee2->SetFillColor(kYellow);
//DYmumu1 histrogram
TH1F *h_lep_eta_DYmumu1 = new TH1F("h_lep_eta_DYmumu1","h_lep_eta_DYmumu1",33,-3,3);
h_lep_eta_DYmumu1->SetFillColor(kBlue);
//DYmumu2 histrogram
TH1F *h_lep_eta_DYmumu2 = new TH1F("h_lep_eta_DYmumu2","h_lep_eta_DYmumu2",33,-3,3);
h_lep_eta_DYmumu2->SetFillColor(kGreen);
//DYtautau1 histrogram
TH1F *h_lep_eta_DYtautau1 = new TH1F("h_lep_eta_DYtautau1","h_lep_eta_DYtautau1",33,-3,3);
h_lep_eta_DYtautau1->SetFillColor(kOrange);
//DYtautau2 histrogram
TH1F *h_lep_eta_DYtautau2 = new TH1F("h_lep_eta_DYtautau2","h_lep_eta_DYtautau2",33,-3,3);
h_lep_eta_DYtautau2->SetFillColor(kPink);
//Merged DY EE,MUMU,TAUTAU histogram
TH1F *h_lep_eta_mergedDY = new TH1F("h_lep_eta_mergedDY","h_lep_eta_mergedDY",33,-3,3);
h_lep_eta_mergedDY->SetFillColor(kGray);
//Stop tchan_top histrogram
TH1F *h_lep_eta_Stchan_top = new TH1F("h_lep_eta_Stchan_top","h_lep_eta_Stchan_top",33,-3,3);
h_lep_eta_Stchan_top->SetFillColor(kRed);
//Stop tchan_antitop histrogram
TH1F *h_lep_eta_Stchan_antitop = new TH1F("h_lep_eta_Stchan_antitop","h_lep_eta_Stchan_antitop",33,-3,3);
h_lep_eta_Stchan_antitop->SetFillColor(kYellow);
//Stop schan histogram
TH1F *h_lep_eta_Sschan = new TH1F("h_lep_eta_Sschan ","h_lep_eta_Sschan",33,-3,3);
h_lep_eta_Sschan->SetFillColor(kBlue);
//Stop wtchan histogram
TH1F *h_lep_eta_Swtchan = new TH1F("h_lep_eta_Swtchan","h_lep_eta_Swtchan",33,-3,3);
h_lep_eta_Swtchan->SetFillColor(kGreen);
//Merged stop
TH1F *h_lep_eta_mergedS = new TH1F("h_lep_eta_mergedS","h_lep_eta_mergedS",33,-3,3);
h_lep_eta_mergedS->SetFillColor(kGray);
//TTBar had histrogram
TH1F *h_lep_eta_TThad = new TH1F("h_lep_eta_TThad","h_lep_eta_TThad",33,-3,3);
h_lep_eta_TThad->SetFillColor(kRed);
//TTBar lep histrogram
TH1F *h_lep_eta_TTlep = new TH1F("h_lep_eta_TTlep","h_lep_eta_TTlep",33,-3,3);
h_lep_eta_TTlep->SetFillColor(kYellow);
//Merged TTBar
TH1F *h_lep_eta_mergedTT = new TH1F("h_lep_eta_mergedTT","h_lep_eta_mergedTT",33,-3,3);
h_lep_eta_mergedTT->SetFillColor(kGray);
//DATA HISTOGRAMS
//Egamma histrogram
TH1F *h_lep_pt_Data = new TH1F("h_lep_pt_Data","h_lep_pt_Data",20,0,200);
h_lep_pt_Data->SetFillColor(kRed);
//Muons histrogram
TH1F *h_lep_pt_Muons = new TH1F("h_lep_pt_Muons","h_lep_pt_Muons",20,0,200);
h_lep_pt_Muons->SetFillColor(kYellow);
//PHI HISTOGRAMS
//Zee histrogram - last cut only
TH1F *h_lep_phi_Zee = new TH1F("h_lep_phi_Zee","h_lep_phi_Zee",40,-4,4);
h_lep_phi_Zee->SetFillColor(kRed);
//Zmumu histrogram - last cut only
TH1F *h_lep_phi_Zmumu = new TH1F("h_lep_phi_Zmumu","h_lep_phi_Zmumu",40,-4,4);
h_lep_phi_Zmumu->SetFillColor(kYellow);
//Ztautau histrogram - last cut only
TH1F *h_lep_phi_Ztautau = new TH1F("h_lep_phi_Ztautau","h_lep_phi_Ztautau",40,-4,4);
h_lep_phi_Ztautau->SetFillColor(kBlue);
//Merged Zee,Zmumu,Ztautau
TH1F *h_lep_phi_mergedZ = new TH1F("h_lep_phi_mergedZ","h_lep_phi_mergedZ",40,-4,4);
h_lep_phi_mergedZ->SetFillColor(kGreen);
//Diboson WW histrogram
TH1F *h_lep_phi_Dww = new TH1F("h_lep_phi_Dww","h_lep_phi_Dww",40,-4,4);
h_lep_phi_Dww->SetFillColor(kRed);
//Diboson ZZ histrogram
TH1F *h_lep_phi_Dzz = new TH1F("h_lep_phi_Dzz","h_lep_phi_Dzz",40,-4,4);
h_lep_phi_Dzz->SetFillColor(kYellow);
//Diboson WZ histrogram
TH1F *h_lep_phi_Dwz = new TH1F("h_lep_phi_Dwz","h_lep_phi_Dwz",40,-4,4);
h_lep_phi_Dwz->SetFillColor(kBlue);
//Merged Diboson WW,ZZ,WZ histogram
TH1F *h_lep_phi_mergedD = new TH1F("h_lep_phi_mergedD","h_lep_phi_mergedD",40,-4,4);
h_lep_phi_mergedD->SetFillColor(kGreen);
//DYee1 histrogram
TH1F *h_lep_phi_DYee1 = new TH1F("h_lep_phi_DYee1","h_lep_phi_DYee1",40,-4,4);
h_lep_phi_DYee1->SetFillColor(kRed);
//DYee2 histrogram
TH1F *h_lep_phi_DYee2 = new TH1F("h_lep_phi_DYee2","h_lep_phi_DYee2",40,-4,4);
h_lep_phi_DYee2->SetFillColor(kYellow);
//DYmumu1 histrogram
TH1F *h_lep_phi_DYmumu1 = new TH1F("h_lep_phi_DYmumu1 ","h_lep_phi_DYmumu1",40,-4,4);
h_lep_phi_DYmumu1->SetFillColor(kBlue);
//DYmumu2 histrogram
TH1F *h_lep_phi_DYmumu2 = new TH1F("h_lep_phi_DYmumu2","h_lep_phi_DYmumu2",40,-4,4);
h_lep_phi_DYmumu2->SetFillColor(kGreen);
//DYtautau1 histrogram
TH1F *h_lep_phi_DYtautau1 = new TH1F("h_lep_phi_DYtautau1","h_lep_phi_DYtautau1",40,-4,4);
h_lep_phi_DYtautau1->SetFillColor(kOrange);
//DYtautau2 histrogram
TH1F *h_lep_phi_DYtautau2 = new TH1F("h_lep_phi_DYtautau2","h_lep_phi_DYtautau2",40,-4,4);
h_lep_phi_DYtautau2->SetFillColor(kPink);
//Merged DY EE,MUMU,TAUTAU histogram
TH1F *h_lep_phi_mergedDY = new TH1F("h_lep_phi_mergedDY","h_lep_phi_mergedDY",40,-4,4);
h_lep_phi_mergedDY->SetFillColor(kGray);
//Stop tchan_top histrogram
TH1F *h_lep_phi_Stchan_top = new TH1F("h_lep_phi_Stchan_top","h_lep_phi_Stchan_top",40,-4,4);
h_lep_phi_Stchan_top->SetFillColor(kRed);
//Stop tchan_antitop histrogram
TH1F *h_lep_phi_Stchan_antitop = new TH1F("h_lep_phi_Stchan_antitop","h_lep_phi_Stchan_antitop",40,-4,4);
h_lep_phi_Stchan_antitop->SetFillColor(kYellow);
//Stop schan histogram
TH1F *h_lep_phi_Sschan = new TH1F("h_lep_phi_Sschan ","h_lep_phi_Sschan",40,-4,4);
h_lep_phi_Sschan->SetFillColor(kBlue);
//Stop wtchan histogram
TH1F *h_lep_phi_Swtchan = new TH1F("h_lep_phi_Swtchan","h_lep_phi_Swtchan",40,-4,4);
h_lep_phi_Swtchan->SetFillColor(kGreen);
//Merged stop
TH1F *h_lep_phi_mergedS = new TH1F("h_lep_phi_mergedS","h_lep_phi_mergedS",40,-4,4);
h_lep_phi_mergedS->SetFillColor(kGray);
//TTBar had histrogram
TH1F *h_lep_phi_TThad = new TH1F("h_lep_phi_TThad","h_lep_phi_TThad",40,-4,4);
h_lep_phi_TThad->SetFillColor(kRed);
//TTBar lep histrogram
TH1F *h_lep_phi_TTlep = new TH1F("h_lep_phi_TTlep","h_lep_phi_TTlep",40,-4,4);
h_lep_phi_TTlep->SetFillColor(kYellow);
//Merged TTBar
TH1F *h_lep_phi_mergedTT = new TH1F("h_lep_phi_mergedTT","h_lep_phi_mergedTT",40,-4,4);
h_lep_phi_mergedTT->SetFillColor(kGray);
//E HISTOGRAMS
//Zee histrogram - last cut only
TH1F *h_lep_E_Zee = new TH1F("h_lep_E_Zee","h_lep_E_Zee",30,0,1000);
h_lep_E_Zee->SetFillColor(kRed);
//Zmumu histrogram - last cut only
TH1F *h_lep_E_Zmumu = new TH1F("h_lep_E_Zmumu","h_lep_E_Zmumu",30,0,1000);
h_lep_E_Zmumu->SetFillColor(kYellow);
//Ztautau histrogram - last cut only
TH1F *h_lep_E_Ztautau = new TH1F("h_lep_E_Ztautau","h_lep_E_Ztautau",30,0,1000);
h_lep_E_Ztautau->SetFillColor(kBlue);
//Merged Zee,Zmumu,Ztautau
TH1F *h_lep_E_mergedZ = new TH1F("h_lep_E_mergedZ","h_lep_E_mergedZ",30,0,1000);
h_lep_E_mergedZ->SetFillColor(kGreen);
//Diboson WW histrogram
TH1F *h_lep_E_Dww = new TH1F("h_lep_E_Dww","h_lep_E_Dww",30,0,1000);
h_lep_E_Dww->SetFillColor(kRed);
//Diboson ZZ histrogram
TH1F *h_lep_E_Dzz = new TH1F("h_lep_E_Dzz","h_lep_E_Dzz",30,0,1000);
h_lep_E_Dzz->SetFillColor(kYellow);
//Diboson WZ histrogram
TH1F *h_lep_E_Dwz = new TH1F("h_lep_E_Dwz","h_lep_E_Dwz",30,0,1000);
h_lep_E_Dwz->SetFillColor(kBlue);
//Merged Diboson WW,ZZ,WZ histogram
TH1F *h_lep_E_mergedD = new TH1F("h_lep_E_mergedD","h_lep_E_mergedD",30,0,1000);
h_lep_E_mergedD->SetFillColor(kGreen);
//DYee1 histrogram
TH1F *h_lep_E_DYee1 = new TH1F("h_lep_E_DYee1","h_lep_E_DYee1",30,0,1000);
h_lep_E_DYee1->SetFillColor(kRed);
//DYee2 histrogram
TH1F *h_lep_E_DYee2 = new TH1F("h_lep_E_DYee2","h_lep_E_DYee2",30,0,1000);
h_lep_E_DYee2->SetFillColor(kYellow);
//DYmumu1 histrogram
TH1F *h_lep_E_DYmumu1 = new TH1F("h_lep_E_DYmumu1 ","h_lep_E_DYmumu1",30,0,1000);
h_lep_E_DYmumu1->SetFillColor(kBlue);
//DYmumu2 histrogram
TH1F *h_lep_E_DYmumu2 = new TH1F("h_lep_E_DYmumu2","h_lep_E_DYmumu2",30,0,1000);
h_lep_E_DYmumu2->SetFillColor(kGreen);
//DYtautau1 histrogram
TH1F *h_lep_E_DYtautau1 = new TH1F("h_lep_E_DYtautau1","h_lep_E_DYtautau1",30,0,1000);
h_lep_E_DYtautau1->SetFillColor(kOrange);
//DYtautau2 histrogram
TH1F *h_lep_E_DYtautau2 = new TH1F("h_lep_E_DYtautau2","h_lep_E_DYtautau2",30,0,1000);
h_lep_E_DYtautau2->SetFillColor(kPink);
//Merged DY EE,MUMU,TAUTAU histogram
TH1F *h_lep_E_mergedDY = new TH1F("h_lep_E_mergedDY","h_lep_E_mergedDY",30,0,1000);
h_lep_E_mergedDY->SetFillColor(kGray);
//Stop tchan_top histrogram
TH1F *h_lep_E_Stchan_top = new TH1F("h_lep_E_Stchan_top","h_lep_E_Stchan_top",30,0,1000);
h_lep_E_Stchan_top->SetFillColor(kRed);
//Stop tchan_antitop histrogram
TH1F *h_lep_E_Stchan_antitop = new TH1F("h_lep_E_Stchan_antitop","h_lep_E_Stchan_antitop",30,0,1000);
h_lep_E_Stchan_antitop->SetFillColor(kYellow);
//Stop schan histogram
TH1F *h_lep_E_Sschan = new TH1F("h_lep_E_Sschan ","h_lep_E_Sschan",30,0,1000);
h_lep_E_Sschan->SetFillColor(kBlue);
//Stop wtchan histogram
TH1F *h_lep_E_Swtchan = new TH1F("h_lep_E_Swtchan","h_lep_E_Swtchan",30,0,1000);
h_lep_E_Swtchan->SetFillColor(kGreen);
//Merged stop
TH1F *h_lep_E_mergedS = new TH1F("h_lep_E_mergedS","h_lep_E_mergedS",30,0,1000);
h_lep_E_mergedS->SetFillColor(kGray);
//TTBar had histrogram
TH1F *h_lep_E_TThad = new TH1F("h_lep_E_TThad","h_lep_E_TThad",30,0,1000);
h_lep_E_TThad->SetFillColor(kRed);
//TTBar lep histrogram
TH1F *h_lep_E_TTlep = new TH1F("h_lep_E_TTlep","h_lep_E_TTlep",30,0,1000);
h_lep_E_TTlep->SetFillColor(kYellow);
//Merged TTBar
TH1F *h_lep_E_mergedTT = new TH1F("h_lep_E_mergedTT","h_lep_E_mergedTT",30,0,1000);
h_lep_E_mergedTT->SetFillColor(kGray);
//PT Histograms
//Zee histrogram
h_lep_pt_Zee->SetMarkerStyle(8);
h_lep_pt_Zee->SetMarkerSize(0.5);
h_lep_pt_Zee->SetLineColor(kRed);
//Zmumu histrogram
h_lep_pt_Zmumu->SetMarkerStyle(8);
h_lep_pt_Zmumu->SetMarkerSize(0.5);
h_lep_pt_Zmumu->SetLineColor(kYellow);
//Ztautau histrogram
h_lep_pt_Ztautau->SetMarkerStyle(8);
h_lep_pt_Ztautau->SetMarkerSize(0.5);
h_lep_pt_Ztautau->SetLineColor(kBlue);
//Merged Zee,Zmumu,Ztautau
h_lep_pt_mergedZ->SetMarkerStyle(8);
h_lep_pt_mergedZ->SetMarkerSize(0.5);
h_lep_pt_mergedZ->SetLineColor(kGreen);
//Diboson WW histogram
h_lep_pt_Dww->SetMarkerStyle(8);
h_lep_pt_Dww->SetMarkerSize(0.5);
h_lep_pt_Dww->SetLineColor(kRed);
//Diboson ZZ histogram
h_lep_pt_Dzz->SetMarkerStyle(8);
h_lep_pt_Dzz->SetMarkerSize(0.5);
h_lep_pt_Dzz->SetLineColor(kYellow);
//Diboson WZ histogram
h_lep_pt_Dwz->SetMarkerStyle(8);
h_lep_pt_Dwz->SetMarkerSize(0.5);
h_lep_pt_Dwz->SetLineColor(kBlue);
//Merged Diboson WW,ZZ,WZ histogram
h_lep_pt_mergedD->SetMarkerStyle(8);
h_lep_pt_mergedD->SetMarkerSize(0.5);
h_lep_pt_mergedD->SetLineColor(kGreen);
//DYee1 histogram
h_lep_pt_DYee1->SetMarkerStyle(8);
h_lep_pt_DYee1->SetMarkerSize(0.5);
h_lep_pt_DYee1->SetLineColor(kRed);
//DYee2 histogram
h_lep_pt_DYee2->SetMarkerStyle(8);
h_lep_pt_DYee2->SetMarkerSize(0.5);
h_lep_pt_DYee2->SetLineColor(kYellow);
//DYmumu1 histrogram
h_lep_pt_DYmumu1->SetMarkerStyle(8);
h_lep_pt_DYmumu1->SetMarkerSize(0.5);
h_lep_pt_DYmumu1->SetLineColor(kBlue);
//DYmumu2 histrogram
h_lep_pt_DYmumu2->SetMarkerStyle(8);
h_lep_pt_DYmumu2->SetMarkerSize(0.5);
h_lep_pt_DYmumu2->SetLineColor(kGreen);
//DYtautau1 histrogram
h_lep_pt_DYtautau1->SetMarkerStyle(8);
h_lep_pt_DYtautau1->SetMarkerSize(0.5);
h_lep_pt_DYtautau1->SetLineColor(kOrange);
//DYtautau2 histrogram
h_lep_pt_DYtautau2->SetMarkerStyle(8);
h_lep_pt_DYtautau2->SetMarkerSize(0.5);
h_lep_pt_DYtautau2->SetLineColor(kPink);
//Merged DY EE,MUMU,TAUTAU histogram
h_lep_pt_mergedDY->SetMarkerStyle(8);
h_lep_pt_mergedDY->SetMarkerSize(0.5);
h_lep_pt_mergedDY->SetLineColor(kGray);
//Stop tchan_top histrogram
h_lep_pt_Stchan_top->SetMarkerStyle(8);
h_lep_pt_Stchan_top->SetMarkerSize(0.5);
h_lep_pt_Stchan_top->SetLineColor(kRed);
//Stop tchan_antitop histrogram
h_lep_pt_Stchan_antitop->SetMarkerStyle(8);
h_lep_pt_Stchan_antitop->SetMarkerSize(0.5);
h_lep_pt_Stchan_antitop->SetLineColor(kYellow);
//Stop schan histogram
h_lep_pt_Sschan->SetMarkerStyle(8);
h_lep_pt_Sschan->SetMarkerSize(0.5);
h_lep_pt_Sschan->SetLineColor(kBlue);
//Stop wtchan histogram
h_lep_pt_Swtchan->SetMarkerStyle(8);
h_lep_pt_Swtchan->SetMarkerSize(0.5);
h_lep_pt_Swtchan->SetLineColor(kGreen);
//Merged stop
h_lep_pt_mergedS->SetMarkerStyle(8);
h_lep_pt_mergedS->SetMarkerSize(0.5);
h_lep_pt_mergedS->SetLineColor(kGray);
//TTBar had histrogram
h_lep_pt_TThad->SetMarkerStyle(8);
h_lep_pt_TThad->SetMarkerSize(0.5);
h_lep_pt_TThad->SetLineColor(kRed);
//TTBar lep histrogram
h_lep_pt_TTlep->SetMarkerStyle(8);
h_lep_pt_TTlep->SetMarkerSize(0.5);
h_lep_pt_TTlep->SetLineColor(kYellow);
//Merged TTBar
h_lep_pt_mergedTT->SetMarkerStyle(8);
h_lep_pt_mergedTT->SetMarkerSize(0.5);
h_lep_pt_mergedTT->SetLineColor(kRed);
//ETA HISTOGRAMS
//Zee histrogram
h_lep_eta_Zee->SetMarkerStyle(8);
h_lep_eta_Zee->SetMarkerSize(0.5);
h_lep_eta_Zee->SetLineColor(kRed);
//Zmumu histrogram
h_lep_eta_Zmumu->SetMarkerStyle(8);
h_lep_eta_Zmumu->SetMarkerSize(0.5);
h_lep_eta_Zmumu->SetLineColor(kYellow);
//Ztautau histrogram
h_lep_eta_Ztautau->SetMarkerStyle(8);
h_lep_eta_Ztautau->SetMarkerSize(0.5);
h_lep_eta_Ztautau->SetLineColor(kBlue);
//Merged Zee,Zmumu,Ztautau
h_lep_eta_mergedZ->SetMarkerStyle(8);
h_lep_eta_mergedZ->SetMarkerSize(0.5);
h_lep_eta_mergedZ->SetLineColor(kGreen);
//Diboson WW histogram
h_lep_eta_Dww->SetMarkerStyle(8);
h_lep_eta_Dww->SetMarkerSize(0.5);
h_lep_eta_Dww->SetLineColor(kRed);
//Diboson ZZ histogram
h_lep_eta_Dzz->SetMarkerStyle(8);
h_lep_eta_Dzz->SetMarkerSize(0.5);
h_lep_eta_Dzz->SetLineColor(kYellow);
//Diboson WZ histogram
h_lep_eta_Dwz->SetMarkerStyle(8);
h_lep_eta_Dwz->SetMarkerSize(0.5);
h_lep_eta_Dwz->SetLineColor(kBlue);
//Merged Diboson WW,ZZ,WZ histogram
h_lep_eta_mergedD->SetMarkerStyle(8);
h_lep_eta_mergedD->SetMarkerSize(0.5);
h_lep_eta_mergedD->SetLineColor(kGreen);
//DYee1 histogram
h_lep_eta_DYee1->SetMarkerStyle(8);
h_lep_eta_DYee1->SetMarkerSize(0.5);
h_lep_eta_DYee1->SetLineColor(kRed);
//DYee2 histogram
h_lep_eta_DYee2->SetMarkerStyle(8);
h_lep_eta_DYee2->SetMarkerSize(0.5);
h_lep_eta_DYee2->SetLineColor(kYellow);
//DYmumu1 histrogram
h_lep_eta_DYmumu1->SetMarkerStyle(8);
h_lep_eta_DYmumu1->SetMarkerSize(0.5);
h_lep_eta_DYmumu1->SetLineColor(kBlue);
//DYmumu2 histrogram
h_lep_eta_DYmumu2->SetMarkerStyle(8);
h_lep_eta_DYmumu2->SetMarkerSize(0.5);
h_lep_eta_DYmumu2->SetLineColor(kGreen);
//DYtautau1 histrogram
h_lep_eta_DYtautau1->SetMarkerStyle(8);
h_lep_eta_DYtautau1->SetMarkerSize(0.5);
h_lep_eta_DYtautau1->SetLineColor(kOrange);
//DYtautau2 histrogram
h_lep_eta_DYtautau2->SetMarkerStyle(8);
h_lep_eta_DYtautau2->SetMarkerSize(0.5);
h_lep_eta_DYtautau2->SetLineColor(kPink);
//Merged DY EE,MUMU,TAUTAU histogram
h_lep_eta_mergedDY->SetMarkerStyle(8);
h_lep_eta_mergedDY->SetMarkerSize(0.5);
h_lep_eta_mergedDY->SetLineColor(kGray);
//Stop tchan_top histrogram
h_lep_eta_Stchan_top->SetMarkerStyle(8);
h_lep_eta_Stchan_top->SetMarkerSize(0.5);
h_lep_eta_Stchan_top->SetLineColor(kRed);
//Stop tchan_antitop histrogram
h_lep_eta_Stchan_antitop->SetMarkerStyle(8);
h_lep_eta_Stchan_antitop->SetMarkerSize(0.5);
h_lep_eta_Stchan_antitop->SetLineColor(kYellow);
//Stop schan histogram
h_lep_eta_Sschan->SetMarkerStyle(8);
h_lep_eta_Sschan->SetMarkerSize(0.5);
h_lep_eta_Sschan->SetLineColor(kBlue);
//Stop wtchan histogram
h_lep_eta_Swtchan->SetMarkerStyle(8);
h_lep_eta_Swtchan->SetMarkerSize(0.5);
h_lep_eta_Swtchan->SetLineColor(kGreen);
//Merged stop
h_lep_eta_mergedS->SetMarkerStyle(8);
h_lep_eta_mergedS->SetMarkerSize(0.5);
h_lep_eta_mergedS->SetLineColor(kGray);
//TTBar had histrogram
h_lep_eta_TThad->SetMarkerStyle(8);
h_lep_eta_TThad->SetMarkerSize(0.5);
h_lep_eta_TThad->SetLineColor(kRed);
//TTBar lep histrogram
h_lep_eta_TTlep->SetMarkerStyle(8);
h_lep_eta_TTlep->SetMarkerSize(0.5);
h_lep_eta_TTlep->SetLineColor(kYellow);
//Merged TTBar
h_lep_eta_mergedTT->SetMarkerStyle(8);
h_lep_eta_mergedTT->SetMarkerSize(0.5);
h_lep_eta_mergedTT->SetLineColor(kRed);
//DATA Histograms
//Egamma
h_lep_pt_Data->SetMarkerStyle(8);
h_lep_pt_Data->SetMarkerSize(0.5);
h_lep_pt_Data->SetLineColor(kRed);
//Muons
h_lep_pt_Muons->SetMarkerStyle(8);
h_lep_pt_Muons->SetMarkerSize(0.5);
h_lep_pt_Muons->SetLineColor(kYellow);
//PHI Histograms
//Zee histrogram
h_lep_phi_Zee->SetMarkerStyle(8);
h_lep_phi_Zee->SetMarkerSize(0.5);
h_lep_phi_Zee->SetLineColor(kRed);
//Zmumu histrogram
h_lep_phi_Zmumu->SetMarkerStyle(8);
h_lep_phi_Zmumu->SetMarkerSize(0.5);
h_lep_phi_Zmumu->SetLineColor(kYellow);
//Ztautau histrogram
h_lep_phi_Ztautau->SetMarkerStyle(8);
h_lep_phi_Ztautau->SetMarkerSize(0.5);
h_lep_phi_Ztautau->SetLineColor(kBlue);
//Merged Zee,Zmumu,Ztautau
h_lep_phi_mergedZ->SetMarkerStyle(8);
h_lep_phi_mergedZ->SetMarkerSize(0.5);
h_lep_phi_mergedZ->SetLineColor(kGreen);
//Diboson WW histogram
h_lep_phi_Dww->SetMarkerStyle(8);
h_lep_phi_Dww->SetMarkerSize(0.5);
h_lep_phi_Dww->SetLineColor(kRed);
//Diboson ZZ histogram
h_lep_phi_Dzz->SetMarkerStyle(8);
h_lep_phi_Dzz->SetMarkerSize(0.5);
h_lep_phi_Dzz->SetLineColor(kYellow);
//Diboson WZ histogram
h_lep_phi_Dwz->SetMarkerStyle(8);
h_lep_phi_Dwz->SetMarkerSize(0.5);
h_lep_phi_Dwz->SetLineColor(kBlue);
//Merged Diboson WW,ZZ,WZ histogram
h_lep_phi_mergedD->SetMarkerStyle(8);
h_lep_phi_mergedD->SetMarkerSize(0.5);
h_lep_phi_mergedD->SetLineColor(kGreen);
//DYee1 histogram
h_lep_phi_DYee1->SetMarkerStyle(8);
h_lep_phi_DYee1->SetMarkerSize(0.5);
h_lep_phi_DYee1->SetLineColor(kRed);
//DYee2 histogram
h_lep_phi_DYee2->SetMarkerStyle(8);
h_lep_phi_DYee2->SetMarkerSize(0.5);
h_lep_phi_DYee2->SetLineColor(kYellow);
//DYmumu1 histrogram
h_lep_phi_DYmumu1->SetMarkerStyle(8);
h_lep_phi_DYmumu1->SetMarkerSize(0.5);
h_lep_phi_DYmumu1->SetLineColor(kBlue);
//DYmumu2 histrogram
h_lep_phi_DYmumu2->SetMarkerStyle(8);
h_lep_phi_DYmumu2->SetMarkerSize(0.5);
h_lep_phi_DYmumu2->SetLineColor(kGreen);
//DYtautau1 histrogram
h_lep_phi_DYtautau1->SetMarkerStyle(8);
h_lep_phi_DYtautau1->SetMarkerSize(0.5);
h_lep_phi_DYtautau1->SetLineColor(kOrange);
//DYtautau2 histrogram
h_lep_phi_DYtautau2->SetMarkerStyle(8);
h_lep_phi_DYtautau2->SetMarkerSize(0.5);
h_lep_phi_DYtautau2->SetLineColor(kPink);
//Merged DY EE,MUMU,TAUTAU histogram
h_lep_phi_mergedDY->SetMarkerStyle(8);
h_lep_phi_mergedDY->SetMarkerSize(0.5);
h_lep_phi_mergedDY->SetLineColor(kGray);
//Stop tchan_top histrogram
h_lep_phi_Stchan_top->SetMarkerStyle(8);
h_lep_phi_Stchan_top->SetMarkerSize(0.5);
h_lep_phi_Stchan_top->SetLineColor(kRed);
//Stop tchan_antitop histrogram
h_lep_phi_Stchan_antitop->SetMarkerStyle(8);
h_lep_phi_Stchan_antitop->SetMarkerSize(0.5);
h_lep_phi_Stchan_antitop->SetLineColor(kYellow);
//Stop schan histogram
h_lep_phi_Sschan->SetMarkerStyle(8);
h_lep_phi_Sschan->SetMarkerSize(0.5);
h_lep_phi_Sschan->SetLineColor(kBlue);
//Stop wtchan histogram
h_lep_phi_Swtchan->SetMarkerStyle(8);
h_lep_phi_Swtchan->SetMarkerSize(0.5);
h_lep_phi_Swtchan->SetLineColor(kGreen);
//Merged stop
h_lep_phi_mergedS->SetMarkerStyle(8);
h_lep_phi_mergedS->SetMarkerSize(0.5);
h_lep_phi_mergedS->SetLineColor(kGray);
//TTBar had histrogram
h_lep_phi_TThad->SetMarkerStyle(8);
h_lep_phi_TThad->SetMarkerSize(0.5);
h_lep_phi_TThad->SetLineColor(kRed);
//TTBar lep histrogram
h_lep_phi_TTlep->SetMarkerStyle(8);
h_lep_phi_TTlep->SetMarkerSize(0.5);
h_lep_phi_TTlep->SetLineColor(kYellow);
//Merged TTBar
h_lep_phi_mergedTT->SetMarkerStyle(8);
h_lep_phi_mergedTT->SetMarkerSize(0.5);
h_lep_phi_mergedTT->SetLineColor(kRed);
//E Histograms
//Zee histrogram
h_lep_E_Zee->SetMarkerStyle(8);
h_lep_E_Zee->SetMarkerSize(0.5);
h_lep_E_Zee->SetLineColor(kRed);
//Zmumu histrogram
h_lep_E_Zmumu->SetMarkerStyle(8);
h_lep_E_Zmumu->SetMarkerSize(0.5);
h_lep_E_Zmumu->SetLineColor(kYellow);
//Ztautau histrogram
h_lep_E_Ztautau->SetMarkerStyle(8);
h_lep_E_Ztautau->SetMarkerSize(0.5);
h_lep_E_Ztautau->SetLineColor(kBlue);
//Merged Zee,Zmumu,Ztautau
h_lep_E_mergedZ->SetMarkerStyle(8);
h_lep_E_mergedZ->SetMarkerSize(0.5);
h_lep_E_mergedZ->SetLineColor(kGreen);
//Diboson WW histogram
h_lep_E_Dww->SetMarkerStyle(8);
h_lep_E_Dww->SetMarkerSize(0.5);
h_lep_E_Dww->SetLineColor(kRed);
//Diboson ZZ histogram
h_lep_E_Dzz->SetMarkerStyle(8);
h_lep_E_Dzz->SetMarkerSize(0.5);
h_lep_E_Dzz->SetLineColor(kYellow);
//Diboson WZ histogram
h_lep_E_Dwz->SetMarkerStyle(8);
h_lep_E_Dwz->SetMarkerSize(0.5);
h_lep_E_Dwz->SetLineColor(kBlue);
//Merged Diboson WW,ZZ,WZ histogram
h_lep_E_mergedD->SetMarkerStyle(8);
h_lep_E_mergedD->SetMarkerSize(0.5);
h_lep_E_mergedD->SetLineColor(kGreen);
//DYee1 histogram
h_lep_E_DYee1->SetMarkerStyle(8);
h_lep_E_DYee1->SetMarkerSize(0.5);
h_lep_E_DYee1->SetLineColor(kRed);
//DYee2 histogram
h_lep_E_DYee2->SetMarkerStyle(8);
h_lep_E_DYee2->SetMarkerSize(0.5);
h_lep_E_DYee2->SetLineColor(kYellow);
//DYmumu1 histrogram
h_lep_E_DYmumu1->SetMarkerStyle(8);
h_lep_E_DYmumu1->SetMarkerSize(0.5);
h_lep_E_DYmumu1->SetLineColor(kBlue);
//DYmumu2 histrogram
h_lep_E_DYmumu2->SetMarkerStyle(8);
h_lep_E_DYmumu2->SetMarkerSize(0.5);
h_lep_E_DYmumu2->SetLineColor(kGreen);
//DYtautau1 histrogram
h_lep_E_DYtautau1->SetMarkerStyle(8);
h_lep_E_DYtautau1->SetMarkerSize(0.5);
h_lep_E_DYtautau1->SetLineColor(kOrange);
//DYtautau2 histrogram
h_lep_E_DYtautau2->SetMarkerStyle(8);
h_lep_E_DYtautau2->SetMarkerSize(0.5);
h_lep_E_DYtautau2->SetLineColor(kPink);
//Merged DY EE,MUMU,TAUTAU histogram
h_lep_E_mergedDY->SetMarkerStyle(8);
h_lep_E_mergedDY->SetMarkerSize(0.5);
h_lep_E_mergedDY->SetLineColor(kGray);
//Stop tchan_top histrogram
h_lep_E_Stchan_top->SetMarkerStyle(8);
h_lep_E_Stchan_top->SetMarkerSize(0.5);
h_lep_E_Stchan_top->SetLineColor(kRed);
//Stop tchan_antitop histrogram
h_lep_E_Stchan_antitop->SetMarkerStyle(8);
h_lep_E_Stchan_antitop->SetMarkerSize(0.5);
h_lep_E_Stchan_antitop->SetLineColor(kYellow);
//Stop schan histogram
h_lep_E_Sschan->SetMarkerStyle(8);
h_lep_E_Sschan->SetMarkerSize(0.5);
h_lep_E_Sschan->SetLineColor(kBlue);
//Stop wtchan histogram
h_lep_E_Swtchan->SetMarkerStyle(8);
h_lep_E_Swtchan->SetMarkerSize(0.5);
h_lep_E_Swtchan->SetLineColor(kGreen);
//Merged stop
h_lep_E_mergedS->SetMarkerStyle(8);
h_lep_E_mergedS->SetMarkerSize(0.5);
h_lep_E_mergedS->SetLineColor(kGray);
//TTBar had histrogram
h_lep_E_TThad->SetMarkerStyle(8);
h_lep_E_TThad->SetMarkerSize(0.5);
h_lep_E_TThad->SetLineColor(kRed);
//TTBar lep histrogram
h_lep_E_TTlep->SetMarkerStyle(8);
h_lep_E_TTlep->SetMarkerSize(0.5);
h_lep_E_TTlep->SetLineColor(kYellow);
//Merged TTBar
h_lep_E_mergedTT->SetMarkerStyle(8);
h_lep_E_mergedTT->SetMarkerSize(0.5);
h_lep_E_mergedTT->SetLineColor(kRed);
int nentriesMC, nbytesMC, i;
nentriesMC = (Int_t)MonteCarlo->GetEntries();
std::cout << nentriesMC;
for (i=0;i<nentriesMC; i++){
nbytesMC = MonteCarlo->GetEntry(i);
//SF
Float_t scaleFactor = MonteCarlo_scaleFactor_ELE*MonteCarlo_scaleFactor_MUON* MonteCarlo_scaleFactor_TRIGGER;
//EventW
Float_t eventWeight = MonteCarlo_mcWeight*MonteCarlo_scaleFactor_PILEUP*MonteCarlo_scaleFactor_ZVERTEX;
//weight = SF * EventW
Float_t weight = scaleFactor*eventWeight;
//first cut
if(MonteCarlo_lep_n>1){
//second cut
Float_t OS = MonteCarlo_lep_charge[0]*MonteCarlo_lep_charge[1];
if(OS == -1.){
//third cut
if(MonteCarlo_lep_type[0] ==MonteCarlo_lep_type[1]){
//forth cut
TLorentzVector Lepton_1 = TLorentzVector();
TLorentzVector Lepton_2 = TLorentzVector();
Lepton_1.SetPtEtaPhiE( MonteCarlo_lep_pt[0]/1000. , MonteCarlo_lep_eta[0] , MonteCarlo_lep_phi[0] , MonteCarlo_lep_E[0]/1000.);
Lepton_2.SetPtEtaPhiE( MonteCarlo_lep_pt[1]/1000. , MonteCarlo_lep_eta[1] , MonteCarlo_lep_phi[1] , MonteCarlo_lep_E[1]/1000.);
Float_t Z_Mass = 91.188;
Float_t ll_mass = TMath::Abs( (Lepton_1 + Lepton_2).M());
if(ll_mass - Z_Mass < 20.){
//drawing histograms
//Zee
if(MonteCarlo_runNumber ==147770){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Zee->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Zee->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Zee->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Zee->Fill(MonteCarlo_lep_E_inGeV,weight);
}
//Zmumu
else if(MonteCarlo_runNumber ==147771){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Zmumu->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Zmumu->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Zmumu->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Zmumu->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//Ztautau
else if(MonteCarlo_runNumber ==147772){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Ztautau->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Ztautau->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Ztautau->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Ztautau->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//Diboson WW
else if(MonteCarlo_runNumber ==105985){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Dww->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Dww->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Dww->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Dww->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//Diboson ZZ
else if(MonteCarlo_runNumber ==105986){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Dzz->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Dzz->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Dzz->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Dzz->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//Diboson WZ
else if(MonteCarlo_runNumber ==105987){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Dwz->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Dwz->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Dwz->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Dwz->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//DY EE1
else if(MonteCarlo_runNumber ==173041){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_DYee1->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_DYee1->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_DYee1->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_DYee1->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//DY EE2
else if(MonteCarlo_runNumber ==173042){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_DYee2->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_DYee2->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_DYee2->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_DYee2->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//DY MUMU1
else if(MonteCarlo_runNumber ==173043){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_DYmumu1->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_DYmumu1->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_DYmumu1->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_DYmumu1->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//DY MUMU2
else if(MonteCarlo_runNumber ==173044){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_DYmumu2->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_DYmumu2->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_DYmumu2->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_DYmumu2->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//DY TAUTAU1
else if(MonteCarlo_runNumber ==173045){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_DYtautau1->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_DYtautau1->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_DYtautau1->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_DYtautau1->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//DY TAUTAU2
else if(MonteCarlo_runNumber ==173046){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_DYtautau2->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_DYtautau2->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_DYtautau2->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_DYtautau2->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//Stop tchan top
else if(MonteCarlo_runNumber ==110090){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Stchan_top->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Stchan_top->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Stchan_top->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Stchan_top->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//Stop tchan antitop
else if(MonteCarlo_runNumber ==110091){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Stchan_antitop->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Stchan_antitop->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Stchan_antitop->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Stchan_antitop->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//Stop schan
else if(MonteCarlo_runNumber ==110119){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Sschan->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Sschan->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Sschan->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Sschan->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//Stop wtchan
else if(MonteCarlo_runNumber ==110140){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Swtchan->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_Swtchan->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_Swtchan->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_Swtchan->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//TTBar had
else if(MonteCarlo_runNumber ==117049){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_TThad->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_TThad->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_TThad->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_TThad->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//TTBar lep
else if(MonteCarlo_runNumber ==117050){
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_TTlep->Fill(MonteCarlo_lep_pt_inGeV, weight);
//ETA
h_lep_eta_TTlep->Fill(MonteCarlo_lep_eta[0], weight);
//PHI
h_lep_phi_TTlep->Fill(MonteCarlo_lep_phi[0], weight);
//E
float MonteCarlo_lep_E_inGeV = MonteCarlo_lep_E[0]/1000.;
h_lep_E_TTlep->Fill(MonteCarlo_lep_E_inGeV, weight);
}
//Data egamma histogram
else{
//PT
float MonteCarlo_lep_pt_inGeV = MonteCarlo_lep_pt[0]/1000.;
h_lep_pt_Data->Fill(MonteCarlo_lep_pt_inGeV, 1);
}
}
}
}
}
}
36288230
//SC Z
scaleFactorZee=(Luminosity*1241.2072)/(203795455568148.0*0.151816307);
scaleFactorZmumu=(Luminosity*1241.2072)/(225316022111048.0*0.124706551);
scaleFactorZtautau= (Luminosity*1240.8988)/(31508540303680.9*0.921502294);
//SC Diboson
scaleFactorDww=(Luminosity*20.90286)/(2489550.66/0.388);
scaleFactorDzz=(Luminosity*1.5376)/(243674.74/0.95088);
scaleFactorDwz=(Luminosity*6.9673)/(996210.5/0.96668);
//SC Drell Yan
scaleFactorDYee1=(Luminosity*92.15)/(4772549.624/0.893071694);
scaleFactorDYee2=(Luminosity*279.19)/(13895470.42/0.945799529);
scaleFactorDYmumu1=(Luminosity*92.08)/(4977032.008/0.960298189);
scaleFactorDYmumu2=(Luminosity*279.2)/(11510844.32/0.999654588);
scaleFactorDYtautau1=(Luminosity*92.12)/(2489620.309/1);
scaleFactorDYtautau2=(Luminosity*279.11)/(13704752.81/1);
//SC Stop
scaleFactorStchan_top=(Luminosity*18.39495)/(0.087147591/0.08947185);
scaleFactorStchan_antitop=(Luminosity*9.95976)/(0.046283642/0.058425933);
scaleFactorSschan=(Luminosity*1.817694)/(5972600.588/0.051100576);
scaleFactorSwtchan=(Luminosity*22.3014)/(996333.734/0.636162836);
//SC TTBar
scaleFactorTThad=(Luminosity*115.518854)/(115.518854/1);
scaleFactorTTlep=(Luminosity*137.29749)/(49761200.21/0.072212854);
//PT Histograms Scale
h_lep_pt_Zee->Scale(scaleFactorZee);
h_lep_pt_Zmumu->Scale(scaleFactorZmumu);
h_lep_pt_Ztautau->Scale(scaleFactorZtautau);
h_lep_pt_Dww->Scale(scaleFactorDww);
h_lep_pt_Dzz->Scale(scaleFactorDzz);
h_lep_pt_Dwz->Scale(scaleFactorDwz);
h_lep_pt_DYee1->Scale(scaleFactorDYee1);
h_lep_pt_DYee2->Scale(scaleFactorDYee2);
h_lep_pt_DYmumu1->Scale(scaleFactorDYmumu1);
h_lep_pt_DYmumu2->Scale(scaleFactorDYmumu2);
h_lep_pt_DYtautau1->Scale(scaleFactorDYtautau1);
h_lep_pt_DYtautau2->Scale(scaleFactorDYtautau2);
h_lep_pt_Stchan_top->Scale(scaleFactorStchan_top);
h_lep_pt_Stchan_antitop->Scale(scaleFactorStchan_antitop);
h_lep_pt_Sschan->Scale(scaleFactorSschan);
h_lep_pt_Swtchan->Scale(scaleFactorSwtchan);
h_lep_pt_TThad->Scale(scaleFactorTThad);
h_lep_pt_TTlep->Scale(scaleFactorTTlep);
//ETA Histograms Scale
h_lep_eta_Zee->Scale(scaleFactorZee);
h_lep_eta_Zmumu->Scale(scaleFactorZmumu);
h_lep_eta_Ztautau->Scale(scaleFactorZtautau);
h_lep_eta_Dww->Scale(scaleFactorDww);
h_lep_eta_Dzz->Scale(scaleFactorDzz);
h_lep_eta_Dwz->Scale(scaleFactorDwz);
h_lep_eta_DYee1->Scale(scaleFactorDYee1);
h_lep_eta_DYee2->Scale(scaleFactorDYee2);
h_lep_eta_DYmumu1->Scale(scaleFactorDYmumu1);
h_lep_eta_DYmumu2->Scale(scaleFactorDYmumu2);
h_lep_eta_DYtautau1->Scale(scaleFactorDYtautau1);
h_lep_eta_DYtautau2->Scale(scaleFactorDYtautau2);
h_lep_eta_Stchan_top->Scale(scaleFactorStchan_top);
h_lep_eta_Stchan_antitop->Scale(scaleFactorStchan_antitop);
h_lep_eta_Sschan->Scale(scaleFactorSschan);
h_lep_eta_Swtchan->Scale(scaleFactorSwtchan);
h_lep_eta_TThad->Scale(scaleFactorTThad);
h_lep_eta_TTlep->Scale(scaleFactorTTlep);
//PHI Histograms Scale
h_lep_phi_Zee->Scale(scaleFactorZee);
h_lep_phi_Zmumu->Scale(scaleFactorZmumu);
h_lep_phi_Ztautau->Scale(scaleFactorZtautau);
h_lep_phi_Dww->Scale(scaleFactorDww);
h_lep_phi_Dzz->Scale(scaleFactorDzz);
h_lep_phi_Dwz->Scale(scaleFactorDwz);
h_lep_phi_DYee1->Scale(scaleFactorDYee1);
h_lep_phi_DYee2->Scale(scaleFactorDYee2);
h_lep_phi_DYmumu1->Scale(scaleFactorDYmumu1);
h_lep_phi_DYmumu2->Scale(scaleFactorDYmumu2);
h_lep_phi_DYtautau1->Scale(scaleFactorDYtautau1);
h_lep_phi_DYtautau2->Scale(scaleFactorDYtautau2);
h_lep_phi_Stchan_top->Scale(scaleFactorStchan_top);
h_lep_phi_Stchan_antitop->Scale(scaleFactorStchan_antitop);
h_lep_phi_Sschan->Scale(scaleFactorSschan);
h_lep_phi_Swtchan->Scale(scaleFactorSwtchan);
h_lep_phi_TThad->Scale(scaleFactorTThad);
h_lep_phi_TTlep->Scale(scaleFactorTTlep);
//E Histograms Scale
h_lep_E_Zee->Scale(scaleFactorZee);
h_lep_E_Zmumu->Scale(scaleFactorZmumu);
h_lep_E_Ztautau->Scale(scaleFactorZtautau);
h_lep_E_Dww->Scale(scaleFactorDww);
h_lep_E_Dzz->Scale(scaleFactorDzz);
h_lep_E_Dwz->Scale(scaleFactorDwz);
h_lep_E_DYee1->Scale(scaleFactorDYee1);
h_lep_E_DYee2->Scale(scaleFactorDYee2);
h_lep_E_DYmumu1->Scale(scaleFactorDYmumu1);
h_lep_E_DYmumu2->Scale(scaleFactorDYmumu2);
h_lep_E_DYtautau1->Scale(scaleFactorDYtautau1);
h_lep_E_DYtautau2->Scale(scaleFactorDYtautau2);
h_lep_E_Stchan_top->Scale(scaleFactorStchan_top);
h_lep_E_Stchan_antitop->Scale(scaleFactorStchan_antitop);
h_lep_E_Sschan->Scale(scaleFactorSschan);
h_lep_E_Swtchan->Scale(scaleFactorSwtchan);
h_lep_E_TThad->Scale(scaleFactorTThad);
h_lep_E_TTlep->Scale(scaleFactorTTlep);
TCanvas *cz = new TCanvas("cz","cz",10,10,700,900);
TText tz; tz.SetTextFont(42); tz.SetTextAlign(21);
cz->Divide(2,2);
cz->cd(1); h_lep_pt_Zee->Draw();
cz->cd(2); h_lep_pt_Zmumu->Draw();
cz->cd(3); h_lep_pt_Ztautau->Draw();
cz->cd(4); h_lep_pt_mergedZ->Draw();
TList *listZ = new TList;
listZ->Add(h_lep_pt_Zee);
listZ->Add(h_lep_pt_Zmumu);
listZ->Add(h_lep_pt_Ztautau);
h_lep_pt_mergedZ->Reset();
h_lep_pt_mergedZ->Merge(listZ);
h_lep_pt_mergedZ->SetFillColor(kBlue);
h_lep_pt_mergedZ->Draw();
cz->Draw();
TCanvas *czEta = new TCanvas("czEta","czEta",10,10,700,900);
TText tzeta; tzeta.SetTextFont(42); tzeta.SetTextAlign(21);
czEta->Divide(2,2);
czEta->cd(1); h_lep_eta_Zee->Draw();
czEta->cd(2); h_lep_eta_Zmumu->Draw();
czEta->cd(3); h_lep_eta_Ztautau->Draw();
czEta->cd(4); h_lep_eta_mergedZ->Draw();
TList *listZEta = new TList;
listZEta->Add(h_lep_eta_Zee);
listZEta->Add(h_lep_eta_Zmumu);
listZEta->Add(h_lep_eta_Ztautau);
h_lep_eta_mergedZ->Reset();
h_lep_eta_mergedZ->Merge(listZEta);
h_lep_eta_mergedZ->SetFillColor(kBlue);
h_lep_eta_mergedZ->Draw();
czEta->Draw();
TCanvas *czPhi = new TCanvas("czPhi","czPhi",10,10,700,900);
TText tzphi; tzphi.SetTextFont(42); tzphi.SetTextAlign(21);
czPhi->Divide(2,2);
czPhi->cd(1); h_lep_phi_Zee->Draw();
czPhi->cd(2); h_lep_phi_Zmumu->Draw();
czPhi->cd(3); h_lep_phi_Ztautau->Draw();
czPhi->cd(4); h_lep_phi_mergedZ->Draw();
TList *listZPhi = new TList;
listZPhi->Add(h_lep_phi_Zee);
listZPhi->Add(h_lep_phi_Zmumu);
listZPhi->Add(h_lep_phi_Ztautau);
h_lep_phi_mergedZ->Reset();
h_lep_phi_mergedZ->Merge(listZPhi);
h_lep_phi_mergedZ->SetFillColor(kBlue);
h_lep_phi_mergedZ->Draw();
czPhi->Draw();
TCanvas *czE = new TCanvas("czE","czE",10,10,700,900);
TText tzE; tzE.SetTextFont(42); tzE.SetTextAlign(21);
czE->Divide(2,2);
czE->cd(1); h_lep_E_Zee->Draw();
czE->cd(2); h_lep_E_Zmumu->Draw();
czE->cd(3); h_lep_E_Ztautau->Draw();
czE->cd(4); h_lep_E_mergedZ->Draw();
TList *listZE = new TList;
listZE->Add(h_lep_E_Zee);
listZE->Add(h_lep_E_Zmumu);
listZE->Add(h_lep_E_Ztautau);
h_lep_E_mergedZ->Reset();
h_lep_E_mergedZ->Merge(listZE);
h_lep_E_mergedZ->SetFillColor(kBlue);
h_lep_E_mergedZ->Draw();
czE->Draw();
TCanvas *cdata = new TCanvas("cdata","cdata",10,10,700,900);
TText tdata; tdata.SetTextFont(42); tdata.SetTextAlign(21);
cdata->Divide(2,2);
cdata->cd(1); h_lep_pt_Data->Draw();
cdata->Draw();
TCanvas *cd = new TCanvas("cd","cd",10,10,700,900);
TText td; td.SetTextFont(42); td.SetTextAlign(21);
cd->Divide(2,2);
cd->cd(1); h_lep_pt_Dww->Draw();
cd->cd(2); h_lep_pt_Dzz->Draw();
cd->cd(3); h_lep_pt_Dwz->Draw();
cd->cd(4); h_lep_pt_mergedD->Draw();
TList *listD = new TList;
listD->Add(h_lep_pt_Dww);
listD->Add(h_lep_pt_Dzz);
listD->Add(h_lep_pt_Dwz);
h_lep_pt_mergedD->Reset();
h_lep_pt_mergedD->Merge(listD);
h_lep_pt_mergedD->SetFillColor(kBlue);
h_lep_pt_mergedD->Draw();
cd->Draw();
TCanvas *cdEta = new TCanvas("cdeta","cdeta",10,10,700,900);
TText tdeta; tdeta.SetTextFont(42); tdeta.SetTextAlign(21);
cdEta->Divide(2,2);
cdEta->cd(1); h_lep_eta_Dww->Draw();
cdEta->cd(2); h_lep_eta_Dzz->Draw();
cdEta->cd(3); h_lep_eta_Dwz->Draw();
cdEta->cd(4); h_lep_eta_mergedD->Draw();
TList *listDEta = new TList;
listDEta->Add(h_lep_eta_Dww);
listDEta->Add(h_lep_eta_Dzz);
listDEta->Add(h_lep_eta_Dwz);
h_lep_eta_mergedD->Reset();
h_lep_eta_mergedD->Merge(listDEta);
h_lep_eta_mergedD->SetFillColor(kBlue);
h_lep_eta_mergedD->Draw();
cdEta->Draw();
TCanvas *cdPhi = new TCanvas("cdPhi","cdPhi",10,10,700,900);
TText tdphi; tdphi.SetTextFont(42); tdphi.SetTextAlign(21);
cdPhi->Divide(2,2);
cdPhi->cd(1); h_lep_phi_Dww->Draw();
cdPhi->cd(2); h_lep_phi_Dzz->Draw();
cdPhi->cd(3); h_lep_phi_Dwz->Draw();
cdPhi->cd(4); h_lep_phi_mergedD->Draw();
TList *listDPhi = new TList;
listDPhi->Add(h_lep_phi_Dww);
listDPhi->Add(h_lep_phi_Dzz);
listDPhi->Add(h_lep_phi_Dwz);
h_lep_phi_mergedD->Reset();
h_lep_phi_mergedD->Merge(listDPhi);
h_lep_phi_mergedD->SetFillColor(kBlue);
h_lep_phi_mergedD->Draw();
cdPhi->Draw();
TCanvas *cdE = new TCanvas("cdE","cdE",10,10,700,900);
TText tdE; tdE.SetTextFont(42); tdE.SetTextAlign(21);
cdE->Divide(2,2);
cdE->cd(1); h_lep_E_Dww->Draw();
cdE->cd(2); h_lep_E_Dzz->Draw();
cdE->cd(3); h_lep_E_Dwz->Draw();
cdE->cd(4); h_lep_E_mergedD->Draw();
TList *listDE = new TList;
listDE->Add(h_lep_E_Dww);
listDE->Add(h_lep_E_Dzz);
listDE->Add(h_lep_E_Dwz);
h_lep_E_mergedD->Reset();
h_lep_E_mergedD->Merge(listDE);
h_lep_E_mergedD->SetFillColor(kBlue);
h_lep_E_mergedD->Draw();
cdE->Draw();
TCanvas *cdy = new TCanvas("cdy","cdy",10,10,900,1000);
TText tdy; tdy.SetTextFont(42); tdy.SetTextAlign(21);
cdy->Divide(4,4);
cdy->cd(1); h_lep_pt_DYee1->Draw();
cdy->cd(2); h_lep_pt_DYee2->Draw();
cdy->cd(3); h_lep_pt_DYmumu1->Draw();
cdy->cd(4); h_lep_pt_DYmumu2->Draw();
cdy->cd(5); h_lep_pt_DYtautau1->Draw();
cdy->cd(6); h_lep_pt_DYtautau2->Draw();
cdy->cd(7); h_lep_pt_mergedDY->Draw();
TList *listDY = new TList;
listDY->Add(h_lep_pt_DYee1);
listDY->Add(h_lep_pt_DYee2);
listDY->Add(h_lep_pt_DYmumu1);
listDY->Add(h_lep_pt_DYmumu2);
listDY->Add(h_lep_pt_DYtautau1);
listDY->Add(h_lep_pt_DYtautau2);
h_lep_pt_mergedDY->Reset();
h_lep_pt_mergedDY->Merge(listDY);
h_lep_pt_mergedDY->SetFillColor(kBlue);
h_lep_pt_mergedDY->Draw();
cdy->Draw();
TCanvas *cdyEta = new TCanvas("cdyEta","cdyEta",10,10,900,1300);
TText tdyEta; tdyEta.SetTextFont(42); tdyEta.SetTextAlign(21);
cdyEta->Divide(4,4);
cdyEta->cd(1); h_lep_eta_DYee1->Draw();
cdyEta->cd(2); h_lep_eta_DYee2->Draw();
cdyEta->cd(3); h_lep_eta_DYmumu1->Draw();
cdyEta->cd(4); h_lep_eta_DYmumu2->Draw();
cdyEta->cd(5); h_lep_eta_DYtautau1->Draw();
cdyEta->cd(6); h_lep_eta_DYtautau2->Draw();
cdyEta->cd(7); h_lep_eta_mergedDY->Draw();
TList *listDYEta = new TList;
listDYEta->Add(h_lep_eta_DYee1);
listDYEta->Add(h_lep_eta_DYee2);
listDYEta->Add(h_lep_eta_DYmumu1);
listDYEta->Add(h_lep_eta_DYmumu2);
listDYEta->Add(h_lep_eta_DYtautau1);
listDYEta->Add(h_lep_eta_DYtautau2);
h_lep_eta_mergedDY->Reset();
h_lep_eta_mergedDY->Merge(listDYEta);
h_lep_eta_mergedDY->SetFillColor(kBlue);
h_lep_eta_mergedDY->Draw();
cdyEta->Draw();
TCanvas *cdyPhi = new TCanvas("cdyPhi","cdyPhi",10,10,900,1300);
TText tdyPhi; tdyPhi.SetTextFont(42); tdyPhi.SetTextAlign(21);
cdyPhi->Divide(4,4);
cdyPhi->cd(1); h_lep_phi_DYee1->Draw();
cdyPhi->cd(2); h_lep_phi_DYee2->Draw();
cdyPhi->cd(3); h_lep_phi_DYmumu1->Draw();
cdyPhi->cd(4); h_lep_phi_DYmumu2->Draw();
cdyPhi->cd(5); h_lep_phi_DYtautau1->Draw();
cdyPhi->cd(6); h_lep_phi_DYtautau2->Draw();
cdyPhi->cd(7); h_lep_phi_mergedDY->Draw();
TList *listDYPhi = new TList;
listDYPhi->Add(h_lep_phi_DYee1);
listDYPhi->Add(h_lep_phi_DYee2);
listDYPhi->Add(h_lep_phi_DYmumu1);
listDYPhi->Add(h_lep_phi_DYmumu2);
listDYPhi->Add(h_lep_phi_DYtautau1);
listDYPhi->Add(h_lep_phi_DYtautau2);
h_lep_phi_mergedDY->Reset();
h_lep_phi_mergedDY->Merge(listDYPhi);
h_lep_phi_mergedDY->SetFillColor(kBlue);
h_lep_phi_mergedDY->Draw();
cdyPhi->Draw();
TCanvas *cdyE = new TCanvas("cdyE","cdyE",10,10,900,1300);
TText tdyE; tdyE.SetTextFont(42); tdyE.SetTextAlign(21);
cdyE->Divide(4,4);
cdyE->cd(1); h_lep_E_DYee1->Draw();
cdyE->cd(2); h_lep_E_DYee2->Draw();
cdyE->cd(3); h_lep_E_DYmumu1->Draw();
cdyE->cd(4); h_lep_E_DYmumu2->Draw();
cdyE->cd(5); h_lep_E_DYtautau1->Draw();
cdyE->cd(6); h_lep_E_DYtautau2->Draw();
cdyE->cd(7); h_lep_E_mergedDY->Draw();
TList *listDYE = new TList;
listDYE->Add(h_lep_E_DYee1);
listDYE->Add(h_lep_E_DYee2);
listDYE->Add(h_lep_E_DYmumu1);
listDYE->Add(h_lep_E_DYmumu2);
listDYE->Add(h_lep_E_DYtautau1);
listDYE->Add(h_lep_E_DYtautau2);
h_lep_E_mergedDY->Reset();
h_lep_E_mergedDY->Merge(listDYE);
h_lep_E_mergedDY->SetFillColor(kBlue);
h_lep_E_mergedDY->Draw();
cdyE->Draw();
TCanvas *cstop = new TCanvas("cstop","cstop",10,10,700,900);
TText tstop; tstop.SetTextFont(42); tstop.SetTextAlign(21);
cstop->Divide(2,2);
cstop->cd(1); h_lep_pt_Stchan_top->Draw();
cstop->cd(2); h_lep_pt_Stchan_antitop->Draw();
cstop->cd(3); h_lep_pt_Sschan->Draw();
cstop->cd(4); h_lep_pt_Swtchan->Draw();
TList *listSt = new TList;
listSt->Add(h_lep_pt_Stchan_top);
listSt->Add(h_lep_pt_Stchan_antitop);
listSt->Add(h_lep_pt_Sschan);
listSt->Add(h_lep_pt_Swtchan);
h_lep_pt_mergedS->Reset();
h_lep_pt_mergedS->Merge(listSt);
h_lep_pt_mergedS->SetFillColor(kBlue);
h_lep_pt_mergedS->Draw();
cstop->Draw();
TCanvas *cstopEta = new TCanvas("cstopEta","cstopEta",10,10,700,900);
TText tstopEta; tstopEta.SetTextFont(42); tstopEta.SetTextAlign(21);
cstopEta->Divide(2,2);
cstopEta->cd(1); h_lep_eta_Stchan_top->Draw();
cstopEta->cd(2); h_lep_eta_Stchan_antitop->Draw();
cstopEta->cd(3); h_lep_eta_Sschan->Draw();
cstopEta->cd(4); h_lep_eta_Swtchan->Draw();
TList *listStEta = new TList;
listStEta->Add(h_lep_eta_Stchan_top);
listStEta->Add(h_lep_eta_Stchan_antitop);
listStEta->Add(h_lep_eta_Sschan);
listStEta->Add(h_lep_eta_Swtchan);
h_lep_eta_mergedS->Reset();
h_lep_eta_mergedS->Merge(listStEta);
h_lep_eta_mergedS->SetFillColor(kBlue);
h_lep_eta_mergedS->Draw();
cstopEta->Draw();
TCanvas *cstopPhi = new TCanvas("cstopPhi","cstopPhi",10,10,700,900);
TText tstopPhi; tstopPhi.SetTextFont(42); tstopPhi.SetTextAlign(21);
cstopPhi->Divide(2,2);
cstopPhi->cd(1); h_lep_phi_Stchan_top->Draw();
cstopPhi->cd(2); h_lep_phi_Stchan_antitop->Draw();
cstopPhi->cd(3); h_lep_phi_Sschan->Draw();
cstopPhi->cd(4); h_lep_phi_Swtchan->Draw();
TList *listStPhi = new TList;
listStPhi->Add(h_lep_phi_Stchan_top);
listStPhi->Add(h_lep_phi_Stchan_antitop);
listStPhi->Add(h_lep_phi_Sschan);
listStPhi->Add(h_lep_phi_Swtchan);
h_lep_phi_mergedS->Reset();
h_lep_phi_mergedS->Merge(listStPhi);
h_lep_phi_mergedS->SetFillColor(kBlue);
h_lep_phi_mergedS->Draw();
cstopPhi->Draw();
TCanvas *cstopE = new TCanvas("cstopE","cstopE",10,10,700,900);
TText tstopE; tstopE.SetTextFont(42); tstopE.SetTextAlign(21);
cstopE->Divide(2,2);
cstopE->cd(1); h_lep_E_Stchan_top->Draw();
cstopE->cd(2); h_lep_E_Stchan_antitop->Draw();
cstopE->cd(3); h_lep_E_Sschan->Draw();
cstopE->cd(4); h_lep_E_Swtchan->Draw();
TList *listStE = new TList;
listStE->Add(h_lep_E_Stchan_top);
listStE->Add(h_lep_E_Stchan_antitop);
listStE->Add(h_lep_E_Sschan);
listStE->Add(h_lep_E_Swtchan);
h_lep_E_mergedS->Reset();
h_lep_E_mergedS->Merge(listStE);
h_lep_E_mergedS->SetFillColor(kBlue);
h_lep_E_mergedS->Draw();
cstopE->Draw();
TCanvas *ctt = new TCanvas("ctt","ctt",10,10,700,900);
TText tt; tt.SetTextFont(42); tt.SetTextAlign(21);
ctt->Divide(2,2);
ctt->cd(1); h_lep_pt_TThad->Draw();
ctt->cd(2); h_lep_pt_TTlep->Draw();
TList *listt = new TList;
listt->Add(h_lep_pt_TThad);
listt->Add(h_lep_pt_TTlep);
h_lep_pt_mergedTT->Reset();
h_lep_pt_mergedTT->Merge(listt);
h_lep_pt_mergedTT->SetFillColor(kBlue);
h_lep_pt_mergedTT->Draw();
ctt->Draw();
TCanvas *cttEta = new TCanvas("cttEta","cttEta",10,10,700,900);
TText ttEta; ttEta.SetTextFont(42); ttEta.SetTextAlign(21);
cttEta->Divide(2,2);
cttEta->cd(1); h_lep_eta_TThad->Draw();
cttEta->cd(2); h_lep_eta_TTlep->Draw();
TList *listtEta = new TList;
listtEta->Add(h_lep_eta_TThad);
listtEta->Add(h_lep_eta_TTlep);
h_lep_eta_mergedTT->Reset();
h_lep_eta_mergedTT->Merge(listtEta);
h_lep_eta_mergedTT->SetFillColor(kBlue);
h_lep_eta_mergedTT->Draw();
cttEta->Draw();
TCanvas *cttPhi = new TCanvas("cttPhi","cttPhi",10,10,700,900);
TText ttPhi; ttPhi.SetTextFont(42); ttPhi.SetTextAlign(21);
cttPhi->Divide(2,2);
cttPhi->cd(1); h_lep_phi_TThad->Draw();
cttPhi->cd(2); h_lep_phi_TTlep->Draw();
TList *listtPhi = new TList;
listtPhi->Add(h_lep_phi_TThad);
listtPhi->Add(h_lep_phi_TTlep);
h_lep_phi_mergedTT->Reset();
h_lep_phi_mergedTT->Merge(listtPhi);
h_lep_phi_mergedTT->SetFillColor(kBlue);
h_lep_phi_mergedTT->Draw();
cttPhi->Draw();
TCanvas *cttE = new TCanvas("cttE","cttE",10,10,700,900);
TText ttE; ttE.SetTextFont(42); ttE.SetTextAlign(21);
cttE->Divide(2,2);
cttE->cd(1); h_lep_E_TThad->Draw();
cttE->cd(2); h_lep_E_TTlep->Draw();
TList *listtE = new TList;
listtE->Add(h_lep_E_TThad);
listtE->Add(h_lep_E_TTlep);
h_lep_E_mergedTT->Reset();
h_lep_E_mergedTT->Merge(listtE);
h_lep_E_mergedTT->SetFillColor(kBlue);
h_lep_E_mergedTT->Draw();
cttE->Draw();