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)