How to access fit parameters after a fit operation ?
The following example illustrate how to
access the fit parameters (eg. CHI2) just after a
HISTO/FIT
or
VECTOR/FIT
command.
Vector/Del *
Vector/Cre x(10) r 1 2 3 4 5 6 7 8 9 10
Vector/Cre y(10) r 0.9 2.1 2.9 4.1 4.9 6.1 6.9 8.1 8.9 10.1
Vector/Cre ey(10) r 10*0.3
Vector/Cre chi2(1) r
*
Vector/Fit x y ey p1 e
*
Application COMIS Quit
SUBROUTINE chiq
* NPFITS Number of points used in the fit
* NFPAR Number of free parameters
* FITCHI Chisquare
* FITPAR Values of parameters
* FITSIG Errors on parameters
COMMON/HCFITS/NCFITS,NPFITS,NFPAR,FITCHI,FITPAR(35),FITSIG(35)
+ ,FITDER(35)
vector chi2
chi2(1)=fitchi
END
Quit
*
Call chiq
Vector/Pri chi2
*
return
|
PAW > exe chi2
**********************************************
* *
* Function minimization by SUBROUTINE HFITV *
* Variable-metric method *
* ID = 0 CHOPT = E *
* *
**********************************************
Convergence when estimated distance to minimum (EDM) .LT. 0.10E+01
FCN= 1.077442 FROM MINOS STATUS=SUCCESSFUL 18 CALLS 58 TOTAL
EDM= 0.38E-11 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER PARABOLIC MINOS ERRORS
NO. NAME VALUE ERROR NEGATIVE POSITIVE
1 P1 -0.33333E-01 0.20507 -0.20494 0.20494
2 P2 1.0061 0.33054E-01 -0.33029E-01 0.33029E-01
CHISQUARE = 0.1347E+00 NPFIT = 10
CHI2(1) = 0.13468
PAW >
|
Paw.Support@cern.ch