[论文解读] Reconstruction of Atmospheric Neutrinos in Antares

A. Heijboer|arXiv (Cornell University)|Aug 6, 2009

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

一句话总结

本文提出了一种先进的重建技术，用于在ANTARES中微子望远镜中对大气中微子进行重建，利用完整似然拟合和在线算法，从深海水中切伦科夫光中精确重建μ子轨迹。该方法实现了亚度级角分辨率，向上行μ子和中微子的数据与模拟结果高度一致，验证了探测器的性能及其对大气中微子的灵敏度。

ABSTRACT

In May 2008, the Antares neutrino telescope was completed at 2.5 km depth in the Mediterranean Sea; data taking has been going on since. A prerequisite for neutrino astronomy is an accurate reconstruction of the neutrino events, as well as a detailed understanding of the atmospheric muon and neutrino backgrounds. Several methods have been developed to confront the challenges of muon reconstruction in the sea water environment, which are posed by e.g. backgrounds due to radioactivity and bioluminescence. I will discuss the techniques that allowed Antares to confidently identify its first neutrino events, as well as recent results on the measurement of atmospheric neutrinos.

研究动机与目标

• 开发适用于在具有高放射性与生物发光背景的复杂水下环境中运行的ANTARES中微子望远镜的稳健μ子轨迹重建算法。
• 通过考虑水中散射光子和次级粒子引起的延迟光子，提升角分辨率与重建效率。
• 通过将重建的大气μ子和中微子通量与模拟预期进行比较，验证探测器性能。
• 通过区分下行μ子与背景事例，实现对上行中微子候选事例的可靠识别。
• 提供一种快速、实时的重建方法（在线算法），适用于初步数据分析，同时保持通量测量所需的足够精度。

提出的方法

• 完整似然拟合使用基于模拟的、详细的、时间残差的概率密度函数（PDF），包含散射光子和次级粒子的影响，并通过脉冲幅度和光学模块（OM）方向对背景事例进行建模。
• 该方法采用解析梯度的迭代最小化似然函数，以提高收敛速度与精度，从多个初始参数猜测出发，避免陷入局部极小值。
• 预拟合阶段使用M-估计器对大残差事例（如背景）进行降权，从而在完整似然拟合前实现稳健的初始轨迹估计。
• 在线算法基于时间残差、幅度和到轨迹的距离，采用类似χ²的拟合方法，并引入一个综合残差与信号强度的质量度量Q，用于剔除质量差的重建结果。
• 对于单线事例，采用四参数拟合确定天顶角，同时引入亮点假设以剔除因电磁簇射而模拟出上行轨迹的事件。
• 两种算法均针对实时处理进行了优化，其中在线方法比完整似然拟合快约10倍，可在单个CPU上实现完整事例的重建。

实验结果

研究问题

• RQ1 在存在高背景放射性与生物发光的深海中微子望远镜中，能否实现μ子轨迹的精确重建？
• RQ2 在真实探测器条件下，完整似然与在线重建算法在角分辨率与鲁棒性方面如何比较？
• RQ3 重建的大气μ子与中微子通量在多大程度上与使用标准通量和相互作用模型的蒙特卡罗模拟结果一致？
• RQ4 在线重建算法能否在保持实时处理能力的同时，为初始通量测量提供足够性能？
• RQ5 时间残差分布在验证重建模型以及识别背景或信号事例方面具有何种重要性？

主要发现

• 完整似然拟合实现了亚度级角分辨率，在严格筛选条件下典型分辨率为1度，显著优于在线算法。
• 在线算法在高能区典型角分辨率为2度，比完整似然拟合差3至6倍，但足以满足初始通量研究需求。
• 2008年多线事例中，重建到582个上行事例，模拟预测为494（13）个来自大气中微子（μ子）的事例，比值为0.87，在系统误差范围内。
• 单线事例中，观测到237个上行事例，模拟预测为190（21）个来自大气中微子（μ子）的事例，比值为0.89，与多线结果一致。
• 重建中微子候选事例周围光点的时间残差分布与模拟结果高度吻合，证实了重建模型的准确性。
• 探测器表现出稳定性能，平均每天检测到约五个高质量上行中微子候选事例，证实其对大气中微子的灵敏度。