2

我正在尝试在 zfit 中执行未合并的 3D 角度拟合，其中输入数据是具有从单独的不变质量峰值拟合分配的每个事件 sWeights 的样本。我想我在角相空间的某些区域遇到了负加权事件的问题，因为 zfit 给出了错误：

Traceback (most recent call last):
  File "unbinned_angular_fit.py", line 282, in <module>
    main()
  File "unbinned_angular_fit.py", line 217, in main
    result = minimizer.minimize(nll)
  File "/home/dhill/miniconda/envs/ana_env/lib/python3.7/site-packages/zfit/minimizers/baseminimizer.py", line 265, in minimize
    return self._hook_minimize(loss=loss, params=params)
  File "/home/dhill/miniconda/envs/ana_env/lib/python3.7/site-packages/zfit/minimizers/baseminimizer.py", line 274, in _hook_minimize
    return self._call_minimize(loss=loss, params=params)
  File "/home/dhill/miniconda/envs/ana_env/lib/python3.7/site-packages/zfit/minimizers/baseminimizer.py", line 278, in _call_minimize
    return self._minimize(loss=loss, params=params)
  File "/home/dhill/miniconda/envs/ana_env/lib/python3.7/site-packages/zfit/minimizers/minimizer_minuit.py", line 179, in _minimize
    result = minimizer.migrad(**minimize_options)
  File "src/iminuit/_libiminuit.pyx", line 859, in iminuit._libiminuit.Minuit.migrad
RuntimeError: exception was raised in user function
User function arguments:
    Hm_amp = +nan
    Hm_phi = +0.000000
    Hp_phi = +0.000000
Original python exception in user function:
RuntimeError: Loss starts already with NaN, cannot minimize.
  File "/home/dhill/miniconda/envs/ana_env/lib/python3.7/site-packages/zfit/minimizers/minimizer_minuit.py", line 121, in func
    values=info_values)
  File "/home/dhill/miniconda/envs/ana_env/lib/python3.7/site-packages/zfit/minimizers/baseminimizer.py", line 47, in minimize_nan
    return self._minimize_nan(loss=loss, params=params, minimizer=minimizer, values=values)
  File "/home/dhill/miniconda/envs/ana_env/lib/python3.7/site-packages/zfit/minimizers/baseminimizer.py", line 107, in _minimize_nan
    raise RuntimeError("Loss starts already with NaN, cannot minimize.")

我可以通过稍微限制一个拟合的可观察范围来避免这个错误，以避免出现少量数据事件的区域，其中一些数据被负权重（信号被 sWeights 稍微过度减去）。但我想知道在 zfit 中是否有另一种方法？

也许 zfit 中的 UnbinnedNLL 方法明确需要正事件，但负加权数据点可以设置为零或一个小的正值？我应该说，与权重的总和相比，负权重的水平似乎很小，并且发生在只有少量数据事件的角度分布之一的边缘。该区域的低数据率是由于实验接受效应。

在测试文件上运行以重现错误的代码在这里： https ://github.com/donalrinho/zfit_3D_unbinned_angular_fit_test

当costheta_X_VV_reco变量的范围限制为 (-0.9, 1.0) 而不是整个范围 (-1.0, 1.0) 时，不会遇到此错误。我相信这是因为它消除了加权数据为负的相空间区域。

4

2 回答 2

1

只是为了关闭这个线程，在某些地方 PDF 是负面的，我认为这是由于接受 PDF。确实h_pst没有使用 ins，但删除它并没有改变任何东西。最后，我刚刚在没有任何负 PDF 值的区域中拟合了数据，这似乎不会影响结果（它只是忽略了costheta_X密度接近零的小区域）。

于 2020-09-25T11:04:33.723 回答
1

从zfit中 NLL 的定义中可以看出，权重只是乘以对数概率，因此负权重应该不是问题。

但是，PDF 似乎对某些值返回负概率，您可以通过简单地获取返回的数组来查看

custom_pdf.pdf(data)

一旦获取日志，这个负概率将变成 NaN。

h_pst由于该变量似乎未使用，因此 PDF 的定义中可能存在拼写错误。

于 2020-08-28T12:50:42.080 回答