[论文解读] Parametrization of the Cosmic Muon Flux for the Generator CMSCGEN

P. Biallass, T. Hebbeker|ArXiv.org|Jul 31, 2009

Astrophysics and Cosmic Phenomena参考文献 7被引用 18

一句话总结

本文为CMSCGEN生成器提出了一种新的宇宙μ子通量参数化方法，基于CORSIKA模拟，采用EPOS和GHEISHA相互作用模型，在3–3000 GeV能量范围内通过多项式拟合描述动量和角度分布。关键成果是获得了一种归一化通量参数化，其中在中心能量区间内不确定性为±7%，并经过实验数据验证，且已针对电荷比一致性进行调整。

ABSTRACT

The cosmic muon generator CMSCGEN is based on a parametrization of the differential muon flux at ground level, as obtained from the air shower simulation program CORSIKA. We present the underlying ansatz for this parameterization and provide an approximation of the momentum and angular distributions in terms of simple polynomials, in the momentum range 3 to 3000 GeV.

研究动机与目标

• 提高CMSCGEN生成器中宇宙μ子通量参数化的精度及其不确定性量化，以适用于对撞机实验。
• 解决先前参数化在高能和低能区域动量与角度分布方面的局限性。
• 引入更新的相互作用模型（EPOS、GHEISHA）和近期实验数据，以提升与实测结果的一致性。
• 确保正负μ子通量比的正确处理，这对探测器模拟至关重要。

提出的方法

• 将微分μ子通量参数化为 $ \frac{d\Phi}{dp\,dc\,d\phi} = C_{\text{norm}} \cdot \frac{1}{p^3} \cdot s(L) \cdot z(c,L) \cdot \frac{1}{2\pi} $，其中 $ L = \log_{10}(p/\text{GeV}) $。
• 使用六次多项式拟合动量谱函数 $ s(L) $，以匹配CORSIKA模拟在3–3000 GeV范围内的结果。
• 通过 $ z(c,L) = b_0(L) + b_1(L)c + b_2(L)c^2 $ 建模天顶角依赖性，其中 $ c = \cos\theta $，并在 $ c \in [-1, -0.1] $ 上归一化为1。
• 利用100 GeV处垂直μ子通量测量值校准 $ C_{\text{norm}} $，结合L3及其他实验数据，得出 $ C_{\text{norm}} = 1.27 \times 10^3 \, \text{GeV}^2 \text{m}^{-2} \text{s}^{-1} \text{sr}^{-1} $。
• 应用1.05的“修正因子”以改善在46和100 GeV处与L3数据的一致性。
• 将电荷比调整为 $ R = 1.280 \pm 0.016 $，其中66.1%的μ子指定为正，43.9%为负，以匹配实验通量。

实验结果

研究问题

• RQ1 多项式参数化在3–3000 GeV范围内对宇宙μ子动量和角度分布的描述精度如何？
• RQ2 使用现代相互作用模型（EPOS、GHEISHA）相比旧模型，对模拟μ子通量有何影响？
• RQ3 新参数化在多大程度上能重现实验μ子通量测量结果，特别是在垂直入射时？
• RQ4 正负μ子的适当电荷比是多少？在生成器中应如何实现以匹配数据？
• RQ5 参数化中的主要不确定性来源是什么，特别是在低能和高能区域？

主要发现

• 动量谱函数 $ s(L) $ 可由六次多项式良好近似：$ s(L) = -1 + 6.2218L - 13.940L^2 + 18.164L^3 - 9.2278L^4 + 1.9923L^5 - 0.15643L^6 $，在3–3000 GeV范围内拟合精度约为5%。
• 100 GeV处垂直μ子的归一化通量为 $ (2.63 \pm 0.06) \times 10^{-3} \, \text{m}^{-2} \text{s}^{-1} \text{GeV}^{-1} \text{sr}^{-1} $，用于校准 $ C_{\text{norm}} $。
• 归一化常数 $ C_{\text{norm}} $ 确定为 $ 1.27 \times 10^3 \, \text{GeV}^2 \text{m}^{-2} \text{s}^{-1} \text{sr}^{-1} $，并应用1.05的修正因子以与L3数据对齐。
• 正负μ子的电荷比设定为 $ R = 1.280 \pm 0.016 $，其中66.1%的μ子为正，43.9%为负，以匹配测量结果。
• 通量不确定性在中心动量区间（10–500 GeV）为±7%，在3 GeV处增至±25%，在3000 GeV处增至±50%，主要由于缺乏实验数据和模型不稳定性。
• 该参数化适用于 $ c \in [-1, -0.1] $，对应 $ \theta \in [0^\circ, 85.2^\circ] $，不应用于3 GeV以下或3000 GeV以上能量。