蒙特卡罗方法模拟137Cs的γ能谱(共6页).doc

精选优质文档-倾情为你奉上蒙特卡罗方法模拟137Cs的能谱问题：蒙特卡罗方法模拟137Cs源的661keV射线在NaI(Tl)闪烁体中的输运过程。一、建立坐标系以闪烁体中轴线为z轴建立柱坐标系，闪烁体位于0z4且r2的区域。Cs电源位于z=-2，r=0点。二、编程思想一个粒子的状态用表示。其中E为粒子能量，为粒子的位置坐标，为粒子的运动方向。其中为粒子运动方向与 z 轴的夹角的余弦，为粒子运动方向在 x y 平面上投影的方位角。 对于一个粒子经历以下过程：1 源抽样：由于是点源，能量和位置的分布均为函数，抽样得到 。运动方向：各向同性分布密度函数为抽样得：。2到达闪烁体：如果，无法到达闪烁体，丢弃，返回源抽样重新产生粒子。粒子进入闪烁体的瞬间，状态为：其他量不变。3输运过程：抽样得到到下次碰撞的距离，根据当前粒子状态中的算出下次碰撞的坐标，如果不在闪烁体区域（0z4且r碰撞有两种可能：光电效应和康普顿散射。根据粒子当前的能量，（由NaI