def draw():
ps = plot_saver(plot_dir('o2t_templates_run2'), size=(600,600))
f = ROOT.TFile(gp.fn)
whiches = [
('multijet', [
(name2isample(f, 'mfv_neu_tau00100um_M0800'), ROOT.kRed, 'multijet M = 800 GeV, #tau = 100 #mum'),
(name2isample(f, 'mfv_neu_tau00300um_M0800'), ROOT.kGreen+2, 'multijet M = 800 GeV, #tau = 300 #mum'),
(name2isample(f, 'mfv_neu_tau01000um_M0800'), ROOT.kBlue, 'multijet M = 800 GeV, #tau = 1 mm'),
(name2isample(f, 'mfv_neu_tau10000um_M0800'), ROOT.kMagenta, 'multijet M = 800 GeV, #tau = 10 mm'),
(name2isample(f, 'mfv_neu_tau30000um_M0800'), ROOT.kOrange+2, 'multijet M = 800 GeV, #tau = 30 mm'),
]),
('dijet', [
(name2isample(f, 'mfv_ddbar_tau00100um_M0800'), ROOT.kRed, 'dijet M = 800 GeV, #tau = 100 #mum'),
(name2isample(f, 'mfv_ddbar_tau00300um_M0800'), ROOT.kGreen+2, 'dijet M = 800 GeV, #tau = 300 #mum'),
(name2isample(f, 'mfv_ddbar_tau01000um_M0800'), ROOT.kBlue, 'dijet M = 800 GeV, #tau = 1 mm'),
(name2isample(f, 'mfv_ddbar_tau10000um_M0800'), ROOT.kMagenta, 'dijet M = 800 GeV, #tau = 10 mm'),
(name2isample(f, 'mfv_ddbar_tau30000um_M0800'), ROOT.kOrange+2, 'dijet M = 800 GeV, #tau = 30 mm'),
]),
]
def fmt(h, name, color, save=[]):
binning = to_array(0., 0.04, 0.07, 0.15)
h = h.Rebin(len(binning)-1, name, binning)
h.Sumw2()
h.SetStats(0)
h.SetLineWidth(3)
h.SetLineColor(color)
h.SetTitle(';d_{VV} (cm);event rate (unit norm.)')
move_overflow_into_last_bin(h)
h.Scale(1./h.Integral(0,h.GetNbinsX()+2))
save.append(h)
return h
for which_name, which in whiches:
hbkg = fmt(f.Get('h_bkg_dvv'), 'bkg', ROOT.kBlack)
hbkg.Draw('hist e')
hbkg.GetYaxis().SetRangeUser(0,1.5)
leg = ROOT.TLegend(0.142, 0.657, 0.702, 0.857)
leg.SetBorderSize(0)
leg.AddEntry(hbkg, 'Background template d_{VV}^{C}', 'LE')
for isample, color, title in which:
h = fmt(f.Get('h_signal_%i_dvv' % isample), title, color)
h.Draw('hist e same')
leg.AddEntry(h, title, 'LE')
leg.Draw()
ps.save(which_name, log=False)
return 6
elif mass == 1200:
return 7
elif mass == 1600:
return 8
elif mass == 3000:
return 9
matched = []
not_matched = []
x = []
y = []
ex = []
ey = []
gs = []
l1 = ROOT.TLegend(0.75,0.1,0.95,0.5)
l2 = ROOT.TLegend(0.75,0.5,0.95,0.9)
for sample in samples:
if ctau not in sample.name:
continue
fn = os.path.join(os.path.expanduser(file_path), '%s.root' % sample.name)
if not os.path.isfile(fn):
print colors.red('%s is missing' % fn) + '\n'
continue
print sample.name
f = ROOT.TFile(fn)
nrec = f.Get('mfvTheoristRecipe%s/h_gen%s_dvv' % (rec_den, match)).GetEntries()
ngen = f.Get('mfvGen%s/h_gen%s_dvv' % (gen_den, match)).GetEntries()