[論文レビュー] Reconstruction accuracy of the surface detector array of the Pierre Auger Observatory

Auger Collaboration, Ave, M.|ArXiv.org|Sep 13, 2007

Astronomical Observations and Instrumentation被引用数 10

ひとこと要約

本論文は、広大な空気シャワーの再構成精度、特に方位角分解能とS(1000)を用いたエネルギー推定に関して、ピエール・オーガー観測所の表面検出器アレイを評価している。イベントごとの再構成と完全なモンテカルロシミュレーションを用いて、3駅イベントでは2°未塔の分解能が達成され、より高い密度では1°未塔にまで向上することが示され、高エネルギー領域ではS(1000)エネルギー推定誤差が8%未塔であり、主にシャワー発達のゆらぎに起因する。

ABSTRACT

The reconstruction of extensive air showers (arrival direction, core position and energy estimation) by the surface detector of the Pierre Auger Observatory is discussed together with the corresponding accuracy. We determine the angular reconstruction accuracy as a function of the station multiplicity by using two different aproaches. We discuss statistical and systematic uncertainties in the determination of the signal at 1000 m from the core, S(1000), which is used to estimate the primary energy.

研究の動機と目的

• 表面検出器アレイの広大な空気シャワーイベントに対する方位角分解能を、駅の多重度とシャワー幾何学的構造に基づいて特定すること。
• 主粒子エネルギー再構成の主要な観測量であるS(1000)エネルギー推定器における統計的および系統的不確実性を定量化すること。
• 完全なモンテカルロシミュレーションと冗長な駅ペア測定を用いて、誤差推定フレームワークの妥当性を検証すること。
• 飽和駅とシャワー発達のゆらぎがS(1000)再構成精度に与える影響を評価すること。

提案手法

• 方位角分解能は、到達方向（θおよびφ）の分散を用いてイベントごとに計算され、空間角分解能関数F(η)はF(η) = 1/2 (V[θ] + sin²(θ)V[φ])として定義され、AR = 1.5√F(η)と表される。
• 時間フィッティング手順では、シャワーの波面を放物線的モデルで近似し、時間測定誤差の影響を最小限に抑える。
• S(1000)は、変動する勾配βを有する修正NKG横方向分布関数を用いて、シャワー核からの1000 m地点での信号サイズを推定する。
• S(1000)における系統的誤差は、パラメータ化されたβ(S(1000))を中心とするガウス分布からβをサンプリングすることで評価され、σβ(S(1000)) = 0.71 × exp(−0.976 × log(S(1000)))と表される。
• 統計的誤差は、隣接する駅ペアを用いて測定された信号サイズの不確実性から導出される。
• 完全なモンテカルロシミュレーションでは、Corsikaで生成されたシャワーを用い、再構成されたS(1000)値と真のS(1000)値を比較することで、S(1000)における総誤差をベンチマークする。

実験結果

リサーチクエスチョン

• RQ1 表面検出器アレイの方位角分解能は、駅の多重度と天頂角の関数としてどのように変化するか？
• RQ2 S(1000)における統計的および系統的不確実性は、信号サイズおよび駅の飽和状態にどのように依存するか？
• RQ3 高エネルギー領域において、S(1000)の不確実性の主な寄与要因がシャワー発達のゆらぎである割合はどの程度か？
• RQ4 S(1000)の経験的誤差推定フレームワークは、完全なモンテカルロシミュレーションとどの程度一致するか？
• RQ5 冗長な駅ペアの使用は、方位角分解能およびS(1000)誤差推定値の一貫性をどのように裏付けるか？

主な発見

• 3駅イベント（E < 4 EeV）では方位角分解能が2°未塔であり、4〜5駅イベント（3 < E < 10 EeV）では1.2°未塔に向上し、6駅以上（E > 10 EeV）では0.9°未塔にまで改善される。
• 飽和駅を含むイベントではS(1000)の総誤差が8%未塔、非飽和イベントでは12%未塔であり、S(1000) < 40 VEMでは統計的誤差が主な寄与要因である。
• S(1000)における系統的誤差は主に横方向分布関数の勾配βのゆらぎに起因し、σβ(S(1000)) = 0.71 × exp(−0.976 × log(S(1000)))としてパラメータライズドされている。
• 信号の非線形性のため、飽和駅を含むイベントではS(1000)の誤差予算が系統的に5%高めに推定される。
• 完全なモンテカルロシミュレーションにより、S(1000)の総誤差推定が確認され、データベースの誤差推定と非常に良好に一致しており、高エネルギー飽和イベントではわずかに4%の過大評価が見られる。
• S(1000)誤差に天頂角依存性は観測されず、S(1000) < 10 VEMの領域でのみわずかなS(1000)再構成バイアスが確認される。