PEGS5の入力データ
2019年10月1日
KEK 波戸、平山

PEGS5入力ファイルの例(sampl5.inpの一部)

ELEM&INP /
FE-RAYLEIGH
FE
FE
ENER
&INP AE=1.5,AP=0.100,UE=2000.0,UP=2000.0 &END
TEST
&INP /PWLF&INP /DECK&INP /PEGS5入力ファイルの例(sampl5.inpの一部)

ELEM&INP /
元素組成等指定部分
FE-RAYLEIGH
FE
FE
ENER
エネルギー範囲&INP AE=1.5,AP=0.100,UE=2000.0,UP=2000.0 &END 指定部分
TEST
&INP /PWLF
定型部分&INP /(通常変更せず)
DECK&INP /

元素組成等指定部分
5種類の例で説明

単体

化合物

混合物

固体、液体鉄

アクリル

鉛ガラス

気体

Xe ガス空気鉄：単体・固体&INP / : オプション入力部「変数名=値」の書式で入力

単体(左詰)ELEM&INP /FE

FEFE
元素記号(左詰)

物質の名称 (ユーザーコードで使用。左詰24文字以内)

密度効果Flag #(31文字目から書く)# : 密度効果表 (表2.3) に列挙された物質を指定する– Sternheimer-Seltezer-Bergerの式 (高精度)– 省略時：Sternheimer-Peierlsの一般式(中精度)

Xe ガス: 単体・気体
1気圧での密度 (g/cm3)“GASP=Gas 圧力”は、「物質を 気体 として指定」

ELEM&INP RHO=5.89E-3, GASP=1.0 /XENON-GAS-STP

XE-GASXE
気体の問題点：密度が圧力に比例。圧力が温度に比例。
RHOとGASPは共通の温度(0℃, 20℃など) に対する値を入力。
参考: 標準状態（STP)=0ºC, 1気圧 ・ 自然状態(NTP)=20℃, 1気圧

アクリル(PMMA): 化合物・固体
元素の数

化合物
(左詰)

密度 (g/cm3)

原子の数の比COMP&INP NE=3, RHO=1.19, PZ=5.0, 8.0, 2.0 /PMMAC
元素記号(A2,1X)
2文字1空白or
1文字2空白

H

O

同じ順番！

鉛ガラス: 混合物・固体混合物
(左詰)

元素の数

密度 (g/cm3)

元素の質量比MIXT&INP NE=5, RHO=3.61, RHOZ=41.8, 21.0, 29.0, 5.0, 2.2 /LEAD GLASSPB SI O

元素記号
(A2,1X)
2文字1空白or
1文字2空白K

NA同じ順番！
<化合物と混合物の使い分け>
原子の数の比で指定したい→化合物
元素の質量比で指定したい→混合物空気 (20 ºC, 1気圧): 混合物・気体
0℃ 1気圧での密度 (g/cm3)

20℃、1気圧を0℃にした
場合の気圧

MIXT

&INP NE=3, RHO=1.2929E-3, GASP=0.93174,
RHOZ=0.75575,0.23143,0.01282 /AIR-20CNO

2文字1空白or1文字2空白AIR-GAS
AR

同じ順番！
RHO=1.205E-3, GASP=1.0でも可

N=78%は容積比であり、間違いエネルギー範囲指定部分
ENER&INP AE=1.5, AP=0.100, UE=2000.0, UP=2000.0 /光子電子• 入射エネルギーEIとの関連– 光子入射

上限

UP

UE

下限

AP

AE

MeV単位
電子のエネルギーは静止質量を含む• AP ≦ EI≦UP, EI≒UP• AE ≦ EI+RM≦UE– 電子入射• AE ≦ EI ≦ UE, EI≒UE• AP ≦ EI-RM ≦ UP– 陽電子入射（電子入射の条件に追加）• AP ≦ EI+RM ≦ UP

• AEとAP,UEとUPの通常の設定– AE-AP=RM, UE-UP=RM

[EOF]の位置に関連するエラー集• DECK&INP /END

[EOF]

[EOF]が最終行にある場合：
正常動作• DECK&INP /END [EOF]

EOF:
ファイルの終わりを示す記号[EOF]が/ENDと同じ行にある場合：異常終了

