東京工業大学 >
大学院理工学研究科 >
基礎物理学専攻 >
中村研究室 >
メンバー >
Takayuki Sako >
PAW/ANAPAW コマンドPAW : 気が向いたものだけ。histogram関係を中心に充実させる予定。 ANAPAW : コマンドをほぼ網羅。 command list
alias list
Tips
元ネタPAWPAW Reference manual
PAW/ANAPAW上でhelpと打っても見られるPAWのコマンドも全てANAPAW上で確認しているため純粋なPAWとは挙動が違う部分あるかもしれない。しかし、PAWでの確認は行っていない。 PAW User's Guide [ps.gz]
PAW tutorial
PAW Frequently Asked Questions
ANAPAWANAPAWのhelpANAPAW上でhelpと打つと見られる。 ANAPAW> help From /... 1: KUIP Command Processor commands. 2: MACRO Macro Processor commands. 3: VECTOR Vector Processor commands. 4: HISTOGRAM Manipulation of histograms. Interface to the HBOOK package. 5: FUNCTION Operations with Functions. Creation and plotting. 6: NTUPLE Ntuple creation and related operations. 7: GRAPHICS Interface to the graphics packages HPLOT and HIGZ. 8: PICTURE Creation and manipulation of HIGZ pictures. 9: ZEBRA Interfaces to the ZEBRA RZ, FZ and DZ packages. 10: FORTRAN Interface to MINUIT, COMIS, SIGMA and FORTRAN Input/Output. 11: NETWORK To access files on remote computers. 12: MLP Multi-Layer Perceptron (MLP). 13: OBSOLETE Obsolete commands 14: ANALYS On(Off)-Line Analysis Software on PAW. Enter a number ('0'=Top, '\'=one level back, 'Q'=command mode):14 Menu /ANALYS 1: /ANALYS On(Off)-Line Analysis Software on PAW. From /ANALYS/... 2: * LOOP This command provides you to see the histograms (1D and 2D-Histograms). For the 2D-Histogram, it is shown by using SCATTER-PLOT. 3: * STATUS Print Current Status. 4: * RDMP Store the Raw Data to a file with a given Gate. 5: * ANTUPLE Store the Raw Data to Ntuple files with a given Gate. 6: * UFLAG Change User Flag defined in ENC_***. 7: * SETNUMBER Set Run Number of a current data file. 8: DEFINITION Define and Book Histograms using ANALYS-CODE. 9: HISTOGRAM Manipulation of Histograms. ______ 0: Top level menu Enter a number ('0'=Top, '\'=one level back, 'Q'=command mode): さらにweb上には Julien氏のページがある。 このページの製作にあたっては上記ページも参考にした。 また、任意のコマンドcommandに対し直接helpを見ることもできる。
ANAPAW> help command
※このページはこれらhelpの翻訳にあらず。あくまでこれらを参考とした備忘録。適宜補足の情報が入っている。
同時に間違いなども含まれている可能性がある。※ANAPAWオフィシャルのhelpには含まれていない項目が追加されているため、helpにあるコマンドの順番とは必ずしも一致しない。 /KUIPshellホストOSのコマンドを実行する。要はLinuxのコマンドが実行出来るということ。shでも可。
PAW> shell [ CMD ]
waitポーズ(一時休止)を発生させる。主としてマクロに用いられる。
PAW> wait [ STRING SEC ]
exitインタラクティブ・セッションを終了する。
PAW> exit
quitインタラクティブ・セッションを終了する。
PAW> quit
versionPAWとパッケージのバージョンを表示する。PAW > version +---------+--------+---------------------+ | Package | Vers.# | Last modified at | +---------+--------+---------------------+ | PAW 2.14/04 12/01/2004 10.45.00 | | KUIP 2.09/00 13/02/2003 15:25.00 | | HBOOK 4.28/01 04/06/2003 17.00.00 | | COMIS 1.24/00 27/10/2004 11:00.00 | | HPLOT 5.24/02 12/11/2002 16.30.00 | | HIGZ 1.29/04 10/06/2004 10:00:00 | | SIGMA 1.10/03 10/02/1994 17.18.05 | +----------------------------------------+ /KUIP/SET_SHOW/host_pagerリードオンリーモードでファイルを読むためのホストコマンドを設定する。
PAW> host_pager [ PAGER ]
例えば、
PAW> host_pager less
でhelpなどをlessで読むことができる。
host_editorエディターを設定する。
PAW> HOST_EDITOR [ editor top left width height dxpad dypad npads ]
デフォルトではviが設定されていたりする。
すると、例えばコマンド'b'などで間違ってバグレポートの編集画面などがでた際にviが起動するのが煩わしい。
とりあえず使い慣れたemacsに設定しておく。
PAW> HOST_EDITOR 'emacs -nw'
/HISTOGRAMhistogram/listヒストグラムのリストを表示する。hi/lでも可
PAW> hi/list [ chopt ]
PAW> hi/l ===> Directory : //PAWC 11 (1) zenith angle distribution 12 (1) azimuthal angle distribution 13 (1) Energy spectrum [GeV] 14 (1) Energy spectrum [GeV] 15 (2) Muon intensity Energy-theta 16 (2) Muon intensity phi-Energy 17 (2) Muon intensity phi-theta PAW> hi/l I ===> Directory : //PAWC 1 .............................................................................................. . . . HBOOK HBOOK CERN VERSION 4.25 HISTOGRAM AND PLOT INDEX 25/05/2011 . . . .............................................................................................. . . . NO TITLE ID B/C ENTRIES DIM NCHA LOWER UPPER ADDRESS LENGTH . . . .............................................................................................. . . . . . 1 zenith angle 11 32 157251 1 X 700 0.000E+00 0.900E+02 9903028 749 . . . . . . 2 azimuthal angle 12 32 157251 1 X 1500 0.000E+00 0.180E+03 9902277 1549 . . . . . . 3 Energy spectrum 13 32 157251 1 X 500 0.000E+00 0.100E+03 9900727 548 . . . . . . 4 Energy spectrum 14 32 157251 1 X 500 0.000E+00 0.100E+02 9900177 548 . . . . . . 5 Energy-theta 15 32 157251 2 X 500 0.000E+00 0.100E+02 9899624 352453 . . Y 700 0.000E+00 0.900E+02 9547199 352416 . . . . 6 phi-Energy 16 32 157251 2 X 1500 0.000E+00 0.180E+03 9547169 754053 . . Y 500 0.000E+00 0.100E+02 8793144 754016 . . . . 7 phi-theta 17 32 157251 2 X 1500 0.000E+00 0.180E+03 87931151054452 . . Y 700 0.000E+00 0.900E+02 77386901054416 . histogram/deleteカレント・ディレクトリ(メモリ)のヒストグラム/Ntupleを消去する。hi/dでも可。
PAW> hi/del ID
histogram/plotヒストグラムをプロット(表示)する。
PAW> hi/pl [ id chopt ]
PAW> hi/pl id(ic1:ic2)
ic1,ic2が整数の時、チャンネルic1からic2の間でプロットする。
PAW> hi/pl id(x1:x2)
x1,x2が実数の時(すなわち.が付いている時)、xの値がx1からx2の間でプロットする。
(左図) PAW> zone 1 2; hi/pl 101 ; hi/pl 101(0:100)
(右図) PAW> zone 1 2; hi/pl 101 ; hi/pl 101(0.:100.) Ic1もしくはx1を省略した場合、最初のチャンネルが使われる。同様に、ic2もしくはx2を省略した場合最後のチャンネルが使われる。
(左図) PAW> zone 1 2; hi/pl 101 ; hi/pl 101(150.:)
(右図) PAW> zone 1 2; hi/pl 101 ; hi/pl 101(25.:) この表示範囲を限定する機能は2Dプロットでも同様に使うことができる。
(左図) PAW> zone 1 2; hi/plo 6000 colz ; hi/plo 6000(20.:50.) colz
(右図) PAW> zone 1 2; hi/plo 6000 colz ; hi/plo 6000(20.:50.,10.:70.) colz 一次元ヒストグラムはLEGOもしくはSURFオプションをつけることが可能である。この場合、角度はTHETA=1,PHI=-1になる。 具体例、
PAW> hi/plo 101 LEGO
PAW> hi/plo 101 SURF
ちなみにカラーインデックス8を白にしている場合(color 8 1 1 1)、当然ながら0に近い値が白になる。
PAW> pal 1 ; color 8 1 1 1 ; hi/plo 101 SURF1
PAW> pal 1 ; color 8 1 1 1 ; hi/plo 101 SURF2
オプション'E'を使っている時、マーカーの種類はSET MTYP、マーカーのサイズはSET KSIZ、マーカーの色はSET PMCIで変更できる。
(左図) PAW> set MTYP ; set KSIZ ; set PMCI ; hi/plo 101 e
(右図) PAW> set MTYP 5 ; set KSIZ 0.5 ; set PMCI 2 ; hi/plo 101 e
オプション'E1'を使っている時、エラーバーの端のマークのサイズはSET KSIZで変更できる。
(左図) PAW> set KSIZ ; hi/plo 101 e1
(右図) PAW> set KSIZ ; hi/plo 101 e1
オプションLOGXとLOGYは等高線プロット(オプションCONT)に使うことはできない。すなわち、等高線プロットは常に線形のx,y軸に対してプロットされる。 ヒストグラムIDのx軸に対するプロジェクション(射影)をプロットするには ヒストグラムIDをy方向に3つにスライスしてプロットするには LEGOおよびSURFプロットに対しては直交座標系、極座標系、円柱座標系、球座標系、、、、、が使える。 ANGLESコマンドにより見る角度を変更することができる。3Dプロットにおいて軸、フロントボックス、バックボックスの表示を消すためにはA,FB,BBオプションを使う。 zoomヒストグラムのチャンネルのうちICMINからICMAXまでをプロットする。すなわち拡大。zでも可。
PAW> zoom [ ID CHPOT ICMIN ICMAX ]
many_plots一つもしくはいくつかの一次元ヒストグラムを同じプロットで表示する。
PAW> many_plots IDLIST
PAW> many 101 102 103 104 105
histogram/copyヒストグラムをコピーする。hi/cでも可。
PAW> hi/copy id1 id2 [ title ]
PAW> hi/c 110 210 NEW_TITLE
また、x軸やy軸の範囲を変更してコピーすることもできる。以下、具体例
PAW> hi/copy 110(0.:1500.) 310 NEW_TITLE
histogram/fit
PAW> hi/fit ID FUNC [ CHOPT NP PAR STEP PMIN PMAN ERRPAR ]
HISTOGRAM/2D PLOT/多様なフォーマットで二次元ヒストグラムをプロットする。lego
PAW> lego [ ID THETA PHI CHOPT ]
PAW > LEGO ID1+ID2+ID3
ヒストグラムIDの上限は10である。それぞれのIDの色は ATT/PALETTEにより与えられる。
PAW > SET HCOL 2
ヒストグラムの色を2 (赤)にする。
PAW > LEGO 20
レゴを線と一緒に表示する。
PAW > LEGO 20 ! ! 1
レゴを異なる光と共に表示する。
PAW > LEGO 20 ! ! 2
レゴを色と共に表示する。
PAW > PALETTE 1 3 2 3 4
PAW > SET HCOL 0.1
リスト1を使った累積レゴプロット
PAW > LEGO 10+20+30
累積レゴプロットを線と共に表示する。
PAW > LEGO 10+20+30 ! ! 1 |
surface
surface [ ID THETA PHI CHOPT ]
PAW > ATT/PAL 1 3 2 3 4
パレット1を3エレメントで定義する。
PAW > SET HCOL 0.1
ヒストグラムの色をリスト1に設定する。
PAW > SET NDVZ 4
Zの分割数を4に設定する。
PAW > SURF id 90 0 2
等高線を描く。
HISTOGRAM/CREATEtitle_globalグローバルタイトルを設定する。
PAW> title_globa [ chtitl chopt ]
contour
contour [ ID NLEVEL CHOPT PARAM ]
PAW> Fun2 2 x*y 40 0 1 40 0 1 ’ ’
二次元ヒストグラムを作る。
PAW> V/Cr PAR(5) R .1 .11 .3 .31 .5
等高線を定義する。
等距離でない等高線プロットを描く。
オプション'3'で等距離でない等高線プロットを描くことはできない。
オプションLOGX,LOGYは等高線プロットに対しては動作しない。すなわち、等高線プロットは常に線形のX,Y軸に対して動作する。
PAW> Contour 2 5 ! PAR /HISTOGRAM/OPERATIONS/ヒストグラムの操作と比較。addヒストグラムの足し算を行う。
PAW> add ID1 ID2 ID3 [ C1 C2 OPTION ]
subtractヒストグラムの引き算を行う。
PAW> subtract ID1 ID2 ID3 [ C1 C2 OPTION ]
multiplyヒストグラムの掛け算を行う。
PAW> multiply ID1 ID2 ID3 [ C1 C2 OPTION ]
divideヒストグラムの割り算を行う。
PAW> divide ID1 ID2 ID3 [ C1 C2 OPTION ]
resetヒストグラムの内容とエラーをリセットする。ビンの定義は変更しない。
PAW> reset ID [ TITLE ]
/HISTOGRAM/SETmaximumy軸の最大値を設定する。maxでも可。
PAW> maximum ID VMAX
minimumy軸の最小値を設定する。minでも可。
PAW> minimum ID VMIN
PAW> min ID 0
こうするとリセットされる。
PAW> max ID 0 maximum,minimumは既存のヒストグラムの設定を変更する。 最大・最小を変更したヒストグラムを作成してすぐに表示する際にはmamiを使う。 /NTUPLEvmemNtupleの解析の間、Ntuple columns の保存に使っている動的メモリバッファのサイズを表示もしくは変更する。デフォルトは10MB。上限は128MB。
PAW> vmem [ MXSIZE ]
PAW> vmem -1 Dynamic Ntuple Memory Buffer Statistics (in bytes): Upper Limit Used Free 10485760 0 10485760 PAW> vmem -2 Variable Lun Ntuple Type Usage Size ----------------------------------------------------------------- No variables in memory. PAW> vmem 128 PAW> vmem -1 Dynamic Ntuple Memory Buffer Statistics (in bytes): Upper Limit Used Free 134217728 0 134217728 /GRAPHICS/グラフィックスパッケージHPLOTとHIGZのインターフェースsetHPLOT設定の詳細を設定する。
PAW> set [ CHATT VALUE ]
STAT : 統計情報
PAW> set stat RMS Mean Entries ID
PAW> set stat 0000
としても
PAW> set stat 1111
と同じになってしまう。統計を消したいときは素直に
PAW> opt3 nsta
とする。
PAW Reference Manualにある情報だけ見ると以上のような解釈になる。だがPAW User's Guide [ps.gz]を見ると正しくは以下のような実装らしい。
PAW> set stat AOURMEI
FIT : フィット情報
PAW> set fit X2/ndf Error Par
FCOL,FWID : 関数の線の色、幅デフォルトではフィットした線が見にくいが以下のようにすると見やすくなる。
ANAPAW> zon 1 2 ; set fcol 1 ; set fwid 1 ; ht 244
ANAPAW> set fcol 2 ; set fwid 7 ; ht 244 HCOL : ヒストグラムの色ヒストグラムの色を設定する。
PAW> set hcol FBI
PAW> set hcol 1203
ならば、
小林氏のページに載っているようなさらに細かい色を使うには
PAW> set ncol 58
でカラーインデックスの数を増やして
PAW> pal 1
適当にパレットを設定し(set ncolの直後はグレースケールになるため)、例えば
PAW> set hcol 10216
とすれば
さらに応用。
PAW> set stat 1110001 *統計情報をID,アンダー/オーバーフロー、カウント数に変更
PAW> opt nbox PAW> set bwid 7 *ボックス(画像の周りの線ではなく軸)の幅 PAW> set ncol 58 PAW> pal 1 PAW> set hcol 10418 PAW> set htyp 2 *ヒストグラム内側のhatch PAW> set hwid 5 *ヒストグラムの線の幅 PAW> hi/pl 101 option
PAW> option [ CHOPTN ]
SQR : 正方形
PAW> size 18 24
とした後表示している
TIT : タイトルHTITとUTIT BOX : ボックス(囲み線)
TIC : 目盛
STA : 統計情報
FIT : フィット情報
GRID : 目盛
DATE : 日付
FILE : ファイル名
P : ページ数
LOGY,LINY : 線形、ログ
metafile※個人的にはこの方法を推奨しない。メタファイルロジカルユニットとメタファイルタイプを設定する。このコマンドはグラフィックスの出力を制御する。
PAW> metafile [ LUN METAFL ]
LUN =0 : 画面にのみ出力する LUN = 10 : ユニット10を画面にもメタファイルにも出力する。 コマンドFORTRAN/FILEで新たなファイルを開き、 FORTRAN/CLOSEで閉じる。
PAW > FOR/FILE 44 test.ps
PAW > NULL 0 10 0 10 PAW > CLOSE 44
-[Format][Nx][Ny][Type] [Format][Format] は0から99までの整数であり、紙のフォーマットを定義している。 例えば、Format=3 はスタンダードなA3 formatである。Format=4 と Format=0 は同じ A4 を定義している。A0 は Format=99である。The US format Letter はFormat=100である。 The US format Legal は Format=200である。The US format Ledger はFormat=300である。[Nx,Ny][Nx, Ny] はx,y軸それぞれのzoneを決定する。 Nx と Ny は1から9までの整数である。[Type]
3: Encapsulated PostScript. : このTypeは他の文章にインクルードされることを許可する。例えば LaTeX である。このTypeでは Nx と Ny は常に1でなければならず、 Format は意味を持たない。画像のサイズはSIZEコマンドを使ってユーザーが指定しなければならない。そこで Encapsulated PostScriptに対する workstation typeは -113となる。例えば、 encapsulated PostScriptファイルのファイル名が example.epsである時、LaTeX にインクルードするためには(LaTeX ファイルの中で)以下のように記述すれば良い。
\begin{figure}
Typeが1,2,4, 5の時、画像はページの中心になり、紙の中で使える領域はA4フォーマットの大きさに比例する。
\epsffile{example.eps} \caption{Example of Encapsulated PostScript in LaTeX.} \label{EXAMPLE} \end{figure} -111 もしくは -4111は A4を分割しない。-6322 は横長 A6 を2行3列に分割する。 +-------+-------+-------+ | 1 | 2 | 3 | +-------+-------+-------+ | 4 | 5 | 6 | +-------+-------+-------+最初の画像は領域1に描かれる。画面を消した後は上図のような順で次の領域に画像が出力される。ページが埋まったら、同じグリッドで新しいページが使われる。紙の節約のために何もないページは印刷されない。 A12より小さいフォーマットは無視される。使用できる PostScript workstation types の合計値は 4x9x9x13+1 = 4213 ! このコマンドは ワークステーションID 2でメタファイルを開く。 個人的な主張PAWには以上の方法で複数ページに渡る、しかも詳細なオプションを指定して画像を出力する方法が存在する。しかし、この方法を使うべきではない。 ロジカルユニットを開いている間にミスをすると残念な画像になるか、もう一度同じ操作を繰り返すことになる。それはANAPAWのpson,psoffを使っても同じこと。ヒストグラムをプロットして、その結果を見て改めてファイル名を指定して画像を保存する方が合理的だ。その方法はpic/priに書いた。こちらを推奨する。workstationワークステーション(ウィンドウ)を作る/消す、ステータスを変える。woでも可。
PAW> workstaton IWKID [ CHOPT IWTYP ]
メタファイルをアクティブにしているならIWKIDは2にすべきではない。コマンドMETAFILEが常にワークステーション2を使うからである。 実際に使ってみる。普通にPAW(ANAPAW)を立ち上げるとHIGZ_01のウィンドウが立ち上がっている。
ANAPAW> work 3 oa
とすると、新たなウィンドウHIGZ_03が立ち上がる。
ANAPAW> work -1 a ; ht 103
とすれば、ワークステーション1のみアクティブにし、それ以外(=3)は相補的なアクションを行っている(この場合ディアクティブにする)。
よってワークステーション1でのみ103をプロット。その後は逆の動作でワークステーション3にのみ6000をプロットしている。
ANAPAW> work -3 a ; ht 6000 colz ワークステーションIDは99までが正しく動作するらしい。100を越えるとIDが**になってしまう。ただ、表示がおかしいだけで正しい動作はしているようだ。 デフォルトのワークステーションの挙動が気に食わないので便利に使えるようにマクロを作った。→ ANAPAW マクロ /GRAPHICS/VIEWING/zone画像を幾つに分割するか設定する。
PAW> zone [ NX NY IFIRST CHOPT ]
PAW> zone 1 1 1
なので、1 x 1 に分割、1番目から始まる。
PAW> zone 2 2
だと、2 x 2に分割、一番目から始まる。
PAW> zone 3 3 5
だと、3 x 3に分割、5番目から始まる。
少し応用。
ANAPAW> zone 1 2
1 x 2に分割、1番目にヒストグラム。
ANAPAW> ht 101
ANAPAW> zone 2 2 3 s
2 x 2に分割。3番目は2行目1列目で左下に来ることになる。
ANAPAW> ht 102
ANAPAW> zone 2 4 6 s
2 x 4に分割。2行目までで4個。3行目が5,6番になるので6番から始めると右下の上の段になる。
ANAPAW> ht 103
ANAPAW> zone 4 4 15 s
4 x 4に分割。全部で16個。最後の1つ手前が15番になる。
ANAPAW> ht 104 コマンドZONEはnormalisation transformation(WN,SVP,SELNTを参照)を定義しない。 size画像のサイズを設定する。
PAW> size [ XSIZE YSIZE ]
selntNormalisation transformation numberを選択する。selでも可。
PAW> sel nt
If ZONE 2 2 is active , then: If ZONE 1 1 is active, then: +------------------------------+ +-----------------------------+ | | | | | +----------+ +---------+ | | +-----------------------+ | | | | | | | | | | | | | NT=10 | | NT=20 | | | | | | | | | | | | | | | | | +----------+ +---------+ | | | | | | | | | NT=10 | | | +----------+ +---------+ | | | | | | | | | | | | | | | | | NT=30 | | NT=40 | | | | | | | | | | | | | | | | | +----------+ +---------+ | | | | | | | | +-----------------------+ | | NT=1 | | NT=1 | +------------------------------+ +-----------------------------+itxで文字を書く際などに利用する。 /GRAPHICS/PRIMITIVESitx文字を書く。
PAW> itx x y text
graphベクターをプロットする(グラフにする)。
PAW> graph n x y [ chopt ]
/GRAPHICS_ATTRIBUTESHIGZの属性を変更する。※PAWのグラフィックスに関しては小林氏 が詳細をまとめている。 color_table色を指定する。coloでも可。
PAW> color_table ICOL [ RED GREEN BLUE ]
paletteパレットを指定する。
PAW> palette PALNB [ NEL LIST ]
A)NCOLが8以下の時、パレットは自動的に8の基本色で埋まる。
ANAPAW> set NCOL 8
このとき、8エレメントは0,5,7,3,6,2,4,1となる。実際にヒストグラムを表示させてみると、
小さい値から順に0(白),5(黄色),7(水色),3(黄緑),6(ピンク),2(赤),4(青),1(黒)になっていることが分かる。
ANAPAW> pale 1 ANAPAW> ht 6000 colz また、
ANAPAW> set NCOL 4
このとき、4エレメントは0,5,7,3となる。実際にヒストグラムを表示させてみると、
小さい値から順に0(白),5(黄色),7(水色),3(黄緑)になっていることが分かる。
ただし、col(z)はエレメント数を8より小さくして表示できないらしい?NCOL 4にしても出力結果は8と同じになってしまった。
ANAPAW> pale 1 ANAPAW> ht 6000 colz B)NCOLが8より大きい時、パレットは青から赤まで連続的に変化するように色が埋まる(グラデーションになる)。 これはジオグラフィカル・パレットと呼ばれる。例えば、
ANAPAW> set NCOL 16
ANAPAW> pale 1 ANAPAW> ht 6000 colz コマンドSET NCOLの後カラーインデックスの8からNCOLまでは階調レベル(グレースケール)設定される。 コマンドPALETTE 1は青から赤までのジオグラフィカル・パレットに設定し直している。 よってSET NCOLの後にPALETTEを設定するのを忘れるとグレースケールで表示される。
ANAPAW> set NCOL 26
ANAPAW> ht 6000 colz /GRAPHICS/HPLOT多様な HPLOT オブジェクト (symbols, errors, key, etc.)を描く。symbols
PAW> symbols X Y N [ ISYMB SSIZE ]
errorsエラーバーを描く。
PAW> gra/hp/errors x y ex ey n [ isymb ssize chopt ]
エラーバーの色は set hcolで変わる。set PLCIではない。
atitle現在プロットしているzoneの軸のタイトルを描く。
PAW> atitle [ XTIT YTIT ZTIT IALGN CHOPT ]
PAW> NUL 0 10 0 10
PAW> NUL 0 100 0 100 S PAW> ATITLE ’End of axis’ ’Middle of axis on the right’ ! 320 R PAW> ATITLE ’Beginning of axis’ ’End of axis’ ! 130 PAW> ATITLE ’Middle of axis on the top’ ’Beginning of axis’ ! 210 T grid
PAW> grid
グリッドを描く。グリッドのパラメーターは line と text attributes により制御される(コマンド SET)。
/PICTUREHIGZ pictureの生成と制御。
PAW> pic/pri [ FILE WIDTH HEIGHT ]
gifファイルに対してはメモリーが使われない。gifファイルの場合はウィンドウナンバー1の直接の「ハードコピー」になる。 WIDTHとHEIGHTはgifの出力ファイルのスケールを許可する。 FILE=HIGZPRINTER もしくは FILE=' 'の時、paw.psが生成され、環境変数HIGZPRINTERにより定義されるOSのコマンドが実行される。 環境変数 HIGZPRINTER は以下のように設定できる。 UNIXならば、
setenv HIGZPRINTER 'lp -dprinter_name paw.ps'
もしくは
export HIGZPRINTER='lp -dprinter_name paw.ps'
である。
※補足予め以下のオプションを設定しておく。
PAW> opt ZFL1
その上でヒストグラムをプロットし、現在表示している画像を保存するときには
PAW> pic/pri filename.eps
とする。すると filename.epsが出力される。
デフォルトではZFLオプションはNZFLになっておりZデータベースを使っていないためにPostScriptを出力することはできない。
PAW> pic/pri filename.ps
となってしまう。ただこの状態でもgifなら出力が可能。
> pic/pri filename.ps *** No picture in memory
PAW> pic/pri filename.gif
何れにしても(ANAPAWなら).analogon.kumacに最初からopt ZFL1を記述しておいて、eps出力を可能にしておくべき。
また、.psと.epsを比較した場合、.epsが最初から他のファイルに挿入されることを前提としているため.epsの方が取り回しが良い(何も考えずにTeXに貼れる)。
その他の(些細な)違いとしては何も考えずに印刷した際、デフォルトでは.psは紙の中央に、.epsは紙の左下に印刷される。このどちらを良しとするかは個人の価値観による。
以下、余談。 PAWで画像を保存するためにはMETAFILEを使う方法が存在し、for/fileで予めファイルを開いておく方法が(数少ない)PAWを説明したWEBサイトで盛んに紹介されているらしい。 また、ANAPAW標準のpsonもこれと全く同じ事を行なっている。その利点として現在表示しているヒストグラムだけでなく複数ページに渡って画像を保存できる点があるのかもしれない。 ただ、非常に使いづらい。今回のpic/priの方が圧倒的に使い易い。 追記PAWのメモリーに余裕がある場合pic/priは非常に使い易い。 一方でPAWのメモリーに余裕がない場合に使うと現在のスクリーン全てが画像として保存されず、複数のヒストグラムがあるならそのうちの一部、もしくはヒストグラムの一部だけが保存される。これを回避するには メモリーサイズを増やす必要がある。 ただ、メモリーサイズを増やすとヒストグラムを表示しただけでPAW自体が落ちることがあるのでその点では使いづらい。メモリーに余裕がない場合は単なるハードコピーである.gifを使うか、metafileを使うことにならざるを得ない。ただ、このような場合は.epsや.psにした際にPAW上ではエラーメッセージが表示されないためでできた画像ファイルを見るか印刷でもしないとそもそもメモリー不足に陥っていることに気付かない。注意して使うことが必要。ANALYS/STATUSANALYS/STATUS現在のステータスを表示する。 Scaler, RDF, ANA, Analyser, その他。ana/stでも可。
ANAPAW> ana/status
ANAPAW-M : Current Status. Online Mode Blocks : 0 Valid Events : 0 RDF : ANA : ana/gcheck.ana Analyzer : 1 Analyzer : 3 ANALYS/LOOPANALYS/LOOPRAW-DATAを読み込んでヒストグラムに入れる。a/lでも可。
ANAPAW> ana/loop [ RDFNAME CHOPT NUM REFTIME ]
以下のような表現でイベント数を指定できる。
ANAPAW> ana/loop filename ! 100000 5.
ana/lopした際のコマンドはEVTLOOPを参照。
rdmpRawデータにゲートをかけて圧縮したファイルを生成する。
ANAPAW> rdmp [ FIO DMPNAME GID ]
ANAPAW> book ANAFILE
ANAPAW> rdmp o DMPNAME 1000
ANAPAW> a/l xxxx.rdf
ANAPAW/EVTLOOP> start
ANAPAW> rdmp c
クローズを忘れると圧縮したデータの一部(最後の数イベント)が保存されないらし
い。
PIされた特定の核種のゲートをかけて以降はその核種についてのみ解析を行う。
このようなとき予めrdmpしておくと余計なイベントを解析しないため効率的。
ANALYS/DEFINITIONbookanaファイルを読み込む。
ANAPAW> book [ FNAME ]
clear全てのヒストグラム定義、ヒストグラム、イベントサマリーを消す。
ANAPAW> clear [ CHOPT ]
anaaddhst1ZONEを定義して、一次元ヒストグラムを一度に表示する。 デフォルト(リターンキーを押す、もしくは0を入れる)では2×2で最初の4つ一次元ヒストグラムが表示される。
ANAPAW> hst1 [ NX NY ]
> i
デフォルトでは
===> Histogram List # HID Kind Title -> 1 101 (2) ID vs 2 102 (1) NEUT101 3 103 (1) NEUT102 4 104 (1) NEUT103 5 105 (1) NEUT104 6 106 (1) NEUT105 7 107 (1) NEUT106 8 108 (1) NEUT107 9 109 (1) NEUT108 10 110 (1) NEUT109 ANAPAW/EVTLOOP> hst1 > hst1 HST1> ( また、hst1 3 3なら ANAPAW/EVTLOOP> hst1 3 3 > hst1 3 3 HST1> ( また、表示するヒストグラムのIDを指定することもできる。 ANAPAW/EVTLOOP> hst1 2 3 > hst1 2 3 HST1> ( hst2hst1の2次元版。デフォルトは1×1。
ANAPAW> hst2 [ NX NY ]
ANALYS/HISTOGRAMhstorePAWC上のヒストグラムを保存する。
ANAPAW> hstore [ HBKFILENAME ]
fetchhbookを読み込む。
ANAPAW> fetch [ HBKFILENAME OFFSET ]
htヒストグラムを表示する。
ANAPAW> ht [ HID CHOPT ]
※htではIDが表示されるが、hi/plではIDは表示されない。
hnt次のヒストグラムを表示する。hnでも可。
ANAPAW> hnt [ CHOPT ]
hbt前のヒストグラムを表示する。hbでも可。
ANAPAW> hbt [ CHOPT ]
blow指定した範囲のヒストグラムを作成。
ANAPAW> blow [ BLXMIN BLXMAX BLYMIN BLYMAX ]
ANAPAW>blow -10 15 500 1500
上段がオリジナル。下段がblow
xblow指定した範囲(x軸)のヒストグラムを作成する(マウス使用)。
ANAPAW> xblow
xyblow指定した範囲(x,y軸)のヒストグラムを作成する(マウス使用)。
ANAPAW> xyblow
prx二次元ヒストグラムのx軸への射影。
ANAPAW> prx
pry二次元ヒストグラムのy軸への射影。
ANAPAW> pry
bnx二次元ヒストグラムのx軸への射影(y軸の範囲を指定)。
ANAPAW> bnx [ YMIN YMAX ]
ANAPAW> bnx 500 2000
bnyx軸の範囲を指定してy軸への射影。
ANAPAW> bny [ XMIN XMAX ]
slx二次元ヒストグラムをスライスしてx軸への射影。
ANAPAW> slx [ SNX ]
slyスライスしてy軸への射影。
ANAPAW> sly [ SNY ]
ANAPAW> sly 66
ANAPAW> ht 117 mamiy軸の最大値、最小値を変更。
ANAPAW> mami [ MRYMIN MRYMAX ]
ANAPAW> mami 200 800
xmamiマウスを用いてy軸の最大値、最小値を変更。
ANAPAW> xmami [ HID ]
ANAPAW> xmami 103
mamiの場合は新たなヒストグラムの#は元のヒストグラムの番号+1になり、HIDは全てのヒストグラムの番号+1になる。xmamiの場合は#は全てのヒストグラムの番号+1、 HIDは元のヒストグラムのIDから+1し続けて、存在しない番号にあたった場合にその番号が振られる。普通は両者のHIDは同一になるが、ヒストグラムを消去したり、FORTRANでHBOOKを作る際に飛び飛びのID設定を 行なっているとHIDの設定のされ方が変わる。また、#の振られ方が異なるため、hn,hbの挙動は明確に異なる。hn,hbはHIDではなく#を基準に動くため。 このような挙動の違いはANPAWのソースそのものに入っているコマンドと.kumacがエイリアスされているだけのものとの間で生じ得る。一般に前者はmamiのような挙動、後者はxmamiのような挙動になる。 ===> Histogram List # HID Kind Title (中略) 27 102 (1) Energy loss NEUT102 28 20231 (1) Mami.(Energy loss NEUT102) -> 29 103 (1) Energy loss NEUT103 30 104 (1) Energy loss NEUT104 31 105 (1) Energy loss NEUT105 (中略) 239 20229 (1) position NEUT229 240 20230 (1) position NEUT230 241 131 (1) Energy loss NEUT103 xstatus現在のヒストグラムのステータスを表示する(マウス使用)。ANAPAW> xstatus Histogram ID : 103 TITLE : NEUT101 Channel : From 509.48 to 1266.10 | Bin : From 41 to 102 Mean Value : 870.75 Standard Deviation : 199.54 Entries: 88610 Integrated content: 46686 Ratio: 52.69 % hstatus現在のヒストグラムの定義を表示する。ANAPAW> hstatus Histogram ID : 103 TITLE : NEUT101 Xbin : 200 Xmin : 0.00 Xmax : 2500.00 Mean Value : 795.04 Standard deviation : 498.19 Entries: 88610 Integrated content: 87340 Ratio: 98.57 % hlistDisplay Formatで全てのヒストグラムのリストを表示する。
ANAPAW> hlist [ CHOPT ]
xfitgガウシアンでフィット。 マウスで最大、最小を指定する。
ANAPAW> xfitg
例えば、
ANAPAW> xfitg MINUIT RELEASE 96.03 INITIALIZED. DIMENSIONS 100/ 50 EPSMAC= 0.89E-15 ********** ** 1 **SET EPS 0.1000E-06 ********** FLOATING-POINT NUMBERS ASSUMED ACCURATE TO 0.100E-06 ********** ** 2 **SET ERR 1.000 ********** 1 ********************************************** * * * Function minimization by SUBROUTINE HFITGA * * Variable-metric method * * ID = 103 CHOPT = TUR * * * ********************************************** Convergence when estimated distance to minimum (EDM) .LT. 0.10E+01 FCN= 124.4512 FROM MIGRAD STATUS=CONVERGED 91 CALLS 92 TOTAL EDM= 0.85E-04 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 Constant 997.82 7.1993 2.5243 -0.67120E-03 2 Mean 852.77 2.5446 2.1574 0.18870E-02 3 Sigma 264.69 4.1004 1.0775 -0.31615E-02 CHISQUARE = 0.2894E+01 NPFIT = 46 figaガウシアンでフィット。figでも可。
ANAPAW> figa
例えば、
ANAPAW> fig MINUIT RELEASE 96.03 INITIALIZED. DIMENSIONS 100/ 50 EPSMAC= 0.89E-15 ********** ** 1 **SET EPS 0.1000E-06 ********** FLOATING-POINT NUMBERS ASSUMED ACCURATE TO 0.100E-06 ********** ** 2 **SET ERR 1.000 ********** 1 ********************************************** * * * Function minimization by SUBROUTINE HFITGA * * Variable-metric method * * ID = 103 CHOPT = TU * * * ********************************************** Convergence when estimated distance to minimum (EDM) .LT. 0.10E+01 FCN= 14911.99 FROM MIGRAD STATUS=CONVERGED 73 CALLS 74 TOTAL EDM= 0.22E-04 STRATEGY= 1 ERROR MATRIX ACCURATE EXT PARAMETER STEP FIRST NO. NAME VALUE ERROR SIZE DERIVATIVE 1 Constant 853.16 4.4955 13.758 -0.38430E-03 2 Mean 834.03 1.7995 7.5618 0.28502E-02 3 Sigma 436.42 1.7651 5.3964 -0.36265E-02 CHISQUARE = 0.7608E+02 NPFIT = 199 hdeleteHID1からHID2までのヒストグラムを消去。0 0にすると現在のヒストグラムを消去。
ANAPAW> hdelete [ HID1 HID2 ]
ANAPAW> i ===> Histogram List # HID Kind Title 1 101 (2) NEUT208 dT vs 2 102 (2) ID vs -> 3 103 (1) NEUT101 4 104 (1) NEUT102 5 105 (1) NEUT103 6 106 (1) NEUT104 7 107 (1) NEUT105 8 108 (1) NEUT106 9 109 (1) NEUT107 10 110 (1) NEUT108 11 111 (1) NEUT109 12 112 (1) multiplicity NEUT100 13 113 (1) multiplicity NEUT200 ANAPAW> hdelete 103 111 ANAPAW> i ===> Histogram List # HID Kind Title 1 101 (2) NEUT208 dT vs 2 102 (2) ID vs -> 3 112 (1) multiplicity NEUT100 4 113 (1) multiplicity NEUT200 eraseヒストグラムの中身とイベントサマリーを消去。 bookで読み込んだanaファイルの情報は消えない。
ANAPAW> erase
anaファイルの情報も含めてすべて消すにはclearを使う。
ANAPAW> clear
確認せずに消去する場合はfを付ける。cle fでも可。
Realy Clear!? > ( ANAPAW-M : All definitions were cleared.
ANAPAW> clear f
xval一元ヒストグラムの値を読み取る。
ANAPAW> xval
例
ANAPAW> xval
ID: 103 | X Value: 813.53 | Contents: 1024 xyval二元ヒストグラムの値を読み取る。
ANAPAW> xyval
例
ANAPAW> xyval
ID: 101 | X Value: 1.12 | Y Value: 908.92 | Contents: 15 cpro二元ヒストグラムをカットして新たなヒストグラムを作る(マウス使用)。最後は右クリックで固定。最低3点必要。
ANAPAW> cpro [ HID COPT ]
avy二元ヒストグラムをプロファイルヒストグラムに射影。
ANAPAW> avy [ HID COPT ]
cdrawHCUTで定義したゲートを表示し、数を数える。 カットファイル名は既に読み込まれている場合省略可能。
ANAPAW> cdraw [ CUTID CUTNAME COPT ]
ANAPAW/EVTLOOP> cdraw 1001 > cdraw 1001 Histogram ID : 101 TITLE : NEUT208 dT vs Xbin : 100 Xmin : -21.00 Xmax : 21.00 Ybin : 250 Ymin : 0.00 Ymax : 2500.00 GID : 1001 KIND : XYGate VAL_X : 1 8 8 11 VAL_Y : 1 8 8 14 Number of Points : 12 Integrated content : 46546 ACCEPTED ENTRIES : 46526 Mean ( 8.0056 , 934.5109 ) STD ( 6.0056 , 208.9780 ) chbinビンを変える。
ANAPAW> chbin [ BF ]
ANAPAW> chbin 5
> chbin 5 Current Histogram ID : 118 gpro二次元ヒストグラムにゲートをかけて表示する。
ANAPAW> gpro [ CUTID CUTNAME COPT ]
ANAPAW/EVTLOOPana/loopした際にEVTLOOPになる。list//PAWC上のヒストグラムのリストを表示する。iとほぼ同じ。ANAPAW/EVTLOOP> list ===> Directory : //PAWC 101 (1) Energy loss NEUT101 102 (1) Energy loss NEUT102 103 (1) Energy loss NEUT103 104 (1) Energy loss NEUT104 105 (1) Energy loss NEUT105 startloopをスタートする。Enterで止める。ANAPAW/EVTLOOP> start 2011年 5月 8日 日曜日 01:48:03 JST ANAPAW-M : Hit CR to stop. ANAPAW-M : Interrupt Event Loop! quitRDFを閉じてイベントループを終了する。qでも可。updateupdateのON/OFFを切り替える。
ANAPAW/EVTLOOP> update
UPDATE> ON UPDATE> OFF
|