Beispiel #1
0
                    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"
Beispiel #2
0
    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 = {}
Beispiel #3
0
    #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:
Beispiel #5
0
    # 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)