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

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)上で２次元の絵（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

ANAPWA command

● ANAPAW command 補足版
・ヒストグラムに色をつける
>set hcol  数字

・ヒストグラムのコピー
>hi/copy h1 h2 "newname"

・プロットのマーカーを変える
PAW>set mtyp 数字
・プロットの大きさを変える
PAW>set mscf 0.5とか

　1・・・(・)
  2・・・ ＋
  3・・・ ＋と×を重ね合わせたやつ。
  4・・・ ○
  5・・・ ×
  8・・・ ●
  20・・・ ●
  21・・・ ■
  22・・・ ▲
  23・・・ ▼
  24・・・ ○
  25・・・ □
  26・・・ △
  27・・・ ◇

 http://nuclear.phys.tohoku.ac.jp/~hirose/PAW/mtyp.htmlを参考にさせていただきました。
・zone上２、下１

>zon 2 2

>h/pl 100

>h/pl 110

>zon 1 2 2 S

>h/pl 130 text x0 y0 'letter' 0.5 (x0,y0）の位置にletterを書きこむ。文字=0.5 title 'Higgs"5#4 leptons' タイトルの位置に Higgs→4 leptonsを書く。 atitle 'xaxsis' 'yaxisis' X軸とY軸に軸タイトルを 書く。 [tt]"5#10^33!m?[gg]! ττ→1033mγγ： [XX]はギリシャ文字指定、~XX!上付、?YY!下付 box x0 x1 y0 y1 箱を(x0-x1,y0-y1)に書く。 igset fais 3でハッチ。 ●参考 http://atlas.kek.jp/physics/nyumon/paw-commands.html http://nuclear.phys.tohoku.ac.jp/~hirose/PAW/bibouroku_PAW.html