マニュアル
Ver2.3Ver2.3

波戸 芳仁、平山 英夫
高エネルギー加速器研究機構〒305-0801 茨城県つくば市大穂 1-1
高村 篤、杉田 武志(有)科学システム研究所〒309-17 茨城県笠間市住吉 1342-6

1

目

次
頁

第 1 章 概要 ・・・・・・・・・・・・・・・・・・・・・・・・・・・

5

第 2 章 動作環境 ・・・・・・・・・・・・・・・・・・・・・・・・・

6

第 3 章 機能概要 ・・・・・・・・・・・・・・・・・・・・・・・・・

7

3.1 体系表示に関する機能・・・・・・・・・・・・・・・・・・・

7

3.2 粒子の飛跡表示に関する機能 ・・・・・・・・・・・・・・・・

9

3.3 マウス操作に関する機能 ・・・・・・・・・・・・・・・・・・

10

3.4 付加機能 ・・・・・・・・・・・・・・・・・・・・・・・・・

11

3.5 体系データ作成機能 ・・・・・・・・・・・・・・・・・・・・

12

3.6 体系の整合性確認機能・・・・・・・・・・・・・・・・・・・

12

3.7 領域体積計算機能 ・・・・・・・・・・・・・・・・・・・・・

12

第 4 章 インストール ・・・・・・・・・・・・・・・・・・・・・・・

13

4.1

WINDOWS 版 ・・・・・・・・・・・・・・・・・・・・・・・・

13

4.2

Linux 版 ・・・・・・・・・・・・・・・・・・・・・・・・・

14

第 5 章 起動方法 ・・・・・・・・・・・・・・・・・・・・・・・・・

15

5.1

WINDOWS 版 ・・・・・・・・・・・・・・・・・・・・・・・・

16

5.2

Linux 版 ・・・・・・・・・・・・・・・・・・・・・・・・・

16

第 6 章 基本操作 ・・・・・・・・・・・・・・・・・・・・・・・・・

17

6.1 データファイルの読み込み ・・・・・・・・・・・・・・・・・・

17

6.2 データファイルの再読込 ・・・・・・・・・・・・・・・・・・・

17

6.3 表示データの保存 ・・・・・・・・・・・・・・・・・・・・・・

18

6.4

DOS 画面を開く ・・・・・・・・・・・・・・・・・・・・・・・

19

6.5 画面の印刷 ・・・・・・・・・・・・・・・・・・・・・・・・・

20

6.6 使用方法のヘルプを表示 ・・・・・・・・・・・・・・・・・・・

20

6.7 アプリケーションの終了方法 ・・・・・・・・・・・・・・・・・

20

第 7 章 表示の設定 ・・・・・・・・・・・・・・・・・・・・・・・・

21

7.1 座標軸の表示・非表示を変更する ・・・・・・・・・・・・・・・

21

7.2 表示する座標軸を変更する ・・・・・・・・・・・・・・・・・・

21

7.3 投影方法を変更する ・・・・・・・・・・・・・・・・・・・・・

22

7.4 体系の表示タイプを変更する ・・・・・・・・・・・・・・・・・

22

7.5 動画表示 ・・・・・・・・・・・・・・・・・・・・・・・・・

23

7.6 体系表示に関する設定を変更する ・・・・・・・・・・・・・・・

24

7.6.1 表示する体系レベルを変更する ・・・・・・・・・・・・・・

25

7.6.2 表示する形状番号、領域名、物質番号を選択する ・・・・・・

25

7.6.3 表示色を変更する ・・・・・・・・・・・・・・・・・・・・

25

2

7.6.4 線幅を変更する ・・・・・・・・・・・・・・・・・・・・・

25

7.6.5 領域作成時の補正オプション ・・・・・・・・・・・・・・・

26

7.7 粒子の表示に関する設定を変更する ・・・・・・・・・・・・・・

27

7.7.1 ヒストリー範囲の指定 ・・・・・・・・・・・・・・・・・・

27

7.7.2 表示する粒子の指定・・・・・・・・・・・・・・・・・・・

28

7.7.3 粒子の線色を変更する ・・・・・・・・・・・・・・・・・・

28

7.7.4 粒子の線種を変更する ・・・・・・・・・・・・・・・・・・

28

7.7.5 粒子の線幅を変更する ・・・・・・・・・・・・・・・・・・

28

7.7.6 エネルギー範囲を指定する ・・・・・・・・・・・・・・・・

29

7.7.7 WEIGHT について ・・・・・・・・・・・・・・・・・・・・・

29

7.7.8 LATCH 機能 ・・・・・・・・・・・・・・・・・・・・・・・

29

7.7.9 時間処理指定 ・・・・・・・・・・・・・・・・・・・・・・・ 29
7.8 画面表示に関する設定を変更する ・・・・・・・・・・・・・・・

30

7.8.1 背景色を変更する ・・・・・・・・・・・・・・・・・・・・

30

7.8.2 文字色を変更する ・・・・・・・・・・・・・・・・・・・・

30

7.8.3 座標軸の設定を行う

・・・・・・・・・・・・・・・・・・・

31

7.8.4 タイトルを表示する

・・・・・・・・・・・・・・・・・・・

31

7.8.5 メモを表示する ・・・・・・・・・・・・・・・・・・・・・

31

7.8.6 粒子の凡例を表示する ・・・・・・・・・・・・・・・・・・

32

7.8.7 ヘルプ表示用ブラウザーを設定する ・・・・・・・・・・・・

32

7.9 部分削除 ・・・・・・・・・・・・・・・・・・・・・・・・・・

33

7.10 設定値を初期値に戻す ・・・・・・・・・・・・・・・・・・・・ 33
第 8 章 マウス操作に関する機能 ・・・・・・・・・・・・・・・・・・

34

8.1 表示体系を回転させる

・・・・・・・・・・・・・・・・・・・・

34

8.2 表示体系を移動させる・・・・・・・・・・・・・・・・・・・・

34

8.3 表示体系を拡大させる

・・・・・・・・・・・・・・・・・・・・

34

8.4 表示体系を縮小させる

・・・・・・・・・・・・・・・・・・・・

34

第 9 章 体系データを作成する ・・・・・・・・・・・・・・・・・・・

35

9.1 体系データ作成画面を開く ・・・・・・・・・・・・・・・・・・

35

9.2 体系データを入力する

・・・・・・・・・・・・・・・・・・・・

36

9.3 作成したデータを保存する ・・・・・・・・・・・・・・・・・・

37

9.4 作成したデータを開く

・・・・・・・・・・・・・・・・・・・・

37

9.5 作成したデータの表示・・・・・・・・・・・・・・・・・・・・

37

9.6

Excel でファイルを編集する ・・・・・・・・・・・・・・・・・

38

9.7

CG 形状一覧 ・・・・・・・・・・・・・・・・・・・・・・・・

39

第 10 章 整合性確認 ・・・・・・・・・・・・・・・・・・・・・・・・

349

第 11 章 領域体積計算 ・・・・・・・・・・・・・・・・・・・・・・・

52

11.1 分割レベルを指定する ・・・・・・・・・・・・・・・・・・・・ 53
11.2 評価位置を指定する・・・・・・・・・・・・・・・・・・・・

53

11.3 収束レベルを指定する ・・・・・・・・・・・・・・・・・・・・ 53
11.4 結果の表示 ・・・・・・・・・・・・・・・・・・・・・・・・・ 53
11.5 計算状況の表示 ・・・・・・・・・・・・・・・・・・・・・・・ 53
11.6 ファイルに保存する・・・・・・・・・・・・・・・・・・・・・ 53

11.7 計算の開始 ・・・・・・・・・・・・・・・・・・・・・・・・・ 54
11.8 閉じる ・・・・・・・・・・・・・・・・・・・・・・・・・・・ 54
11.9 ヘルプの表示 ・・・・・・・・・・・・・・・・・・・・・・・・ 54
第 12 章 ファイル形式 ・・・・・・・・・・・・・・・・・・・・・・・

55

12.1 軌跡データの構成 ・・・・・・・・・・・・・・・・・・・・・・ 55
12.2 軌跡データ(フリーフォーマット)の形式 ・・・・・・・・・・・・ 58
12.3 軌跡データ(CSV フォーマット)の形式 ・・・・・・・・・・・・・

58

12.4 軌跡データ(固定フォーマット)の形式 ・・・・・・・・・・・・・ 58
12.5 体系ファイルの構成

・・・・・・・・・・・・・・・・・・・・・ 59

12.6 体系ファイル(フリーフォーマット)の形式 ・・・・・・・・・・・ 60
12.7 体系ファイル(CSV フォーマット)の形式 ・・・・・・・・・・・・
第 13 章

EGS ユーザーコード ・・・・・・・・・・・・・・・・・・・・

60
61

13.1 メインプログラム ・・・・・・・・・・・・・・・・・・・・・・ 61
13.2 SUBROUTINE AUSGAB ・・・・・・・・・・・・・・・・・・・・・
13.3 SUBROUTINE PLOTXYZ
第 14 章63

・・・・・・・・・・・・・・・・・・・・・ 63

CG を用いた EGS 用計算体系作成例 ・・・・・・・・・・・・・・

4

64

