Re: [ROOT] Getting status of a fit.

From: Drew Alton (dalton@fnal.gov)
Date: Wed Apr 14 2004 - 15:20:10 MEST


thanks...fCstatu gives me what I am after...but it's not clear to me
what fStatus gives.  It seems to always be 0 regardless of the fit
results. 
                                          -Drew




On Wed, 14 Apr 2004, Rene Brun wrote:

> Hi Drew,
> 
> You can access all TMinuit members (they are public) after the fit
> root > hist.Fit("gaus");
> root > gMinuit->fStatus
> root > gMinuit->fCstatu
> 
> Rene Brun
> 
> Drew Alton wrote:
> > 
> > Hello,
> >    I am doing fits in Root.  When I do them interactively I see something
> > like this:
> >  FCN=1.26222e-15 FROM MINOS     STATUS=SUCCESSFUL   3927 CALLS        4804
> > TOTAL
> >                      EDM=2.52464e-15    STRATEGY= 1      ERROR MATRIX
> > ACCURATE
> >   EXT PARAMETER                  PARABOLIC         MINOS ERRORS
> >   NO.   NAME      VALUE            ERROR      NEGATIVE      POSITIVE
> >    1  p0           2.06813e-01   7.75185e-02  -6.69639e-02   8.88198e-02
> >    2  p1           7.93187e-01   7.86323e-02  -8.98656e-02   6.83877e-02
> >    3  p2           6.99663e-01   1.64727e-01  -1.58526e-01   1.73633e-01
> >    4  p3           3.00337e-01   1.86156e-01  -1.93731e-01   1.85260e-01
> >    5  p4           7.35839e-01   1.21140e-01  -1.08260e-01   1.40031e-01
> >    6  p5           1.16457e-01   3.56975e-02  -3.99971e-02   3.22014e-02
> >    7  p6           6.41985e-01   9.11096e-02  -8.25384e-02   1.01704e-01
> >    8  p7           1.50805e-01   3.03124e-02  -3.35608e-02   2.75935e-02
> > 
> > The top line says STATUS=SUCCESSFUL, however sometimes I have FAILED
> > or PROBLEMS.  I assume there is some way to get this information but
> > I cant find any access methods that sound right.
> > 
> > I am basically defining a histogram and then a fit fuction and then
> > I do:
> > histo->Fit("System8","rnem");
> > 
> > I figured out how to get the parameters, the errors, the chi2, but
> > the success/failed flag escapes me.
> > 
> > Anyone know how to get this information?
> >                                           -Drew
> 



This archive was generated by hypermail 2b29 : Sun Jan 02 2005 - 05:50:07 MET