[論文レビュー] 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ジェネレータ用に、EPOSおよびGHEISHA衝突モデルを用いたCORSIKAシミュレーションを用いて、3–3000 GeVのエネルギー範囲で運動量および方位角分布を多項式フィットによりモデル化することで、宇宙線ミューオンフラックスの新しいパrametrizationを提示する。主な結果は、中央領域で±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}) $ である。
• 3–3000 GeVの範囲でCORSIKAシミュレーション結果に一致するように、運動量スペクトル関数 $ s(L) $ を6次多項式でフィットする。
• 仰角依存性を $ 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} $ を導出する。
• L3データとの一致を向上させるために、1.05の「フィッジファクタ」を適用する。
• 電荷比を $ R = 1.280 \pm 0.016 $ に調整し、陽性ミューオンに66.1%、陰性ミューオンに43.9%を割り当てて、実験測定値と一致させる。

実験結果

リサーチクエスチョン

• RQ1 3–3000 GeVの範囲で、多項式ベースのパラメータ化は、宇宙線ミューオンの運動量および方位角分布をどれほど正確に記述できるか？
• RQ2 最新の衝突モデル（EPOS、GHEISHA）を用いることで、古いモデルと比較してシミュレートされたミューオンフラックスにどのような影響を与えるか？
• RQ3 新しいパラメータ化は、特に垂直入射時において、実験的ミューオンフラックス測定値をどれほどよく再現できるか？
• RQ4 陽性ミューオン対陰性ミューオンの適切な電荷比は何か、そしてこの比をジェネレータに実装する際の適切な方法は何か？
• RQ5 特に低エネルギーおよび高エネルギー領域において、パラメータ化の主な不確実性要因は何か？

主な発見

• 運動量スペクトル関数 $ s(L) $ は6次多項式でよく近似され、$ 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} $ に決定され、L3データとの一致を高めるために1.05のフィッジファクタが適用されている。
• 陽性ミューオン対陰性ミューオンの電荷比は $ 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を超える領域では使用すべきではない。