将 boost-histogram 写入在 uproot 中创建的根文件

问题描述

我目前正在尝试在根树中读取 root 树，使用 boost-histogram 预订和填充直方图，然后将这些写回新的 ROOT 文件。 这对于 bin 内容似乎非常有效，但我似乎无法在最终 ROOT 文件中获得直方图的正确 bin 错误。当我阅读这些直方图时，h.GetBinError(b) 似乎等于 h.GetBinContent(b)，而我想要熟悉的 sumW2 错误，这些错误通常是为 ROOT 直方图产生的。 我的工作流程如下，使用python 3.8：

import boost_histogram as bh
import uproot3 as uproot 
bins = np.array( [1,2,3,4],dtype='d' )
hist = bh.Histogram( bh.axis.Variable(bins))

然后我使用 uproot4.iterate 遍历输入树，为每个“块”生成一个笨拙的数组，这里用 e 表示：

hist.fill(e["pt"],weight=(e["weight"]))
with uproot.recreate(outfilename) as outfile:
   outfile["hist_pt"] = hist.to_numpy()

我不确定如何在构建 boost-histogram 时获得 sumw2 权重，但根据文档，我将直方图定义更改为：

hist = bh.Histogram( bh.axis.Variable(bins),storage=bh.storage.Weight() )

但是，当我运行我的代码时，这会导致 .to_numpy() 行出现问题。错误说： “what = uproot3_methods.convert.towriteable(what)” “TypeError: cannot perform reduce with flexible type”

是否有其他方法可以将 sumW2 样式的错误从 boost-histogram 传播到在 uproot 中创建的 ROOT 文件？或者我对这个 storage= 变量的理解有误.... 任何帮助将不胜感激。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）