在RooUnfold框架中初始化2D分布的响应矩阵
我正在尝试使用T.Adye开发的RooUnfold软件包来展开2D分布。一维案例效果很好,我已经成功地使它成功了。现在,对于2D情况,我面临着问题。我想初始化二维分布的响应矩阵。 RooUnfold网站推荐的格式为
RooUnfoldResponse响应(nbins_measured,x_lo_measured,x_hi_measured,nbins_true,x_lo_true,x_hi_true);
或者,如果需要不同的装箱,
RooUnfoldResponse响应(hist_measured,hist_truth);
在最后一种情况下,例如使用hist_measured和hist_truth来指定分布的维数(此处不使用直方图的内容)。用于2D或3D分布或非均匀装仓。
Response_ProtonPion = new RooUnfoldResponse (301, -150.5, 150.5, 301, -150.5, 150.5);
但这不起作用,并且报告的错误是这样的:
virtual Int_t RooUnfoldResponse::Fill(Double_t, Double_t, Double_t, Double_t, Double_t): Assertion `_mdim==2 && _tdim==2' failed.