Abstract:SURA is a self-triggered radio array on the roof of physics faculty at Semnan university in Iran. It is designed to detect radio emissions from air showers produced by ultra-high energy (UHE) cosmic rays with energies exceeding 1017 eV. The array consists of 4 LPDA radio antennas operating in the 40 MHz to 80 MHz range. In this study, we present a method that compares the signal intensities of simulated and experimental data. Specifically, we use a simulated dense array with a large number of antennas as a reference. By comparing the experimental signal intensity of each antenna to that of the corresponding antenna in the reference array, we can reconstruct the cosmic ray air shower core location. We first validate our method on simulated events to estimate the associated error. Afterward, we apply the technique to the cosmic ray candidates detected by the SURA array. Our results show that the core location can be reconstructed with a minimum error of about 6 m. However, when the characteristics of the shower being reconstructed differ significantly from the reference array, the error increases. Finally, we propose optimizations to improve reconstruction accuracy and reduce simulation time.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: how to reconstruct the core position of the extensive air shower (EAS) generated by ultra - high - energy cosmic rays (UHECRs) in the atmosphere by comparing the radio intensity observed in experiments with that in simulation data. Specifically, the author proposes a method. Using a dense array as a reference, the core position of the air shower is reconstructed by comparing the signal intensity received by the SURA experimental array with that of the corresponding antennas in the simulation array. ### Problem Background When ultra - high - energy cosmic rays (with energy exceeding \(10^{17}\) eV) pass through the Earth's atmosphere, they will trigger extensive air showers (EAS), which contain a large number of secondary particles. Studying these showers can help scientists understand the origin, composition, and energy distribution of cosmic rays. However, because the flux of ultra - high - energy cosmic rays drops sharply as the energy increases, direct detection becomes very difficult. Therefore, researchers usually adopt indirect methods, such as detecting secondary particles generated by EAS and analyzing their lateral distribution functions (LDF) to study these rays. ### Research Objectives The main objective of this paper is to develop a method that can accurately reconstruct the core position of cosmic - ray air showers by comparing experimental data and simulation data. Specifically, the author uses a dense simulation array as a reference and determines the core position of the air shower by comparing the actually observed signal intensity with that of the corresponding antennas in the simulation array. ### Method Overview 1. **Simulating a Dense Array**: The author first simulates a dense radio antenna array, which contains a large number of antennas and is used to generate reference signals. 2. **Comparing Signal Intensities**: The signal intensity received by the SURA experimental array is compared with that of the corresponding antennas in the simulation array. 3. **Core Position Reconstruction**: By minimizing the \(\chi^2\) statistic, the experimental signal closest to the simulation signal is found, so as to determine the core position of the air shower. The formula is as follows: \[ \chi^2_{ijk}=\chi^2(x_i, y_j)=\frac{1}{n}\sum_{k = 1}^{n}\left(\frac{E_{\text{sim},ijk}-E_{\text{exp},k}}{E_{\text{exp},k}}\right)^2 \] where: - \((x_i, y_j)\) is the tested core position; - \(k\) is the antenna number (1 to 4 for SURA); - \(E_{\text{exp},k}\) is the experimental signal received by the \(k\) - th antenna; - \(E_{\text{sim},ijk}\) is the signal received by the simulated virtual antenna; - \(C_{ij}\) is a scaling factor used to adjust the intensity difference between the simulation signal and the experimental signal. 4. **Optimization and Verification**: The author also explores how to optimize the method by reducing the scale of the dense array and evaluates the influence of the primary energy and the incident direction on the reconstruction accuracy of the core position. ### Results and Conclusions By this method, the author successfully reduces the reconstruction error of the core position of the air shower to about 6 meters. However, when the characteristics of the shower are significantly different from those of the reference array (for example, when there are large differences in primary energy or incident direction), the error will increase. For this reason, the author proposes further optimization methods to improve the reconstruction accuracy and reduce the simulation time. In general, this research provides an effective method for reconstructing the core position of ultra - high - energy cosmic - ray air showers, which is of great significance for understanding the properties of cosmic rays.

Reconstruction of cosmic ray air shower core location at SURA experiment

Estimation of aperture of the Tunka-Rex radio array for cosmic-ray air-shower measurements

Reconstruction of air shower muon lateral distribution functions using integrator and binary modes of underground muon detectors

Radio wavefront of very inclined extensive air-showers: a simulation study for extended and sparse radio arrays

A reconstruction procedure for near horizon extensive air showers based on radio signals

Towards a Cosmic-Ray Energy Scale with the Auger Engineering Radio Array

A new potential method for the $X_{\rm max}$ measurement of extensive air showers based on backtracking radio signals

A high-precision interpolation method for pulsed radio signals from cosmic-ray air showers

Analytical planar wavefront reconstruction and error estimates for radio detection of extensive air showers

Particle content of very inclined air showers for radio signal modeling

The Auger Engineering Radio Array and multi-hybrid cosmic ray detection (TAUP 2015)

Energy reconstruction of hadron-initiated showers of ultra-high energy cosmic rays

Results of a self-triggered prototype system for radio-detection of extensive air showers at the Pierre Auger Observatory

Reconstruction of air-shower measurements with AERA in the presence of pulsed radio-frequency interference

Validation of straight-line signal propagation for radio signal of very inclined cosmic ray air showers

Cosmic ray mass independent energy reconstruction method using Cherenkov light and muon content in LHAASO

Detection of extensive cosmic ray air showers by measuring radio emission

Background identification algorithm for future self-triggered air-shower radio arrays

Modelling uncertainty of the radiation energy emitted by extensive air showers

UHECR primary identification using the lateral profile of muons in EAS