How to Get the Value of the Width of the Error Band from VisualizeError

Roofit and RooStats
,
1

I use the following code to draw the error band of a fitting curve:

RooRealVar x("x", "x", pd[i], pu[i]);
RooDataHist data("data", "data", x, hvar[i]);
RooRealVar mg1("mg1", "mg1", pmean[i], pd[i], pu[i]);
RooRealVar sg1("sg1", "sg1", insg1[i], lsg1[i], hsg1[i]);
RooGaussian gauss1("gauss1","gauss1", x, mg1, sg1);
RooRealVar mg2("mg2", "mg2", pmean[i], pd[i], pu[i]);
RooRealVar sg2("sg2", "sg2", insg2[i], lsg2[i], hsg2[i]);
RooGaussian gauss2("gauss2","gauss2", x, mg2, sg2 );
RooRealVar frac("frac","frac",infrac[i],lfrac[i],hfrac[i]);
x.setRange("mini", pd[i], pu[i]);
RooAddPdf gauss("gauss","gauss",RooArgList(gauss1,gauss2),RooArgList(frac));

RooRealVar co1("co1","coefficienct #1", inco1[i], dco1[i], uco1[i]);
RooRealVar co2("co2","coefficienct #2", inco2[i], dco2[i], uco2[i]);
RooChebychev  poly("poly", "poly PDF", x, RooArgList(co1,co2));
RooRealVar f1("Sig", "fra1", innsig[i], lnsig[i], hnsig[i]);
RooRealVar f2("Bkg", "fra2", innbkg[i], lnbkg[i], hnbkg[i]);
RooAddPdf sum("sum", "sig+bak", RooArgList(gauss,poly), RooArgList(f1,f2));
RooFitResult *fitresult;
RooPlot *xframe=x.frame();
data.plotOn(xframe);
sum.plotOn(xframe);
sum.plotOn(xframe, VisualizeError(*fitresult), Components("poly"), FillColor(kSpring+10));
sum.plotOn(xframe, VisualizeError(*fitresult), Components("gauss1"), FillColor(kRed-9));

However, the error band seems abnormal. I want to obtain the width of the error band.

a
a956×1174 48.2 KB

Would you like to add more details or context to this? For example, any specific issues you’re encountering or the expected output?

@moneta

2

Hi,

I am sorry to read you are encountering this issue.
Since you already pinged @moneta , I allow myself to ask @couet if he sees an obvious issue with the graphics settings.

Cheers,
D

3

Hi, do you have some minimal C++ macro reprducing the problem you are experimenting ?

4

fithp1.cxx (15.1 KB)

5

I think another file is required:

Processing fithp1.cxx...
In file included from input_line_9:1:
/Users/couet/Downloads/fithp1.cxx:1:10: fatal error: '/Volumes/D/Work/BESIII/bes3plotstyle.c' file not found
#include "/Volumes/D/Work/BESIII/bes3plotstyle.c"
         ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
6

cms.txt (565 Bytes)
fith.txt (162 Bytes)
bes3plotstyle.c (12.2 KB)
bes3plotstyle.h (5.9 KB)
mc_2232S_02a.root (11.4 KB)

7

Thanks. Now I get this:

root [0] 
Processing fithp1.cxx...
Error in <TFile::TFile>: file /Users/couet/Downloads/mc_S_02.root does not exist
8

PID1.zip (1.4 MB)
Please use fithp1.cxx

9

same thing:

root [0] .x fithp1.cxx
Processing fitp.cxx...
Error in <TFile::TFile>: file /Users/couet/Downloads/PID1/mc_S_02.root does not exist
10

What does Line127 shows?

image
image1284×324 82.7 KB

11

Screenshot 2024-06-20 at 13.28.56
Screenshot 2024-06-20 at 13.28.56856×172 20.1 KB

12

PID1.zip (1.5 MB)

13

Hi, the script runs now. However, I had to modify the first two lines of the script to find the included files.

When I run your macro, several EPS and PDF files are generated. They look like the following picture (if I understand correctly, the first plot you posted is a close-up zoom on the top of the peak). What’s wrong with this picture?

Screenshot 2024-06-20 at 13.47.04
Screenshot 2024-06-20 at 13.47.042368×1698 368 KB

14

Please zoom in. I want read the value of the width of the red error band. This plot just an example.
image

15

The width is defined in pixels in the PDF file, and the PDF viewer renders the line accordingly. I’m not sure what value you aim to achieve with this significant zoom-in. It doesn’t seem to make much sense to me to zoom in so much.

16

I you want to make a zoom you should make in at the C level not in the pdf file.

17

Hi

This error band is from VisualizeError
I just wish to know how to draw the errorband from VisualizeError(*fitresult), not means this plot have problem.

Cheers

WX20240620-214131@2x
WX20240620-214131@2x1010×836 27.5 KB

18

The error bands in ROOT are drawn using the option E3. I guess (but I’m not sure) it is the same in RooPlot? (@moneta?)

19

This error bands is from Rootfit
sum.plotOn(xframe, VisualizeError(*fitresult), Components(“gauss1”), FillColor(kRed-9));
But I cannot read the value.

20

The error bands are drawn by connecting the extremities of the error bars. So, at a given point (x,y), the top of the band is equal to y + erry and the bottom is y - erry.