BEAMnrc参数设置_MainInputs（3）

2019独角兽企业重金招聘Python工程师标准>>>

      全局截止能量, ECUT:

      当电子的总能量低于ECUT时，追踪结束，其能量沉积在当前区域。计算的时间很大地依赖于ECUT的值，所以尽可能地高。

      一般的规律是，ECUT的选取应该使得：能量为ECUT的电子的射程<13(the smallest dimension in a dose scoring region)

      全局的ECUT和各个部件的ECUT都应该>=材料的截止能量（也就是521icru或者700icru所设定的截至能量），否则是无效的。

一般情况都取0.7。

 

      光子的全局截止能量，要求不必精确，因为低能量的光子也不会花很多时间，通常取0.01 。

 

        电子的射程排斥（射程否决）：

        三个选项："off" (0), "on with varying ECUTRR" (1), and "on with set ECUTRR" (2).

        射程排斥用来节省时间。基本方法是计算带电粒子的射程，如果它不能离开当前的区域，则追踪结束，能量沉积在当前的点上。

        ECUTRR是带电粒子能够离开当前区域并且以大于ECUT的能量到达加速器底部所具备的最小能量。

 

        BEAM对于各种材料，计算了射程随电子能量（截至ECUT）变化的列表，射程的计算用了有限的阻止本领，给出了最长的可能射程。

        如果射程否决变量IREJCT_GLOBAL=2，射程否决将应用在区域上，当（截至ECUT的）射程<到最近的区域边界的距离时，进程结束，能量沉积的当前区域。

        如果IREJCT_GLOBAL=1，射程否决的应用将考虑带电粒子是否能够到达加速器的底部，并且仍然具有能量ECUT。这种情况下，在每个区域都会计算ECUTRR，当（截至ECUTRR的）射程<到最近的区域边界的距离时，进程结束，能量沉积的当前区域。

        IREJCT_GLOBAL=1相比IREJCT_GLOBAL=2节省时间，但是只有在记录平面只有一个且位于加速器底部时可以使用，其它情况下你可以对不同区域仔细选择ECUT来近似IREJCT_GLOBAL=1.

        射程否决是一种近似，因为在停止粒子进程、沉积能量在当前区域上，它假设了任何该粒子可能产生的韧致辐射光子都不能离开当前区域。不过你有办法减少这种近似带来的误差：使用输入变量ESAVE_GLOBAL，它定义了对带电粒子使用射程否决的最大能量（MeV）。ESAVE_GLOBAL的选取取决于入射射束的能量和介质是何种材料。ESAVE is treated internally on a region by region basis, 但是目前你只能对CM SLABS的不同区域（通过ESAVEIN）指定不同的值，因为CM SLABS是用作韧致辐射的靶，开发者认为这里可能需要更好地控制。

 

        能量大于该值的电子不必考虑上面所说的射程排斥，因为有足够的能量可以离开当前区域。

 

        强制光子与特定的CM发生作用。

        当光子作用很少的时候有利于提高散射光子的统计，例如薄板材料或者材料的密度很小，执行这个选项的目的之一是研究光子束中的污染电子的产生。

 

 

        通量（fluence）记录平面的个数：

        记录平面通常位于某一个CM下方，需要选择平面上的记录区（记录带）是方形、环形还是格子，记录区仅仅只是为了输出.egslst文件。

           

        Square： 

        或者grid：

        方形、环形需要指定区的数目，如果设为0，BEAM会自动赋予最大值，且每个区的面积相等。

 

 

转载于:https://my.oschina.net/u/1180521/blog/150869