MCNP使用教程.docx

《MCNP使用教程.docx》由会员分享，可在线阅读，更多相关《MCNP使用教程.docx（18页珍藏版）》请在三一文库上搜索。

1、页眉内容第1章MCNP 既述1.1 MCN时算过程MCN(P Monte Carlo N-Particle Transport code)是计算粒子输运过程的一套蒙特卡罗模拟计算程序。 这个程序需要用户通过输入文件给出计算模型。 计算模型中需要提供源的属性、 感兴区内各种物体的属性、 记录粒子信息的方法等。例如，若想计算一个1MeV勺X射线透过2cm铁的概率是多少，我们可以通过下面的模型进行计算，如图 1 所示。图 1 计算模型在上面的计算模型中， 感兴区是一个球的内部， 其中包含X 射线源、铁块和记录面， 而其他位置均为真空。 由于当粒子被输运到感兴区外时， 它将肯定不会再对记录结果产生贡献

2、， 所以程序会自动停止这个粒子的输运过程， 这也正是设定感兴区的原因。源的属性主要包括位置、 能量、 出射方向、 粒子种类等。 图 1 的计算模型中，源的能量为单能1MeV方向为单向垂直于铁块的左表面，粒子种类为光子(Pho ton )。感兴区内物体的属性包括几何尺寸、 材料成分、 密度等。 图 1 中使用了一块铁块，它的厚度为2cm,其他方向的尺寸对我们的计算结果没有影响，但要保证 铁块完整地包含于感兴区内。记录方法有多种，其中包括通过某个面的特定种类粒子的个数。在图 1 中，我们可以利用MCNPE录通过“记录面”的能量为1MeV勺光子个数。计算图1的模型时，MCN至首先根据源的属性描述，抽

3、样出一个起始粒子。图 1 中的源为单能且单向的点源， 所以每次抽样出的粒子都是能量、 方向、 种类 相同的粒子。这个粒子会沿着它的出射方向（垂直于铁块左表面）飞行，当它入射到铁块里时， 会有一定的概率发生康普顿散射、 电子对效应和光电效应。 发生三种反应的概率由MCNP勺截面库中的微观截面数据、输入文件中铁的密度以及抽样得到的随机数共同决定。若X射线发生了康普顿散射，原来的X射线将被具 有新属性的X射线取代，它将有不同的出射方向、能量。 MCN除继续输运这个 新产生的 X 射线直到它发生下一次反应或者飞出感兴区；X 射线还会有一定的概率不发生任何反应，直接透过铁块。当基于一个起始粒子的输运过程

4、结束后,MCN会重复上述的过程。由于每次抽样的得到的随机数会不同， 所以每个起始粒子的输运过程都是独立而且不同的。多次的输运过程会使得 MCNP勺记录结果越来越接近它的期望值，重复的次 数越多，统计涨落就会越小。尽管MCN唆记录所有通过记录面的光子。但由于光子一旦发生反应，它的 能量就会变化，所以若只观察能量为1MeV勺光子，我们就可以排除那些经过多次散射而打到记录面的光子了。由于每次输运均是基于一个起始源粒子， 所以最后记录得到的结果也都是相对于一个起始粒子的概率。对于图 1 来说，记录结果就是1 个 X 射线透过 2cm的铁后打到记录面的概率。1.2 MCNP文件结构MCNPJ输入文件名不

5、可以多于七个字符,所以不建议对输入文件增加扩展名；MCN的输出文件白扩展名为.o； MCNPJ部计算过程的记录文件的扩展名为r 。若一次计算结果就已经达到满意的统计结果，那么 .r 文件就没有什么价值保存下来； 但是若一次计算结果达不到满意的统计可信度， 那么我们可以在已得到的计算结果的基础上继续模拟输运过程，来减小统计涨落。而为了这么做， MCNP 需要上一次计算结果中使用的 .r 文件。为了这个目的，有时候.r 文件还是值得保留的。在基于上一次的计算结果继续模拟计算时,MCN招不再读入输入文件，而单单读入 .r 文件的信息。所以我们想继续计算时，需要保证 .r 文件就是我们所设想的模型的前

