ResolutionModel=tres.model(),
Bins=bin_spec,
Relative=rel_spec,
)
pdf = Single_Exponent_Time(Name="pdf", time=t, resolutionModel=res_model)
pdf = pdf.pdf()
gen_observables = [t, hlt1_excl_biased, hlt2_unbiased, hlt2_biased]
## Get proto data
proto_data = None
if options.protodata:
proto_file = TFile.Open(options.protodata)
if not proto_file.IsOpen():
raise OSError
for key in proto_file.GetListOfKeys():
if key.GetClassName() == "RooDataSet":
proto_data = key.ReadObj()
break
else:
raise RunTimeError
if not proto_data:
raise RunTimeError
## run the toy
toy.set_fit_opts(**dict(Verbose=False))
toy.run(Observables=gen_observables, Pdf=pdf, GenPdf=pdf, ProtoData=proto_data)
toy.write_output()
