#
# 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')