6、一次模拟计算中得到的 .r 文件。除了上面三种常用的文件类型，MCN还有两种不长用的输出文件。这两种文件只有用户通过输入文件要求MCN输出这些文件时才会得到。其中一个为.w文件，这个文件是我们想通过MCN输出surface source的时候才会产生；另外一个为.p文件，这个文件时我们想通过MCN聆踪粒子的输运过程的时候才会得到。上面的几种文件是用户常用到的文件类型，下面将介绍MCN程序内部使用的若干的文件。1. MCNP5.exe： MCN的应用程序，这个程序是在 dos下运行的程序。所以用户不要试图双击这个程序来运行它。2. vised.exe:这个程序可以让用户观察 MCNPt入的输入文

7、件中的计算模型。 我们可以利用这个程序观察我们所设想的计算模型与实际书写出来的计算模型是否有偏差。止匕外，这个程序还可以方便我们找出MCN输入文件中潜在的一些错误，例如：两个物体有相交的部分等。3. 目录 Documents: 这个文件夹内包含LANL(Los Alamos NationalLaboratory)提供的MCNP勺官方使用说明、林谦老师的蒙卡课程讲义和西安 交通大学翻译的中文MCNP勺使用手册等4. 目录MCNPData:这个目录里存放的是 MCN的截面库。5. xs52: 这个文件内说明了各个核素的属性，包括质量，截面库内各种反应类型对应的截面的存放位置等。在初始使用MCN时，

8、我们有可能要修改这个文件中的一部分内容，具体的修改方法将在下面介绍。6. X11.dll : MCN程序的动态链接库，没有这个文件 MCNPI无法计算，所以不要把这个文件弄丢了。除了上面的文件外，为了在dos下使用MCN程序我们编写了几个批处理文 件。1) ccmd.bat: 只要双击这个程序就可以弹出一个dos 界面。2) g5.bat: 这个批处理将是用户用到最多的批处理程序。我们可以利用这个批处理直接调用MCNP勺主程序。3) gogo.bat :这个批出文件可以使MCNP1续的计算若干个输入文件，具体 的使用方法将在后面介绍。4) resume.bat:这个批处理文件可以使得MCN而上

9、一次计算结果的基础上 继续进行计算，具体使用方法将在下面介绍。上面所提到的文件就是用户在使用 MCN时最频繁接触到的文件。1.3 MCNP使用方法在使用MCN时，我建议用户安装一个比较方便的文本编辑器，如： Ultrae dit 。为了使用MCNP用户需要首先进行一些文件路径的指明等工作，具体如下：1）建议用户把MCN的整个目录（名为LANL放在一个硬盘分区的根目录下。2）修改g5.bat中文件路径：我们用Ultraedit打开g5.bat可以发现， 这个批处理的内容为前5句的作用是删除与输入文件相关的已经存在的输出文件，这样做的目的是，如果MCN程序发现已有输入文件对应的输出文件存在，那么程

10、序就认为这 个输入文件已经得到了模拟计算并且没有必要进行重复的计算。批处理中的最后一句就是调用 MCN的主程序一一MCNP5.exe并指明管理截 面库的文件（xs52）的路径。由于我已经把这个文件放在了 LANL的目录内，所以 用户不用修改这个文件的路径了。3）修改xs52文件：利用Ultraedit 打开xs52文件，可以发现文件的第一 行为。为了让MCN找到截面库，用户需要正确地修改它的路径。因为我把LANL放在了 E盘的根目录下，所以这里的datapath分区的根目录下，所以在这里修改时，用户只需把e更；经过上面的修改后， 用户就可以使用MCNP?序了； A,那么我们可以通过下面的步骤进

