21x = ROOT.RooRealVar(
"x",
"x", 0, 10)
25mean = ROOT.RooRealVar(
"mean",
"mean of gaussians", 5)
26sigma1 = ROOT.RooRealVar(
"sigma1",
"width of gaussians", 0.5)
27sigma2 = ROOT.RooRealVar(
"sigma2",
"width of gaussians", 1)
29sig1 = ROOT.RooGaussian(
"sig1",
"Signal component 1", x, mean, sigma1)
30sig2 = ROOT.RooGaussian(
"sig2",
"Signal component 2", x, mean, sigma2)
33a0 = ROOT.RooRealVar(
"a0",
"a0", 0.5, 0.0, 1.0)
34a1 = ROOT.RooRealVar(
"a1",
"a1", 0.2, 0.0, 1.0)
35bkg = ROOT.RooChebychev(
"bkg",
"Background", x, [a0, a1])
38sig1frac = ROOT.RooRealVar(
"sig1frac",
"fraction of component 1 in signal", 0.8, 0.0, 1.0)
39sig = ROOT.RooAddPdf(
"sig",
"Signal", [sig1, sig2], [sig1frac])
42bkgfrac = ROOT.RooRealVar(
"bkgfrac",
"fraction of background", 0.5, 0.0, 1.0)
43model = ROOT.RooAddPdf(
"model",
"g1+g2+a", [bkg, sig], [bkgfrac])
48d = model.generate({x}, 10000)
56ll = ROOT.RooLinkedList()
57model.chi2FitTo(dh, ll)
65dsmall = d.reduce(ROOT.RooFit.EventRange(1, 100))
66dhsmall = dsmall.binnedClone()
67chi2_lowstat = model.createChi2(dhsmall)
68print(chi2_lowstat.getVal())