[論文レビュー] Gammapy - A prototype for the CTA science tools

Christoph Deil, R. Zanin+29 more|arXiv (Cornell University)|Jan 1, 2017

Gamma-ray bursts and supernovae参考文献 9被引用数 11

ひとこと要約

Gammapy は、Numpy、Scipy、Astropy を基盤とするオープンソースの Python パッケージであり、高エネルギーガンマ線データ解析のためのもので、Cherenkov Telescope Array (CTA) の科学ツールのプロトタイプを提供する。イベントリストと機器応答関数（IRF）からスカイ画像、スペクトル、光度曲線の生成を可能にし、H.E.S.S.、Fermi-LAT、およびシミュレートされた CTA データに対する妥当性検証がなされた。CTA のオープンかつ協働的な科学フレームワークに適した状態にあることが示された。

ABSTRACT

Gammapy is a Python package for high-level gamma-ray data analysis built on Numpy, Scipy and Astropy. It enables us to analyze gamma-ray data and to create sky images, spectra and lightcurves, from event lists and instrument response information, and to determine the position, morphology and spectra of gamma-ray sources. So far Gammapy has mostly been used to analyze data from H.E.S.S. and Fermi-LAT, and is now being used for the simulation and analysis of observations from the Cherenkov Telescope Array (CTA). We have proposed Gammapy as a prototype for the CTA science tools. This contribution gives an overview of the Gammapy package and project and shows an analysis application example with simulated CTA data.

研究の動機と目的

• 将来の Cherenkov Telescope Array (CTA) 観測所に対応する、柔軟でオープンソースのソフトウェアフレームワークを構築すること。
• 科学的 Python を基盤とし、Astropy や Numpy などの既存の天文学ソフトウェアスタックを活用することで、共通でコミュニティ主導の解析パイプラインを確立すること。
• 現在および将来のガンマ線実験間で、データフォーマットと解析手法の相互運用性と標準化を実現すること。
• CTA に適したスケーラブルで拡張可能かつ保守可能な科学ツールスタックのプロトタイプを提供し、1次元および3次元の解析手法をサポートすること。
• Astropy などの広く採用されているツールを基盤とし、アクセスしやすく、十分にドキュメント化されたソフトウェアを提供することで、ガンマ線天文学コミュニティにおける協働を促進すること。

提案手法

• 高性能で読みやすく保守可能なコードを実現するため、科学的 Python スタック（Numpy、Scipy、Astropy）を基盤として利用する。
• イベントリストと機器応答関数（IRF）には FITS 形式を採用し、スカイ画像およびスペクトル生成のためのエネルギーおよび空間のビンングをサポートする。
• リングモデルやデータ駆動型テンプレートを用いたバックグラウンド推定を実装し、イベントエネルギーおよび方向情報から露出マップを計算する。
• 2次元および3次元（キューブ）解析手法をサポートし、後者により、源の放射の空間的およびスペクトル的同時モデリングが可能になる。
• モデリングおよびフィッティングには Sherpa を、可視化には Matplotlib を、空間および天球座標演算には regions/reproject/healpy をオプションで統合する。
• EventList、SkyImage、Background2D、SourceSpectralModel などのクラスを備えたモジュラーでオブジェクト指向の設計を採用し、再利用可能で組み合わせ可能な解析ワークフローを実現する。

実験結果

リサーチクエスチョン

• RQ1 Astropy や Numpy を基盤とする Python ベースのオープンソースパッケージが、CTA 科学ツールの実用的で拡張可能なプロトタイプとして成立するか。
• RQ2 Gammapy は、源検出、形態フィッティング、スペクトル抽出、光度曲線生成といった高レベルなガンマ線解析タスクをどれほど効果的にサポートするか。
• RQ3 Gammapy は、現在の IACT（例：H.E.S.S.、MAGIC）および Fermi-LAT からの既存のデータフォーマットやパイプラインと、どの程度統合可能か。
• RQ4 シミュレートされた CTA データに対して、Gammapy の3次元（キューブ）解析アプローチは、従来の1次元スペクトル解析と比較してどれほど効果的か。
• RQ5 Gammapy のアーキテクチャおよび開発モデルは、CTA の30年間の運用期間にわたり、長期的な進化とコミュニティの貢献を支えることができるか。

主な発見

• Gammapy は、イベントリストと機器応答関数（IRF）を用いて、スカイ画像、スペクトル、光度曲線の生成を含む高レベルなガンマ線データ解析を成功裏に実現した。
• 本パッケージは、H.E.S.S. や Fermi-LAT の実データに対して妥当性検証が行われ、現在は CTA データのシミュレーションおよび解析に活用されている。
• Gammapy は 2次元および 3次元解析手法を両方サポートしており、3次元キューブ解析は複雑な源領域における源の特徴抽出の観点で有望であることが示された。
• Astropy および Numpy の使用により、大規模なイベントリストやスカイ画像における効率的な計算が可能であり、Numpy および Astropy の C 拡張により性能が向上している。
• Gammapy は PyPI、conda、多数の Linux パッケージマネージャー経由で配布されており、科学コミュニティにおける広範なアクセス性とインストールの容易さを実現している。
• 約 30 名のガンマ線天文学者が貢献しており、コミュニティによる採用が進んでおり、協働開発の可能性が明確に示された。