11、行计算；1） 双击ccmd.bat,系统会弹出下面的窗口； 2）在窗口内输入g5?Xray7并接回车， MCN计算时MCN脍占据单个CPU勺所有资源；B,关于停止一次计算，可以通过 在输入文件内限定计；在这个界面下，输入 s并按回车，程分区的根目录下， 所以在这里修改时， 用户只需把e 更换为实际存放的硬盘分区的盘符就可以了。经过上面的修改后，用户就可以使用MCN程序了。例如我们已经书写好了一个文件名为 Xray7 的输入文件，A,那么我们可以通过下面的步骤进行计算。1） 双击 ccmd.bat ，系统会弹出下面的窗口。2）在窗口内输入g5?Xray7并接回车，MCN就开始进行计算了。界面如下

12、图所示。计算时MCN强占据单个CPU勺所有资源。若用户的计算机是单核的，那么CPUS用量是100%若用户的计算机是双核的，那么CPU勺使用量是50%B, 关于停止一次计算，可以通过在输入文件内限定计算时间或输运的起始粒子个数。 此外， 还可以通过dos 下强制停止计算过程的方法。 若想强制停止计算过程，可以使用 Ctrl+C 键停止计算。同时按下这两键时，屏幕会出现下面的界面。在这个界面下，输入 s 并按回车，程序会输出已经消耗的计算时间；按下m是调用MCNP勺画图软件，因为有了 vised.exe程序，所以不建议使用这里的画图程序；按下q是停止计算程序，并输出结算结果。有时候，MCN脍等待正

13、在模拟的一个粒子的输运过程结束后在结束程序，但并不会拖延太长时间；按下 k 是直接停止计算程序，而并不会输出计算结果。C, 继续计算一个已计算完的模型。用户需要保留初始计算过程中产生的 .r文件。还以上面的计算情况为例，在上一次计算中我们得到了 Xray7.o （输出文件）和Xray7.r （计算过程的记录文件）。在 dos命令行下输入resume?Xray7 并按下回车，可以得到下面的界面。在继续运算的情况，由于MCN不会再次读入输入文件，而仅仅读入.r文件, 所以这次计算的终止需要我们通过强制停止的方法。D, 连续计算多个软件这里只给出了计算一个文件的命令，若我们想连续地计算多个文件（如D

14、1，D2然后在dos命令行下，键入gogo并按下回车，MCN就可以依次计算D1,D2 和 D3 了。第2章MCNP俞入文件书写方法2.1 MCNP输入文件规范MCNPJ输入文件在结构上基本上可以分为三块，分别为cell块、surface块和其他内容。块内除了注释外，不得有空行，而两个块之间只能有一个空行。通常情况下， 每一行中只能有一个主导性的助记符， 这个助记符可以有很多参数和附属的次要助记符。MCNPI方手册中称这卞的助记符为卡（card）。每一列不得超过 80 个字符， 若书写不下时， 可以在行尾写入 &符号， 并从下行继续书写。输入文件的第一行不论写了什么都会被 MCN眼略掉，这一行的

15、作用是方便 用户给这个输入文件做一些备注。第二行后写计算模型，我的习惯是先写 cell 块，然后是surface 块，最后是其他。输入文件中有两种注释方法，助记符分别为 c 和$符号。在第一行写 c 就是把这一行都注释掉；在某一行的某列上写$就是注释该行该列后的内容。但是即使是注释，请用户也不要超过每一列 80 个字符（包括空格）的要求。还有三个常用的助记符在这里介绍个用户，分别为 r, i和j。r代表的是r epeat,例如1 3r就等于写了 1 1 1 1 ; i代表线性插入点，例如1 2i 4 就等 于写了 1 2 3 4 ; j代表的是jump,表示该位置使用默认值，例如1 j 2就等

