''' 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(