[论文解读] A Simple Parameterization of the Cosmic-Ray Muon Momentum Spectra at the Surface as a Function of Zenith Angle

D. Reyna|ArXiv.org|Apr 17, 2006

Neutrino Physics Research参考文献 4被引用 45

一句话总结

本文提出了一种简单且通用的参数化方法，用于描述地球表面宇宙射线μ子动量谱随天顶角的变化，采用缩放变量 $\zeta = p_\mu \cos\theta$。通过将垂直谱按 $I(p_\mu, \theta) = \cos^3\theta \cdot I_\text{v}(p_\mu \cos\theta)$ 进行缩放，该模型在所有天顶角和动量范围（1–2000 GeV/$\cos\theta$）内准确预测了微分μ子强度，实现了在204个自由度下 $\chi^2/\text{ndf} = 165.6$ 的拟合效果，且系数经过优化。

ABSTRACT

The designs of many neutrino experiments rely on calculations of the background rates arising from cosmic-ray muons at shallow depths. Understanding the angular dependence of low momentum cosmic-ray muons at the surface is necessary for these calculations. Heuristically, from examination of the data, a simple parameterization is proposed, based on a straighforward scaling variable. This in turn, allows a universal calculation of the differential muon intensity at the surface for all zenith angles and essentially all momenta.

研究动机与目标

• 开发一种适用于所有天顶角的简单、通用的地球表面宇宙射线μ子动量谱参数化方法，尤其适用于浅深度中微子实验。
• 解决浅深度（< 几百米水当量）下低动量μ子（3–20 GeV）缺乏准确模型的问题，此时μ子本底占主导。
• 提供一种可扩展、数据驱动的方法，将角度和动量相关的μ子强度与已知的垂直谱关联起来。
• 通过拟合涵盖多个天顶角和动量范围的全面表面μ子数据集，改进现有模型。

提出的方法

• 作者引入缩放变量 $\zeta = p_\mu \cos\theta$，表示μ子动量垂直于地表的分量。
• 提出函数形式 $I(p_\mu, \theta) = \cos^3\theta \cdot I_\text{v}(p_\mu \cos\theta)$，其中 $I_\text{v}$ 为垂直微分强度谱。
• 该方法使用来自六个地表μ子实验的数据，覆盖天顶角从0°到89°，动量范围从1到2000 GeV/$\cos\theta$。
• 系统误差与统计误差通过平方和开方合并，对存在重叠角响应的数据进行了仔细处理，以避免偏差。
• 通过在完整数据集上优化Bugaev等人函数形式中 $I_\text{v}(p_\mu)$ 的系数，获得最佳拟合参数化。
• 通过 $\chi^2$ 比较验证了模型，包括在有无角缩放假设下的表现。

实验结果

研究问题

• RQ1 是否能通过一种简单、通用的参数化方法，准确描述所有天顶角下地球表面宇宙射线μ子动量谱？
• RQ2 缩放关系 $I(p_\mu, \theta) = \cos^3\theta \cdot I_\text{v}(p_\mu \cos\theta)$ 是否能为地表μ子数据提供稳健且精确的拟合？
• RQ3 与Tang、Bogdanova和Bugaev等人的现有参数化方法相比，该模型的精度如何？
• RQ4 通过将所有地表数据（包括角度依赖性）联合拟合，能否改进Bugaev等人函数形式？

主要发现

• 所提出的参数化 $I(p_\mu, \theta) = \cos^3\theta \cdot I_\text{v}(p_\mu \cos\theta)$ 在所有天顶角和1至2000 GeV/$\cos\theta$ 动量范围内，对地表μ子数据提供了极佳的拟合效果。
• 通过优化Bugaev函数形式中的系数获得的最佳拟合模型，在204个自由度下达到 $\chi^2/\text{ndf} = 165.6$，相比现有模型有适度改进。
• 垂直谱的优化系数为 $c_1 = 0.00253$，$c_2 = 0.2455$，$c_3 = 1.288$，$c_4 = -0.2555$，$c_5 = 0.0209$。
• 该模型优于原始的Bugaev参数化方法，并且在通过角缩放关系应用时，其拟合效果优于Tang和Bogdanova模型。
• 该模型推荐用于μ子能量为几十GeV（小于100 m.w.e.）的模拟中，此时低能μ子具有显著贡献。
• 在更深的探测器中，当低能μ子可忽略时，建议使用原始Bugaev系数，因其在高能区具有更平滑的外推特性。