16、于 写了 1默认值2 。由于MCNPt输入文件的规范要求苛刻，所以请用户谨记这样的规范。2.1 7MeV初致辐射X射线算例下面将通过书写计算7MeV勺物致辐射X射线的角分布、能谱的方法介绍MC NP输入文件的书写方法。在书写MCN输入文件前，请用户首先想好要计算的模型，包括坐标系以及 各个物体在这个坐标系内的相对位置。我的计算模型如下图所示。图2计算模型；在图2的计算模型中感兴区是半径为100cm的球内；为了书 写上面的计算模型，我们首先写 surface ; 1,源所在面的书写方法：1PX-2；直 于X轴的面；-2代表的是该面在X轴上的轴距；2,同理Au的左底面的书写方 法：2PXQ 3,同

17、理Au的右底面的书写方法：3PX1; 4, Au的侧面的书写方法： 4CX1;到类似地,还会有CY和CZ的助记符）；5,圈图2计算模型在图2的计算模型中感兴区是半径为100cm的球内，源是能量为7MeV位 置在（-2,0,0）的单向电子束；X射线靶是半径（Y或者Z方向）1cm,高度（X方 向）1cm的金；其他位置均为真空。我将记录面设置在圈定感兴区的球面上。为了书写上面的计算模型，我们首先写 surface 模块。这里需要定义的 sur face 有源所在面； Au 的三个表面（ 两个底面和一个侧面） ，圈定感兴区的球面；划分感兴区表面的记录用面。1，源所在面的书写方法： 1 PX -2 。

18、这里的第一个1 代表这个面的序号，这个序号将在书写 cell 块时用到； px 为助记符表示该面是垂直于X轴的面；-2代表的是该面在X轴上的轴距。2，同理Au 的左底面的书写方法：2 PX 0 。3，同理Au 的右底面的书写方法：3 PX 1 。（可以猜到类似地，还会有PY和 PZ 的助记符）。4, Au的侧面的书写方法：4 CX 1。这里的4代表面的序号，cx代表这个面 是轴心是X 轴的圆柱面，后面的1 代表圆柱面的半径。（可以猜到类似地，还会有CY和CZ的助记符）。5，圈定的感兴区的面的书写方法：100 SO 100。第一个100 仍然代表面的序号。定义面的需要时，并不需要连续地定义，但要

19、保证定义的顺序是单调递增的。第二个so 代表的是圆心在原点的球面，之后的100就是圆的半径。6，除了这些面，为了记录时把记录面分割为若干个小面，我们需要额外定义一些分割用面。为了把记录面（ 球面 ） 分割为若干个球带，我们可以定义若干个圆锥体，如 30 kx 0 0.1 1 ，其中 30 代表这个面的序号；kx 代表这个面是轴心是X 轴的圆锥面；0 代表圆锥的顶点在原点；0.1 代表圆锥体顶角的tan 值的平方值；我们可以想象，这样定义的圆锥体应有两叶，在 Z0 的一叶。所以最后的 1 代表的是，我们想要的圆锥面是Z0 的一叶。我们可以类似地定义很多这样的分割面。有了这些 surface 的定

20、义后，我们就可以开始书写 cell 块了。1） 首先定义 Au 块 （cell） ，它由三个面围成，它们是序号分别为 2， 3， 4。Au块由2号面的向右方向（或X轴的正方向），3号面的向左方向（或X轴的负方 向）以及 4 号面的内部围成。所以这个Au 块（ cell ）的书写方法是： 1 1 -19.3 2 -3 -4 。第一个 1 代表的是这个cell 的序号；第二个1 代表的是这个cell所对应的材料种类为 1，关于 1 号材料的定义将在后面介绍；第三个-19.3 中，负号代表的是实际密度，若没有负号则代表的是原子密度。所以 -19.3 代表的是Au块（cell ）的密度是19.3g/c

