[論文レビュー] Instrument simulation for the analysis of cosmic ray electron with the Fermi LAT

C. Sgrò, J. Bregeon+1 more|arXiv (Cornell University)|Jul 2, 2009

Dark Matter and Cosmic Phenomena参考文献 5被引用数 3

ひとこと要約

本論文では、宇宙線電子（CRE）解析のためのフェルミLAT機器応答をモデル化するために、GEANT4を用いた詳細なモンテカルロシミュレーションを提示している。これにより、正確なイベント選別と高精度なスペクトル再構成が可能になった。シミュレーションは20 GeVまでで10%未満のエネルギー分解能を達成し、系統的不確実性を制御しており、20 GeVから1 TeVまでの範囲で、これまでで最も精度の高いCREスペクトル測定を実現した。

ABSTRACT

The Fermi LAT collaboration has built up a detailed Monte Carlo simulation to characterize the instrument response and tune its performance. The simulation code is built around the widely used GEANT4 toolkit and was carefully validated against beam test and flight data. This poster shows how the full LAT simulation is used to develop the event selection for the Cosmic-Ray Electron (\emph{CRE}) analysis so as to optimize the instrument performance. In particular, we will show in detail the determination of the geometry factor and the residual hadron contamination. The very accurate MC simulation proved to be fundamental to control the systematic uncertainties on the CRE spectrum measured by the Fermi LAT.

研究の動機と目的

• フェルミ大型宇宙望遠鏡（LAT）を用いた宇宙線電子（CRE）スペクトル測定のための高忠実度の機器シミュレーションを開発すること。
• 電磁シャワーのトポロジーをモデル化し、バリオン的バックグラウンドと区別することで、CREのイベント選別基準を最適化すること。
• スペクトルのデコンボリューションに用いるため、幾何因子およびエネルギー分解能を含む機器応答関数を評価・定量化すること。
• ビームテストおよび軌道上データとの比較を通じて、シミュレーションを検証し、CREスペクトルにおける系統的不確実性を最小限に抑えること。

提案手法

• シミュレーションは、GEANT4ツールキットを用いて、すべてのLAT検出器部品（活性および非活性材料を含む）における粒子の伝播、エネルギー損失、多次散乱、エネルギー損失をモデル化している。
• デジタイゼーションアルゴリズムは、ゲインのばらつき、ノイズ、非線形性、センサのずれといった現実的な検出器効果を組み込み、実データ処理を模倣している。
• 軌道上での粒子環境は、宇宙線、地球アルベドガンマ線、二次粒子を含む詳細なモデルでシミュレートされており、地磁気カットオフおよび東西効果も含む。
• 機器応答関数の独立評価のため、エネルギースペクトルがべき乗則（E⁻¹）に従う等方的電子フラックスをシミュレートしている。
• イベント選別は、MCデータでトレーニングされた分類ツリーに基づき、シャワーのトポロジーと応答特徴を用いて、電子とバリオンを区別する。
• 機器応答関数（幾何因子およびエネルギー分解能）はMCシミュレーションから導出され、測定されたイベントレートを真のCREスペクトルにデコンボリューションするために使用されている。

実験結果

リサーチクエスチョン

• RQ1 フェルミLATは、GEANT4を用いて、20 GeVから1 TeVのエネルギー範囲にわたり、宇宙線電子の応答をどの程度正確にシミュレートできるか？
• RQ2 LATにおける電子誘発シャワーとバリオン的シャワーを区別するための、主な検出器応答特徴は何か？
• RQ3 シミュレーションに基づくイベント選別は、CREサンプルにおける残存バリオン汚染をどの程度低減するか？
• RQ4 MCシミュレーションは、再構成されたCREスペクトルにおける系統的不確実性をどの程度制御できるか？
• RQ5 真の電子エネルギーの関数として、LATのエネルギー分解能および幾何因子は、シミュレーションからどのように導出されるか？

主な発見

• ヒット多重度分布の比較により、100 GeV/cおよび200 GeV/cの電子ビームに対するビームテストデータと、モンテカルロシミュレーションの間に良好な一致が得られている。
• 30 GeV以上で幾何因子が急激に増加しているのは、搭載されたイベント選別が20 GeV未満のエネルギー損失を持つイベントを除外するためであり、これはバックグラウンドの排除を強化している。
• エネルギー分解能は20 GeVで10%未満であり、1 TeVまで30%未満で維持されており、再構成エネルギーと真のエネルギーの68%含有窓はほぼ1の周辺に集中している。
• 再構成エネルギーの最頻値を真のエネルギーで割った値は、全エネルギー範囲にわたり1と整合しており、1 TeVまでにかけても偏りのない応答であることが確認された。
• ビームおよび軌道上データとの比較による検証を通じて、シミュレーション最適化済みのイベント選別により、残存バリオン汚染が最小限に抑えられている。
• 高忠実度のMCモデルにより、20 GeVから1 TeVの範囲で、宇宙線電子スペクトルの最も精度の高い測定が可能になった。