limitHistos = rt.TFile.Open('T2tt-plots-0.root')
limit = limitHistos.Get('Limits')
susyNLO = limit.Clone('susyNLO')
susyNLO.Reset()
crossSec = rt.TFile.Open('referenceXSecs.root')
stop = crossSec.Get('stop')
#rt.gROOT.Macro('/Users/wreece/Downloads/SMS_30_04_12/INCL/T2tt/sigma_limit_HAD.C')
#rt.gROOT.Macro('/Users/wreece/Downloads/SMS_30_04_12/INCL/T2tt/sigma_limit_combined.C')
#rt.gROOT.Macro('/Users/wreece/Downloads/SMS_30_04_12/BJET/T2tt/sigma_limit_combined.C')
rt.gROOT.Macro('/Users/wreece/Downloads/SMS_30_04_12/BJET/T2tt/sigma_limit_HAD.C')
incHad = rt.gDirectory.Get('plot_me_13')
for point in sorted(norms.keys()):
#we exclude if xsec is less than th
xs = getBinContent(limit, point[0], point[1])
inc = getBinContent(incHad, point[0], point[1])
th = getBinContent1D(stop,point[0])
if xs < inc:
print point, xs, inc, th
setBinContent(susyNLO, point[0], point[1], xs/inc)
else:
setBinContent(susyNLO, point[0], point[1], 1e-4)
limitFile = rt.TFile.Open('T2tt-Inc-ratio.root','recreate')
susyNLO.Write()
limitFile.Close()
def changePlot(plot, t, x):
excluded = (t >= x)
if True:
setBinContent(plot, point[0], point[1], t/x)