#cuts Fulltag = "(topcandtau2/topcandtau1>0.1&(lepcut2Drel>25.||lepcut2Ddr>0.5)&heavytopcandmass>250.)&(topcandtau3/topcandtau2<0.55&topcandmass<250&topcandmass>140)&isLoose>0." fout = TFile("Zprime_Theta_Feed.root", "UPDATE") # fout.cd() # Plots: mZP15 = TH1F("MU__zp1500", "", 30, 500, 3500) eZP15 = TH1F("EL__zp1500", "", 30, 500, 3500) mZP2 = TH1F("MU__zp2000", "", 30, 500, 3500) eZP2 = TH1F("EL__zp2000", "", 30, 500, 3500) mZP3 = TH1F("MU__zp3000", "", 30, 500, 3500) eZP3 = TH1F("EL__zp3000", "", 30, 500, 3500) # Fill Plots: writeplot(z15FileName, lumi*z15xs/z15n, mZP15, "EventMass", "("+Fulltag+"&isMuon>0.&isGTt>0.)", "(1.0)") writeplot(z15FileName, lumi*z15xs/z15n, eZP15, "EventMass", "("+Fulltag+"&isElec>0.&isGTt>0.)", "(1.0)") writeplot(z2FileName, lumi*z2xs/z2n, mZP2, "EventMass", "("+Fulltag+"&isMuon>0.&isGTt>0.)", "(1.0)") writeplot(z2FileName, lumi*z2xs/z2n, eZP2, "EventMass", "("+Fulltag+"&isElec>0.&isGTt>0.)", "(1.0)") writeplot(z3FileName, lumi*z3xs/z3n, mZP3, "EventMass", "("+Fulltag+"&isMuon>0.&isGTt>0.)", "(1.0)") writeplot(z3FileName, lumi*z3xs/z3n, eZP3, "EventMass", "("+Fulltag+"&isElec>0.&isGTt>0.)", "(1.0)") fout.Write() fout.Save() fout.Close()
edZPe = TH1F("EL__NT", "", 25, 0, 2500) # data: mdZP = TH1F("MU__DATA", "", 25, 0, 2500) edZP = TH1F("EL__DATA", "", 25, 0, 2500) # MC measurements: mtZPm = TH1F("MU__tt", "", 25, 0, 2500) msZPm = TH1F("MU__st", "", 25, 0, 2500) etZPm = TH1F("EL__tt", "", 25, 0, 2500) esZPm = TH1F("EL__st", "", 25, 0, 2500) # Now we fill these all up: #Subtraction: for i in range(len(tFileName)): # ttbar writeplot(tFilePrefix+tFileName[i], lumi*txs[i]/tn[i], etZPs, "EventMass", "("+PreSel+"&"+AntiTag+"&"+El+")", "("+ntW_el+"*"+TW+")") writeplot(tFilePrefix+tFileName[i], lumi*txs[i]/tn[i], mtZPs, "EventMass", "("+PreSel+"&"+AntiTag+"&"+Mu+")", "("+ntW_mu+"*"+TW+")") for i in range(len(sFileName)): # signle top writeplot(sFilePrefix+sFileName[i], lumi*sxs[i]/sn[i], esZPs, "EventMass", "("+PreSel+"&"+AntiTag+"&"+El+")", "("+ntW_el+")") writeplot(sFilePrefix+sFileName[i], lumi*sxs[i]/sn[i], msZPs, "EventMass", "("+PreSel+"&"+AntiTag+"&"+Mu+")", "("+ntW_mu+")") #Estimates: writeplot(dFileNameE, 1.0, edZPe, "EventMass", "("+PreSel+"&"+AntiTag+"&"+El+")", "("+ntW_el+")") writeplot(dFileNameM, 1.0, mdZPe, "EventMass", "("+PreSel+"&"+AntiTag+"&"+Mu+")", "("+ntW_mu+")") #data: writeplot(dFileNameE, 1.0, edZP, "EventMass", "("+PreSel+"&"+TopTag+"&"+El+")", "1.0") writeplot(dFileNameM, 1.0, mdZP, "EventMass", "("+PreSel+"&"+TopTag+"&"+Mu+")", "1.0") # Fill MC: for i in range(len(tFileName)): # All versions of the ttbar writeplot(tFilePrefix+tFileName[i], lumi*txs[i]/tn[i], etZPm, "EventMass", "("+PreSel+"&"+TopTag+"&"+El+")", "("+TW+")") writeplot(tFilePrefix+tFileName[i], lumi*txs[i]/tn[i], mtZPm, "EventMass", "("+PreSel+"&"+TopTag+"&"+Mu+")", "("+TW+")")
mZPn_Nup = TH1F("MU__NT__N__up", "", 30, 500, 3500) eZPn_Nup = TH1F("EL__NT__N__up", "", 30, 500, 3500) mZPt_Ndn = TH1F("MU__TT__N__down", "", 30, 500, 3500) eZPt_Ndn = TH1F("EL__TT__N__down", "", 30, 500, 3500) mZPn_Ndn = TH1F("MU__NT__N__down", "", 30, 500, 3500) eZPn_Ndn = TH1F("EL__NT__N__down", "", 30, 500, 3500) # Errors from MC: # Now we fill them: #cuts Fulltag = "(topcandtau2/topcandtau1>0.1&(lepcut2Drel>25.||lepcut2Ddr>0.5)&heavytopcandmass>250.)&(topcandtau3/topcandtau2<0.55&topcandmass<250&topcandmass>140)&isLoose>0." Antitag = "(topcandtau2/topcandtau1>0.1&(lepcut2Drel>25.||lepcut2Ddr>0.5)&heavytopcandmass>250.)&(topcandtau3/topcandtau2>0.55&topcandmass<250&topcandmass>140)&isLoose>0." # Subtractions: # ttbar: for i in range(len(tFileName)): writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], mtZPs, "EventMass", "("+Antitag+"&isMuon>0.)", "("+ntW+"*"+TW+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], mtZPs_aup, "EventMass", "("+Antitag+"&isMuon>0.)", "("+ntW+"*"+TW_aup+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], mtZPs_adn, "EventMass", "("+Antitag+"&isMuon>0.)", "("+ntW+"*"+TW_adn+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], mtZPs_Nup, "EventMass", "("+Antitag+"&isMuon>0.)", "("+ntW+"*"+TW_Nup+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], mtZPs_Ndn, "EventMass", "("+Antitag+"&isMuon>0.)", "("+ntW+"*"+TW_Ndn+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], mtZPsU, "EventMass", "("+Antitag+"&isMuon>0.)", "("+ntWu+"*"+TW+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], mtZPsD, "EventMass", "("+Antitag+"&isMuon>0.)", "("+ntWd+"*"+TW+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], etZPs, "EventMass", "("+Antitag+"&isElec>0.)", "("+ntW+"*"+TW+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], etZPs_aup, "EventMass", "("+Antitag+"&isElec>0.)", "("+ntW+"*"+TW_aup+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], etZPs_adn, "EventMass", "("+Antitag+"&isElec>0.)", "("+ntW+"*"+TW_adn+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], etZPs_Nup, "EventMass", "("+Antitag+"&isElec>0.)", "("+ntW+"*"+TW_Nup+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], etZPs_Ndn, "EventMass", "("+Antitag+"&isElec>0.)", "("+ntW+"*"+TW_Ndn+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], etZPsU, "EventMass", "("+Antitag+"&isElec>0.)", "("+ntWu+"*"+TW+")") writeplot(tFilePrefix+tFileName[i]+'.root', lumi*txs[i]/tn[i], etZPsD, "EventMass", "("+Antitag+"&isElec>0.)", "("+ntWd+"*"+TW+")") # single top for i in range(len(sFileName)):