0

我有一些数据正在通过 ROOT 的命令行界面在 TH2F 上绘制。我有一个名为 goodtree 的 TTree*,我在检测器中绘制事件的 XY 位置,如下所示:

root [1] TCanvas *can = new TCanvas("can","can",800,800)
root [2] goodtree->Draw("y:x>>h1(400,-200,200,400,-200,200)","r<200","colz")

我还想通过查看 r^2 与 theta 来制作归一化面积图:

root [3] goodtree->Draw("r*r:t>>h2(400,-3.14,3.14,400,0,41000)","r<200","colz")

这部分很好。我接下来要做的是将地图叠加到 XY 图上,并让它自动变形到 R^2T 图上的正确位置。我的意思是,这是一个粒子探测器,使用在 XY 方向具有圆形横截面的光电倍增管 (PMT)。我希望能够将地图叠加到 h1 上,显示这些 PMT 的轮廓(呈蜂窝状)。我也可以用脚本很快地做到这一点。

困难的一点是我希望能够在 XY 中定义此地图,并将其绘制在 R^2T 数据点的上方。有没有一种方法可以轻松做到这一点,而无需手动计算所有这些变形椭圆的位置、宽度和高度?

4

0 回答 0