[论文解读] MUPAGE: a fast atmospheric MUon GEnerator for neutrino telescopes based on PArametric formulas

G. Carminati, M. Bazzotti+6 more|ArXiv.org|Jul 31, 2009

Astrophysics and Cosmic Phenomena参考文献 2被引用 10

一句话总结

MUPAGE 是一种用于水下/冰下中微子望远镜中大气μ子簇射的快速蒙特卡洛生成器，采用参数化公式模拟μ子通量、能谱和簇射多重性，计算效率极高。它将模拟时间减少了数个数量级，同时准确建模簇射特性，使 ANTARES、NEMO 和 KM3NeT 等探测器的大规模本底和校准研究成为可能。

ABSTRACT

Neutrino telescopes are opening new opportunities in observational high energy astrophysics. In these detectors, atmospheric muons from primary cosmic ray interactions in the atmosphere play an important role. They provide the most abundant source of events for calibration and for testing the reconstruction algorithms. On the other hand, they represent the major background source. The simulation of a statistically significant number of muons in large volume neutrino telescopes requires a big effort in terms of computing time. Some parameterizations are currently available, but they do not explicitly take into account the arrival of muons in bundles. A fast Monte Carlo generator (MUPAGE) was developed to generate single and multiple atmospheric muon events in underwater/ice neutrino telescopes. The code reduces the computing time for the simulation of atmospheric muons significantly. The event kinematics is produced on the surface of a user-defined cylinder, virtually surrounding the detector volume. The flux of muon bundles at different depths and zenith angles, the lateral spread and the energy spectrum of the muons in the bundles are based on parametric formulas. These formulas were obtained according to a specific primary cosmic ray flux model and constrained by the measurements of the muon flux in the MACRO underground experiment. Some MUPAGE applications are presented.

研究动机与目标

• 开发一种计算高效的替代方案，用于大型体积中微子望远镜中的大气μ子模拟，替代全蒙特卡洛模拟。
• 解决现有参数化模型中缺乏对大气μ子簇射多重性的考虑的问题。
• 实现探测器校准、本底估计和重建算法测试的大规模蒙特卡洛模拟。
• 通过提供快速、精确的μ子事例生成，支持 ANTARES、NEMO 和 KM3NeT 等中微子望远镜的设计与性能优化。
• 促进对μ子诱导本底效应的研究，包括在触发窗口内多个μ子引起的耦合事例。

提出的方法

• MUPAGE 使用从 CORSIKA 模拟中导出的参数化公式，并通过 MACRO 实验数据进行约束，以建模μ子通量、能谱和方位分布作为深度和天顶角的函数。
• 该生成器在用户定义的圆柱形“探测器罩”表面生成μ子运动学参数（位置、方向、能量），该探测器罩包围探测体积，深度范围为 1.5 至 5.0 千米水当量。
• 显式建模μ子簇射多重性，包括μ子能量随与簇射轴距离和簇射大小的变化关系。
• 代码使用 C++ 实现，结合 ROOT 库进行随机数生成，并输出格式化的 ASCII 表格，供后续模拟切伦科夫光和 PMT 信号使用。
• 该生成器支持可配置参数，如随机种子、运行编号和生成事例数，且模块化设计，包含独立的头文件、源代码和配置脚本目录。
• 可模拟来自同一宇宙射线相互作用的单个μ子和多个μ子，这对建模相关本底事例至关重要。

实验结果

研究问题

• RQ1 如何在不牺牲精度的前提下，加速大型体积中微子望远镜中大气μ子簇射的模拟？
• RQ2 参数化公式在多大程度上能准确再现深水下/冰中环境中μ子簇射的通量、能谱和多重性？
• RQ3 μ子簇射多重性对中微子望远镜中向上运动轨迹重建和本底率有何影响？
• RQ4 与全蒙特卡洛模拟（如 CORSIKA）相比，MUPAGE 在生成大规模事例样本方面的效率如何？
• RQ5 MUPAGE 是否能有效用于模拟多个μ子在短时间内（短于触发窗口）到达所引起的耦合事例，从而影响触发和重建性能？

主要发现

• 与全蒙特卡洛模拟相比，MUPAGE 将μ子模拟的计算时间减少了 10 倍以上，使大规模研究成为可能。
• 在 2×Intel Xeon 四核系统上，单个包含 10^7 个事例的 MUPAGE 文件生成时间少于 2 小时，证明了其极高的计算效率。
• 对于 ANTARES 5 线配置，MUPAGE 生成了一个月运行时间的数据集（分为 322 个文件），用于验证轨迹重建效率和角度残差。
• MUPAGE 模拟支持开发了一种去卷积方法，以提取作为深度函数的真实大气μ子通量，结果与已发表值一致。
• ‘耦合事例’率（即两个μ子在 4 μs 内到达，可能被重建为单个向上运动轨迹）测量为 7×10⁻⁴ Hz，90% 置信水平上限为 9×10⁻⁷ Hz。
• 在高能中微子本底研究中（E > 100 TeV），采用每簇射中每个μ子的保守能量阈值 3 TeV，一年运行时间的模拟在 232 小时内完成。