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

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

相关推荐


Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其他元素将获得点击?
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。)
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbcDriver发生异常。为什么?
这是用Java进行XML解析的最佳库。
Java的PriorityQueue的内置迭代器不会以任何特定顺序遍历数据结构。为什么?
如何在Java中聆听按键时移动图像。
Java“Program to an interface”。这是什么意思?
Java在半透明框架/面板/组件上重新绘画。
Java“ Class.forName()”和“ Class.forName()。newInstance()”之间有什么区别?
在此环境中不提供编译器。也许是在JRE而不是JDK上运行?
Java用相同的方法在一个类中实现两个接口。哪种接口方法被覆盖?
Java 什么是Runtime.getRuntime()。totalMemory()和freeMemory()?
java.library.path中的java.lang.UnsatisfiedLinkError否*****。dll
JavaFX“位置是必需的。” 即使在同一包装中
Java 导入两个具有相同名称的类。怎么处理?
Java 是否应该在HttpServletResponse.getOutputStream()/。getWriter()上调用.close()?
Java RegEx元字符(。)和普通点?