第 1 章 概要
軌跡とジオメトリー表示プログラム CGVIEW は電子・光子・陽電子の 3 次元軌跡と計算ジ
オメトリーを PC 上で表示するためのシステムである。粒子の軌跡は EGS コードシステムで計算される。
EGS の軌跡表示システムは、最初に IBM 大型機の時代に専用グラフィック端末用のものが開発され、次に UNIX 上でグラフィックシステムの一つである PHIGS を利用した「EGS4_Windows_2」が作成された。これらのシステムは電磁カスケードシャワーの理解に大変有用なツールであったが、専用のハードまたはソフトが必要でありどこでも使用できる状況ではなかった。PC 用の EGS 軌跡表示システムは 1994 年に MS-DOS 用が高エネルギー
加速器研究機構(以後、KEK)で開発され、Windows 対応などのバージョンアップを重ねながら広く使用されてきた。しかしながらこのシステムで表示可能なジオメトリーは円筒と平板のみに限られていた。
一方、複雑なジオメトリーを比較的少ない労力で記述できる Combinatorial Geometry (以後、CG)と呼ばれるジオメトリーパッケージの EGS 計算コードでの利用が進められてきた。
CGVIEW は、粒子の軌跡を CG で記述されたジオメトリーとともに 3 次元表示するためのものであり、CGVIEW を用いると、意図したとおりのジオメトリーが入力されているかどうかを
3 次元表示で確認することができるため、計算準備の効率が大幅に上がることが期待できる。
このマニュアルの 2 章から 10 章では CGVIEW の動作環境、インストール手順、使用方法などについて述べる。ユーザーはこの部分を読みながら、サンプルデータを用いて CGVIEW
を動かすことによりその機能を知ることができる。11 章では EGS から CGVIEW に送る軌跡・
ジオメトリーデータの構造について、12 章ではそのデータを作成するための EGS のユーザ
ーコードについて述べる。

5

第 2 章 動作環境
ハードウエア ：PC/AT 互換機

ＯＳ

Windows2000/XP、Linux + X-Window システム(推奨 RedHat 8.0)：

ビデオカード ：OpenGL 対応のもの
但し、ドライバーのバージョンが古い場合 OpenGL に対応していない可能性があります。うまく動作しない場合は使用しているビデオカー
ドのメーカーから最新版をダウンロードしてください。

画面の色

：

32 ビット(TRUE COLOR)を推奨。
但し、ビデオカードによっては正しく表示できない場合も有るのでその場合は、16 ビットに変更してみる。

OpenGL：

Linux の場合
ディストリビューションによっては収録されていない場合があります。その場合は OpenGL 類似の Mesa ライブラリをインストールしてください。
ダウンロード先：http://www.mesa3d.org/

6

第 3 章 機能概要
3.1 体系表示に関する機能
1)表示できる体系データは以下の 4 通りとする。
・同心球
・同心円筒
・3 次元平板・CG(Combinatory Geometry)2)表示方法は以下の 2 通りとする。(デフォルト値は面表示)
・面表示
・線表示
3)表示レベル
表示レベルは以下の 3 通りとし、表示するレベルを選択できる。
・形状 ・・ 入力した形状毎に表示の ON/OFF と色指定が可能
・領域 ・・ 領域(Zone)毎に表示の ON/OFF と色指定が可能
・物質 ・・ 物質毎に表示の ON/OFF と色指定が可能
4)部分切取表示
仮想 BOX を入力することで仮想 BOX と重なる部分を非表示にすることができる。
入力するデータは X 方向の最小値と最大値、Y 方向の最小値と最大値、Z 方向の最小値と最大値。空白の場合表示しているデータから自動で最大値、最小値を決定。
5)表示範囲表示範囲は入力されたデータから判断する。
6)表示する座標軸の変更
表示する座標軸以下の 4 通りから選択できる。但し、座標系は右手系とする。
・Z-Y ・・ 右方向に Z 軸プラス、上方向に Y 軸プラス(デフォルト)・Y-Z ・・ 右方向に Y 軸プラス、上方向に Z 軸プラス
・X-Z ・・ 右方向に X 軸プラス、上方向に Z 軸プラス
・X-Y ・・ 右方向に X 軸プラス、上方向に Y 軸プラス7)背景色の変更
背景色の変更が可能(デフォルトは黒)8)画面文字色の変更
座標軸、タイトル、メモ等の画面上に表示される文字の色を変更可能(デフォルトは白)
9)線幅の変更
表示する線の太さを変更する。線表示の場合のみ。
10)領域作成時の補正オプション
領域作成の際、コンピュータ精度の問題から、正しくデータが表示されない場合が7

ある。特にジオメトリーの論理差をとる場合。そのため、論理差を行う場合のみ差をとる方のジオメトリーを大きく補正する。補正量は以下の項目と以下の式より算出する。
・形状補正率・・寸法に対しての割合で補正量を決める(デフォルト 1E-6)・形状補正値・・直接補正量を指定(デフォルト 0.005cm)
補正後の寸法=補正前寸法+補正前寸法*形状補正率+形状補正値

8

3.2 粒子の軌跡表示に関する機能
1)表示できる粒子の種類は以下の 3 通りとする。
・光子(γ、x)・電子(ｅ-)
・陽電子(ｅ+)2)軌跡の表示に条件を指定することができる。
・ヒストリー範囲・・ ヒストリーの範囲の上限値と下限値を入力することで、その範囲外のヒストリーを非表示にする。
空白の場合は下限値、上限値ともにデータから最小ヒストリ
ー番号、最大ヒストリー番号が割り当てられる。また、最後のヒストリー番号からの個数を入力することで最後のヒス
トリー番号から入力した個数分の粒子のみ表示できる。

・線色

・線種

・・ 粒子毎に設定。デフォルト値を以下に示す。
光子

－

黄

電子

－

赤

陽電子

－

青

・・ 粒子毎に以下の 4 つから選択。デフォルト値は全て実線［実線、破線、一点鎖線、二点鎖線］

・線幅

・・ 粒子毎に線幅を指定できる。デフォルト値は 1

・表示指定・・ 粒子毎に ON/OFF を設定できる。

・エネルギー範囲

・・ 粒子の運動エネルギー範囲を指定することで指定範囲外の粒子を非表示にする。空白の場合下限値に 0、上限値に1E+15 が割り当てられる。

3)粒子の WEIGHT
各粒子の WEIGHT 情報から表示する粒子の軌跡の太さを変化させることが出来る。表示
する軌跡の太さは各粒子で指定した線幅×WEIGHT になる。WEIGHT 情報の無いデータを読み込んだ場合は全て 1 が割り当てられる。
4)LATCH
入力データでヒストリー毎にフラグを与え、表示の ON/OFF を制御できる。また、LATCH機能を使用するかどうかの設定も出来る。
5)動画表示粒子の生存時間で分割した動画表示が出来る。

9

3.3 マウス操作に関する機能
マウス操作による機能は以下の 4 通りがある。
1)拡大縮小・・ マウス右ボタンを押したままマウスを移動させる。
ボタンを押したポイントから下方向で拡大、ボタンを押したポイントから上方向で縮小。

2)移動・・ マウス右左ボタンの同時押し、もしくは中ボタンを押したままマウスを移動させる。

3)回転・・ マウス左ボタンを押したままマウスを移動させる。

4)範囲指定して拡大・・ Shift キーとマウス右ボタンを押したままマウスを移動させる。移動させると白い枠線が表示されるので範囲を
指定する。

10

3.4 付加機能
1)軸表示に関する機能
・軸表示の ON/OFF
・目盛り単位

・・ 表示する座標軸の目盛り単位。(デフォルト値=cm)・目盛り間隔・・ 目盛線の間隔。(デフォルト値=1)・目盛りサイズ・・ 軸に表示される目盛線のサイズ。(デフォルト値=0.25)

・矢印サイズ・・ 軸に表示される矢印の大きさ。(デフォルト値=1)

・文字サイズ・・ 目盛りラベル文字の大きさ。(デフォルト値=1)・文字間隔

・・ 目盛りラベル文字を表示する間隔。(デフォルト値=5)2)タイトルボックス表示に関する機能
表示した体系図にタイトルボックスを表示できる。タイトルボックスの内容は各自入力する。また、タイトルボックスの位置はマウスで操作することができる。
(複数行入力可)3)メモボックス表示に関する機能
表示した体系図にメモボックスを表示できる。メモボックスの内容は各自入力する。
また、メモボックスの位置はマウスで操作することができる。(複数行入力可)
4)凡例表示に関する機能
粒子の凡例を表示できる。表示する内容は以下の通り。
・表示しているヒストリー範囲
・各粒子の線種と線色
5)画面印刷に関する機能
現在表示している画面を印刷することができる。(Windows 版のみ)
6)画面図形表示
現在表示している画面をビットマップ形式か JPEG 形式にて保存できる。
7)DOS 画面の起動
メイン画面から DOS のウインドウを起動することができる。Linux 版の場合は X 端末を
起動する。
8)再読込
表示しているデータファイルを再読込できる。EGS の計算終了後に再度読み込み、変化を確認できる。
9)設定の初期化
表示のために設定した値を初期値に戻すことが出来る。

11

3.5 体系データ作成機能
EGS コードの計算に使用する計算体系の作成を補助する。使用する体系データは CG とする。EGS コードで使用する体系データを入力し、画面上に表示することができる。また、作成したデータをファイルに保存し、EGS コードの入力データに使用することができる。
3.6 体系の整合性確認機能
反応の無い粒子を飛ばして通過した領域の定義が正しく行われているかをチェックする。チェック用粒子の発生箇所と方向を指定し、発生点と最初の境界との中点及び、形状の境界と境界の中点をチェック点として領域の定義が正しいかどうか確認する。検出
するエラーはチェック点での領域の未定義や重複定義、一番外側で物質番号が 0 以外の場合とする。

図 3.6 整合性確認方法
3.7 領域体積計算機能
各領域の体積を計算することができる。計算は以下の方法を用いて行う。
２のべき数で等間隔に分割した yz 面、xz 面、xy 面の中心から垂直に粒子を飛ばし、領域を横切る長さと分割した面の面積を掛け算し、それらを集計することで領域の体積を求めます。値の確かさは、yz 面、xz 面、xy 面の３つ値の収束レベル、２のべき数の分割レベル間の収束レベルで判断します。２のべき数での等間隔の分割は、べき数が 0 から指定した数値まで行います。

12

第 4 章 インストール
4.1

WINDOWS 版

1) インストール媒体内の cgview-***.exe(***はバージョン番号)をダブルクリックする。
2) 図 4.1 のような自己解凍ソフトが起動するので解凍先を指定し実行します。デフォル
トでは c:\program files\cgview に解凍されます。

図 4.1 解凍画面

13

4.2

Linux 版

