MCNP中如何记录剂量值啊？

Evan

琢磨了下手册，对剂量(dose)和通量(flux)间关系没有搞明白，请问谁能告诉我下？
另外，查了不少文献，就是没有找到能够明明白白、简简单单地告诉我如何由中子和伽马通量转化为剂量值。
请帮忙啊！


DEn 剂量能量卡
DFn 剂量函数卡

格式En A E1……Ek
DFn B F1……Fk
n ＝计数序号
Ei ＝能量(MeV)
Fi ＝剂量函数的相应值
A ＝能量表的LOG 或LIN 插值方法。
Ｂ ＝剂量函数表的LOG 或LIN 插值方法。

缺省：如果没有输入Ａ或Ｂ，则选择LOG 插值方法。
应用：建议使用计数注释卡。
这个特征允许用户输入一个点值响应函数（如通量剂量转换因子）作为能量的函数，用以修
改常规计数。这两个卡必须有相同的输入数据项，并且他们是按能量单调增加的。当粒子的能量
位于此卡定义的能量范围之外时，则使用此卡的最高或最低能量。
按省缺处理，对能点之间的响应函数值使用双对数值插值法给出，而不是象EMn 卡那样使用
直方图函数。在DEn 卡上指定用的能量点不必等于对Fn 计数指定的En 卡上的能量间隔。与使用
EMn 卡不同，在DEn 和DFn 卡上可以有许多能量点，但可能只在很少几个能量间隔内记录这个响
应，比如一个无边界的能量间隔（即一个总的间隔）．
如果在两个卡上的n 是零，这个函数将用于没有给出DEn 和DFn 卡的所有计数。
对这两个表每一个都可以单独选择 LIN 或LOG。这样，就可能有任何组合的插值方式（如果
对数――对数，线性――线性，线性――对数或对数――线性）。缺省的对数――对数插值方式
适用于ANSI／ANS 通量――剂量率转换因子（见［1］的附录H）、空气、水以及细胞组织的比释
动能因子、能量吸收系数。
例：DE5 E1 E2 E3 E4．．．．．．Ek
DF5 LIN F1 F2 F3 F4．．．．．．Fk
这个例子将按照剂量函数F（E）修改点探测器计数，该剂量函数在能量表上使用对数插值，
在剂量函数表上使用线性插值




我认为有个概念需要明确，这也是很多人的误区。
fluence，粒子注量，单位1/cm2，mcnp里对应的是dose equivalent剂量当量（ICRP60以后叫equivalent dose当量剂量）单位是Sv（或是rem）。
fluence rate （就是以前的叫法flux，我老板说flux已经作废，现在叫fluence rate，他说只是很多人还在用flux而已），粒子的注量率，单位是1/(cm2*s)，mcnp里对应的才是剂量率（Sv/s）。
mcnp里提供的注量－剂量转换因子（也可以叫通量－剂量率转换因子）就是把tally计数的注量转换成需要的当量剂量了。当然我们一般做光子的时候，吸收剂量和当量剂量相等，就无所谓了。




转自52mc