振幅分析的系统误差和统计误差笔记

日边孤舟明

已于 2023-05-10 20:41:05 修改

阅读量332

点赞数

分类专栏：粒子物理 BES学习文章标签：笔记

于 2023-04-26 11:23:39 首次发布

本文链接：https://blog.csdn.net/Ares1121/article/details/130380143

版权

粒子物理同时被 2 个专栏收录

14 篇文章 4 订阅

订阅专栏

BES学习

12 篇文章 0 订阅

订阅专栏

1.FF的统计误差

To get the statistical uncertainties of FF, we float the parameters in fit results randomly 500 times according to the covariant matrix and fit uncertainty and calculate the FFs. Then, we fit the distribution of each FF with a Gaussian and the width of Gaussian is defined as the uncertainty of the FF.

为了获得FF的统计不确定性，我们根据协变矩阵和拟合不确定性，将拟合结果中的参数随机浮动500次，并计算FF。然后，我们用高斯拟合每个FF的分布，并将高斯的宽度定义为FF的不确定性。

体现在代码如下：

关闭拟合
随机FF500次

然后用高斯函数拟合该结果，所得的σ值即是该比份的统计误差，中心值应该是在比份的最佳解.

2.振幅的统计误差

拟合结果的Error就是振幅的统计误差

===============================系统误差===================================

我们做FF和 $\phi$ 的系统误差，是以统计误差为单位

振幅模型的系统误差Amplitude mode

根据PDG的值改变固定的振幅的质量和宽度 $\pm 1\sigma$ ，每个参数的最大相对变化的二次和被视为系统的不确定性。

计算公式 FF或者 $\phi$ 系统误差 = （改变后的值-最佳解初值）/ 统计误差

改变一个共振态的质量和宽度其他的共振态也会有影响，所以我们要做出上面的表格.

有效比率R Effective ratius

Variation of effective radii of Blatt-Weisskopf Barrier within the range $[2.0,4.0] Gev/c^{-1}$ for intermediate resonances and $[4.0,6.0] Gev$ for D mesons.

例如取R=2，D=6和R=4，D=4，两个极端去拟合两次，选出最大的变化那个

拟合偏差 Fit bias

也就是输入输出分为 $\phi$ 和FF两个的输入和输出

需要300份由振幅分析得到产生子的产生的Signal MC 作为假的data，相关程序在
```
/scratchfs/bes/zhous/share/02_readIO
```
余下的600份作为信号MC积分，但是算比份还是需要用PHSP MC积分
做 $\rho\phi$ 的输入输出需要将拟合纯度设为1.0，然后data的root替换为300份假data的root，用于做积分的余下600份sigMC 需要加入权重程序如下

获得权重的时候需要关闭拟合。DPFPWAEVA.cu文件里面添加
然后拟合300次获得结果代码如下
```
/scratchfs/bes/zhous/share/IO_rhophi
```

日边孤舟明

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
振幅分析的系统误差和统计误差笔记

为了获得FF的统计不确定性，我们根据协变矩阵和拟合不确定性，将拟合结果中的参数随机浮动500次，并计算FF。的输入输出需要将拟合纯度设为1.0，然后data的root替换为300份假data的root，用于做积分的余下600份sigMC 需要加入权重程序如下。===============================系统误差===================================然后用高斯函数拟合该结果，所得的σ值即是该比份的统计误差，中心值应该是在比份的最佳解.和FF两个的输入和输出。
复制链接

扫一扫