ここでは CD-ROM から/usr/local/cgview にインストールするものとして説明を行います。
1)インストールするディレクトリを作成し移動する。
[ / ]# mkdir /usr/local/cgview (ディレクトリを作成)[ / ]# cd /usr/local/cgview (ディレクトリに移動)
2)CD-ROM 内の cgview-***.tar.gz(***はバージョン番号)をコピーし解凍します。
[ /usr/local/cgview ]# mount /mnt/cdrom (CDROM をマウント、環境によって異なる)[ /usr/local/cgview ]# cp /mnt/cdrom/cgview-***.tar.gz ./ (ファイルコピー)[ /usr/local/cgview ]# tar zxvf cgview-***.tar.gz (解凍実行)
3)root 権限が有る場合は 4)を、無い場合は 5)、6)を実行します。
4)共有ライブラリを/usr/lib に作成し、シンボリックリンクを作成します。
[ /usr/local/cgview ]# cp ./libborqt-6.9.0-qt2.3.so /usr/lib[ /usr/local/cgview ]# cd /usr/lib[ /usr/lib ]# ln -s libborqt-6.9.0-qt2.3.so libborqt-6.9-qt2.3.so5)同一ディレクトリ内にシンボリックを作成します。
[ /usr/local/cgview ]# ln –s libborqt-6.9.0-qt2.3.so libborqt-6.9-qt2.3.so6)LD_LIBRARY_PATH 変数を Linux 起動時に設定するようにホームディレクトリ(~)内の.bashrc もしくは.cshrc ファイルの最後に以下の行を追加します。
(.bashrc の場合) export(.cshrc の場合)

LD_LIBRARY_PATH=$LD_LIBRARY_PATH:./setenv LD_LIBRARY_PATH '$LD_LIBRARY_PATH':./

注)Linux の起動時に設定しない場合は CGVIEW 起動時に毎回上記コマンドを実行する必要があります。

14

第 5 章 起動方法
本ソフトが正常に起動すると図 5.1 のようなメイン画面が表示されます。

図 5.1 メイン画面

15

5.1

WINDOWS 版
マイコンピュータを開き、インストール先フォルダ内の cgview.exe をダブルクリック

します。
5.2

Linux 版
インストールしたディレクトリ内にある cgview を実行します。(/usr/local/cgview に

インストールされている場合)[ /usr/local/cgview ]# ./cgview &(注)CGVIEW は文字コードが euc で作成されているため、X-Window のデフォルトが euc 以
外の場合(FedoraCore 等)は起動時に変換処理を行わないと文字化けします。したがって起動コマンドの前に以下のコマンドを追加してください。
[ /usr/local/cgview ]# LANG=ja_JP.eucJP ./cgview &16

第 6 章 基本操作
6.1 データファイルの読み込み
データを読み込むにはメインメニューの[[ファイル]
ファイル]-[体系・体系・飛跡データ
飛跡データ読込
データ読込]
読込]を選択します。メニューを選択するとファイルオープンダイログが表示されるのでファイルダイ
アログの中から読み込むファイルを選択します。読み込むファイルの形式は第 11 章ファ
イル形式を参照してください。
注)読み込める軌跡データの数には制限が設けられており指定したヒストリー番号まで
しか読み込むことができません。デフォルトは 1000 になっていますが変更する場合は同一フォルダ内にある CGVIEW.ini を編集してください。
編集箇所) [Particle]以下の ReadMax=1000 <-この値を増減させます。(10 にした場合は11 から上のヒストリー番号を持つ粒子は無視されます)

図 6.1 [体系・飛跡データ読込]メニュー
6.2 データファイルの再読込
再計算を行った場合などで現在読み込んでいるデータに変更が加えられた場合、再度読み込むことが出来ます。再度読み込みなおすにはメインメニューの[[ファイル]
ファイル]-[再読込]を選択します。

図 6.2 [再読込]メニュー

176.3 表示データの保存表示されたデータの画面をキャプチャーして保存することができます。保存する形式
はビットマップ(拡張子 bmp)か JPEG(拡張子 jpg)になります。メインメニューの[[ファイル]-[名前を名前を付けて保存
けて保存]
保存]を選択します。選択するとファイルを保存するためのファイル
ダイアログが表示されます。そのファイルダイアログでファイル名を指定すると、表示
画面が保存されます。6.5 節の[[ファイル]
ファイル]-[印刷]印刷]では画像を印刷し、この処理では画像のファイル保存をします。

図 6.3 [名前を付けて保存]メニュー

18

6.4

DOS 画面を開く(X 端末を開く)
DOS 画面上で EGS 計算などを実行させたい場合にメインメニューの[[ファイル]
ファイル]-[DOS 画

面を開く(Linux は X 端末を端末を開く)]を選択すると図 6.4-2 のようなコマンドプロンプトウ
ィンドウを呼び出します。

図 6.4-1 [DOS 画面を開く]メニュー

図 6.4-2 DOS 画面
ここで cygwin.bat を実行すると、DOS 窓が cygwin 窓に変わり、cygwin 上にインストー
ルされた egs5 を実行することができます。cygwin.bat が存在するディレクトリにパスが通っていない場合には、例えば、
c:\cygwin\cygwin.batと、cygwin のインストールディレクトリに対応した入力を行ってください。

19

6.5 画面の印刷
現在表示している画面を印刷します。メインメニューの[[ファイル]
ファイル]-[印刷]印刷]を選択します。プリンター選択ダイアログが表示されますので、出力先のプリンターを選択すると、印刷が実行されます。(WINDOWS 版のみ有効な機能です)

図 6.5 [印刷]メニュー
6.6 使用方法のヘルプを表示
使用方法についてのヘルプを表示したい場合は、メインメニューの[[ヘルプ]ヘルプ]-[使用方法]を選択します。各項目についてのヘルプが表示されます。必要に応じて閲覧・印刷してください。

図 6.6 [ヘルプ]メニュー
6.7 アプリケーションの終了方法
メインメニューの[[ファイル]
ファイル]-[終了]終了]を選択します。

図 6.7 [終了]メニュー20

第 7 章 表示の表示の設定
7.1 座標軸の表示、非表示を変更する
・座標軸を表示にする→ メインメニューの[[座標軸表示]
座標軸表示]-[ON]を選択します。
[ON]・座標軸を非表示にする → メインメニューの[[座標軸表示]
座標軸表示]-[OFF]を選択します。
[OFF]

図 7.1 [座標軸表示]メニュー
7.2 表示する座標軸を変更する
座標軸を指定することで視点を移動し、表示する座標の方向を変更します。
・ウインドウに対して右側を Z 軸プラス、上側を Y 軸プラスにする。
→

メインメニューの[[座標軸]座標軸]-[Z[Z-Y]を選択します。
Y]・ウインドウに対して右側を Y 軸プラス、上側を Z 軸プラスにする。
→

メインメニューの[[座標軸]座標軸]-[Y[Y-Z]を選択します。
Z]・ウインドウに対して右側を X 軸プラス、上側を Z 軸プラスにする。
→

メインメニューの[[座標軸]座標軸]-[X[X-Z]を選択します。
Z]・ウインドウに対して右側を X 軸プラス、上側を Y 軸プラスにする。
→

メインメニューの[[座標軸]座標軸]-[X[X-Y]を選択します。
Y]

図 7.2 [座標軸]メニュー

21

7.3 投影方法を変更する
投影方法には以下の 2 種類があります。
・平行投影

→

投影面に向かう投射線は平行であり、遠くの形状も近くの形状も同じ
大きさで表示されます。平行投影にするにはメインメニューの[[投影方法]-[平行投影]
平行投影]を選択します。

・透視投影

→

投影物から視点に向かって投射線が集中する方式。透視投影は人が物を見る場合と同様に、遠くの物は小さく表示され、近くの物ほど大きく表示され、遠近感を出すことが出来ます。透視投影にするにはメイ
ンメニューの[[投影方法]
投影方法]-[透視投影]
透視投影]を選択する。

図 7.3 [投影方法]メニュー
7.4 体系の表示タイプを変更する
表示する体系の表示の仕方を変更する。
・面表示にする → メインメニューの[[表示タイプ
表示タイプ]
タイプ]-[面表示]
面表示]を選択する。
・線表示にする → メインメニューの[[表示タイプ
表示タイプ]
タイプ]-[線表示]線表示]を選択する。

図 7.4 [表示タイプ]メニュー

22

7.5 動画表示
この機能は、時間情報を含む飛跡ファイルを読み込んだ場合に使用可能です。例えば、
cgsamp-free-time.pic を読み込んで下さい。
粒子を動画表示したい場合は、メインメニューの[[動画表示]
動画表示]を選択します。

図 7.5-1 [動画表示]メニュー
粒子の動画表示を中断したい場合は、描画領域の左上のコメント領域を左ダブルクリックします。

図 7.5-2 動画表示画面

23

7.6 体系表示に関する設定を変更する
本ソフトでは表示する体系のレベルが以下の三つから選択できます。
・形状 ・・ 形状毎に表示する。
・領域 ・・ 領域毎に表示する。(形状を組み合わせたもの)
・物質 ・・ 物質毎に表示する。(同物質を使用する領域をまとめたもの)
また、レベル毎に表示の ON/OFF や色の指定ができます。
変更を行うにはメインメニューの[[設定]
設定]-[表示設定]
表示設定]を選択します。選択すると図 7.6-2のような画面が表示されますので、上部にあるタブの中から[[体系表示設定]
体系表示設定]を選択します。

図 7.6-1 [表示設定]メニュー

図 7.6-2 体系表示設定画面

24

7.6.1 表示する体系レベルを変更する
図 7.6-2 の[[表示する
表示する体系情報する体系情報]体系情報]枠内の項目内から表示したい体系レベルを選択します。
7.6.2 表示する形状番号、領域名、物質番号を選択する
各表示項目の下にあるリスト内に現在読み込まれたデータの一覧が表示されます。各
形状番号、領域名、物質番号の前にあるチェックの有無で表示/非表示を切り換えられます。ただし、物質番号が 0 のものについては常に非表示になります。
チェック有り － 表示
チェック無し － 非表示
また、全表示、全非表示はリスト内でマウスの右クリックをすると、図 7.6-3 のようなポップアップメニューが表示されますので、全て表示にするならメニュー内の[[全て表示][全て非表示]表示]を、全て非表示にするなら[非表示]を選択してください。

