%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;iGetEntry(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); } } } } } } //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();