2010年5月28日金曜日

pawの2次元fit(2次元ヒストグラムのフィッティング)

c for fitting by x** function @ ANAPAW                                            
      REAL FUNCTION FUNC(X,Y) !for a 2-Dim histogram                              
      COMMON/PAWPAR/PAR(2)
      FUNC=-Y+PAR(1)*X+PAR(2)
      END
こういうテキストファイルを作り、名前をfunc.fとする。
anapaw(paw)上で2次元の絵(histID=100)を出したら、
ANAPAW> His/fit 100 func.f ! 2 par
と打つと、以下のようにフィッティング結果が出る。
 *** Unknown vector
 PAR                                                                            
1

     **********************************************
     *                                            *
     * Function minimization by SUBROUTINE HFITH  *
     * Variable-metric method                     *
     * ID =        105  CHOPT =                   *
     *                                            *
     **********************************************
 Convergence when estimated distance to minimum (EDM) .LT.  0.10E+01

 FCN=   58617.65     FROM MIGRAD    STATUS=CONVERGED     20 CALLS       21 TOTAL
                     EDM=  0.42E-08  STRATEGY=1  ERROR MATRIX UNCERTAINTY=  0.3%

  EXT PARAMETER                                   STEP         FIRST  
  NO.   NAME        VALUE          ERROR          SIZE      DERIVATIVE
   1      P1      -0.23138E-01   0.74368E-03  -0.96424E-07  -0.13032   
   2      P2        139.58       0.67839E-01   0.10123E-05  -0.52769E-03

 CHISQUARE = 0.4187E+02  NPFIT =  1402

0 件のコメント:

コメントを投稿