-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathplotFR.py
More file actions
83 lines (67 loc) · 2.65 KB
/
Copy pathplotFR.py
File metadata and controls
83 lines (67 loc) · 2.65 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
import ROOT,os
# ______________________________________________________________________________ ||
#inputElPath = "Data/FakeRate_Vukasin_190326/FR_WZremoved_17.root"
#inputMuPath = "Data/FakeRate_Vukasin_190326/FR_WZremoved_17.root"
#
#histNameMuEB = "Data_FRmu_EB"
#histNameMuEE = "Data_FRmu_EE"
#histNameElEB = "Data_FRel_EB"
#histNameElEE = "Data_FRel_EE"
#
#outputDir = "/home/lucien/public_html/Higgs/DarkZ/FakeRate/20190327_plotFR_2017/"
#factor = 1.2
#legPos = [0.1,0.75,0.35,0.9]
#tag = "_2017"
inputElPath = "Data/FakeRate_Vukasin_190301/fakeRates_el_2016.root"
inputMuPath = "Data/FakeRate_Vukasin_190301/fakeRates_mu_2016.root"
histNameMuEB = "h1D_FRmu_EB"
histNameMuEE = "h1D_FRmu_EE"
histNameElEB = "h1D_FRel_EB"
histNameElEE = "h1D_FRel_EE"
outputDir = "/home/lucien/public_html/Higgs/DarkZ/FakeRate/20190327_plotFR_2016/"
factor = 1.2
legPos = [0.1,0.75,0.35,0.9]
tag = "_2016"
# ______________________________________________________________________________ ||
ROOT.gROOT.SetBatch(ROOT.kTRUE)
# ______________________________________________________________________________ ||
inputElFile = ROOT.TFile(inputElPath)
inputMuFile = ROOT.TFile(inputMuPath)
histElEB = inputElFile.Get(histNameElEB)
histElEE = inputElFile.Get(histNameElEE)
histMuEB = inputMuFile.Get(histNameMuEB)
histMuEE = inputMuFile.Get(histNameMuEE)
# ______________________________________________________________________________ ||
histMuEB.SetLineColor(ROOT.kRed)
histMuEE.SetLineColor(ROOT.kBlue)
histMuEB.SetStats(0)
histMuEE.SetStats(0)
histMuEB.SetTitle("Muon Fake Rate")
histMuEE.SetTitle("Muon Fake Rate")
maximum = max([histMuEB.GetMaximum(),histMuEE.GetMaximum(),])
histMuEB.GetYaxis().SetRangeUser(0.,1.2*maximum)
c = ROOT.TCanvas()
leg = ROOT.TLegend(*legPos)
leg.AddEntry(histMuEB,"Barrel")
leg.AddEntry(histMuEE,"Endcap")
histMuEB.Draw()
histMuEE.Draw("same")
leg.Draw("same")
c.SaveAs(outputDir+"FR_mu%s.pdf"%tag)
# ______________________________________________________________________________ ||
histElEB.SetLineColor(ROOT.kRed)
histElEE.SetLineColor(ROOT.kBlue)
histElEB.SetStats(0)
histElEE.SetStats(0)
histElEB.SetTitle("Electron Fake Rate")
histElEE.SetTitle("Electron Fake Rate")
maximum = max([histElEB.GetMaximum(),histElEE.GetMaximum(),])
histElEB.GetYaxis().SetRangeUser(0.,1.2*maximum)
c = ROOT.TCanvas()
leg = ROOT.TLegend(*legPos)
leg.AddEntry(histElEB,"Barrel")
leg.AddEntry(histElEE,"Endcap")
histElEB.Draw()
histElEE.Draw("same")
leg.Draw("same")
c.SaveAs(outputDir+"FR_el%s.pdf"%tag)