Gamma-Hadron Separation Methods for the VERITAS Array of Four Imaging Atmospheric Cherenkov Telescopes

Ground-based arrays of imaging atmospheric Cherenkov telescopes have emerged as the most sensitive gamma-ray detectors in the energy range of about 100 GeV and above. The strengths of these arrays are a very large effective collection area on the order of 100,000 square meter, combined with excellent single photon angular and energy resolutions. The sensitivity of such detectors is limited by statistical fluctuations in the number of Cosmic Ray initiated air showers that resemble gamma-ray air showers in many ways. In this paper, we study the performance of simple event reconstruction methods when applied to simulated data of the Very Energetic Radiation Imaging Telescope Array System (VERITAS) experiment. We review methods for reconstructing the arrival direction and the energy of the primary photons, and examine means to improve on their performance. For a software threshold energy of 300 GeV (100 GeV), the methods achieve point source angular and energy resolutions of sigma[63%]= 0.1 degree (0.2 degree) and sigma[68%]= 15% (22%), respectively. The main emphasis of the paper is the discussion of gamma-hadron separation methods for the VERITAS experiment. We find that the information from several methods can be combined based on a likelihood ratio approach and the resulting algorithm achieves a gamma-hadron suppression with a quality factor that is substantially higher than that achieved with the standard methods used so far.