[论文解读] HESS II Data Analysis with ImPACT

R. D. Parsons, M. Gajdus+1 more|arXiv (Cornell University)|Sep 17, 2015

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

一句话总结

本文提出了一种针对H.E.S.S. II阵列的ImPACT事件重建算法的改进版本，实现了使用混合望远镜设计的超高能伽马射线 shower 的高性能立体与单台重建。该方法在子100 GeV能量范围内实现了改进的角分辨率与能量分辨率，达到了亚100 GeV灵敏度，通过96σ的蟹状星云探测（最低至100 GeV）以及与先前结果一致的对数抛物谱拟合得到验证。

ABSTRACT

The High Energy Stereoscopic System (H.E.S.S.) very high energy gamma-ray telescope array has added a fifth telescope of 600 m$^{2}$ mirror area to the centre of the 4 existing telescopes, lowering its energy threshold to the sub-100 GeV range and becoming the first operational IACT array using multiple telescope designs. In order to properly access this low-energy range however, some adaptation must be made to the existing event analysis. We present an adaptation of the high-performance event reconstruction algorithm, Image Pixelwise fit for Atmospheric Cherenkov Telescopes (ImPACT), for performing mono and stereo event reconstruction with the H.E.S.S. II array. The reconstruction algorithm is based around the likelihood fitting of camera pixel amplitudes to an expected image template, directly generated from Monte Carlo simulations. This advanced reconstruction is combined with a multi variate analysis based background rejection scheme to provide a sensitive and stable analysis scheme in the sub-100 GeV gamma-ray energy range. We will present the latest results of the ImPACT analysis on both simulated and real H.E.S.S. II data, demonstrating the behaviour of the ImPACT analysis at the lowest energies.

研究动机与目标

• 为了实现在子100 GeV能量范围内的H.E.S.S. II数据的敏感且稳定的分析，这对探测低能伽马射线源至关重要。
• 为了将高性能的ImPACT重建算法适配于新型600 m²中心望远镜（CT5），其可将能量阈值降低至50 GeV以下。
• 通过量身定制的重建与背景抑制技术，解决仅在CT5中检测到的单台事件（mono events）和在CT5及至少另一台望远镜中检测到的立体事件（stereo events）的挑战。
• 通过模拟数据和蟹状星云的真实观测，验证新的分析链，确认其在低能区的性能与灵敏度。

提出的方法

• ImPACT算法对相机像素幅值进行最大似然拟合，匹配基于蒙特卡洛模拟生成的图像模板，模板分别针对CT5和原始H.E.S.S. I望远镜预先计算。
• 对于立体事件，使用MINUIT进行多维拟合，优化源位置、撞击点、能量和簇射最大深度，模板基于固定初级粒子参数的库进行插值。
• 对于单台事件，采用两步拟合程序：首先进行自由拟合，Xmax固定在预期值；随后固定源位置以解决信息有限导致的退化问题。
• 背景抑制通过多变量分析（MVA）和基于Hillas参数训练的神经网络实现，为单台事件优化了不同等级的截断阈值（安全、标准、宽松）。
• 分析基于像素幅值与模板比较的似然函数，即使在低信号条件下也能实现精确重建。
• 应用了对均值缩放宽度/长度的盒形截断作为初步背景抑制步骤，未来计划通过高级MVA技术进一步优化。

实验结果

研究问题

• RQ1 ImPACT算法能否成功适配于H.E.S.S. II阵列，该阵列使用两种不同的望远镜设计（CT5与H.E.S.S. I望远镜）？
• RQ2 在100 GeV以下能量范围内，单台事件与立体事件重建在角分辨率与能量分辨率方面的表现如何比较？
• RQ3 与以往基于H.E.S.S. I的分析方法相比，新分析链在灵敏度与能量阈值方面提升了多少？
• RQ4 ImPACT重建在真实数据上的表现如何，特别是在强而著名的源（如蟹状星云）上的表现？

主要发现

• ImPACT分析利用7.2小时的H.E.S.S. II数据，实现了对蟹状星云的96σ显著性探测，产生超过7,700个多余伽马样事例，信号与背景比为2.3。
• 该分析成功探测到100 GeV的能量范围，相比此前300–400 GeV的能量阈值实现了显著提升。
• 使用对数抛物模型对蟹状星云数据进行谱拟合，得到最优参数：φ₀ = 3.951 ± 0.09 × 10⁻¹¹ cm⁻² s⁻¹ TeV⁻¹，α = 2.49 ± 0.03，β = 0.203 ± 0.02，E₀ = 1 TeV。
• 立体重建在全能量范围内（尤其在300 GeV以上）相比单台重建提供了显著更优的角分辨率与能量分辨率，得益于多台望远镜图像的使用。
• 单台事件重建虽然约束较少，但实现了稳定的性能，每事件拟合时间仅0.2秒，且通过神经网络截断实现了有效的背景抑制。
• 由于引入了4台H.E.S.S. I望远镜，立体分析的有效面积在300 GeV以上显著提升，证实了多望远镜重建的优势。