21、m3;之后的2 -3 -4描述的是Au块是由2号面的正方向， 3 号面的负方向， 4 号面的内部围成的。2） 源不是一个物体，所以不需要在cell 块中定义。而在感兴趣内，除了 Au 块（ cell ），其他位置均为真空。尽管是真空，我们仍需要在输入文件内把真空定义为一个物体。定义的方法如下： 2 0 -100 #1 。其中第一个2 代表 cell的序号；第二个0 代表定义的块内什么东西都没有；后面的 -100 #1 限定了这个cell 的几何空间，它表示这个块是在100号面的内部，（#助记符代表扣除的意思）除了 1 号 cell 的所有空间。3） 我们同样需要定义感兴区外的空间， 定义的方法

22、如下： 100 0 100 。 第一个 100是这个 cell 的序号；第二个0代表真空；第三个100代表第 100 号面的外部。2号 cell 和 100号 cell 均为真空， 区分哪一个是感兴区的方法是利用另外一个助记符imp,它代表的是importance ,也就是粒子在该cell内的重要性。这个助记符可以在每个cell 的定义的末端写上，也可以在第三块中写入。在cell 末端写入的方法如： 2 0 -100 #1 imp:e,p=1 ， 它代表了 2 号 cell 中电子 (e)及光子 (p) 的重要性为 1，即它是感兴区内的物体； 100 0 100 imp:e,p=0 ，它代表了

23、 100 号 cell 的电子及光子的重要性为0，即它在感兴区外。注： 感兴区与非感兴区并不需要一定是位置划分分明的两块区域， 我们可以在感兴区内设定一个非感兴区， 只要粒子被输运到该区域就会被杀死 (类似于黑洞)。最后我们需要书写第三个模块， 在这个模块内将指明模拟计算中所涉及的反应类型，记录方法，运算时间限定等。1 . 模式卡： mode e p 。这个卡说明了计算中将涉及两种粒子，分别为电子及光子。(源粒子是电子，它轰击Au 块会产生轫致辐射X 射线)。2 .材料卡：ml 79197 -1.0。其中ml代表第一号材料，同理会有 m2 m3等材料； 79197的格式是ZZAAA， 即 79

24、是原子序数(Au) ， 197是质量数；-1.0 中负号代表的是原子个数的比重，若正好则代表的是质量比重。材料卡的各种核素的比重的标识中并不要求归一化，MCN强自动给用户归一化。也就是说ml79197 -2.0 79198 -2.0 与 ml 79197 -0.5 79198 -0.5的写法对于 MCN而言是等效的，都说明1号材料中197Au与198Au各占一半。3. 源定义： sdef par=3 sur=1 pos=-2 0 0 vec=1 0 0 dir=1 erg=7 。 其中 sdef 是源定义的助记符代表source definition ； par 代表源粒子的种类， 1 代表中

25、子， 2 代表光子， 3 代表电子； sur 代表源所在的面，上面定义中 sur=1 说明源在序号为 1 的面上； pos 代表源所在的位置， 按照上面的定义， 源的位置为( -2,0,0 )。请谨记若sur 和 pos 同时被用上，请保证pos 的点是在 sur 的面上，不然可能导致计算结果有误； vec 代表参考向量，向量的起始点是原点，而终止点是 (1,0,0) 。这说明了我们设定的参考向量就是X 轴的正方向； dir 代表粒子出射方向与参考向量的夹角的余弦，上面定义中为 1，说明了源粒子沿着X 轴正方向出射；erg代表源粒子的能量，按照上面的定义，源粒子的能量为 7MeV关 于源的定义

26、还有许多其他的助记符以及各种分布情况，请参考MCNP勺官方使用手册。4. 粒子重要性表示： imp:e 1 1 0 。 其中 imp:e 代表了这张卡描述的是各个cell 中电子的重要性， 1 1 0 代表了在 cell 块中写的各个cell (以书写顺序)的重要性分别为 1、 1、 0。这说明了cell 块中第三个描述的 cell ( cell 的序号可能不是三，只是代表它是第三个被定义的 cell 。)是电子的黑洞，只要电子被输运到这个cell 就会被杀死，停止输运。同理我们还可以定义imp:p 11 0 。在前面已经描述过，这里的 imp 卡可以在 cell 块中写入如：5 .能量截断卡

