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