egs5job.out: PEGS5-call comes next
pgs5job.pgs5lst: Stopped in pegs5 because namelist/INP/ data was missing.

•DECK&INP /END
(空白行)
[EOF]

・<空白> [EOF][EOF]の前に空白行がある場合：
異常終了[EOF]が最終行にあるものの、[EOF]のまえに空白などがある場合：
異常終了

CALL オプション
ELEM.

&INP /
PB
PB

フィッティング前の量を出力するオプションCALL&INP XP(1)=49.99 /49.99 MeV 光子に対する鉛の平均自由行程を出力

GMFP
<出力>
OPT=CALLFUNCTIONCALL: 1.95522 = GMFP OF 49.99

放射長単位低エネルギー光子輸送関連のフラッグ• IRAYL=1 (レイリー散乱断面積)• IBOUND =1 (束縛電子コンプトン断面積)• INCOH=1 (束縛コンプトン散乱角度分布)• ICPROF=-3 (ドップラー広がり)• IMPACT=1-6 (K殻電子衝突電離)• すべて、0で無視（既定値）
他の多くのオプション、機能については
pegs_user_manual.pdf を参照

pegs5の初歩的な練習問題• uc_examin.f を egs5run を用いて走らせよ. egs5job.outの各列の数値のエネルギー依存を説明せよ。 (グラフ化を推奨） egs5/extra_ucodes/uc_examin/ucexamin.f参照• uc_examin.inpの2行目のIMPACT=1,の後にIUNRST=2を追加し、非制限全阻止能を出力せよ。それを文献値(例えば、NISTのestarと比較せよ。）

pegs5の中級練習問題• Feに対する pegs5 入力ファイルを作成し、 uc_examin.f を走らせよ。元の uc_examin.inp と比べて何が変わったか?
• 水に対する pegs5 入力ファイルを作成し、 uc_examin.f を走らせよ。 外部の参考書などの数値と比較するなどして数値の妥当性を検討せよ。

uc_examin.fを変更なしに走らせる場合の入力• Material Identifier: AL<Enter>• Title after: <Enter>• Electron and photons (0), only electrons (1) orphotons (2): 0<Enter>• Include Rayleigh (coherent) scattering (1) or not (0):1<Enter>• Tables (0) or individual energies (1) : 0<Enter>• Output terminal (0), disk file (1), terminal&plot (2),disk&plot (3) : 0<Enter>• 結果はegs5job.outに出力される。

PEGS5入力とユーザーコードの関係
PEGS5入力*

PEGS5

物質データ
ユーザーコード*EGS5* egs5runで指定計算がうまく実行されない場合には、これが作成されているかどうかを調べる。数100行以上作成されていればOK。空白または数行なら失敗。

改訂記録• 22JUL2004 EGS5用記述• 05AUG2005 charDに対応してEFRACはオプション化• 2006-05-01 EFRACを削除。空気組成の誤った比率（容積比）の例を追加。
• 2006-06-21 Efrac廃止を追加。練習問題を追加。[EOF]の位置に関連するエラーを追加• 2008-08-01 efracに関する記述を削除。
• 2009-03-09 EOFの位置に関連するエラーを1個追加• 2009-07-28 uc_examin.fの入力を追加• 2009-7-29 「pegs5の初歩的な問題」を改訂• 2010-7-21 吹き出し記号を利用して記述を簡略化• 2012-6-19 空気の説明文で、RHO=1.0 → GASP=1.0と訂正• 2019.10.1 ucexamin.fの場所を記載

以下は予備のスライド

CO2 ガス(20 ºC, 1気圧) : 化合物・気体
COMP&INP NE=2, RHO=1.977E-3, GASP=0.93174, PZ=1.,2. /CO2-20C

CCO2-GASO• GASP: 0.93174 atm (=273ºC/293ºC).
• この圧力は、 20ºC、1気圧のガスを体積一定のまま0ºCに冷却したときに得られる。
• CO2-GAS：密度効果表でのCO2ガスの名称

H2 Gas (STP): 化合物, 気体
COMP&INP NE=2, RHO=8.99E-5, GASP=1.0, PZ=1.,1. /H2-GAS-STP

HH2-GASH• 単元素の分子気体 (例. H2) は化合物として扱う• NE=1 はエラーになる• H2-GAS (31 列目) : 密度効果表での水素ガスの名称