[论文解读] Studying the mass sensitivity of air-shower observables using simulated cosmic rays

Benjamin Flaggs, Alan Coleman+1 more|arXiv (Cornell University)|Jun 22, 2023

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

一句话总结

本研究利用 CORSIKA 模拟评估了在南极点和阿根廷马尔尔古埃地区，针对质子、氦、氧和铁初生粒子，空气簇射可观测量（如 $X_{\text{max}}$、μ子含量 ($N_{\nu}$) 和纵向剖面参数）的质量敏感性。结果表明，结合可观测量可实现强大的质子-铁分离，即使存在重建不确定性，但中等质量核素的区分仍具挑战性。

ABSTRACT

Using CORSIKA simulations, we investigate the mass sensitivity of cosmic-ray air-shower observables for sites at the South Pole and Malargüe, Argentina, the respective locations of the IceCube Neutrino Observatory and the Pierre Auger Observatory. Exact knowledge of observables from air-shower simulations was used to study the event-by-event mass separation between proton, helium, oxygen, and iron primary cosmic rays with a Fisher linear discriminant analysis. Dependencies on the observation site as well as the energy and zenith angle of the primary particle were studied in the ranges from $10^{16.0}-10^{18.5}\,$eV and $0^\circ$ to $60^\circ$: they are mostly weak and do not change the qualitative results. Promising proton-iron mass separation is achieved using combined knowledge of all studied observables, also when typical reconstruction uncertainties are accounted for. However, even with exact measurements, event-by-event separation of intermediate-mass nuclei is challenging and better methods than the Fisher discriminant and/or the inclusion of additional observables will be needed. As an individual observable, high-energy muons ($> 500\,$GeV) provide the best event-by-event mass discrimination, but the combination of muons of any energy and $X_{ ext{max}}$ provides already a high event-by-event separation between proton-iron primaries at both sites. We also confirm that the asymmetry and width parameters of the air-shower longitudinal profile, $R$ and $L$, are mass sensitive. Only $R$ seems to be suitable for event-by-event mass separation, but $L$ can potentially be used to statistically determine the proton-helium ratio. Overall, our results motivate the coincident measurement of several air-shower observables, including at least $X_{ ext{max}}$ and the sizes of the muonic and electromagnetic shower components, for the next generation of air-shower experiments.

研究动机与目标

• 评估空气簇射可观测量在事件级宇宙射线质量识别中的内在潜力。
• 研究 $X_{\text{max}}$、$N_{\nu}$ 和纵向剖面形状 ($R$、$L$) 等可观测量如何随初生粒子质量、能量和天顶角变化。
• 评估观测地点（南极点 vs. 马尔尔古埃）以及重建不确定性对质量分辨率的影响。
• 确定结合多个可观测量是否可克服单一可观测量在区分中等质量核素方面的局限性。

提出的方法

• 在能量范围 $10^{16.0}$ 至 $10^{18.5}$ eV 和天顶角 $0^\text{°}$ 至 $60^\text{°}$ 范围内，对质子、氦、氧、铁初生宇宙射线进行了 CORSIKA 模拟。
• 从模拟中精确提取了 IceCube（南极点）和 Pierre Auger（马尔尔古埃）观测点的 $X_{\text{max}}$、$N_{\nu}$、电子-μ子比值以及纵向剖面参数 ($R$、$L$)。
• 应用费雪线性判别分析以量化不同初生粒子种类之间的事件级质量分离能力。
• 采用参数化 Gaisser-Hillas 拟合方法建模纵向剖面，并施加质量筛选以确保拟合稳定性。
• 使用区分度量（FOM）以优化在两个观测点均实现质量敏感性的 μ子环带（800–850 m）的选择。
• 通过向可观测量施加系统性误差并评估其对分离能力的影响，对重建不确定性进行了建模。

实验结果

研究问题

• RQ1 在南极点和马尔尔古埃观测点，$X_{\text{max}}$、$N_{\nu}$ 和纵向剖面参数的质量敏感性如何随初生粒子质量变化？
• RQ2 能量和天顶角对质子、氦、氧和铁初生粒子之间事件级分离能力有何影响？
• RQ3 结合多个可观测量（如 $X_{\text{max}}$ 和 $N_{\nu}$）是否能显著提升质子-铁分离度，相比单一可观测量？
• RQ4 典型重建不确定性如何影响真实实验中可实现的质量分辨率？
• RQ5 纵向剖面的 $R$ 和 $L$ 参数在统计分析或事件级质量组成分析中可被利用到何种程度？

主要发现

• 在南极点和马尔尔古埃观测点，$X_{\text{max}}$ 和 μ子含量 ($N_{\nu}$) 的组合均能实现对质子和铁初生粒子的高事件级分离能力。
• 即使测量精确，使用费雪判别法对中等质量核素（如氦与氧）进行事件级分离仍具挑战性。
• 仅靠高能 μ子（>500 GeV）即可实现最佳个体质量区分，但结合低能 μ子与 $X_{\text{max}}$ 已能实现强大的分离能力。
• 纵向剖面的不对称性参数 $R$ 适用于事件级质量分离，而宽度参数 $L$ 更适合用于统计确定质子-氦比例。
• 参数化 Gaisser-Hillas 拟合比平移版本更稳定，因此在纵向剖面分析中更受青睐。
• 重建不确定性会降低分离能力，但即使在真实误差条件下，$X_{\text{max}}$ 和 $N_{\nu}$ 的联合使用仍能保持优异性能。