[論文レビュー] Photon Reconstruction for H.E.S.S. Using a Semi-Analytical Shower Model

M. Holler, A. Balzer+3 more|arXiv (Cornell University)|Sep 9, 2015

Astrophysics and Cosmic Phenomena参考文献 3被引用数 19

ひとこと要約

本論文では、H.E.S.S. IIの5台の望遠鏡アレイにおけるガンマ線再構築のための半解析的シャワー・モデルの適応が提示され、モノスコピックおよびステレオスコピック解析を可能にし、エネルギー閾値と感度が向上した。コンビネッド解析モードは、モノスコピックおよびステレオスコピック再構築を統合することで、エネルギー範囲全体にわたる最も広いエネルギーカバレッジを達成し、効果的断面積および感度の面で個々のモードを上回った。

ABSTRACT

The High Energy Stereoscopic System (H.E.S.S.) is an array of five Imaging Atmospheric Cherenkov Telescopes (IACTs) designed to detect cosmogenic gamma-rays with very high energies. Originally consisting of just four identical IACTs (CT1-4) with an effective mirror diameter of 12$\,$m each, it was expanded with a fifth IACT (CT5) with a mirror diameter of 28$\,$m in 2012. Being the largest IACT worldwide, CT5 allows to lower the energy threshold of H.E.S.S., making the array sensitive at energies where space-based detectors run out of statistics. Events can be analysed either monoscopically (i.e. using only information of CT5) or stereoscopically (requiring at least two triggered telescopes per event). To achieve a good performance, a sophisticated event reconstruction and analysis framework is indispensable. This is particularly important for H.E.S.S. since it is now the first IACT array that consists of different telescope types. An advanced reconstruction method is based on a semi-analytical model of electromagnetic particle showers in the atmosphere (model analysis). The properties of the primary particle are reconstructed by comparing the image recorded by each triggered telescope with the Cherenkov emission from the shower model using a log-likelihood maximisation. Due to its performance, this method has become one of the standard analysis techniques applied to CT1-4 data. Now it has been modified for use with the five-telescope array. We present the adapted model analysis and its performance in both monoscopic and stereoscopic analysis mode.

研究の動機と目的

• H.E.S.S. Iの4台の12 m望遠鏡で既に使用されていた半解析的モデル解析を、新しい5台の望遠鏡からなるH.E.S.S. IIアレイに拡張すること。
• 異なる望遠鏡タイプおよびサイズを有する非均質なアレイにおけるガンマ線シャワー再構築の課題に対処し、統一された再構築フレームワークを要求すること。
• モノスコピック（CT5のみ）およびステレオスコピック（複数台の望遠鏡）再構築を活用する新しいコンビネッド解析モードの開発および検証を行い、感度およびエネルギーカバレッジを最大化すること。
• モンテカルロシミュレーションおよびカブスネビュラのデータを用いて、モノスコピック、ステレオスコピック、およびコンビネッド解析モードの性能を評価し、効果的断面積および微分感度に注目すること。

提案手法

• 半解析的電磁シャワー・モデルは、一次エネルギー（E）、シャワー深さ（T）、およびインパクト距離（R）の関数として大気中のチェレンコフ光分布を計算し、さまざまなパラメータに対して事前にテンプレートを計算している。
• 各望遠鏡のカメラ画像に対して、測定されたピクセル強度（s）とモデル予測値（μ）を比較することで対数尤度関数を計算する。この際、ピエドストルノイズ（σp）、光電子ピーク幅（σγ）、およびキャリブレーション不確実性（σc）を考慮し、尤度式：P(s|μ,σp,σγ,σc) = ∑n [μ^n * e^(-μ) / (n! * √(2π(σp² + nσγ² + n²σc²))) * exp( - (s - n)² / (2(σp² + nσγ² + n²σc²)) )] を用いる。
• 対数尤度をLevenberg–Marquardtアルゴリズムを用いて最大化し、一次粒子パラメータ（方向（Dir）、エネルギー（E）、シャワー深さ（T）、インパクト距離（R））を再構築する。
• イベント単位の不確実性は、対数尤度の2階微分から推定され、例えばΔDir = 1 / √(∂²lnL/∂Dir²) のように表される。
• 3つの解析モードを実装した：モノ（CT5のみ）、ステレオ（CT1–CT5）、およびコンビネッド（各イベントについて方向不確実性が低い方の結果を選択）。
• 性能評価は、点源を18°の天の赤道角および0.5°のウッブルオフセットで設定したモンテカルロシミュレーションと、実データからのバックグラウンドイベントを用いて行い、現実的な条件を確保した。

実験結果

リサーチクエスチョン

• RQ1 H.E.S.S. IIの非均質な望遠鏡構成（特に28 mのCT5望遠鏡を含む）に対応するため、半解析的モデル解析をどのように適応できるか？
• RQ2 効果的断面積および微分感度の観点から、モノスコピック、ステレオスコピック、およびコンビネッド解析モードの相対的性能はいかほどか？
• RQ3 モノスコピックおよびステレオスコピック再構築を統合するコンビネッド解析モードは、単独のモードよりも優れた全体的なエネルギーカバレッジを達成できるか？
• RQ4 新しいモデル解析のエネルギー閾値および感度は、元のH.E.S.S. I（CT1–CT4）ステレオスコピックモードと比較してどの程度か？
• RQ5 特にθ²カットが、コンビネッドモードにおけるバックグラウンド抑制および感度に与える影響はいかほどか？

主な発見

• コンビネッド解析モードは、モノスコピック再構築（CT5のみ）の低エネルギー閾値と、高エネルギー域におけるステレオスコピック再構築の高い検出効率を統合することで、最も広いエネルギーカバレッジを達成した。
• 300 GeV未満の低エネルギー域では、コンビネッドモードの効果的断面積がモノスコピックおよびステレオスコピックモードを上回り、図4に示すように高エネルギー域でも高い性能を維持した。
• 50時間の観測において、コンビネッドモードはモノスコピックモードと同等以上の微分感度を達成しており、中・高エネルギー域ではそれを上回ったが、現在のところ最適でないカットのため、ステレオスコピックモードに劣っている。
• ステレオスコピックモードは最高の信号対バックグラウンド（S/B）比を達成したが、コンビネッドモードは最高の過剰率および有意水準を示し、全体的な感度が優れていることを示した。
• モデル解析は、大きなθ²値に対してもシミュレーションと良好な一致を示しており、源検出およびバックグラウンド抑制の信頼性を裏付けた。
• 本手法は、異なるサイズの望遠鏡を有するH.E.S.S. IIアレイにおいても高感度再構築を可能にし、初めてのIACTシステムとして、複雑な機器構成への強靭性および適応性を示した。