蒙特卡罗综合习题报告问题：蒙特卡罗方法模拟137Cs 源的661keV γ射线在NaI(Tl)闪烁体中的输运过程。
一、建立坐标系以闪烁体中轴线为z 轴建立柱坐标系，闪烁体位于0<z<4且r<2的区域。
Cs 电源位于z=-2，r=0点。
二、编程思想一个粒子的状态用(,,)E r Ω表示。
其中E 为粒子能量，(,,)r r z θ= 为粒子的位置坐标，(,)μϕΩ=为粒子的运动方向。
其中μ为粒子运动方向与 z 轴的夹角的余弦，ϕ为粒子运动方向在 x y 平面上投影的方位角。
对于一个粒子经历以下过程：1． 源抽样：由于是点源，能量和位置的分布均为δ函数，抽样得到 661,0,2E KeV r z ===-。
运动方向：各向同性分布密度函数为12121()()()411(),()22f f f f f μϕπμϕπ=⋅===Ω抽样得：1221,2μξϕπξ=⋅-=⋅。
2．到达闪烁体：如果μ≤，无法到达闪烁体，丢弃，返回源抽样重新产生粒子。
粒子进入闪烁体的瞬间，状态为：,0r z θϕ===其他量不变。
3．输运过程：1ln ()t m L E ξ-=-∑抽样得到到下次碰撞的距离，根据当前粒子状态中的,r Ω 算出下次碰撞的坐标(,,)r r z θ=，如果不在闪烁体区域（0<z<4且r<2），输运过程结束。
如果在，继续－>碰撞有两种可能：光电效应和康普顿散射。
根据粒子当前的能量，（由NaI(Tl)闪烁体宏观界面数据）线性插值确定它的光电效应截面和康普顿散射界面。
抽样得到本次反应的类型。
如果光电效应，E=0，输运过程结束。
如果康普顿效应，抽样获取碰撞后的能量和运动方向（康普顿散射的能量分布密度函数知道，具体抽样方法参考讲义。
）如果E<1KeV 输运过程结束，反之，重复本过程直到输运过程结束。
4．记录与统计：记录末态能量f E ，计算沉积能量0d f E E E =-，考虑到测量系统分辨率，多道记录能量为沉积能量的高斯展宽。
