egs5の線源形状の変化（等方線源編）

本記事の紹介

本記事は、放射線輸送計算コードegs5の公式より配布されている『ucnaicgv』のシミュレーションにおいて、線源の形状を変化させ、デフォルトのペンシルビームから等方線源（全方向(４π)に同一の確率で放射線を出す線源）にする方法を紹介します。（egs5公式はこちら）

手順

ucnaicgv.fのSTEP8のdo文内のSelect incident angleの下に

        call randomset(rnnow)
        win = 2.0*rnnow-1.0
        call randomset(rnnow)
        phai = PI*(2.0*rnnow-1.0)
        uin = dsqrt(1.0-win*win)*cos(phai)
        vin = dsqrt(1.0-win*win)*sin(phai)

と入力します。call randomset(rnnow)はrnnowに0~1の値をランダムで格納するための呼び出しです。

次にSTEP1に上記で使用したphaiという変数を定義するために、zi0のあとに

,phai

と入力します。

実行結果（CGVIEWの表示：maxpict=200）

等方線源になっていることがわかります。

見やすいようにmaxpictを50から200に変更しています。