#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)):
