analysis/fsmp/batch.py: 修复在 nu_lc 的计算中 likelihood 公式的错误。 (!9) · Merge requests · TJJK / FSMP

Benda Xu requested to merge poisson-point-process into master May 25, 2024

xmass/waveform!82 (merged) 没有被充分 review，在 q(z|t_0, mu_0) 的绝对数值计算中引入了两个错误。Poisson point process 的 pseudo likelihood 的形式为

[ \nu_{lc}=\mathrm{e}^{-\mu_0}\prod_i^N{\mu_0\phi(t_i-t_0)} ]

连乘项为有序时间，转化为 Poisson point process 的采样后为无序时间，消去了 N! 项。delta log likelihood 的计算中，-log(NPE+1) 和 -log(NPE) 是随机取一个 PE 消灭掉的概率，不是 N!/(N-1)! 得到的。
不应该有 NPE * np.log(mu_t)，因为 v_lc 已经是 light_curve，带了 mu_t 项。

analysis/fsmp/batch.py: 修复在 nu_lc 的计算中 likelihood 公式的错误。