EGS4程序介绍与安装.docVIP

EGS4程序介绍与安装 核技术应用 2011082709 李福龙 1、EGS程序介绍 EGS是英文Electron-Gamma Shower 的缩写，它是一个利用蒙特卡罗方法模拟在任意几何中，能量从几个KeV到几个TeV的电子光子射过程的通用程序包。EGS4作为利用蒙特卡罗方法进行模拟计算的一个主要程序已经公布很长一段时间。在这段时间里中，有很多人为EGS4的发展贡献了自己的力量。尤其是在辐射治疗方面的进展突飞猛进，对于精确定位治疗的模拟计算成为当今科技发展中一个热门的领域。伴随计算机速度的发展，也同时由于网络的发展，模拟计算成为解决问题的一种主要方法和手段。 EGS4程序的改进： 对EGS4物理模型的修改：为改善剂量精度，改进了物理模型，提供了一些特殊算法（如PRESTA，EDGSET等）和宏，并考虑了电磁场对正负电子的影响。PRESTA（the Parameter Reduced electron-Step Transport Algorithm）改善了电子路程的修正算法（PLC）、侧面相关算法（LCA）和边界穿越算法（BCA），解决了电子步长和计算精度的矛盾，在保持计算精度的情况下，增大了电子步长，提高了计算速度；EDGSET使EGS4考虑了物质K线和L线的荧光发射。 工具和技术的改进：修正了截面数据库，并且为了提高计算精度，为PEGS4增加APRIM选项，以使电子在单位步长内能量损失（Electron radiative stopping power）符合ICRU-37标准。还开发了如绘图工具，PEGS工具，随机数产生器等。 系统和其它的支持：提供了几何库和一些对计算结果进行数据处理及图形显示的工具，如EGS-Windows（NRCC1991）和EGS4PICT（KEK1994）等。 EGS4程序的特点元素和介质材料齐全：可以模拟在任何元素、化合物、或者混合物中的电子（正电子和负电子）、或者光子辐射传输过程。数据预处理包——PEGS4可以产生1到100号的元素的反应截面表供EGS4使用。 粒子输运方式是随机的：带电粒子和光子的输运采用随机游动方式（random）进行，而非离散（discrete）步长。 能量范围广泛：带电粒子动能范围从几十KeV到几千GeV。尽管更高的能量范围可能有效，但是其物理有效性待检验。对于更低的能量，需要对反应截面数据和步长等进行特殊处理。光子的能量范围从1KeV到几千GeV。 开放式结构：允许用户自己编写程序完成对源粒子特征的定义、几何关系处理、结果的记录和对EGS4系统程序的调用。 涉及的物理过程比较全面，包括： 轫致辐射（不包括对低能的Elwert修正）。 正电子的飞行和静止湮灭。 Molière多重散射（例如：原子核的Coulomb散射），减小的角度以连续分布而不是离散分布方式进行抽样，通过任意步长，随机选定来实现。避免过大或者过小而与理论冲突。 Mller（e-e-）和Bhabha（e+e-）散射。使用了精确的而不是渐进的公式。 在离散反应之间对带电粒子使用了连续能量损失跟踪急速。 对产生。 康普顿散射。 可以通过选项将瑞利散射包括进来。 光电效应 PEGS4是一个包含了12个子程序和85个函数的离线处理进程。输出数据直接可以被EGS4使用。 PEGS4在大量能量间隔之间进行了精确的线性拟和。 一般来说，用户只需使用一次PEGS4来获得EGS4所需的媒质数据文件。 EGS4是一个具有灵活使用者接口的增加了块（block）数据的子程序包。 不需要用户对其内部代码非常熟悉。 引入了宏，避免用户代码对EGS4代码构成bug。 任何给定问题的几何本质在用户写的子程序HOWFAR中给定，在HOWFAR中可以使用附加的子程序。 EGS4代码系统为不想写平面（plane）、圆柱体（cylinder）、圆锥（cone）、球（sphere）等标准几何形状的用户代码的用户提供了EGS4代码子程序。 上面子程序的宏在EGS4MAC文件可以得到，使用这些宏可以获得更快的模拟运行速度 可将MORSE-CG程序的组合几何包引入到子程序HOWFAR中使用（见UCSAMPCG文件）。但是，经验表明，这样做可能模拟速度很慢（大约慢四个量级）。 可以通过在HOWFAR中写一个特殊的子程序模拟在磁场中的传输。 使用者的记录和输出信息可以在用户写的子程序AUSGAB中记录。 实用程序ECNSV1用于监视和保存能量变化轨迹。 实用程序NTALLY用于保存事件。 实用程序WATCH用于按事件或按步模拟。 EGS4允许对重要抽样和其它变量衰减技术进行补充。（例如，主粒子偏移抽样，俄国轮盘赌与分裂技巧，步长偏移抽样等。） 辐射传输的初始化： 可以初始化对-0衰减的光子。（在CALL S