先固定一份谁也不能改写的计数
设一台探测器在不同方向接收散射中子。每个角度只覆盖一个很窄的接收范围,束流、靶厚度、接收立体角与探测效率合并成有效曝光量。给定核势,就能从微分截面算出每个角度的预期计数。
例题固定质心能量 ,有效曝光量 。在二十、三十、四十、五十、六十度处,保存的合成计数依次是 。它们由上一章的势阱产生,伪随机种子为三七零七。生成后的计数保持不动,每个候选都解释同一份记录。
预期计数可以是小数,单次观测只能是整数。比如二十度的预期约为九十点六,保存值为一百零一;这种偏离就是计数涨落。记录还包含五个较大角度,起初把它们留作模型检验。
来源:A. M. Moro
把曲线像不像变成可计算的问题
泊松分布给出取得整数计数 的概率。假设各角度独立,把这些概率相乘就是似然。不同核势改变预期值,也改变这份记录出现的概率。拟合就是寻找似然较大的参数。
直接比较很小的概率不方便,可以取对数,再与“每个预期都恰好等于所见计数”的最大值比较。两倍对数差叫泊松偏差。它越小,所选模型越贴近同一组计数。这个统计量在零计数处也有定义。
误差棒用于快速看图。例题画出的误差棒取计数加一的平方根后除以曝光量,只是一个可见的涨落尺度;在小计数处,它不等于精确的泊松置信区间。拟合直接使用整数计数的泊松概率。
完整例题:三口势阱的前向图样
固定入射能量,搜索势深度、吸收深度和半径。一个局部拟合给出 ,前两项以兆电子伏计,最后一项以飞米计。五角度偏差为二点七九八。另一口很浅的势阱 ,偏差为二点八零四;一口很深的势阱 ,偏差也只有二点八二六。
这三组参数差别很大,拟合好坏却几乎相同。前向振幅里混合了多个分波的相位与大小,改变深度以后,可以借助半径和吸收补偿。有限角度的模平方又丢失了一部分相位信息,反推内部势便容易留下多种解释。
这里展示的是网格搜索加局部细化找到的几个解。候选之间相差多少,能说明多解有多严重;要给参数报一个标准差,还需要研究整个似然形状、实验涨落与模型误差。把候选范围直接写成置信区间会混淆这些问题。
来源:David Tong
把探测器转到候选分开的地方
原来留出的角度为 ,计数依次是 。把这些点加入同一个偏差统计量,上面三口势阱的偏差分别升到约三十四、九十六与四十。浅势候选在前向表现很好,却在新角度付出了很大代价。
再用全部十个角度重新拟合,找到 ,偏差约为八点零。这个解更接近生成记录的势阱,计数涨落仍会把最优参数推离生成值。扩大角度并未证明势阱唯一:另一些重新调整后的浅势解仍能给出较好的十角度拟合。
实验设计因此有一个具体目标:先列出尚能解释数据的候选,再寻找它们预言相差大、仪器又能可靠测量的位置。补测一个恰当角度,可能比在原来的角度上多积累很多相似事件更有用。增加入射能量也能提供独立信息,但本例把能量固定,以便只比较角度的作用。
沿着一个参数看剩下多少自由
若最关心势深度 ,可以每次固定它,再让半径和吸收自由寻找最佳值。这样得到一条深度剖面。半径与吸收在这个问题里叫干扰参数:它们必须参与拟合,却不是这次想单独报告的量。
例题在同一个有限网格上比较五角度与十角度的剖面,每条曲线都减去自己的最低值。谷底代表仍然合适的深度;多个谷底表示尚有不同解释。较宽的测量范围会改变谷底,但不会自动把似然变成单峰。
网格间距与局部极值会影响数值。尤其记录中有零计数,不能未经检验就把某个固定偏差阈值当作置信区间。若要报告区间,应细化优化,并按同样的采集过程重复生成记录,检查区间覆盖生成参数的频率。