'''
outputtexfile.write(tail)
outputtexfile.close()
#gr_med = PyROOTUtils.Graph(masses, medians)
gr_obs = PyROOTUtils.Graph(masses, lim_obs)
gr_obs.SetMarkerStyle(20)
gr_med = PyROOTUtils.Graph(masses, lim_med)
#gr_med.SetMarkerStyle(7)
gr_med.SetLineStyle(R.kDashed)
gr_med.SetLineWidth(2)
#band_1sigma = PyROOTUtils.Band(masses, minus1sigmas, plus1sigmas, fillColor=R.kGreen)
#band_2sigma = PyROOTUtils.Band(masses, minus2sigmas, plus2sigmas, fillColor=R.kYellow)
band_1sigma = PyROOTUtils.Band(masses, minus1sig, plus1sig, fillColor=R.kGreen)
band_2sigma = PyROOTUtils.Band(masses,
minus2sig,
plus2sig,
fillColor=R.kYellow)
# Prepare frame
hframe = R.TH1I("frame", "frame", 1, min(masses), max(masses))
hframe.GetXaxis().SetTitle("m_{H} [GeV]")
#hframe.GetXaxis().SetRangeUser(300,1000)
#hframe.GetYaxis().SetTitle("95% Limit on #sigma/#sigma_{SM}")
if (mode == 'VBF'):
hframe.GetYaxis().SetTitle(
"95% Limit on #sigma_{VBF} * BR(H #rightarrow ZZ) [fb]")
else:
hframe.GetYaxis().SetTitle(
