JUNO offline 电子学模拟所得的 MCP-PMT ENF 计算

豪森，

这里王宇逸(@strawberry_str) 用 JUNO offline 上观察计算 MCP-PMT ENF 的步骤，你看看是否能够复现。

使用最新版本 J22.2.0-rc2，具体的探测器模拟和电子学模拟命令如下：

$ python3 /cvmfs/juno.ihep.ac.cn/centos7_amd64_gcc1120/Pre-Release/J22.2.0-rc2/junosw/Examples/Tutorial/share/tut_detsim.py --evtmax 10000 --anamgr-normal-hit gun --particles e+ --momentums 0.0 --positions 0 0 -15000
$ python3 /cvmfs/juno.ihep.ac.cn/centos7_amd64_gcc1120/Pre-Release/J22.2.0-rc2/junosw/Examples/Tutorial/share/tut_det2elec.py --evtmax 10000 --rate 1.0

具体来说是模拟了动量为0的正电子，在 z=-15m 的点源。事例数单个文件为10000，我模拟了10个文件，总共十万个事件。

从文件 /cvmfs/juno.ihep.ac.cn/centos7_amd64_gcc830/Pre-Release/J21v1r0-Pre1/data/Simulation/ElecSim/PmtData_Lpmt.root 可以得知 PMT 的种类， 0号 PMT 是一个 MCP PMT 是 Hamamatsu PMT 。查阅 offline 的几何可知，0号的位置为 theta=2.7702, phi=186.857（都是角度），这是一个远离点源的 PMT，选取它是为了降低 occupancy。最终这个 PMT 的 occupancy 为 2.7%。

从 elecsim 输出的数据直接拿到每个 PE 的 charge，然后只取 PMT 0，拿到每个波形的 charge，得到附件1的直方图。

可以看到，offline 的长尾非常小。ENF 约为 1.09。 这是 Hamamatsu PMT 的结果。

本达

1.pdf