在是用仿真软件(如Geant4 MCNP 等等)进行核仿真,获取的能量数据,会与现实测量不一致,这是因为,现实的探测器在测量时会存在一个展宽,电子学获取数据存在误差,展宽是必然存在的,有以下能谱展宽半高宽公式。
其中FWHM是半高宽,a、b、c是探测器展宽的三个系数。需要与实际的相结合。
使用matlab进行编程对原始的能量数据(非谱数据)进行展宽有以下函数代码。
代码块
function wide_data = wide(data,a,b,c)
%WIDE FWHM a + b * sqrt(E + c * E ^ 2) 为展宽公式
% data 原始能量数据 a b c 三个系数
% 对能谱数据进行展宽
% randn 产生一个正态分布随机数
for i = 1:length(data)
wide_data(i) = data(i) + randn()*(a + b * sqrt(data(i) + c *data(i)));
end
end
复制成功
使用展宽函数代码进行展宽得到下图。

展宽和未展宽的比较
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删