27、：cut:e j 1。其中cut代表能量截断卡的助记符；e代表电子；j代表第一个参数选用默认值；1代表截断能量为1MeV这个卡的作用是当 输运过程中，电子的能量小于1MeV就杀死这个电子。这样做的目的是，若我 们只关心高能（1MeV的X射线，而能量小于1MeV勺电子已经不会再产生能量 大于1MeV勺X射线了，所以我们可以停止这些低能电子输运来节省计算机机时 而且不会影响计算结果。同理我们还可以写 cut:p j 1。6 .记录方法描述卡1: F1:p 100。其中F1代表的是第一种记录模式，也就是通过某个面的粒子的个数（相对于一个源粒子）； p 代表记录的是光子；100代表记录面是序号为 10

28、0 的面。 故这个记录卡的作用是记录所有通过100号面的光子。7 .记录方法描述2: F2:p 100。其中F2代表的是第二种记录模式，也就是通过某个面的粒子的注量，它等于通过面的粒子个数（相对于一个源粒子）除以该面的面积。8 . 记录面分割： Fs2 -30 -31 -32 -33 -34 -35 -36 -37 -38 -39 40 4142 43 44 。其中Fs2代表的是分割第二号记录卡的记录面（也就是分割记录第100号面）；后面的-30等代表分割方法，输出时 MCN至给出1） 100 号面上，从30号面的内部（圆锥体右叶的内部）出射的粒子2） 100 号面上，从30 号面的外部且40

29、 号面的内部出射的粒子3） 100 号面上， 从 30 号面的外部且 31 号面的外部且 32 号面的内部出射的 粒子4） 100 号面上，从30 号面的外部且31 号面的外部且32号面的外部且33号面的内部出射的粒子。5） 以此类推9） 能量分割：我们可以对任意一个记录卡分别设置能量分割卡或所有记录卡同意设置相同的能量分割卡。 这里我们将对2 号记录卡设置一个能量分割卡：E2 1 19i 7。其中E2代表这是对2号记录卡的能量分割卡，后面的1 19i7代表的是分割方法，也就是在1MeV 7Med甜S入19个点。（若把E2写成E0则代表对所有记录卡使用相同的能量分割方法。）10） 计算时间设定

30、卡：ctme 10。这个卡说明了要求 MCN针算10分钟。可以想象， 计算的时间越长， 模拟的源粒子的个数就会越多， 模拟结果就会越接近它的期望值， 统计涨落就会越小。 所以计算时间需要用户根据自己希望得到的统计可信度，设定计算时间。除了 ctme 还有利用模拟源粒子的个数进行时间限制的方法，如NPS 10000,它说明了需要MCNPU1 10000个源粒子的输运过程。在这个文件中只使用了两个记录方法，分别为F1和F2。MCN呻一共有8中记录方法， 用户请参考使用手册， 选用适合计算模型的记录卡。 此外一种记录卡可以有多个记录面，只要保证记录卡的个位数是8 种记录卡的一种，例如 F11和 F2

31、1 都是第一种记录类型， 即通过某个记录面的粒子个数。 另外请谨记， 在使用 F1 和 F2 的记录卡时，记录面一定要是构成某个cell 的一个面，不然MCNP不会给出记录结果，也就是说如果用户随意定义了一个面（ 如上面定义的某个分割面），而且这个面没有被使用在定义某个cell上，那么我们不能期望MCN阿以记录通过这个面的粒子信息。使用分割面时，其实还是记录通过那个记录面的粒子个数，分割面只是把那个记录面分割为若个块而已。总结上面的书写过程，我们最后的输入文件如下。文件名为Xray7。我们可以利用vised.exe观察这个已经写好的；误，我们需要修改输入文件； 关闭；若打开后，观察到明确的红线