図 7.6-3 ポップアップメニュー
7.6.3 表示色を変更する
各形状番号、領域名、物質番号毎の色を変更します。まず、リスト内の項目をマウスの左クリックで選択します。選択すると図 7.6-3 のように項目の色が反転します。続
いてマウスの右ボタンをクリックしてポップアップメニューを表示します。CTRL キー
を押したまま選択すると複数の項目が選択できます。
注）このとき先のとき先に項目を
項目を左クリックで
クリックで選択しておかないと選択しておかないと正常しておかないと正常に正常に動作しないので動作しないので注意
しないので注意。
注意。
ポップアップメニューの[[色変更]
色変更]を選択すると、色選択用ダイアログボックスが表示されますので変更する色を選択してください。正常に選択されるとリスト内の色識別
コードが選択したものに変更されます。
7.6.4 線幅を変更する
線表示した場合に体系を描画する線の太さを決定します。線幅は 1 ピクセル×線幅 となります。(デフォルトは 1.0)25

7.6.5 領域作成時の補正オプション
領域を作成する際、下図のように重なった面が削除されない場合がある。そこで表示用に削除する側の形状を少し大きく補正し正しい領域を表示できるようにする。(計算用のインプットファイルでは不要)。

本来ならば削除されるのだが、削除する面と削除される面が重なっている
ため正しく削除されない。

図 7.6-4 補正前画面
設定する項目は図 7.6-2 にある以下の二つである。
・形状補正率 ・・ 寸法に対しての割合で補正量を決める(デフォルト 1E-6)・形状補正値 ・・ 直接補正量を指定(デフォルト 0.005cm)
以下の式で削除する側の寸法を設定する。
補正後の寸法=補正前寸法+補正前寸法*形状補正率+形状補正値
補正後の領域は下図のように表示される補正値 0.005 を指定した結果、正しく表示される。

図 7.6-5 補正後画面

26

7.7 粒子の表示に関する設定を変更する
ヒストリー範囲の選択や、粒子毎の線色、線種、エネルギー範囲の指定を行います。
変更を行うにはメインメニューの[[設定]
設定]-[表示設定]
表示設定]を選択します。選択すると図 7.7-1のような画面が表示されますので、上部にあるタブの中から[[粒子表示設定]
粒子表示設定]を選択しま
す。

①

②
⑨

④
⑤

⑧
⑥

③

⑦

図 7.7-1 粒子表示設定画面
7.7.1 ヒストリー範囲の指定
ヒストリー範囲を指定することで範囲内の粒子のみ表示することができます。指定方
法には以下の二通りがあり,どちらか片方のみ選択できます。
・ヒストリーの範囲で指定する
図 7.7-1 内の[[ヒストリー範囲指定[指定]
ヒストリー範囲指定]
範囲指定]枠内にある[指定]をチェックする→ ヒストリーの範囲で指定する

チェックしない → ヒストリーの範囲で指定しないで設定できます。範囲指定しない場合は全ての粒子が表示されます。
範囲指定するには[指定]チェックボックスの横にあるテキストボックスに範囲の下限値(①のボックス)と上限値(②のボックス)を入力します。
ボックスが空白の場合は下限値、上限値ともにデータから最小ヒストリー数、最大ヒストリー数が割り当てら
れます。

27

・最終ヒストリーからの表示個数で指定する
図 7.7-1 内の[[ヒストリー範囲指定[最終ヒストリー
ヒストリー範囲指定]
範囲指定]枠内にある[最終ヒストリーからのヒストリーからの表示個数からの表示個数]
表示個数]をチェックする→ 最終ヒストリーからの表示個数で指定する

チェックしない → 最終ヒストリーからの表示個数で指定しないで設定できます。範囲指定しない場合は全ての粒子が表示されます。表示する個数は⑨のボックスに入力します。（デフォルトは 1）
7.7.2 表示する粒子の指定
本ソフトで表示する粒子の種類は以下に示す三つです。
・光子(γ、x)・電子(ｅ-)
・陽電子(ｅ+)③のリストボックス内に表示されている各項目のチェックをすると、チェックされた粒子を表示します。非表示にする場合は、非表示にする粒子のチェックをはずします。
7.7.3 粒子の線色を変更する
各粒子の線色を変更します。光子の線色を変更するには④のボタンをクリックします。
クリックすると色選択用ダイアログボックスが表示されますので変更する色を選択してください。正常に選択されるとボックス内の色識別コードが選択したものに変更されます。電子、陽電子についても各枠内のボタンで同様に行えます。
7.7.4 粒子の線種を変更する
各粒子の線種を変更できます。変更できる線種の種類は以下の 4 通りです。
・実線
・破線
・一点鎖線
・二点鎖線
⑤のプルダウンメニューから変更する線種を選択します。
7.7.5 粒子の線幅を指定する
基本になる粒子の線幅を設定する。デフォルトは 1。通常の線は 1 ピクセルで表示しているので線幅を変更すると 1×線幅ピクセルで表示される。⑧のテキストボックス内に正の整数で入力。

28

7.7.6 エネルギー範囲を指定する
図 7.7-1 内の各粒子の枠内にある[[エネルギー範囲指定
エネルギー範囲指定]
範囲指定]をチェックすると、運動エネ
ルギー範囲を指定できます。指定すると運動エネルギー範囲内の粒子のみ表示されます。チェックをはずすと運動エネルギー範囲を指定しません。
範囲指定するには⑥のテキストボックスに範囲の下限値を⑦のテキストボックスに上限値を入力します。ボックスが空白の場合は下限値に 0、上限値に 1e15 が割り当てられます。単位は[MeV]。
7.7.7 WEIGHT について
粒子の WEIGHT(入力データより読込)によって表示する線の太さが変わる。表示する線の太さは 7.7.5 節で指定した線幅×WEIGHT になる。
7.7.8 LATCH 機能
図 7.7-1 内の[[ヒストリー範囲指定[LATCH 機能をヒストリー範囲指定]
範囲指定]枠内にある[LATCH機能を使用する
使用する]
する]をチェックすると LATCH 機能が有効になり、入力データで LACTH を-1 に設定したヒストリーの粒子を表示しないようにします。チェックをはずすと LATCH 機能は無効になります。
7.7.9 時間処理指定
図 7.7-1 内の[[時間処理指定[指定]
時間処理指定]指定]枠内にある[
指定]をチェックすると、TIME オプションが設定されているファイルであれば、動画表示が有効になります。動画の開始時間と終了時間は[[時間範囲][刻み幅]の絶対値で 1 画面に表示する粒子
時間範囲]で指定します。また、[飛跡の時間幅を設定します。ただし、[[刻み幅]が負の場合には開始時間は[[時間範囲]
時間範囲]の開始生存時間に固定され、時間と共に飛跡が伸びていくように表示されます。単位は、TIME オプションが設定されているファイルの単位と同じとなります。

29

7.8 画面表示に関する設定を変更する
背景色やタイトル、座標軸に関する設定を行います。変更を行うにはメインメニューの[[設定]
設定]-[表示設定]
表示設定]を選択します。選択すると図 7.8-1 のような画面が表示されますので、上部にあるタブの中から[[画面表示設定]
画面表示設定]を選択します。

①②

図 7.8-1 画面表示設定画面
7.8.1 背景色を変更する
画面の背景色を変更します。①のボタンをクリックします。クリックすると色選択用
ダイアログボックスが表示されますので変更する色を選択してください。正常に選択されるとボックス内の色識別コードが選択したものに変更されます。
7.8.2 文字色を変更する
座標軸、タイトル、メモ、凡例の色を変更します。②のボタンをクリックします。クリックすると色選択用ダイアログボックスが表示されますので変更する色を選択してください。正常に選択されるとボックス内の色識別コードが選択したものに変更されます。

30

7.8.3 座標軸の設定を行う
以下の項目について変更できます。[[座標軸設定]座標軸設定]枠内の各変更項目の後ろにあるテキストボックスに値を入力します。
・目盛り間隔・・ 目盛線の間隔。

・目盛りサイズ・・ 軸に表示される目盛線のサイズ。

・矢印サイズ・・ 軸に表示される矢印の大きさ。

・目盛り単位

・・ 表示する座標軸の目盛り単位。

・文字サイズ・・ 目盛りラベル文字の大きさ。

・文字間隔

・・ 目盛りラベル文字を表示する間隔。

7.8.4 タイトルを表示する
画面上にタイトルを表示します。設定画面内の[[タイトルを
タイトルを表示する
表示する]
する]をチェックします。チェックをすると画面上に下図のようなタイトルが表示されます。タイトルの内容は後ろに続くエディター内に入力します。また、画面内で左クリックしたまま移動
させることができます。

図 7.8-2 タイトル画面
7.8.5 メモを表示する
画面上にメモを表示します。設定画面内の[[メモをメモを表示する
表示する]
する]をチェックします。チェックをすると画面上に下図のようなタイトルが表示されます。タイトルの内容は後ろに続くエディター内に入力します。また、画面内で左クリックしたまま移動させることができます。

図 7.8-3 メモ画面

31

7.8.6 粒子の凡例を表示する
画面上に下図のような凡例を表示します。表示される内容はヒストリー範囲の番号、各粒子の色と線種です。また、枠内の上下ボタンで表示するヒストリーを変更できます。このボタンは範囲指定で粒子を表示しているときのみ有効になります。移動量は指定範囲を 1 ブロックとして次の粒子を表示します。
(例)範囲が 11-11(1 粒子のみ表示)の場合は 12-12、13-13...と移動します。 範囲が11-15(5 粒子を 1 ブロックとする)の場合は 16-20、21-25...と移動します。

図 7.8-4 凡例画面
7.8.7 ヘルプ表示用ブラウザーを設定する
オンラインヘルプを表示するための html ブラウザー名を設定する。(Linux 版のみ)
ヘルプ表示用ブラウザーの[..]ボタンを押下するとファイル選択ダイアログが表示されるので表示するための html ブラウザーを選択する。デフォルトは mozilla を使用する。

32

7.9 部分削除
表示体系の断面図や内部を見たい場合に使用し、体系の一部を削除して表示出来ます。
設定を行うにはメインメニューの[[設定]
設定]-[表示設定]
表示設定]を選択します。選択すると図 7.9-1のような画面が表示されますので、上部にあるタブの中から[[部分削除]
部分削除]を選択します。

