import ROOT
"""
Generate binned Asimov dataset for a continuous pdf.
One should in principle be able to use
pdf.generateBinned(x, n_events, RooFit::ExpectedData()).
Unfortunately it has a problem: it also has the bin bias that this tutorial
demonstrates, to if we would use it, the biases would cancel out.
"""
return data_h
"""
Force numeric integration and do this numeric integration with the
RooBinIntegrator, which sums the function values at the bin centers.
"""
"""
Reset the integrator config to disable the RooBinIntegrator.
"""
fit1 =
expo.fitTo(expo_data, Save=
True, PrintLevel=-1, SumW2Error=
False)
fit2 =
expo.fitTo(expo_data, Save=
True, PrintLevel=-1, SumW2Error=
False)
fit3 =
powerlaw.fitTo(powerlaw_data, Save=
True, PrintLevel=-1, SumW2Error=
False)
fit4 =
powerlaw.fitTo(powerlaw_data, Save=
True, PrintLevel=-1, SumW2Error=
False)
fit5 =
powerlaw.fitTo(powerlaw_data, IntegrateBins=1e-3, Save=
True, PrintLevel=-1, SumW2Error=
False)
fit6 =
model.fitTo(model_data, Save=
True, PrintLevel=-1, SumW2Error=
False)
fit7 =
model.fitTo(model_data, Offset=
"bin", Save=
True, PrintLevel=-1, SumW2Error=
False)
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
[#1] INFO:Eval -- RooRealVar::setRange(x) new range named 'bin' created with bounds [0.1,0.6]
RooFitResult: minimized FCN value: 4754.37, estimated distance to minimum: 3.09852e-09
covariance matrix quality: Full, accurate covariance matrix
Status : MINIMIZE=0 HESSE=0
Floating Parameter FinalValue +/- Error
-------------------- --------------------------
c -1.6862e+00 +/- 1.70e-02
[#0] WARNING:Integration -- RooBinIntegrator::RooBinIntegrator WARNING: integrand provide no binning definition observable #0 substituting default binning of 10 bins
[#1] INFO:NumericIntegration -- RooRealIntegral::init(expo_Int[x]) using numeric integrator RooBinIntegrator to calculate Int(x)
RooFitResult: minimized FCN value: 4440.6, estimated distance to minimum: 5.599e-07
covariance matrix quality: Full, accurate covariance matrix
Status : MINIMIZE=0 HESSE=0
Floating Parameter FinalValue +/- Error
-------------------- --------------------------
c -1.8000e+00 +/- 1.87e-02
RooFitResult: minimized FCN value: 15816.4, estimated distance to minimum: 4.97037e-07
covariance matrix quality: Full, accurate covariance matrix
Status : MINIMIZE=0 HESSE=0
Floating Parameter FinalValue +/- Error
-------------------- --------------------------
a -2.6106e-01 +/- 1.06e-02
[#0] WARNING:Integration -- RooBinIntegrator::RooBinIntegrator WARNING: integrand provide no binning definition observable #0 substituting default binning of 10 bins
[#1] INFO:NumericIntegration -- RooRealIntegral::init(powerlaw_Int[x]) using numeric integrator RooBinIntegrator to calculate Int(x)
RooFitResult: minimized FCN value: 15739.9, estimated distance to minimum: 4.99474e-07
covariance matrix quality: Full, accurate covariance matrix
Status : MINIMIZE=0 HESSE=0
Floating Parameter FinalValue +/- Error
-------------------- --------------------------
a -3.1481e-01 +/- 1.15e-02
RooFitResult: minimized FCN value: 15739.6, estimated distance to minimum: 3.92419e-05
covariance matrix quality: Full, accurate covariance matrix
Status : MINIMIZE=0 HESSE=0
Floating Parameter FinalValue +/- Error
-------------------- --------------------------
a -3.0010e-01 +/- 1.07e-02
[#0] PROGRESS:Generation -- RooAbsPdf::generateBinned(model) Performing costly accept/reject sampling. If this takes too long, use extended mode to speed up the process.
RooFitResult: minimized FCN value: -1.47174e+08, estimated distance to minimum: 0.162058
covariance matrix quality: Full, accurate covariance matrix
Status : MINIMIZE=-1 HESSE=3
Floating Parameter FinalValue +/- Error
-------------------- --------------------------
c -1.7972e+00 +/- 7.39e-04
mu 2.9756e+00 +/- 3.90e-02
nbkg 1.0001e+07 +/- 3.25e+03
nsig 9.4264e+03 +/- 7.36e+02
sigma 4.6849e-01 +/- 2.75e-02
RooFitResult: minimized FCN value: 3416.14, estimated distance to minimum: 0.000238216
covariance matrix quality: Full, accurate covariance matrix
Status : MINIMIZE=0 HESSE=0
Floating Parameter FinalValue +/- Error
-------------------- --------------------------
c -1.7971e+00 +/- 7.26e-04
mu 2.9939e+00 +/- 3.63e-02
nbkg 1.0001e+07 +/- 3.24e+03
nsig 9.2425e+03 +/- 6.91e+02
sigma 4.5747e-01 +/- 2.59e-02