[論文レビュー] A Monte Carlo template-based analysis for very high definition imaging atmospheric Cherenkov telescopes as applied to the VERITAS telescope array

Stéphane Vincent|arXiv (Cornell University)|Sep 7, 2015

Astrophysics and Cosmic Phenomena参考文献 4被引用数 11

ひとこと要約

本論文は、VERITAS超高エネルギーγ線望遠鏡アレイ向けに、モンテカルロテンプレートに基づく尤度再構成アルゴリズムを提示する。観測されたカメラ画像を、ピクセル単位でシミュレートされたテンプレートと比較することで、エネルギーおよび方向再構成の精度を向上させる。この手法は、従来のモーメントベース手法に比べ、特に高い仰角および低エネルギー領域で優れた角分解能およびエネルギー分解能を達成し、モンテカルロシミュレーションを用いてその性能が実証されている。

ABSTRACT

We present a sophisticated likelihood reconstruction algorithm for shower-image analysis of imaging Cherenkov telescopes. The reconstruction algorithm is based on the comparison of the camera pixel amplitudes with the predictions from a Monte Carlo based model. Shower parameters are determined by a maximisation of a likelihood function. Maximisation of the likelihood as a function of shower fit parameters is performed using a numerical non-linear optimisation technique. A related reconstruction technique has already been developed by the CAT and the H.E.S.S. experiments, and provides a more precise direction and energy reconstruction of the photon induced shower compared to the second moment of the camera image analysis. Examples are shown of the performance of the analysis on simulated gamma-ray data from the VERITAS array.

研究の動機と目的

• イメージング大気チェレンコフ望遠鏡（VERITASなど）のより高精度な再構成手法を開発することを目的とし、詳細なモンテカルロシミュレーションを活用する。
• 従来の2次モーメント画像解析に課題があること、特に高仰角および低エネルギー領域で性能が低下することを克服することを目的とする。
• 観測データとシミュレートされたシャワー・テンプレートとのピクセルレベルでの尤度比較を用いて、角分解能およびエネルギー分解能を向上させることを目的とする。
• 尤度関数の数値最適化により、主粒子パラメータ（エネルギー、方向、到達距離）の堅牢な再構成を可能とすることを目的とする。
• テンプレート適用範囲を拡張することで、鉄核などの非ガンマ粒子の再構成の可能性を検討することを目的とする。

提案手法

• 本手法は、エネルギー、仰角、到達距離、衝突深さの異なる条件で実行されたCORSIKAを用いた電磁気的空気シャワーのモンテカルロシミュレーションから、多次元テンプレートライブラリを生成する。
• 大気密度、吸収、検出器応答（例：量子効率、ミラー集光面積）を組み込んだレイトレーシングを用いて、チェレンコフ光子の経路を大気および望遠鏡光学系を通って追跡する。
• 各テンプレートについて、カメラ内のピクセルの強度をモデル化し、0.15°ピクセル分解能とデイビス・コットン光学系を有する焦点面に画像を投影する。
• 観測データと予測されたテンプレート強度とのピクセル単位比較に基づき、ポisson尤度モデルを用いて尤度関数を定義する。
• Levenberg–Marquardtアルゴリズムを用いて尤度関数を最大化し、パrameterまたはコスト関数の安定性の閾値によって収束を判定する。
• 局所的最小値のリスクを低減し、初期パrameter値への感受性を軽減するために、微分進化（differential evolution）を用いたグローバル最適化スキームを採用する。

実験結果

リサーチクエスチョン

• RQ1 イメージング大気チェレンコフ望遠鏡において、従来の2次モーメント解析に比べ、テンプレートベース尤度再構成手法がエネルギー分解能および角分解能を向上させられるか。
• RQ2 エネルギー、仰角、到達パラメータに応じて、テンプレート手法の性能はどのように変化するか。特に、モーメント解析が劣化する領域（高仰角・低エネルギー）での性能は。
• RQ3 低エネルギーおよび高仰角領域において、エネルギーバイアスをどの程度低減でき、エネルギー分解能を改善できるか。
• RQ4 同じテンプレートフレームワークを用いて、鉄核などの非ガンマ粒子の再構成が一般化可能か。
• RQ5 グローバル最適化（微分進化）と局所最小化（Levenberg–Marquardt）の組み合わせが、局所的最小値への収束をどの程度効果的に回避できるか。

主な発見

• テンプレートベース手法は、画像の伸びが顕著になる高仰角領域で、モーメントベース解析に比べて顕著に優れた角分解能を達成する。
• エネルギー分解能は全エネルギー範囲で向上し、エネルギーバイアス分布の68%含有範囲が一貫した性能を示すが、低エネルギーおよび高仰角領域でバイアスが増加する傾向を示す。
• 低仰角領域では、全エネルギー範囲でエネルギーバイアスが低減されるが、最大バイアスは低エネルギーおよび大規模な仰角領域で発生する。これはエネルギー依存補正の必要性を示唆する。
• 尤度ベース手法は、2次モーメントのみに依存するモーメント解析に比べ、画像全体の形態を活用しているため、シャワーのパラメータ再構成において優れた性能を発揮する。
• 微分進化によるグローバル最小化ステップの導入により、収束のロバスト性が向上し、初期パラメータ値への感受性が低下する。
• テンプレートライブラリに強子的シャワーのモデルを組み込むことで、鉄核などの他の粒子種の再構成にも応用可能である。