28t = ROOT.RooRealVar(
"t",
"t", 0, 5)
29tmin = ROOT.RooRealVar(
"tmin",
"tmin", 0, 0, 5)
32t.setRange(tmin, ROOT.RooFit.RooConst(t.getMax()))
35tau = ROOT.RooRealVar(
"tau",
"tau", -1.54, -10, -0.1)
36model = ROOT.RooExponential(
"model",
"model", t, tau)
42dall = model.generate(ROOT.RooArgSet(t), 10000)
45tmp = ROOT.RooGaussian(
"gmin",
"gmin", tmin, ROOT.RooFit.RooConst(
46 0), ROOT.RooFit.RooConst(0.5)).generate(ROOT.RooArgSet(tmin), 5000)
49dacc = model.generate(ROOT.RooArgSet(t), ROOT.RooFit.ProtoData(tmp))
54r = model.fitTo(dacc, ROOT.RooFit.Save())
60frame = t.frame(ROOT.RooFit.Title(
"Fit to data with per-event acceptance"))
61dall.plotOn(frame, ROOT.RooFit.MarkerColor(ROOT.kRed),
62 ROOT.RooFit.LineColor(ROOT.kRed))
69c = ROOT.TCanvas(
"rf314_paramranges",
"rf314_paramranges", 600, 600)
70ROOT.gPad.SetLeftMargin(0.15)
71frame.GetYaxis().SetTitleOffset(1.6)
74c.SaveAs(
"rf314_paramranges.png")