示例#1
0
    hsd[float(lxy)] = h
hs = [hsd[k] for k in sorted(hsd)]
leg = Legend(hs,
             pad=c,
             leftmargin=0.5,
             margin=0.1,
             entryheight=0.02,
             textsize=12)
draw(hs, pad=c, xlimits=(1e-1, 500), logy=True, logx=True)
leg.Draw()
t = LuminosityLabel('XX#rightarrow2A#rightarrow4#mu ({}, {}) GeV'.format(
    mxx, ma))
t.draw()
c.SaveAs(os.path.join(outdir, 'ch4mu_dplxy.pdf'.format(outdir)))

c.clear()

massdir = getattr(f.ch2mu2e, masstag)
lxytags = [x.name for x in massdir.keys()]  # lxy-0p3_ctau-0p0025

hsd = {}
for tag in lxytags:
    lxy_, ctau_ = tag.split('_')
    lxy = lxy_.split('-')[1].replace('p', '.')
    ctau = ctau_.split('-')[1].replace('p', '.')
    h = getattr(massdir, tag).dplxy
    h.linewidth = 2
    h.drawstyle = 'PLC hist'
    h.scale(1. / h.integral())
    h.legendstyle = 'L'
    h.title = 'lxy: {}cm c#tau: {}mm'.format(lxy, ctau)
示例#2
0
    h.linewidth=2
    h.title=sigtag+' (norm. 30fb)'
    h.legendstyle='L'
    sighs.append(h)

draw([bkgh]+sighs, ylimits=(1e-4,1e3), logy=True)
leg = Legend([bkgh]+sighs, pad=canvas,
            margin=0.25, leftmargin=0.45, topmargin=0.02,
            entrysep=0.01, entryheight=0.02, textsize=10)
leg.Draw()
title = TitleAsLatex('[4#mu] muon-type lepton-jet isolation')
title.Draw()


canvas.SaveAs('{}/ch4mu_muljiso.pdf'.format(outdir))
canvas.clear()


signalsigs = []
for h in sighs:
    h_ = h.clone()
    for i in range(1, h.nbins()+1):
        b = math.sqrt(bkgh.integral(xbin1=i, overflow=True))
        s = h.integral(xbin1=i, overflow=True)
        if s==0 or b==0: h_[i]=h_[i-1]
        else: h_[i] = calculate_simpsig(s, b)
        h_[i].error = 0
    h_.scale(1/h_.integral())
    signalsigs.append(h_)

ROOT.gPad.SetLogy(0)