①
図 7.9-1 部分削除設定画面[ジオメトリーのジオメトリーの部分削除を
部分削除を行う]にチェックすると、体系の一部を削除して表示します。削除する部分の設定は[[削除範囲][削除範囲]
削除範囲]枠内にて行います。[
削除範囲]枠内にある項目に数値を入力すると仮想のボックスを作成し、その範囲に入る形状部分を非表示にします。
入力する項目の意味は以下の通りです。単位は[cm]・X 最小値・・ X 方向の最小値(空白の場合はデータから最小値を決定します)・X 最大値・・ X 方向の最大値(空白の場合はデータから最大値を決定します)・Y 最小値・・ Y 方向の最小値(空白の場合はデータから最小値を決定します)・Y 最大値・・ Y 方向の最大値(空白の場合はデータから最大値を決定します)・Z 最小値・・ Z 方向の最小値(空白の場合はデータから最小値を決定します)・Z 最大値・・ Z 方向の最大値(空白の場合はデータから最大値を決定します)

7.10 設定値を初期値に戻す
図 7.9-1①[初期値に戻す]ボタンを押下すると、設定値が初期値に戻ります。

33

第 8 章 マウス操作
マウス操作
体系表示の際にデータが多い場合は、動作が遅れることがあります。その場合は、マウスによる操作を行う前に代表的なデータ以外を非表示にしてから操作し、視点が確定してから残りのデータを表示すると効率的に動作できます。
8.1 表示体系を回転させる
表示した図形を回転回転させるにはマウスの左左ボタンをクリックしたままマウスを移動さ
回転
ボタン
せます。
8.2 表示体系を移動させる
表示した図形を移動移動させるにはマウスの左
左と右ボタンを同時にクリックしたままマウ
移動ボタン
スを移動させます。3 ボタンマウスを使用している場合は、中中ボタンでも移動できます。
ボタン
8.3 表示体系を拡大させる
表示した図形を拡大
拡大させるにはマウスの右右ボタンをクリックしたまま、クリックした
拡大
ボタン
ポイントより下方向にマウスを移動させます。また Shift キーを押しながらマウスの左ボタンをクリックすると下図に示すような枠を表示して拡大範囲を指定できます。

図 8.3 拡大範囲
8.4 表示体系を縮小させる
表示した図形を縮小縮小させるにはマウスの右右ボタンをクリックしたまま、クリックした
縮小
ボタン
ポイントより上方向にマウスを移動させます。

34

第 9 章 体系データ
体系データ作成
データ作成
EGS の体系データを入力しながら作成中のデータを表示できます。また、作成したデー
タを EGS の入力データとして使用することもできます。
9.1 体系データ作成画面を開く
体系データ作成画面を開くにはメインメニューの[[ファイル]
ファイル]-[体系データ
体系データ作成
データ作成]
作成] を選択します。選択すると、図 9.1-2 のような画面が表示されます。

図 9.1-1 [体系データ作成]メニュー

図 9.1-2 体系データ作成画面

35

9.2 体系データを入力する
体系データは CG(Combinatory Geometry)方式を用いて作成します。入力が完了したら[OK]ボタンをクリックすると、入力したデータを基に図形が表示されます。
[OK]
1) CG 形状を入力する
形状データは CG 形状を示す記号、CG 形状の番号及び CG 形状を定義するパラメータで構成され、形状番号は 1～99999 の間で指定します。EGS で使用できる CG 形状のパラメータは 9.6 節 CG 形状一覧を参照してください。
2) CG 領域を入力する[領域番号]
領域番号は Z で始まる 5 文字の識別記号です。番号は自動で割り振られます。
[物質番号]EGS で定義する物質番号を指定します。
[その他の項目]定義した CG 形状を組み合わせて一つの領域を作成します。CG 形状の組み合わせは+,-,OR の 3 種類の記号と形状番号で成り立ちます。+1 は形状番号 1 の内側、-1 は形状
番号の外側を意味します。+1 -2 の様に形状番号が続く場合は+1 と-2 の両方の領域に入ることを意味します。+1 OR +2 のように形状番号の間に OR が入る場合は+1 と+2 のどちらかの領域に入ることを意味します。

図 9.2 CG 領域

36

9.3 作成したデータを保存する[ファイルに
ファイルに保存]
保存]ボタンをクリックすると保存するためのファイルダイアログが表示されるのでファイル名を入力してください。デフォルトのファイルの拡張子は GeometryFile(*.geo)になります。
拡張子が(*.csv)以外の場合はフリーフォーマット形式で保存します。
拡張子が(*.csv)の場合は CSV 形式で保存します。
デフォルトの egs5run では*.data を読み込むように処理されるため、拡張子(geo,csv)
ファイル名を.data に変更する必要があります。また、egs5run を変更して.geo もしく
は.csv を読み込むように修正すれば拡張子を変更する必要がなくなります。
9.4 作成したデータを開く[ファイルを
ファイルを開く]ボタンをクリックするとファイルオープンダイアログが表示されるのでファイルを選択してください。デフォルトのファイルの拡張子は GeometryFile(*.geo)になります。
拡張子が(*.csv)以外の場合はフリーフォーマット形式で読み込みます。
拡張子が(*.csv)の場合は CSV 形式で読み込みます。
9.5 作成したデータの表示
CG Data を作成後、[[表示]表示]ボタンをクリックすると、画面上にジオメトリーが表示され、ジオメトリーのチェックを行うことができます。

37

9.6

Excel でファイルを編集する
拡張子が(*.csv)の形式を使用すると Excel で同様の作業を行えます。Excel で入力する場合の例を以下に示します。Excel で作成後*.csv でファイルを保存することで CGVIEW
上に同一データを再現できます。領域番号は CGVIEW で読み込むときに自動で割り振られるので任意の文字でかまいません。

図 9.6 Excel 入力例

38

9.7

CG 形状一覧

EGS で使用できる CG 体系は以下に示す 14 種類です。
1)四角形(RPP)XYZ について指定した範囲を有する四角形を作成します。

図 9.7-1 RPP[ 形式 ] Xmin(1) Xmax(2) Ymin(3) Ymax(4) Zmin(5) Zmax(6)(注) ( )内の数字は体系データ入力画面のパラメータ入力番号を示します。(1)ならば
表中 1 番の行に入力します。(以降同じ)
Xmin,Xmax : X 方向の最小値と最大値
Ymin,Ymax : Y 方向の最小値と最大値
Zmin,Zmax : Z 方向の最小値と最大値2)円柱(RCC)
基準点 P を中心に半径 R で作成される円をベクトル H 方向に沿って移動した跡にできる円柱を作成します。

図 9.7-2 RCC[ 形式 ] Px(1) Py(2) Pz(3) Hx(4) Hy(5) Hz(6) R(7)Px,Py,Pz : 基準点となる円の中心 X,Y,Z 座標
Hx,Hy,Hz : 円柱の高さ方向の X,Y,Z ベクトル

39

(マイナスも入力可。円柱の高さは√(Hx×Hx+Hy×Hy+Hz×Hz)となる)
R : 円の半径3)球(SPH)
基準点 P を中心に半径 R で作成される球を作成します。

図 9.7-3 SPH[ 形式 ] Px(1) Py(2) Pz(3) R(4)Px,Py,Pz : 基準点となる球の中心 X,Y,Z 座標
R : 球の半径4)円錐(TRC)
基準点 P を中心に半径 R1 で作成される円をベクトル H 方向に沿って半径 R2 になるよ
うに移動した跡にできる円錐を作成します。

図 9.7-4 TRC[ 形式 ] Px(1) Py(2) Pz(3) Hx(4) Hy(5) Hz(6) R1(7) R2(8)Px,Py,Pz : 基準点となる円の中心 X,Y,Z 座標40

Hx,Hy,Hz : 円錐の高さ方向の X,Y,Z ベクトル（マイナスも入力可。円錐の高さは√(Hx×Hx+Hy×Hy+Hz×Hz)となる)
R1 : 円錐基準面の半径R2 : 円錐ベクトル面の半径5)トーラス(TOR)
基準点 P から半径 R1 の円に沿って半径 R2 の円を移動した跡にできるトーラスを作成します。トーラスは必ず XYZ のいずれかの軸に平行であり、th1,th2 で有効範囲を指定することができます。

図 9.7-5 TOR[ 形式 ] Px(1) Py(2) Pz(3) R1(4) R2(5) th1(6) th2(7) ax(8)Px,Py,Pz : 基準点となるトーラスの中心 X,Y,Z 座標
R1 : トーラスの輪の半径R2 : トーラスの管の半径
th1,th2 : トーラスの開始角度、終了角度(度(Degree)で入力)(0～360 の場合は切れ目のないトーラスとなる)
ax : トーラスの輪の中心軸(1=X 軸 2=Y 軸 3=Z 軸)

416)四角柱(BOX)
基準点 P とベクトル A1 と A2 の辺で形成される四角形をベクトル A3 方向に沿って移動した跡にできる四角柱を作成します。ただし、A1、A2 及び A3 のベクトルは直行していること。

