getHiggsMasswcorr(higgsjets[0][0], higgsjets[0][1],
higgsjets[0][2], higgsjets[0][3],
higgsjets[0][4], higgsjets[1][0],
higgsjets[1][1], higgsjets[1][2],
higgsjets[1][3], higgsjets[1][4]), iw + 1)
if iw == 0:
bbmass.FillnormHisto(bbmass_noreg[0], 0)
#else:
# if bbmass_fromTree[00] > 0:
# print 0,"|", bbmass_fromTree[0]
errors.printerr()
errvals = errors.getvalues()
plot_errabs.addPoint(errvals["err_abs_loss"])
plot_errsqu.addPoint(errvals["err_squared_loss"])
targetvoutput.WriteTwoDplot(c1, "ttHbb", True, None, pdfout, False, False,
[0.6, 0.7, 1.6, 1.5])
del tree, regression, errors, targetvoutput
regcorr.WriteHisto(c1, "ttHbb", False, False, pdfout)
regcorr_parton.WriteHisto(c1, "ttHbb", False, False, pdfout)
regcorr_genjet.WriteHisto(c1, "ttHbb", False, False, pdfout)
regcorr_targets.WriteHisto(c1, "ttHbb", False, False, pdfout)
bbmass.WriteHisto(c1, "ttHbb", False, False, pdfout, True)
plot_errabs.WritePointPlot(c1, "ttHbb", pdfout)
plot_errsqu.WritePointPlot(c1, "ttHbb", pdfout)
pdfout.closePDF()
ratio_nEnergyFractions.projecttoHisto(isel, tree, "RegJet_nHEFrac / RegJet_nEmEFrac",selection)
ratio_JetvParton_b_regVnoreg.projecttoHisto(1,tree,"RegJet_Pt / RegJet_MatchedPartonPt",samplesel+"abs(Jet_Flav) == 5")
ratio_JetvParton_b_regVnoreg.projecttoHisto(0,tree,"RegJet_preregPt / RegJet_MatchedPartonPt",samplesel+"abs(Jet_Flav) == 5")
for plot in plotlist:
plot.WriteHisto(c1,sampletype[isample],False,True,pdfout)
ratio_JetvParton.fitGauss(0,0.7,1.2)
ratio_JetvParton.fitGauss(1,0.7,1.25)
ratio_JetvParton.setmanualegendsize("right",0.5,0.58,0.90,0.88)
ratio_JetvParton.WriteHisto(c1,sampletype[isample],False,True,pdfout, False, None, False,False,None,True,[1,2,ROOT.kBlack,2])
ratio_JetvPartonpostreg.fitGauss(0,0.7,1.25)
ratio_JetvPartonpostreg.fitGauss(1,0.7,1.25)
ratio_JetvPartonpostreg.setmanualegendsize("right",0.5,0.58,0.90,0.88)
ratio_JetvPartonpostreg.WriteHisto(c1,sampletype[isample],False,True,pdfout, False, None, False,False,None,True,[1,2,ROOT.kBlack,2])
ratio_JetvParton_b_regVnoreg.WriteHisto(c1,sampletype[isample],False,True,pdfout, False, None, False,False,None,False,[1,2,ROOT.kBlack,2])
ratio_JetvParton_b_regVnoreg.fitGauss(0,0.7,1.2)
ratio_JetvParton_b_regVnoreg.fitGauss(1,0.7,1.25)
ratio_JetvParton_b_regVnoreg.setmanualegendsize("right",0.5,0.58,0.90,0.88)
ratio_JetvParton_b_regVnoreg.WriteHisto(c1,sampletype[isample],False,True,pdfout, False, None, False,False,None,True,[1,2,ROOT.kBlack,2])
pdfout.closePDF()
del targetvsBDTGHnames
del targetvsBDTGregcorrHnames, targetvsBDTGPTJetPartonHnames,filenames,names
#raw_input("Press ret")
"""
pdfout = PDFPrinting("outname_configcomp")
outputfile = ROOT.TFile("outname_configcomp"+".root","RECREATE")
outputfile.cd()
c1 = ROOT.TCanvas()
for plots in compplots:
for plot in plots:
plot.WriteHisto(c1,"ttbar",False,False,pdfout)
pdfout.closePDF()
"""
pdfout_norm = PDFPrinting(outputname)
outputfile_norm = ROOT.TFile(outputname+".root","RECREATE")
outputfile_norm.cd()
c1_norm = ROOT.TCanvas()
for plots in compplots:
for plot in plots:
plot.WriteHisto(c1_norm,"ttbar",False,True,pdfout_norm)
pdfout_norm.closePDF()
del targetvsBDTGregcorrH
del targetvsBDTGPTJetPartonH
del inputfile, tree
del targetvsBDTGHnames
del targetvsBDTGregcorrHnames, targetvsBDTGPTJetPartonHnames, filenames, names
#raw_input("Press ret")
"""
pdfout = PDFPrinting("outname_configcomp")
outputfile = ROOT.TFile("outname_configcomp"+".root","RECREATE")
outputfile.cd()
c1 = ROOT.TCanvas()
for plots in compplots:
for plot in plots:
plot.WriteHisto(c1,"ttbar",False,False,pdfout)
pdfout.closePDF()
"""
pdfout_norm = PDFPrinting(outputname)
outputfile_norm = ROOT.TFile(outputname + ".root", "RECREATE")
outputfile_norm.cd()
c1_norm = ROOT.TCanvas()
for plots in compplots:
for plot in plots:
plot.WriteHisto(c1_norm, "ttbar", False, True, pdfout_norm)
pdfout_norm.closePDF()
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()
pdfout_scatter.addCanvastoPDF(c1)
del htmp_data, htmp_ttbar
pdfout_data.closePDF()
pdfout_ttbar.closePDF()
pdfout_scatter.closePDF()
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()
pdfout_scatter.addCanvastoPDF(c1)
del htmp_data, htmp_ttbar
pdfout_data.closePDF()
pdfout_ttbar.closePDF()
pdfout_scatter.closePDF()
