beamEGS参数设置

BEAMnrc基本参数设置1)Title为输入文件的标题，建议写得详细。当输入文件很多时，标题有助于查找合适的输入文件。另外如果输入文件的参数修改了，最好把标题也相应更新一下，以免误导。2)Medium是指环绕在加速器模型各元件模块周围的媒质，一般可选择空气。3)IWATCHOutput一般按默认的none就可以了。如果想结合EGS_Windows察看粒子径迹，可选择为forEGS_Windows（便于立体形象地显示你做的模型和结果，包括粒子的轨迹）。不过如果选择了forEGS_Windows参数，要注意设置模拟的历史数为一个较小值。可惜在windows下没有这个功能。4)OutputOptions可选择是否生成相空间文件，就是保存记录平面的各种参数的文件，个人认为很重要，一般都要选择，这样才可以做后期的数据图像处理。5)Numberofhistories表示要模拟的历史数（入射粒子数，一般为了得到误差小的结果，需要设置为105或者更高，不然画出的图像就会出现很到的误差bar，不过也可以使用bremstrahlungsplitting这个选项，下面说明）。6)MaximumCPUhoursallowed用于设置一个模拟计算程序最多可用的小时数。如果模拟过程执行的时间已经超过这个值，模拟过程将自动终止。一般可取一个较大的值，比如100。7)bremstrahlungsplitting这个选项好像是为了减少统计误差而增加的选项。在同样入射粒子数目下，结果的误差更小，建议使用。8)Incidentparticle表示入射粒子的类型。9)Sourcenumber表示选择何种粒子源。10)ECUT表示全局范围的电子截止能量，一般取默认值。11)PCUT表示全局范围的光子截止能量，一般取默认值。12)Numberofscoringplanes用于设置记录面个数、位置和尺寸。假设有两个相空间的文件，文件名相同以.egsphsp1，.egsphsp2不同的后缀结尾，其中数字1表示第1个记录面处产生的相空间文件，数字2表示第2个记录面处产生的相空间文件。EGSnrc基本参数设置1、Mediumstepsize：可以改变为更小。2、Boundarycrossingalgorithm（边界算法）：EXACT适合低能粒子，PRESTA-I对高能粒子有效。3、Electron-stepalgorithm：PRESTA-II更精确，PRESTA-I高能粒子时候使用。4、Bremsangularsampling：Simple只是使用Koch和Motz的2BS方程的主要部分，适合高能粒子，时间长。KM选项使用完整的2BS方程适合低能粒子。5、Bresmcrosssections：BH和NIST在粒子能量大于10MeV时候没有差别，在keV量级时候，NIST效果显著。（Bethe-Heitler是对50MeV以上电子轫致辐射截面公式的库仑修正）6、BoundComptonscattering：on选项特别费时，仅仅在能量低于1MeV时候或者使用Raleighscatting时候使用。7、Raleighscatting：仅仅在能量低于1MeV时候使用。8、AtomicRelaxations：低能时候使用。9、其余，我一般都取默认值。EGS中ctcreate的说明管理提醒：本帖被fada执行置顶操作(2008-06-02)希望下面的内容对大家有所帮助，需要提醒的一点是，注意执行文件的目录。说明书Pirs794讲了ctcreate，它可以读取CT片子，生成egsphant文件给dosxyz_nrc用作模体文件。所谓模体phantom，可以这样理解，一个立方体分成很多个小格子，每个小格子里面填充一种物质。ctcreate的执行顺序为：1.读入CT文件2.读入CTheader参数3.读入二进制CT数据4.选择CT数据的子集5.重新安排CT数据，得到满足要求的相应的需要记录剂量的区域6.将CT数据转换成物质和密度，对应每个小体元7.生成egsphant文件，给dosxyz¬_nrc用。Ctcreate的功能就是：1.选择原CT片子的一个子集作为输出文件的整体轮廓。2.重新细分网格小体元，重新给小体元分配物质（密度），用作后续记录剂量的小区域。以dicom格式的CT片子为例，例题中的slice_names包含了10个CT文件的路径名，将里面的$OMEGA_HOME改为实际的目录，这里为/home/egs/HEN_HOUSE/omega。CT_create_DICOM.inp为ctcreate的输入文件，里面记录了需要的数据。第一行，CT片子类型，这里为DICOM。第二行，CT文件的路径存放文件（slices_names）。第三行，新的模体文件的xyz三个方向的边界。第四行，重组模体的小体元的尺寸大小第五行，新模体的物质密度选项，默认为0。Ctcreate的运行有两种方式，读入输入文件×××.inp或是交互式输入。1.终端中输入ctcreate×××.inp，回车。执行结束，生成了×××.egsphant文件，给dosxyz_nrc用。2.输入ctcreate，回车。（1）选择CT类型，1pinnacle2CADPLAN3DICOM（2）输入DICOM文件名，这里为slice_names（3）给定新模体xyz三个方向的边界，单位cm，中间用逗号分隔。（4）输入新的小体元尺寸，单位cm。（5）输入media号，最大为7，如果使用hard-wired抽样函数，则输入0，一般输入0.（6）运行结束，生成了×××.egsphant文件，给dosxyz_nrc用。