図 9.7-6 BOX[ 形式 ] Px(1) Py(2) Pz(3) A1x(4) A1y(5) A1z(6) A2x(7) A2y(8) A2z(9)A3x(10) A3y(11) A3z(12)Px,Py,Pz : 基準点となる頂点の X,Y,Z 座標
A1x,A1y,A1z : 基準点から延ばした辺 1 の X,Y,Z ベクトル
A2x,A2y,A2z : 基準点から延ばした辺 2 の X,Y,Z ベクトル
A3x,A3y,A3z : 基準点から延ばした辺 3 の X,Y,Z ベクトル7)楕円柱(REC)
基準点 P を中心に短軸 R1 と長軸 R2 で作成される楕円をベクトル H 方向に沿って移動した跡にできる楕円柱を作成します。

図 9.7-7 REC

42

[ 形式 ] Px(1) Py(2) Pz(3) Hx(4) Hy(5) Hz(6) R1x(7) R1y(8) R1z(9)R2x(10) R2y(11) R3z(12)Px,Py,Pz : 基準点となる楕円の中心 X,Y,Z 座標
Hx,Hy,Hz : 楕円柱の高さ方向の X,Y,Z ベクトル (マイナスも入力可。円柱の高さは√(Hx×Hx+Hy×Hy+Hz×Hz)となる)
R1x,R1y,R1z : 楕円の短軸を示す X,Y,Z ベクトル
R2x,R2y,R2z : 楕円の長軸を示す X,Y,Z ベクトル8)楕円錐(TEC)
基準点 P を中心に短軸 R1 と長軸 R2 で作成される楕円をベクトル H 方向に沿って比率 wの楕円に移動した跡にできる楕円錐を作成します。

図 9.7-8 TEC[ 形式 ] Px(1) Py(2) Pz(3) Hx(4) Hy(5) Hz(6) R1x(7) R1y(8) R1z(9)R2x(10) R2y(11) R3z(12) w(13)Px,Py,Pz : 基準点となる楕円の中心 X,Y,Z 座標
Hx,Hy,Hz : 楕円柱の高さ方向の X,Y,Z ベクトル (マイナスも入力可。円柱の高さは√(Hx×Hx+Hy×Hy+Hz×Hz)となる)
R1x,R1y,R1z : 楕円の短軸を示す X,Y,Z ベクトル
R2x,R2y,R2z : 楕円の長軸を示す X,Y,Z ベクトルw : 楕円の比率

439)楕円球(ELL)
焦点座標 P1 と P2 主軸 R の長さが一定の距離を保つ軌跡でできる楕円球を作成します。

図 9.7-9 ELL[ 形式 ] P1x(1) P1y(2) P1z(3) P2x(4) P2y(5) P2z(6) R(7)P1x,P1y,P1z : 1 つ目の焦点座標(X,Y,Z)
P2x,P2y,P2z : 2 つ目の焦点座標(X,Y,Z)
R : 主軸の長さ10)三角柱(WED)
基準点 P とベクトル A1 と A2 の辺で形成される三角形をベクトル A3 方向に沿って移動した跡にできる三角柱を作成します。ただし、A1、A2 及び A3 のベクトルは直行していること。

図 9.7-10 WED

44[ 形式 ] Px(1) Py(2) Pz(3) A1x(4) A1y(5) A1z(6) A2x(7) A2y(8) A2z(9)A3x(10) A3y(11) A3z(12)Px,Py,Pz : 基準点となる頂点の X,Y,Z 座標
A1x,A1y,A1z : 基準点から延ばした三角形辺 1 の X,Y,Z ベクトル
A2x,A2y,A2z : 基準点から延ばした三角形辺 2 の X,Y,Z ベクトル
A3x,A3y,A3z : 基準点から延ばした高さの X,Y,Z ベクトル
11)多面体(ARB)
任意の 4,5 または 6 面体。頂点の XYZ 座標 P1～P8 を指定し、面データ S1～S6 を十進数で表した 4 桁の整数で指定する。4 桁の整数は各桁が面の頂点番号を示し、順序は任意とする。4 点より少ない頂点で構成される面は余りの頂点番号に 0 を入力する。使用しない頂点座標は XYZ に 0 を入力する。

図 9.7-11 ARB[ 形式 ] P1x(1) P1y(2) P1z(3) P2x(4) P2y(5) P2z(6) P3x(7) P3y(8) P3z(9) P4x(10)P4y(11) P4z(12) P5x(13) P5y(14) P5z(15) P6x(16) P6y(17) P6z(18) P7x(19)P7y(20) P7z(21) P8x(22) P8y(23) P8z(24) S1(25) S2(26) S3(27) S4(28)S5(29) S6(30)P1x,P1y,P1z : 頂点 1 の X,Y,Z 座標
P2x,P2y,P2z : 頂点 2 の X,Y,Z 座標
P3x,P3y,P3z : 頂点 3 の X,Y,Z 座標
P4x,P4y,P4z : 頂点 4 の X,Y,Z 座標
P5x,P5y,P5z : 頂点 5 の X,Y,Z 座標
P6x,P6y,P6z : 頂点 6 の X,Y,Z 座標
P7x,P7y,P7z : 頂点 7 の X,Y,Z 座標
P8x,P8y,P8z : 頂点 8 の X,Y,Z 座標
S1,S2,S3,S4,S5,S6 : 6 つの面を記述する 4 桁の整数。4 桁の整数は各桁が面の頂点番号を示し、順序は任意とする。4 点より少ない頂点で構成される面は余りの頂点番号に0 を入力する。

4512)平面(HAF)
指定したベクトル N 方向を領域とする半空間を作成する。

図 9.7-12 HAF[ 形式 ] Nx(1) Ny(2) Nz(3) L(4)Nx,Ny,Nz : 内側領域を示すベクトル X,Y,Z (マイナスも入力可)
L : ベクトル N に沿った原点(0,0,0)から平面までの距離 (マイナスも入力可)

4613)六角形(HEX)
基準点 P を中心に幅 W で側面がベクトル A を向く六角形をベクトル H 方向に沿って移動した跡にできる六角形を作成します。

図 9.7-13 HEX[ 形式 ] Px(1) Py(2) Pz(3) Hx(4) Hy(5) Hz(6) W(7) Ax(8) Ay(9) Az(10)Px,Py,Pz : 基準点となる円の中心 X,Y,Z 座標
Hx,Hy,Hz : 円柱の高さ方向の X,Y,Z ベクトル (マイナスも入力可。円柱の高さは√(Hx×Hx+Hy×Hy+Hz×Hz)となる)
W : 六角形の幅
Ax,Ay,Az : 側面の一つに垂直なベクトル A (長さは任意)

4714)一般楕円体(GEL)
基準点 P を中心に A1、A2、A3 ベクトルを軸に持つ楕円を回転させてできる楕円体を作成します。ただし、A1、A2 及び A3 のベクトルは直行していること。

図 9.7-14 GEL[ 形式 ] Px(1) Py(2) Pz(3) A1x(4) A1y(5) A1z(6)A2x(7) A2y(8) A2z(9) A3x(10) A3y(11) A3z(12)Px,Py,Pz : 基準点となる楕円体の中心 X,Y,Z 座標
A1x,A1y,A1z : 楕円の軸方向を示す X,Y,Z ベクトル
A2x,A2y,A2z : 楕円の軸方向を示す X,Y,Z ベクトル
A3x,A3y,A3z : 楕円の軸方向を示す X,Y,Z ベクトル48

第 10 章 整合性確認
9 章で作成した体系データの整合性を確認するにはメインメニューの[[設定]
設定]-[体系整
合性確認]
合性確認]を選択します。選択すると、図 10-2 のような画面が表示されます。

図 10-1[体系整合成確認]メニュー

図 10-2 整合性確認画面
検索を行うには以下の値を入力して[検索開始]ボタンを押下する。
開始位置 : 粒子を発生させる位置を指定する。(デフォルトは中心)幅を持たせた場合は範囲内からランダムに決定される。
検索角度 : 粒子を飛ばす方向を角度(°)で指定する。(デフォルトは全方向)幅を持たせた場合は範囲内からランダムに決定される。

49

theta : 0 ～ 360 (X 軸から Y 軸方向への角度)
phi: -90 ～ 90 (X 軸から Z 軸方向への角度)

図 10-3 検索方向
検索回数 : 飛ばす粒子の数を指定する。(デフォルトは 1000)[検索開始]ボタンは実行中[検索中止]となり、クリックすると検索を中止することができます。検索中にエラーが発生すると[結果出力]欄にエラーが発生した XYZ 座標と指定デー
タに従って検索を実行します。実行中は[検索中止]となり、クリックすると検索を中止することができます。検索の結果として、体系表示画面では、エラーが発生した領域内は線で表示されます。また、結果出力にはエラー点を包含する形状のリストが表示されます。
また、エラーの数が 10 を超えると強制的に検索が中止されます。

50

図 10-4 結果出力

図 10-5 画面表示

51

第 11 章 領域体積計算
9 章で作成した体系データの領域の体積を計算するにはメインメニューの[[設定]
設定]-[領域
体積計算]
体積計算]を選択します。選択すると、図 11-2 のような画面が表示されます。各領域の体積を計算することができます。
計算は以下の方法を用いて行います。
２のべき数で等間隔に分割した yz 面、xz 面、xy 面の中心から垂直に粒子を飛ばし、領域を横切る長さと分割した面の面積を掛け算し、それらを集計することで領域の体積を求めます。値の確かさは、yz 面、xz 面、xy 面の３つ値の収束レベル、２のべき数の分割レベル間の収束レベルで判断します。２のべき数での等間隔の分割は、べき数が 0 から指定した数値まで行います。

図 11-1[領域体積計算]メニュー

図 11-2 領域体積計算画面

52

