def draw_canvas_1d(name): f_mc = TFile("plots/mc/%s_mcc7.root" % name) print(f_mc.ls()) print(name) h_mc = gDirectory.Get("h_%s_reco" % name) h_mc.SetName("h_%s_reco" % name) h_mc.SetLineColor(kRed + 1) h_mc.SetFillColor(kRed + 1) h_mc.SetFillStyle(0) h_mc.GetYaxis().SetTitle("Efficiency") h_mc.GetYaxis().SetRangeUser(0.61, 1.04) f = TFile("plots/data/e_%s_pandoraCosmic.root" % name) h = gDirectory.Get("h_%s" % name) h.SetMarkerStyle(20) h.SetLineColor(1) f_sys = TFile("plots/data/e_%s_sys_pandoraCosmic.root" % name) h_sys = gDirectory.Get("h_%s_sys" % name) h_ratio = h_sys.Clone() h_ratio.Divide(h_mc) x_minbin = h_mc.FindFirstBinAbove() low = h_mc.GetXaxis().GetBinLowEdge(x_minbin) x_maxbin = h_mc.FindLastBinAbove() high = h_mc.GetXaxis().GetBinUpEdge(x_maxbin) h_mc.GetXaxis().SetRangeUser(low, high) h_mc.GetYaxis().SetTitleSize(0.08) h_mc.GetYaxis().SetLabelSize(0.065) h_mc.GetYaxis().SetTitleOffset(0.59) h_ratio.GetXaxis().SetRangeUser(low, high) g_stat = histo2graph(h) g_sys = histo2graph(h_sys) g_sys.SetMarkerStyle(20) g_sys.SetLineWidth(2) g_stat.SetFillStyle(3001) g_stat.SetFillColor(kGray + 2) g_stat.SetLineWidth(2) g_stat.SetLineColor(kGray + 2) pt = TPaveText(0.09, 0.855, 0.32, 0.98, "ndc") pt.AddText("MicroBooNE") pt.SetFillColor(0) pt.SetBorderSize(0) pt.SetShadowColor(0) leg = TLegend(0.55, 0.09, 0.86, 0.4) leg.AddEntry(g_sys, "Data - stat. #oplus sys.", "ep") leg.AddEntry(g_stat, "Data - stat. only", "f") leg.AddEntry(h_mc, "Monte Carlo", "f") canvas = draw_ratio("c_%s" % name, [h_mc, g_stat, g_sys, pt, leg], [h_ratio], ["hist", "2", "p", "", ""], ["ep"]) canvas.SaveAs("plots/%s.pdf" % name) canvas.SaveAs("plots/%s.C" % name)
def draw_canvas_2d(name): f_mc = TFile("plots/mc/%s_mcc7.root" % name) h_mc = gDirectory.Get("h_%s_reco" % name) f = TFile("plots/data/e_%s_pandoraCosmic.root" % name) h = gDirectory.Get("h_%s" % name) h.Divide(h_mc) h.GetZaxis().SetRangeUser(0.5, 1.5) #h.GetZaxis().SetTitle("Data/Monte Carlo") #h.GetZaxis().RotateTitle() h.SetMarkerSize(2) x_minbin = h.FindFirstBinAbove(0, 1) low_x = h.GetXaxis().GetBinLowEdge(x_minbin) x_maxbin = h.FindLastBinAbove(0, 1) high_x = h.GetXaxis().GetBinUpEdge(x_maxbin) y_minbin = h.FindFirstBinAbove(0, 2) low_y = h.GetYaxis().GetBinLowEdge(y_minbin) y_maxbin = h.FindLastBinAbove(0, 2) high_y = h.GetYaxis().GetBinUpEdge(y_maxbin) h.GetXaxis().SetRangeUser(low_x, high_x) h.GetYaxis().SetRangeUser(low_y, high_y) c = TCanvas("c_%s" % name) h.SetMarkerSize(2.5) gPad.SetBottomMargin(0.17) gPad.SetLeftMargin(0.13) gPad.SetTopMargin(0.15) gPad.SetRightMargin(0.15) h.GetYaxis().SetTitleSize(0.07) h.GetXaxis().SetTitleSize(0.07) h.GetYaxis().SetTitleOffset(0.8) pt = TPaveText(0.13, 0.855, 0.42, 0.98, "ndc") pt.AddText("MicroBooNE") pt.SetFillColor(0) pt.SetBorderSize(0) pt.SetShadowColor(0) h.Draw("colz texte") pt.Draw() c.Update() c.SaveAs("plots/%s.pdf" % name) c.SaveAs("plots/%s.p" % name) return c
g_reco_eff_data.SetMarkerStyle(30) g_reco_eff_data.SetMarkerSize(1.5) g_reco_eff_data.SetLineStyle(2) l = TLegend(0.46, 0.16, 0.85, 0.425) l.SetBorderSize(0) l.SetShadowColor(0) l.AddEntry(g_purity, "Purity #it{P}", "lp") l.AddEntry(g_acceptance, "Acceptance #it{A}", "lp") l.AddEntry(g_eff, "MuCS-MC tagging efficiency #it{#epsilon}_{tag}^{MuCS-MC}", "lp") l.AddEntry(g_data_eff, "Data tagging efficiency #it{#epsilon}_{tag}^{data}", "lp") l.AddEntry(g_reco_eff, "MuCS-MC reconstruction efficiency #it{#epsilon}_{MuCS-MC}", "lp") l.AddEntry(g_reco_eff_data, "Data reconstruction efficiency #it{#epsilon}_{data}", "lp") l.Draw() pt = TPaveText(0.09, 0.91, 0.34, 0.98, "ndc") pt.AddText("MicroBooNE") pt.SetFillColor(0) pt.SetBorderSize(0) pt.SetShadowColor(0) pt.Draw() c.Update() input()
h_theta_res_xy_corr.SetLineColor(kBlack) h_theta_res_xy.SetMarkerStyle(4) h_theta_res_xy.SetLineColor(kBlack) print("no_corr", h_theta_res_xy.GetMean() - h_theta_res_xy_mc.GetMean()) print("corr", h_theta_res_xy_corr.GetMean() - h_theta_res_xy_mc.GetMean()) l.AddEntry(h_theta_res_xy_mc, "Monte Carlo", "l") l.AddEntry(h_theta_res_xy, "Data", "ep") l.AddEntry(h_theta_res_xy_corr, "Data - SCE corrected", "ep") l.SetTextFont(43) l.SetShadowColor(0) l.SetBorderSize(0) l.SetTextSize(14) l.Draw() preliminary2 = TPaveText(0.13, 0.78, 0.46, 0.86, "NDC") preliminary2.SetShadowColor(0) preliminary2.SetBorderSize(0) preliminary2.SetFillStyle(0) preliminary2.AddText("MicroBooNE Preliminary") preliminary2.Draw() c_res_xy.Update() c_res_xy.SaveAs(output + "/xy_" + algo + "_res.pdf") c_res_yz = TCanvas("c_res_yz") f_eff_yz = TFile("yz_mc_pandoraCosmic_eff.root") eff_yz = gDirectory.Get("e_yz_mc") h_theta_res_yz_corr.Scale(h_theta_res_yz_mc.Integral() / h_theta_res_yz.Integral()) h_theta_res_yz.Scale(h_theta_res_yz_mc.Integral() / h_theta_res_yz.Integral()) h_theta_res_yz_mc.SetLineColor(kRed + 1) h_theta_res_yz_mc.Draw("hist")