예제 #1
0


    histolist = [resprereg,respostreg,outputdif_full,outputdif_zoom]



    tree.Project("resprereg","(Jet_MatchedPartonPt - Jet_Pt)/Jet_Pt:Jet_Pt","","") 
    tree.Project("respostreg","(Jet_MatchedPartonPt - ("+regpt+"))/("+regpt+"):("+regpt+")","","") 
    tree.Project("outputdif_full","( BDTG - "+target+") : "+target,"","")
    tree.Project("outputdif_zoom","( BDTG - "+target+") : "+target,"","")



    pdfout = PDFPrinting(outputname)
    outputfile = ROOT.TFile(outputname+".root","RECREATE")
    outputfile.cd()
    c1 = ROOT.TCanvas()    
    c1.cd()


    for histo in histolist:
        histo.Draw("colz")
        c1.Update
        pdfout.addCanvastoPDF(c1)


    pdfout.closePDF()

    del histolist,resprereg,respostreg,outputdif_full,outputdif_zoom,inputfile,tree,pdfout
        inputbinning[var0][0],
        inputbinning[var0][1],
        inputbinning[var0][2],
        inputbinning[var1][0],
        inputbinning[var1][1],
        inputbinning[var1][2],
    )

    tree_data.Project(str(rannum_data), var1 + ":" + var0, "")
    tree_ttbar.Project(str(rannum_ttbar), var1 + ":" + var0, mcweight)

    htmp_data.GetXaxis().SetTitle(var0)
    htmp_data.GetYaxis().SetTitle(var1)
    htmp_data.Draw("colz")
    c1.Update()
    pdfout_data.addCanvastoPDF(c1)
    outputfile_data.cd()
    htmp_data.Write()

    htmp_ttbar.GetXaxis().SetTitle(var0)
    htmp_ttbar.GetYaxis().SetTitle(var1)
    htmp_ttbar.Draw("colz")
    c1.Update()
    pdfout_ttbar.addCanvastoPDF(c1)
    outputfile_ttbar.cd()
    htmp_ttbar.Write()

    htmp_data.SetMarkerColor(ROOT.kRed)
    htmp_data.Draw()
    htmp_ttbar.Draw("same")
    c1.Update()
예제 #3
0
                               binning_outputdiv_zoom_div[i][0],
                               binning_outputdiv_zoom_div[i][1],
                               binning_outputdiv_zoom_div[i][2])

    histolist = [resprereg, respostreg, outputdif_full, outputdif_zoom]

    tree.Project("resprereg", "(Jet_MatchedPartonPt - Jet_Pt)/Jet_Pt:Jet_Pt",
                 "", "")
    tree.Project(
        "respostreg", "(Jet_MatchedPartonPt - (" + regpt + "))/(" + regpt +
        "):(" + regpt + ")", "", "")
    tree.Project("outputdif_full", "( BDTG - " + target + ") : " + target, "",
                 "")
    tree.Project("outputdif_zoom", "( BDTG - " + target + ") : " + target, "",
                 "")

    pdfout = PDFPrinting(outputname)
    outputfile = ROOT.TFile(outputname + ".root", "RECREATE")
    outputfile.cd()
    c1 = ROOT.TCanvas()
    c1.cd()

    for histo in histolist:
        histo.Draw("colz")
        c1.Update
        pdfout.addCanvastoPDF(c1)

    pdfout.closePDF()

    del histolist, resprereg, respostreg, outputdif_full, outputdif_zoom, inputfile, tree, pdfout
                          inputbinning[var0][2], inputbinning[var1][0],
                          inputbinning[var1][1], inputbinning[var1][2])
    htmp_ttbar = ROOT.TH2F(str(rannum_ttbar),
                           var0 + "_" + var1 + str(rannum_ttbar),
                           inputbinning[var0][0], inputbinning[var0][1],
                           inputbinning[var0][2], inputbinning[var1][0],
                           inputbinning[var1][1], inputbinning[var1][2])

    tree_data.Project(str(rannum_data), var1 + ":" + var0, "")
    tree_ttbar.Project(str(rannum_ttbar), var1 + ":" + var0, mcweight)

    htmp_data.GetXaxis().SetTitle(var0)
    htmp_data.GetYaxis().SetTitle(var1)
    htmp_data.Draw("colz")
    c1.Update()
    pdfout_data.addCanvastoPDF(c1)
    outputfile_data.cd()
    htmp_data.Write()

    htmp_ttbar.GetXaxis().SetTitle(var0)
    htmp_ttbar.GetYaxis().SetTitle(var1)
    htmp_ttbar.Draw("colz")
    c1.Update()
    pdfout_ttbar.addCanvastoPDF(c1)
    outputfile_ttbar.cd()
    htmp_ttbar.Write()

    htmp_data.SetMarkerColor(ROOT.kRed)
    htmp_data.Draw()
    htmp_ttbar.Draw("same")
    c1.Update()