11.1 分割レベルを指定する[分割レベル
分割レベルのレベルの指定]指定]では、計算を行う分割レベルを２のべき数で指定します。最大は 14です。(デフォルトは 12)
11.2 評価位置を指定する[評価位置]
評価位置]では、体積を計算するための面の範囲を指定します。求めたい領域を包含する範囲で指定します。単位は cm です。
11.3 収束レベルを指定する[収束レベル収束レベル]
レベル]では、計算の打ち切りを行う収束レベルを指定します。
３方向の収束：yz 面、xz 面、xy 面の３つの値間の収束値で判断します。
分割レベルの収束：分割レベル（２のべき数）が一つ前で求めた値との収束を yz 面、
xz 面、xy 面の全てについて判断します。
11.4 結果の表示
結果は、領域番号（Z）、計算された体積（yz 面、xz 面、xy 面の順）および分割レベル（２のべき数）が一つ前で計算された値との収束値（yz 面、xz 面、xy 面の順）が表示されます。なお、□をチェックした領域の数値が収束レベルでの計算打ち切りの判断の対象領域になります。
11.5 計算状況の表示
計算の進捗状況が、率（切り上げ表示）、分割レベル（nlp）、分割個数（nlpcnt)と計算点（外のループ i(ilp)、内のループ（jlp)）で表示されます。
11.6 ファイルに保存する[ファイルに
ファイルに保存]
保存]ボタンをクリックすると、計算結果をファイルに保存できます。保存
されるデータは、領域番号（Z）および計算された体積（平均値、yz 面、xz 面、xy 面の順）です。
Zone_No.,,average, x_direction, y_direction, z_direction(cm3),(cm3),(cm3),(cm3)

Z1, 1.0395E+002, 1.0395E+002, 1.0395E+002, 1.0396E+002
Z2, 1.3347E+001, 1.3356E+001, 1.3356E+001, 1.3329E+001
Z3, 1.2529E+003, 1.2529E+003, 1.2529E+003, 1.2530E+003
Z4, 6.9723E+001, 6.9787E+001, 6.9787E+001, 6.9595E+001
Z5, 7.7361E+003, 7.7363E+003, 7.7363E+003, 7.7357E+003

53

Z6, 2.4868E+002, 2.4844E+002, 2.4844E+002, 2.4915E+002
Z7, 8.5841E+002, 1.2876E+003, 1.2876E+003, 0.0000E+000
11.7 計算の開始[計算開始]
計算開始]ボタンをクリックすると、指定データに従って計算を実行します。実行中は[計算中止]となり、クリックすると計算を中止することができます。
11.8 閉じる[閉じる]じる]ボタンをクリックすると、画面を閉じてメイン画面に戻ります。
11.9 ヘルプの表示[ヘルプ]ヘルプ]ボタンをクリックすると、ヘルプを表示します。

54

第 12 章 ファイル形式
ファイル形式CGVIEW で読み込むための軌跡データのファイル形式及び体系データ作成画面で使用するファイル形式について述べる。
12.1 軌跡データの構成
軌跡ファイルは以下のデータ順で構成されます。
a)体系・領域データ[CG 形式の場合](1)CSTACG データの入力開始とファイルの形式を意味する記号を入力。CSTA：固定フォー
マット CSTA-FREE：フリーフォーマット CSTA-CSV：CSV フォーマット
時間設定の場合[CSTA-TIME：固定フォーマット CSTA-FREE-TIME：フリーフォーマ
ット CSTA-CSV-TIME：CSV フォーマット](2)CSYM
CG 形状記号(RPP,RCC,SPH 等)を入力 (2)～(4)は 1 セットのデータとし、複数デー
タの記述ができる。
(3)CNUM
CG 形状番号(5 桁の正の整数)を入力 (2)～(4)は 1 セットのデータとし、複数デー
タの記述ができる。
(4)CPRM[ ]形状パラメータ(入力数は形状記号に異なる) (2)～(4)は 1 セットのデータとし、
複数データの記述ができる。
(5)END : 形状データの終わりを意味する記号 END を入力(6)ZNUM
領域番号(Z で始まる 5 文字の記号 Z0001 等) (6)～(7)は 1 セットのデータとし、
複数データの記述ができる。
(7)ZPRM[ ]領域パラメータ(OR と形状番号の組み合わせ) (6)～(7)は 1 セットのデータとし、
複数データの記述ができる。
(8)END : 領域データの終わりを意味する記号 END を入力(9)CEND : CG データの入力終了を意味する記号 CEND を入力[同心球の場合](1)GSTACG 以外の入力開始とファイルの形式を意味する記号を入力。GSTA：固定フォーマット GSTA-FREE：フリーフォーマット GSTA-CSV：CSV フォーマット

55

時間設定の場合[GSTA-TIME：固定フォーマット GSTA-FREE-TIME：フリーフォーマ
ット GSTA-CSV-TIME：CSV フォーマット](2)SPHR : 同心球データの入力開始を意味する記号 SPHR を入力(3)GCNT : 同心球の数(1 以上の整数)(4)GPRM[ ] : 同心球の半径を小さいほうから GCNT 分入力(5)GEND : CG 以外の入力終了を意味する記号 GEND を入力[同心円筒の場合](1)GSTACG 以外の入力開始とファイルの形式を意味する記号を入力。GSTA：固定フォーマット GSTA-FREE：フリーフォーマット GSTA-CSV：CSV フォーマット
時間設定の場合[GSTA-TIME：固定フォーマット GSTA-FREE-TIME：フリーフォーマ
ット GSTA-CSV-TIME：CSV フォーマット](2)CYLS : 同心円筒データの入力開始を意味する記号 CYLS を入力(3)GCNT1 : 同心円筒の数(1 以上の整数)(4)GCNT2 : Z 軸方向の断面数(1 以上の整数)(5)GPRM1[ ] : 同心円筒の半径を小さいほうから GCNT1 分入力(6)GPRM2[ ] : Z 軸方向の断面位置(Z 寸法)を小さいほうから GCNT2 分入力(7)GEND : CG 以外の入力終了を意味する記号 GEND を入力[平板の場合](1)GSTACG 以外の入力開始とファイルの形式を意味する記号を入力。GSTA：固定フォーマット GSTA-FREE：フリーフォーマット GSTA-CSV：CSV フォーマット
時間設定の場合[GSTA-TIME：固定フォーマット GSTA-FREE-TIME：フリーフォーマ
ット GSTA-CSV-TIME：CSV フォーマット](2)SLAB : 平板データの入力開始を意味する記号 SLAB を入力(3)GCNT1 : X 軸方向の断面数(1 以上の整数)(4)GCNT2 : Y 軸方向の断面数(1 以上の整数)(5)GCNT3 : Z 軸方向の断面数(1 以上の整数)(6)GPRM1[ ] : X 軸方向の断面位置(X 寸法)を小さいほうから GCNT1 分入力(7)GPRM2[ ] : Y 軸方向の断面位置(Y 寸法)を小さいほうから GCNT2 分入力(8)GPRM3[ ] : Z 軸方向の断面位置(Z 寸法)を小さいほうから GCNT3 分入力(9)GEND : CG 以外の入力終了を意味する記号 GEND を入力
b)物質データ(1)MSTA : 物質データの入力開始を意味する記号 MSTA を入力(2)MCNT : 入力する物質データの数(領域数の等しくなければならない)(3)MMAT[ ] : 領域に対応する物質番号(0 以上の整数)

56

(4)MEND : 物質データの入力終了を意味する記号 MEND を入力
c)粒子の軌跡データ(1)0 : 1 粒子のヒストリーデータ開始を意味する数値 0 を入力(2)PHIS : ヒストリー番号を入力(1 以上の整数)(3)PTYP
粒子の種類を入力 ( 4:光子 5:電子 6:陽電子) (3)～(11)は 1 セットのデータとし、
複数データの記述ができる。
(4)PX粒子の反応点の X 座標 (3)～(11)は 1 セットのデータとし、複数データの記述ができる。
(5)PY
粒子の反応点の Y 座標 (3)～(11)は 1 セットのデータとし、複数データの記述ができる。
(6)PZ
粒子の反応点の Z 座標 (3)～(11)は 1 セットのデータとし、複数データの記述ができる。
(7)PENR
粒子の反応エネルギー (3)～(11)は 1 セットのデータとし、複数データの記述ができる。
(8)PZON
粒子の反応領域の領域番号(1 以上の整数) (3)～(11)は 1 セットのデータとし、複数データの記述ができる。
(9)PWGT
粒子のウエイト(通常は 1、数値に合わせて粒子の描画線の太さが変わる) (3)～(11)は 1 セットのデータとし、複数データの記述ができる。
(10)PTIMEPTIME :粒子の生存時間（TIME オプション設定の場合のみ記述する) (3)～(11)は 1
セットのデータとし、複数データの記述ができる。
(11)PSTP
次の粒子データが起点となる場合は-1 を入力する。そのまま軌跡が続く場合は入力しない。(3)～(11)は 1 セットのデータとし、複数データの記述ができる。
(12)9 : 1 粒子のヒストリーデータ終了を意味する数値 9 を入力(13)PLAT
直前に入力したヒストリーの LATCH。(通常は 1、このヒストリーの粒子を非表示にする場合は-1)57

12.2 軌跡データ(フリーフォーマット)の形式
フリーフォーマットの場合は軌跡データの区切りを空白としてファイルを作成します。
空白の数は複数でも 1 つとして判断されます。
12.3 軌跡データ(CSV フォーマット)の形式
CSV フォーマットの場合は軌跡データの区切りをカンマ( , )としてファイルを作成します。
12.4 軌跡データ(固定フォーマット)の形式
固定フォーマットは旧バージョンで使用された方法で軌跡データは決められたフォーマットで入力する必要があります。
a)体系・領域データ[CG 形式の場合](1)Format('CSTA')(2)Format(1h, 2x, a3, i5, 1p8e12.4) CSYM, CNUM, CPRM[ ]Format(1h, 2x, 3x, 5x, 1p8e12.4) CSYM, CNUM, CPRM[ ] (CPRM が 8 より多い場合)(3)Format(1h, 'END')(4)Format(1h ,2x,a3,5x,1x,10(a2,i5,1x)) ZNUM, ZPRM[ ](a2 は OR か空白、i5 は正負の形状番号)
Format(1h ,2x,3x,5x,1x,10(a2,i5,1x)) ZPRM[ ] (ZPRM が 10 より多い場合)(5)Format('END')(6)Format('CEND')[同心球の場合](1)Format('GSTA')(2)Format('SPHR')(3)Format(3i6) GCNT(4)Format(4E15.7) GPRM[ ](5)Format('GEND')[同心円筒の場合](1)Format('GSTA')(2)Format('CYLS')(3)Format(3i6) GCNT1, GCNT2(4)Format(4E15.7) GPRM1[ ](5)Format(4E15.7) GPRM2[ ](6)Format('GEND')[平板の場合]58

