yp.colorDict[samp + "_Rcs_SB"] = yp.kRed hRcsSB = yp.makeSampHisto(yds, samp, "Rcs_SB", samp + "_Rcs_SB") #hRcsSB.SetTitle("R_{CS} (SB)") hRcsSB.SetTitle( "R_{CS}^{MC}(n_{b-tag}, n_{jet}#in [4,5], EWK)") # hRcsSB.SetTitle("R_{CS}^{MC}(n_{b-tag}, n_{jet}#in [4], EWK)") # Kappa yp.colorDict[samp + "_Kappa"] = yp.kBlack hKappa = yp.makeSampHisto(yds, samp, "Kappa", samp + "_Kappa") hKappa.SetTitle("#kappa_{MC}") yp.prepKappaHist(hKappa) yp.prepRatio(hKappa) # remove marker #hKappa.SetMarkerStyle(0) hKappa.SetLineWidth(2) #hKappa.SetLineColor(yp.kViolet-5) hRcsMB.GetYaxis().SetRangeUser(0, 0.14) width = 1200 height = 600 legPos = "TM" if doSquare == True: width = 600 height = 600 legPos = "Square"
var = "SR_MB" # canvs canvs = [] allhists = [] for mass in masses: hists = [] # read in central value signame = samp + "_" + mass hCentral = yp.makeSampHisto(yds, signame, var, signame + "_central") #print hCentral yp.prepRatio(hCentral, True) #hCentral.GetYaxis().SetRangeUser() ## Define flat values (i.e. for PU) flats = {} #flats["PU"] = 0.05 ### Add flat systs (lumi, lepSF, triggEff, etc.) #flats["stat"] = 0.1 flats["lepSF"] = 0.02 flats["lumi"] = 0.026 flats["iso"] = 0.04 flats["trig"] = 0.02 flats["PU"] = 0.1 #flats = {}
#samp = "EWK" samps = ["EWK", "TTJets", "WJets", "SingleTop", "DY", "TTV"] #samps = ['T_tWch','TToLeptons_tch','TBar_tWch', 'TToLeptons_sch',"EWK"] samp = samps[0] var = "Kappa" #var = "SR_SB" # canvs and hists hists = [] canvs = [] # read in central value hCentral = yp.makeSampHisto(yds, samp, var) print "hCentral done" yp.prepRatio(hCentral) for i, syst in enumerate(systs): yp.colorDict[syst + "_syst"] = sysCols[i] sname = samp + "_" + syst + "_syst" print "Making hist for", sname hist = yp.makeSampHisto(yds, sname, var, syst + "_syst") if syst in systNames: hist.SetTitle(systNames[syst]) else: hist.SetTitle(syst) hist.GetYaxis().SetTitle("Relative uncertainty") hist.GetYaxis().SetTitleSize(0.04) hist.GetYaxis().SetTitleOffset(0.8)
yp.colorDict[samp+"_Rcs_MB"] = yp.kBlue hRcsMB = yp.makeSampHisto(yds,samp,"Rcs_MB",samp+"_Rcs_MB") hRcsMB.SetTitle("R_{CS} (MB)") # RCS SB yp.colorDict[samp+"_Rcs_SB"] = yp.kRed hRcsSB = yp.makeSampHisto(yds,samp,"Rcs_SB",samp+"_Rcs_SB") hRcsSB.SetTitle("R_{CS} (SB)") # Kappa yp.colorDict[samp+"_Kappa"] = yp.kBlack hKappa = yp.makeSampHisto(yds,samp,"Kappa",samp+"_Kappa") hKappa.SetTitle("#kappa") yp.prepKappaHist(hKappa) yp.prepRatio(hKappa) canv = yp.plotHists(samp+"_RcsKappa_",[hRcsMB,hRcsSB],hKappa) canvs.append(canv) if not yp._batchMode: raw_input("Enter any key to exit") # Save canvases exts = [".pdf",".png"] #exts = [".pdf"] odir = "BinPlots/Rcs/test/" for canv in canvs: for ext in exts:
# Sample and variable samp = "EWK" #samps = ["TTJets","WJets","SingleTop","DY","TTV"] #samps = ['T_tWch','TToLeptons_tch','TBar_tWch', 'TToLeptons_sch',"EWK"] #samp = samps[4] var = "Kappa" #var = "SR_SB" # canvs and hists hists = [] canvs = [] # read in central value hCentral = yp.makeSampHisto(yds,samp,var) yp.prepRatio(hCentral) for i,syst in enumerate(systs): yp.colorDict[syst+"_syst"] = sysCols[i] sname = samp+"_"+syst+"_syst" print "Making hist for", sname hist = yp.makeSampHisto(yds,sname,var,syst+"_syst") if syst in systNames: hist.SetTitle(systNames[syst]) else: hist.SetTitle(syst) hist.GetYaxis().SetTitle("Relative uncertainty") hist.GetYaxis().SetTitleSize(0.04) hist.GetYaxis().SetTitleOffset(0.8)