leg1.AddEntry(
profiles[4], "400-500 GeV, slope = %.3f #pm %.3f" %
(frRes4.Parameter(1), frRes4.ParError(1)), "LEP")
profiles[5].Draw("same")
leg1.AddEntry(
profiles[5], "1000-1100 GeV, slope = %.3f #pm %.3f" %
(frRes5.Parameter(1), frRes5.ParError(1)), "LEP")
profiles[7].Draw("same")
leg1.AddEntry(
profiles[7], "All jets, slope = %.4f #pm %.4f" %
(frRes7.Parameter(1), frRes7.ParError(1)), "LEP")
leg1.Draw("")
cmslabel_sim(canvas1, '2016', 11)
canvas1.Update()
canvas1.SaveAs("rho_" + file.replace(".root", ".pdf"))
sleep(10)
QCDPtPlotters = getPlotters("QCD_Pt_", False)
QCD = MergedPlotter(QCDPtPlotters)
# QCDherwigPlotters = getPlotters('QCD_Pt-',False)
# QCDherwig = MergedPlotter(QCDherwigPlotters)
QCD.setLineProperties(1, 434, 2)
QCD.setFillProperties(1001, 0)
file1 = doDDT(QCD, cuts['common'], "rho_pythia.root")
file2 = doDDT(QCD, cuts['commonTOT'], "rho_pythia_FullSel.root")
doFit("rho_pythia.root", False)
doFit("rho_pythia_FullSel.root", True)
QCDherwigPlotters = getPlotters('QCD_Pt-', False)
QCD = MergedPlotter(QCDPtPlotters)
QCDht = MergedPlotter(QCDhtPlotters)
QCDherwig = MergedPlotter(QCDherwigPlotters)
#
# # DATAPlotters = getPlotters('SingleJet',True)
# # data=MergedPlotter(DATAPlotters)
#
# SigPlotters = getPlotters('QstarToQW_2000',False)
# sig = MergedPlotter(SigPlotters)
#
#
QCD.setFillProperties(1001, 921)
QCDht.setLineProperties(1, 633, 2)
QCDht.setFillProperties(1001, 0)
#
QCDherwig.setLineProperties(1, 434, 2)
QCDherwig.setFillProperties(1001, 0)
#
# sig.setFillProperties(3001,ROOT.kGreen+2)
#Stack
# jjStack = StackPlotter()
# jjStack.addPlotter(QCD,"QCD","QCD multijet","background")
# jjStack.addPlotter(data,"data_obs","Data","data")
for p in purities:
canvs = []
print "Plotting variables for category %s" % (p)