(1)Format('GSTA')(2)Format('SLAB')(3)Format(3i6) GCNT1, GCNT2, GCNT3(4)Format(4E15.7) GPRM1[ ](5)Format(4E15.7) GPRM2[ ](6)Format(4E15.7) GPRM3[ ](7)Format('GEND')b)物質データ(1)Format('MSTA')(2)Format(i4) MCNT(3)Format(15i4) MMAT[ ](4)Format('MEND')c)粒子の軌跡データ(1)Format('0',i5) PHIS(2)Format(i1,1p3E13.6,1pE10.3,i4,1p2E10.3)
PTYP, PX, PY, PZ, PENR, PZONO, PWGT, PTIME(3)Format('-1')(4)Format('9',i5) PLAT
12.5 体系ファイルの構成
体系作成画面で入出力されるデータファイルは以下の順で構成されます。
(1)CSYM
CG 形状記号(RPP,RCC,SPH 等)を入力 (1)～(3)は 1 セットのデータとし、複数データの記述ができる。
(2)CNUM
CG 形状番号(5 桁の正の整数)を入力 (1)～(3)は 1 セットのデータとし、複数データの記述ができる。
(3)CPRM[ ]形状パラメータ(入力数は形状記号に異なる) (1)～(3)は 1 セットのデータとし、複数データの記述ができる。
(4)END : 形状データの終わりを意味する記号 END を入力(5)ZNUM
領域番号(Z で始まる 5 文字の記号 Z0001 等) (5)～(6)は 1 セットのデータとし、複数データの記述ができる。
(6)ZPRM[ ]59領域パラメータ(OR と形状番号の組み合わせ) (5)～(6)は 1 セットのデータとし、複数データの記述ができる。
(7)END : 領域データの終わりを意味する記号 END を入力(8)MMAT[ ] : 領域に対応する物質番号(0 以上の整数、領域数分入力)
12.6 体系ファイル(フリーフォーマット)の形式
フリーフォーマットの場合は体系ファイルの区切りを空白としてファイルを作成します。空白の数は複数でも 1 つとして判断されます。拡張子を.data に変更することで EGSのインプットファイルにすることができます。
12.7 体系ファイル(CSV フォーマット)の形式
CSV フォーマットの場合は体系ファイルの区切りをカンマ( , )としてファイルを作成します。拡張子を.data に変更することで EGS のインプットファイルにすることができます。ただし、フリーフォーマットと異なり 1 形状のパラメータに改行を使用しないため 1
行のデータが 250 カラムを超えるとデフォルトの EGS ではエラーとなります。
フリーフォーマット例
RCC

6

0

0

7.9

1.9

6.1

0

0 <-改行有り

CSV 例
RCC,6,0,0,6.1,0,0,7.9,1.9<-改行無し60

第 13 章 EGS ユーザーコード
軌跡の表示のためには、粒子の種類、位置、エネルギーをデータとして出力しておく必要がある。同じ種類の一個の粒子のみを計算している場合にはその座標を単につないでいくだけでよい。EGS の計算では反応により通常 1 個または 2 個の粒子が生成し、そのうちエネルギーの小さい方の粒子を先に追跡する。このため、新しい粒子が発生した場
合や粒子の種類が変わった場合に線で結ぶ反応点のセットを適切に出力する必要がある。
このために SUBROUTINE PLOTXYZ を用いる。SUBROUTINE GEOMGT を用いて CG 情報読み込
む際にエコーバックが出力される。これを体系・軌跡データファイルに出力するによって CG 情報を利用する。各領域に対して指定される物質の番号も体系・軌跡データファイ
ルに出力する。円筒・平板形状の情報はこれまでに作成してある SUBROUTINE GEOMOUT を
そのまま利用して出力する。以下、軌跡ファイル(フリーフォーマット)出力用ユーザー
コードの作り方を紹介する。
13.1 メインプログラム
粒子情報の書式と出力される体系情報は表 13.1 に示すように NPRECI の値によって設定される。円筒平板では、XMIN-ZMAX を用いて各座標の最小値と最大値を指定する。CG では XMIN-ZMAX は使用しない。FNORM は常に使用しない。
表 13.1 NPRECI の値とジオメトリー、プログラムの対応
NPRECI

円筒CG

Medium平板

粒子情報

対応プログラム

書式

0

○

×

×

I5

EGS4WIN32

1

○

×

×

I8

EGS4WIN32

2

○

○

○

1PE13.6

CGVIEW(固定フォーマット)

3

○

○

○

1PE14.6

CGVIEW(フリーフォーマット)

軌跡ファイル(フリーフォーマット)用には NPRECI=3;を指定する。
NPRECI=3;
体系・軌跡データのファイル名を指定して書き込み状態にする。
open(UNIT=39,FILE='egs5job.pic',STATUS='unknown')
ifto = 3961

CG 情報を得るために CALL GEOMGT の前後に次の行を追加する。
write(ifto,fmt="('CSTA-FREE')")
write(ifto,110)
110

FORMAT('CEND')

各リージョンに割り当てられる物質の番号を出力する。
write(ifto,310)
310

FORMAT('MSTA')
write(ifto,320) nreg

320

FORMAT(I4)write(ifto,330) (med(i),i=1,nreg)

330

FORMAT(15I4)
write(ifto,340)

340

FORMAT('MEND')

ヒストリーの始まりを示す記号としての 0 の後に、ヒストリー番号を出力する。
if(npreci.eq.3) thenif (i.le.imax) write(ifto,197) i
endif
197

format('0',I12)

Call Shower の後に全ての粒子が終了したことを示すサブルーチン PLOTXYZ を Call する。
if (i.le.imax) call plotxyz(99,0,0,0.D0,0.D0,0.D0,0.D0,0,0.D0)
ヒストリーの終わりを示す記号 9 と LACTH を出力する。
if(npreci.eq.3) thenif (i.le.imax) write(ifto,207)
endif
207

FORMAT('91')

シャワー表示のためのデータを出力するファイルを閉じる。
close(UNIT=ifto)6213.2 SUBROUTINE AUSGABAUSGAB では、粒子のスタック番号(NP)、タイプ(IQ)、位置(X, Y, Z)、エネルギー(E)、
現在いる領域(IR)、ウェイト(WT)を記録するために、PLOTXYZ を Call する必要がある。
if (i.le.imax) thencall plotxyz(iarg,np,iq(np),x(np),y(np),z(np),e(np),ir(np),wt(np))
end if13.3 SUBROUTINE PLOTXYZ
最初の SUBROUTINE PLOTXYZ は SLAC で Unified Graphics System を使用した 2 次元軌
跡表示のために開発された。次にこれを改良し DOS および windows 用の 3 次元軌跡表示
システム EGS4WIN32 に対応する PLOTXYZ が開発されている。このときにはファイルサイ
ズをなるべく小さくするために粒子の情報を整数で出力していた。CGVIEW では実数で出力するように変更されている。これらの出力形式は表 13.1 に示したように制御変数
NPRECI で制御される。

63

第 14 章 CG を用いた EGS 用計算体系作成例
EGS で CG を使用するための入力ファイル(*.data)を作成する例を図 14-1 のようなγ線検出器を用いて説明します。作成時の条件は以下に従うものとします。(このデータは sample
フォルダ内の g-detector.geo と同じです)
・ A 地点を基準点(0,0,0)とし、上方向を Z プラスとする。
・ 作成する範囲は A 地点より上とする。
・ 物質番号は AIR が 1、ACRILIC が 2 とする。

図 14-1 γ線検出器例
以下に定義する領域と必要になる形状データを示します。また、領域の範囲を図 14-2 に示します。
領域 1 ・・ 検出部(ACRYLIC)
検出部は円柱を用いて表現する。
必要な形状
RCC 1 0.0 0.0 0.0 0.0 0.0 17.0 0.235
領域定義
Z0001

2 +1

64

領域 2 ・・ 検出器内部領域(AIR)検出器の内部領域は円柱と球の組み合わせで表現する。このとき領域 1 の分をマイナスする必要がある。
必要な形状
RCC 2 0.0 0.0 2.81

0.0 0.0

SPH 3 0.0 0.0 2.81

2.67

SPH 4 0.0 0.0 14.79

2.67

11.98 2.67

領域定義
Z0002

1 +2 -1

OR

+3 -1

OR

+4

-1

領域 3 ・・ 検出器外周部領域(ACRYLIC)
内部領域同様は円柱と球の組み合わせで表現する。このとき領域 1,2 の分をマイナスする必要がある。
必要な形状
RCC 5 0.0 0.0 2.81

0.0 0.0

SPH 6 0.0 0.0 2.81

2.81

SPH 7 0.0 0.0 14.79

2.81

11.98 2.81

領域定義
Z0003

2 +5 -2

OR+6 -5 -3 -1

OR +7 -5 -4

領域 4 ・・ 検出器外領域(AIR)検出器より外側の空気層を作成する必要があります。空気層の範囲は基準点から 100.0の球とします。
必要な形状
SPH 8 0.0 0.0 0.0 100.0
領域定義
Z0004

1 +8 -5 -6 -7 -1

領域 5 ・・ 計算終了領域
EGS+CG では必ず計算終了領域が必要になります。計算終了領域は物質番号を 0 とし、この領域に粒子が到達すると計算を終了するようになっています。終了領域は領域 4 から 1
大きい球を想定します。
必要な形状
SPH 9 0.0 0.0 0.0 101.0
領域定義
Z0005

0 +9 -865

図 14-2 領域範囲図

66