# # raise Exception # myeff = ROOT.TEfficiencyBugFixed(pass_histo, all_histo) myeff.SetStatisticOption(0) # 0 means ClopperPearson ROOT.SetOwnership(myeff, False) efficiency = parse_formula(args.efficiency, args.parameters) efficiency.SetName('efficiency') efficiency.SetTitle('efficiency') ROOT.SetOwnership(efficiency, False) myeff.Fit(efficiency, "LMI") # fitta con l'esponenziale ma crasha if args.plot: canvas = ROOT.TCanvas("asdf", "asdf", 800, 600) canvas.SetLogz(True) myeff.Draw("LEGO") efficiency.Draw("SURFSAME") # import pdb; pdb.set_trace() canvas.SaveAs(args.output) plot_name = args.output.replace('.root', '.png') # log.info("Saving fit plot in %s", plot_name) canvas.SaveAs(plot_name) canvas.SaveAs(plot_name.replace('.png', '.pdf')) canvas.SetLogz(False) canvas.SetLogy(True) graph_proj_x = asrootpy( myeff.Projection(ROOT.TEfficiencyBugFixed.xaxis)) graph_proj_x.SetMarkerStyle(20)
file_dataA = ROOT.TFile('results/newNtuple_5Nov/LFVHETauAnalyzerMVA/ZetauEmbedded_Run2012A.root') file_dataB = ROOT.TFile('results/newNtuple_5Nov/LFVHETauAnalyzerMVA/ZetauEmbedded_Run2012B.root') file_dataC = ROOT.TFile('results/newNtuple_5Nov/LFVHETauAnalyzerMVA/ZetauEmbedded_Run2012C.root') file_dataD = ROOT.TFile('results/newNtuple_5Nov/LFVHETauAnalyzerMVA/ZetauEmbedded_Run2012D.root') dataA=file_dataA.Get('os/gg/ept30/h_collmass_pfmet') dataB=file_dataB.Get('os/gg/ept30/h_collmass_pfmet') dataC=file_dataC.Get('os/gg/ept30/h_collmass_pfmet') dataD=file_dataD.Get('os/gg/ept30/h_collmass_pfmet') data=dataC.Clone() data.Add(dataB) data.Add(dataA) data.Add(dataD) c= ROOT.TCanvas("c","c", 800, 1000) c.Draw() c.SetGridx(1) c.SetGridy(1) njets=[0,1,2,3,4] file_MC0 = ROOT.TFile('results/newNtuple_5Nov/LFVHETauAnalyzerMVA/Z0jets_M50_skimmedTT.root') file_MC1 = ROOT.TFile('results/newNtuple_5Nov/LFVHETauAnalyzerMVA/Z1jets_M50_skimmedTT.root') file_MC2 = ROOT.TFile('results/newNtuple_5Nov/LFVHETauAnalyzerMVA/Z2jets_M50_skimmedTT.root') file_MC3 = ROOT.TFile('results/newNtuple_5Nov/LFVHETauAnalyzerMVA/Z3jets_M50_skimmedTT.root') file_MC4 = ROOT.TFile('results/newNtuple_5Nov/LFVHETauAnalyzerMVA/Z4jets_M50_skimmedTT.root') mc0 = file_MC0.Get('os/gg/ept30/h_collmass_pfmet') mc1 = file_MC1.Get('os/gg/ept30/h_collmass_pfmet')