[論文レビュー] Towards a common analysis framework for gamma-ray astronomy

J. Knödlseder, M. Mayer+5 more|arXiv (Cornell University)|Jul 21, 2013

Astrophysics and Cosmic Phenomena被引用数 4

ひとこと要約

本論文は、高エネルギー伽マ線データ解析のためのオープンソースで機器に依存しないフレームワーク、GammaLib および ctools を提示する。このフレームワークにより、複数の望遠鏡を統合してスペクトル的・形態的解析が可能となり、1 MeV から 10 TeV のエネルギー範囲で、Fermi-LAT、H.E.S.S.、COMPTEL のデータを用いてカブールネビュラのマルチインストルメント同時フィッティングを実施。7つのエネルギーデケードにわたり一貫性のあるスペクトルモデリングを達成した。

ABSTRACT

Thanks to the success of current gamma-ray telescopes (Fermi, H.E.S.S., MAGIC, VERITAS), and in view of the prospects of planned observatories such as the Cherenkov Telescope Array (CTA) or the High-Altitude Water Cherenkov Observatory (HAWC), gamma-ray astronomy is becoming an integral part of modern astrophysical research. Analysis today relies on a large diversity of tools and software frameworks that were specifically and independently developed for each instrument. With the aim of unifying the analysis of gamma-ray data, we are currently developing GammaLib (http://sourceforge.net/projects/gammalib), a C++ library interfaced to Python that provides a framework for an instrument independent analysis of gamma-ray data. On top of GammaLib we have created ctools (http://cta.irap.omp.eu/ctools), a set of analysis executables that is being developed as one of the prototypes for the CTA high-level science analysis framework, but which is equally suited for the analysis of gamma-ray data from the existing Fermi-LAT telescope and current Cherenkov telescope arrays. In particular, ctools and GammaLib provide the novel opportunity of a simultaneous multi-instrument analysis. We present the status of the software development, and illustrate its capabilities with a spectral analysis of the Crab nebula emission over seven decades in energy (1 MeV to 10 TeV) using multi-instrument (COMPTEL, Fermi-LAT, H.E.S.S.) gamma-ray observations as well as a simulation of a CTA observation of the supernova remnant RX J1713.7-3946.

研究の動機と目的

• ガンマ線天文学において、現在は特許権を有する機器固有のツールに依存しているが、標準的で機器に依存しない分析フレームワークの欠如に対処する。
• 人工衛星搭載（例：Fermi-LAT）および地上望遠鏡（例：H.E.S.S.、MAGIC）のガンマ線望遠鏡から得られるデータの分析を、1つの相互運用可能なソフトウェアスタックに統合する。
• パrametricスペクトルモデル（例：対数放物線、累乗則）を用いて、複数のインストルメントのデータセットを同時に統合的にモデリングすることで、パラメータの制約と不確実性の評価を向上させる。
• FITS および IRAF 標準に準拠した公開型でコミュニティ主導の科学ツールを提供することで、近い将来に稼働予定の Cherenkov Telescope Array (CTA) をオープン観測所として支援する。
• 共通のソフトウェア基盤を提供することで、今後の高エネルギー天体物理学研究におけるマルチインストルメントの検証と共同モデリングを促進する。

提案手法

• イベントデータ、機器応答関数、源および背景モデル、天球マップ、FITS/XML I/O を機器に依存しない抽象化で実装した C++ ライブラリである GammaLib を開発する。
• ctobssim、ctselect、ctbin、ctlike などの FTOOLS を模したモジュラーでコマンドライン指向のツール群（例：ctobssim、ctselect、ctbin、ctlike）として ctools を実装し、イベントのシミュレーション、選択、ビニング、最尤フィッティングを実行する。
• レイヤードソフトウェアアーキテクチャを採用：上位レイヤーは機器に依存しない分析、2番目のレイヤーは機器固有のインターフェース（COMPTEL、Fermi-LAT、CTA）、3番目のレイヤーは数値計算および最適化サービス、4番目のレイヤーは FITS/XML および IRAF パrameter インターフェースのサポート。
• ctlike を用いたビンありおよびビンなしの最尤フィッティングにより、パラメトリックスペクトルモデル（例：対数放物線、累乗則）を用いて複数インストルメントのデータセットを同時にモデリングする。
• 超新星残光のような拡大源を、位置、半径、厚さ、スペクトル指数を自由パラメータとするシェル型空間モデルでモデリングする。
• カブールネビュラおよび RX J1713.7−3946 をテストケースとして、実データ（Fermi-LAT、H.E.S.S.）およびシミュレーテッド CTA 観測データを用いてフレームワークの妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 1つのオープンソースソフトウェアフレームワークが、人工衛星望遠鏡および地上望遠鏡の両方をカバーするガンマ線データの一貫性があり機器に依存しない分析を可能にするか？
• RQ2 1台のインストルメントでの分析と比較して、マルチインストルメント同時スペクトルフィッティングが、パラメータの制約と不確実性推定をどの程度改善できるか？
• RQ3 本フレームワークは、現実的な空間的およびスペクトル的モデルを用いて、超新星残光のような拡大ガンマ線源を正確にモデリングできるか？
• RQ4 特に拡大源の研究において、非常に高エネルギーのガンマ線天文学における複雑な粒子背景を、本フレームワークはどの程度適切に処理できるか？
• RQ5 本フレームワークは、近い将来に稼働予定の CTA 観測所のオープンデータポリシーを踏まえ、実用的でコミュニティ主導の公開科学ツールとして機能できるか？

主な発見

• 本フレームワークは、Fermi-LAT、H.E.S.S.、COMPTEL のデータを用いて、1 MeV から 10 TeV の7つのエネルギーデケードにわたるカブールネビュラの統合的スペクトルフィッティングに成功。広いエネルギー範囲と一貫性のあるモデリングを実証した。
• 対数放物線モデルを用いた統合的フィッティングにより、Fermi-LAT および H.E.S.S. の個別データにのみフィットさせた累乗則フィッティングよりも、より一貫性のあるスペクトル記述が得られ、パラメータの共分散と不確実性評価が向上した。
• RX J1713.7−3946 の5時間にわたる CTA シミュレーション観測を、ctlike を用いた非ビンあり最尤フィッティングで再現。入力されたシェル型の形態とスペクトルパラメータが正確に回復された。
• RX J1713.7−3946 のフィットモデルは、有効面積の変動に起因する非対称的特徴を含めて、シミュレーテッドデータと優れた一致を示し、空間モデリングの能力を裏付けた。
• 本フレームワークは、標準的なデータフォーマット（FITS）および分析パラダイム（例：IRAF パrameter インターフェース）と互換性を有しており、既存の高エネルギー分野のソフトウェアエコシステムとの相互運用性を確保した。
• GammaLib および ctools のオープンソース性により、コミュニティ主導の開発と相互検証が可能となり、実際の H.E.S.S. および MAGIC データを用いた粒子背景モデリングの精錬のための活発なテストが継続されている。