如何解决将 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 举报,一经查实,本站将立刻删除。