[Paper Review] The Disp Method for Analysing Large Zenith Angle Gamma-Ray Data

G. D. Şentürk|arXiv (Cornell University)|Sep 27, 2011

Astrophysics and Cosmic Phenomena3 citations

TL;DR

This paper presents a multi-dimensional Disp method for improving angular resolution in very-high-energy gamma-ray astronomy using ground-based Cherenkov telescopes, particularly for large-zenith-angle (LZA) observations. By using separate width and length parameters in multidimensional lookup tables from Monte Carlo simulations, the method reduces reconstruction errors and improves angular resolution and signal-to-noise ratio, especially for sources like the Crab Nebula and Sgr A* at high zenith angles.

ABSTRACT

The Disp method is an algorithm that is used for reconstruction of primary gamma ray direction in ground- based atmospheric Cherenkov telescope experiments -measuring very-high-energy (VHE) gamma rays in the energy range between 100GeV and 30 TeV. In general terms, the geometric information obtained from one single shower image is sufficient for the algorithm to find the sky location of the primary. Various versions of the Disp method were implemented and used in the past. In this study, we present a multi-dimensional implementation of the Disp method for the VERITAS instrument and show (using Monte Carlo simulations and the Crab Nebula observations) that it significantly improves the angular resolution for large-zenith-angle (LZA) observations. We also applied the disp method to VERITAS data taken from the galactic center region which is detected by VERITAS.

Motivation & Objective

• To address the degradation in angular resolution for large-zenith-angle (LZA) observations in ground-based Cherenkov telescopes due to projection effects and near-parallel shower image axes.
• To improve direction reconstruction accuracy for VHE gamma rays in the 100 GeV to 30 TeV energy range using a novel implementation of the Disp method.
• To enable high-sensitivity detection of faint and extended sources, such as the Galactic center, by enhancing angular resolution and signal-to-noise ratio in LZA data.
• To validate the method using Monte Carlo simulations and real observations of the Crab Nebula and Sgr A*.

Proposed method

• The method uses multidimensional lookup tables derived from Monte Carlo simulated air shower data, storing Hillas parameters (size, width, length, disp) as functions of zenith angle, azimuth angle, telescope ID, and noise level.
• For each observed shower image, the algorithm calculates size, width, and length from the image, then queries the lookup table to retrieve the corresponding disp value, representing the angular offset from the image centroid along the major axis.
• To resolve the head-tail ambiguity, the algorithm selects the closest cluster of points from each ellipse's two possible directions, estimating the true arrival direction per telescope.
• The final source direction is computed as a weighted average of individual telescope reconstructions, improving robustness and accuracy.
• The key innovation lies in treating width and length as separate dimensions rather than combining them into a width/length ratio, which enhances resolution for LZA events.
• The method is applied to VERITAS data, including LZA Crab Nebula observations and Galactic center observations at ~62° zenith angle.

Experimental results

Research questions

• RQ1 Can the Disp method significantly improve angular resolution for large-zenith-angle (LZA) gamma-ray observations compared to the default geometric reconstruction method?
• RQ2 How does the separation of width and length parameters in the lookup table affect angular resolution and reconstruction accuracy in LZA conditions?
• RQ3 To what extent does the Disp method enhance the signal-to-noise ratio for faint sources (e.g., 1% Crab Nebula flux) at high zenith angles?
• RQ4 Can the Disp method enable high-significance detection of extended or weak sources like Sgr A* in LZA observations?

Key findings

• The Disp method maintains stable angular resolution across all zenith angles, with a 68% containment radius of 0.13° for Sgr A* observations, significantly improving over the default method.
• For the Crab Nebula, the method improves the signal-to-noise ratio by approximately 20% for a source with Crab Nebula flux and by 30% for a source with 1% of Crab Nebula flux at zenith angles above 50°.
• The default geometric reconstruction method shows a sharp decline in angular resolution with increasing zenith angle, while the Disp method exhibits minimal degradation, demonstrating robustness for LZA data.
• The first science application of the Disp method with VERITAS achieved a detection significance of over 11 standard deviations for Sgr A* in 15 hours of observation at an average zenith angle of 62°.
• The use of separate width and length parameters in the lookup table leads to better angular resolution than the previous width/length ratio approach, particularly in LZA conditions.
• The method effectively resolves the head-tail ambiguity by selecting the closest cluster of points from each ellipse, enabling accurate direction reconstruction even with highly flattened images.