My unbinned fit is set over variable
RooRealVar sbpQvalue(“sbQvalue”, “Q [GeV]”, 0.003, 0.210);
RooDataSet* datSbp defined over sbpQvalue is subjected to fit with model RooAddPdf* modelTot . After the successful fit I plot the profile of the fit results on the binned distribution (histogram):
const Double_t xLo = 0.003;
const Double_t xHi = 0.210;
const Int_t numBins = 69;
RooPlot* qframe6 = sbpQvalue.frame(xLo, xHi, numBins);
qframe6 ->SetYTitle(“candidates per 3 MeV”);
modelTot ->plotOn(qframe6, LineStyle(kSolid), LineColor(kRed+2), LineWidth(2),
I have skipped from the fit the first effective bin to avoid the
non-physical cases where fitter could get into trouble.
The plot is nice and neat but without first bin: I want to have this data
bin (which was out of the fit) still to be shown.
Question: how do I plot the fitted pdf defined over (0.003, 0.210)
and resulted from the fit defined also over (0.003, 0.210) but
over the histogram filled with data over (0.0, 0.210) to show the
first skipped bin ? Important caveat: I do not want to set the fit over
(0., 0.210) and run the fitter over sub-range or “signal range” of (0.003, 0.210)
imposing limits to MINUIT as this latter case is not actually equivalent
to the original case as MINUIT will supplement the fitter
with some functions regularizing the left limit.
I will appreciate if somebody could give me a hint or a couple of
lines of codes…