[論文レビュー] 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.

研究の動機と目的

• 全天頂角における地上における宇宙線ミューオン運動量スペクトルの簡単で普遍的なパラメータ化を開発すること。特に、浅い深さのニュートリノ実験に向けたものである。
• 浅い深さ（数100 m.w.e.未満）における低運動量ミューオン（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$ の運動量をカバーする、6つの表面ミューオン実験からの実験データを用いた。
• 系統的誤差と統計的誤差を平方和で結合し、重複する角度領域を有するデータはバイアスを避けるために慎重に処理した。
• 全データセットを用いて、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 らの既存のパラメータ化と比較して、どの程度正確か？
• 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 モデルよりも良好なフィットを示した。
• 数10 GeVのミューオンエネルギー（100 m.w.e.未満）を含むシミュレーションには、本モデルの使用が推奨される。
• 低エネルギーミューオンが無視できる深さでは、高エネルギー領域への滑らかな外挿が可能な、元のBugaev係数が好ましい。