32、，代表两个物体分享了一；中是不允许的； 与我们想像的模型是一致的；入文件的正确性，不然计算出来的结果就很有可能 是错；Vised.exe是windows下的应用程序，；好的输入文件，就可以观察写好 的模型了，如下图所示；图3宏观形状；图4Au块的形状；利用前我们可以利用vised.exe观察这个已经写好的输入文件的模型，如果书写有误，我们需要修改输入文件。有错误时会导致vised.exe无法打开文件并被强制关闭。出现这种情况时，请看输入文件的书写规范，如空行是否多了或少了若打开后，观察到明确的红线，代表两个物体分享了一部分空间，这在MCNP中是不允许的。在书写复杂的模型时，vised.exe可

33、以帮助我们确定写出来的模型与我们想像的模型是一致的。建议用户反复地观察各个方向的视图，从而确 保输入文件的正确性，不然计算出来的结果就很有可能是错误Vised.exe 是 windows 下的应用程序， 我们双击它并用它选中我们已经书写好的输入文件，就可以观察写好的模型了，如下图所示。图 3 宏观形状图 4 Au 块的形状利用前面介绍的调用MCNP勺方法就可以计算这个模型了。下一章帮助各位解读输出文件的内容。第3章MCNP俞出文件解读利用MCN针算上面的输入文件可以得到输出文件为Xray7.o。利用Ultraedit 可以看到计算结果，其中关于记录信息的结果对我们最有用，如 下。其中 nps

34、代表一共运行了 2958959个源粒子的输运过程。Tally type 1 代表记录类型为 1， 也就是 number of particles crossing a surface 。Tally for photons 代表的是记录的是光子。Surface 100 代表记录面是第100 号面。记录结果是2.90042E-01 0.0010。这个结果表示当有1个7MeV勺电子轰击1cm的Au时，会有0.29个能量大于1MeV（我彳门在cut中限定了能量）0.001光子透过第 100 号记录面，而这个数据的相对误差为第二号记录结果如下。其中 nps 代表一共运行了 2958959个源粒子的输运过

35、程。Tally type 1 代表记录类型为2，也就是particle flux averaged overa surface ，请注意它的单位是1/cm2。Tally for photons 代表的是记录的是光子。Surface 100 代表记录面是第100 号面。Segment代表的是利用分割面分割后，100号面的各个子面的面积。由于我们使用的等立体角分割，所以各个子面的面积均相同。其中 surface 100 segment: -30代表统计的是从记录面100 且从分割面30 的内部（圆锥右叶的内部）通过的光子。同理surface 100 segment: 30 -31代表的是统计的是从

36、记录面100 且从分割面30 的外部且从分割面31 的内部通过的光子。最左边一列代表的是能量， 第二列代表的是概率， 第三列代表的是第二列的相对误差。1.0000E+00 0.00000E+00 0.00001.3000E+00 1.58985E-06 0.00581.6000E+00 1.25669E-06 0.0065对于这三行统计结果的解读如下，在相应的记录面内，1）能量为01MeV勺光子的注量为0,这是因为我们把能量小于1MeV勺光子都杀死了（利用 cut 卡）。2）能量为11.3MeV的光子的注量为1.58985E-6,单位是1/cm2,且是相对于一个电子。3）能量为1.31.6MeV的光子的注量为1.25669E-06,单位是1/cm2,且是相对于一个电子。4） 以此类推。到此为止，我相信各位初学者已经初步了解了MCN效用程序的调用方法、输入文件的书写格式、 输出文件的解读方法等内容。 有了这些知识后， 各位用户就可以较容易地通过阅读官方使用手册来提高自己在编写MCN输入文件的台匕能力了。 希望这篇文档对各位初学者在入门中给予了一些帮助。祝各位好运！