[Paper Review] CORSIKA Simulation of the Telescope Array Surface Detector

Telescope Array Collaboration|arXiv (Cornell University)|Mar 4, 2014

Astronomical Observations and Instrumentation11 citations

TL;DR

This paper presents a high-accuracy CORSIKA Monte Carlo simulation framework for the Telescope Array surface detector, using detailed shower generation and validation against real data to calculate detector acceptance with 27% energy scale correction. The method achieves excellent agreement between simulation and data, enabling precise cosmic ray energy spectrum measurements with a total systematic uncertainty of 37% due to energy scale.

ABSTRACT

The Telescope Array is the largest experiment studying ultra-high energy cosmic rays in the northern hemisphere. The detection area of the experiment consists of an array of 507 surface detectors, and a fluorescence detector divided into three sites at the periphery. The viewing directions of the 38 fluorescence telescopes point over the air space above the surface array. In this paper, we describe a technique that we have developed for simulating the response of the array of surface detectors of the Telescope Array experiment. The two primary components of this method are (a) the generation of a detailed CORSIKA Monte Carlo simulation with all known characteristics of the data, and (b) the validation of the simulation by a direct comparison with the Telescope Array surface detector data. This technique allows us to make a very accurate calculation of the acceptance of the array. We also describe a study of systematic uncertainties in this acceptance calculation.

Motivation & Objective

• To develop a high-fidelity CORSIKA Monte Carlo simulation of the Telescope Array surface detector for accurate acceptance calculations.
• To validate the simulation against real Telescope Array surface detector data to ensure reliability in energy spectrum measurements.
• To quantify systematic uncertainties in the acceptance and energy scale for cosmic ray flux measurements.
• To correct the simulated energy scale using fluorescence detector-calibrated air shower data, improving spectral analysis accuracy.
• To enable precise exposure calculations for cosmic ray anisotropy and flux studies.

Proposed method

• The study uses the CORSIKA Monte Carlo code with the QGSJET-II-03 interaction model to simulate extensive air showers for proton primaries.
• Simulations include full detector geometry, scintillator counter response, and realistic atmospheric conditions to model surface detector response.
• Detector response is validated by comparing simulated and real data in counter pulse height, time, and lateral distribution functions.
• A 27% energy scale correction is derived by comparing simulated air showers with fluorescence detector-calibrated calorimetric measurements.
• Systematic uncertainties are evaluated using bin-by-bin energy resolution corrections and spectral index consistency checks.
• The aperture calculation accounts for energy resolution and spectral shape, including the GZK cutoff, to minimize resolution-induced bias.

Experimental results

Research questions

• RQ1 How accurately can CORSIKA simulations model the response of the Telescope Array surface detector to ultra-high energy cosmic rays?
• RQ2 What is the magnitude and origin of systematic uncertainties in the detector acceptance calculation?
• RQ3 To what extent does the energy scale correction derived from fluorescence data improve the accuracy of the surface detector simulation?
• RQ4 How well do simulated lateral distributions and counter pulse heights match real data across energy and zenith angle ranges?
• RQ5 What is the impact of energy resolution and spectral shape on the final flux measurement uncertainty?

Key findings

• The dethinned CORSIKA + QGSJET-II-03 simulation accurately models the Telescope Array surface detector response for energies >10^18.2 eV and zenith angles <45°.
• Simulated counter pulse heights and arrival times show excellent agreement with real data, confirming simulation fidelity.
• A 27% energy scale correction is derived from fluorescence detector-calibrated air shower data, significantly improving simulation accuracy.
• The systematic uncertainty in the flux due to energy scale is 37%, dominated by the 22% uncertainty in the fluorescence energy scale.
• All other systematic contributions are <1%, confirming that energy scale uncertainty is the dominant source of flux uncertainty.
• The simulation enables precise aperture calculations for cosmic ray energy spectrum and anisotropy studies with minimal resolution-induced bias.