Study of Angular Resolution Using Imaging Atmospheric Cherenkov Technique

: Angular resolution is crucial for the detailed study of gamma-ray sources and current Cherenkov telescopes (e.g., HESS, MAGIC, and VERITAS) that operate below tens of TeV. Several gamma-ray sources with a photon energy larger than 100 TeV have been revealed by the LHAASO in recent years; the angular resolution of the LHAASO is around 0.3 ◦ . A gamma-ray detector with an angular resolution of less than 0.1 ◦ operating beyond 100 TeV is needed to study the detailed morphology of ultra-high-energy gamma-ray sources further. The cost-effectiveness is crucial for such large-area detectors. In this paper, the impact of telescope aperture, field of view, pixel size, optical point spread function, and signal integration time window on angular resolution is studied. These results can provide essential elements for the design of telescope arrays.


Introduction
Ground-based gamma-ray detectors detect primary gamma rays by measuring the air shower in the Earth's atmosphere.For gamma rays with an energy below ∼10 TeV, only a few secondaries reach ground level.However, there are Cherenkov photons reaching ground level, and these are emitted by high-energy charged secondary particles during the longitudinal development of an air shower [1].The imaging atmospheric Cherenkov technique (referred to as IACT hereafter) is an instrument used in ground-based gamma-ray astronomy to detect and study high-energy gamma rays from astrophysical sources [2].It involves using large optical telescopes equipped with fast photomultiplier tubes or silicon photomultipliers to detect Cherenkov photons in an air shower.The telescopes are typically positioned in an array, with multiple telescopes covering a large area to increase detection efficiency.By measuring the arrival time, intensity, and spatial distribution of the Cherenkov photons, the energy, direction, and nature (gamma ray or background particle) of the original gamma ray are reconstructed [3].An angular resolution of less than 0.1 degree above 10 TeV was achieved using this technique [4][5][6][7][8].
Angular resolution of the gamma-ray detector is essential for the detailed study of gamma-ray sources.For example, angular resolution is a key element in increasing the sensitivity of instruments because the signal/background ratio is proportional to the square of the angular resolution times the rejection efficiency for cosmic rays for the point source, and angular resolution is also critical for the morphology study of the gamma-ray source; the detection of fine structures in morphologies would help to identify the origin of the gamma rays.The inverse Compton scattering of electrons should show narrow structures, which are governed by the rapid cooling of the radiating electrons [9,10], whereas hadronic interactions are expected to generate much smoother structures [11,12].
Current IACT experiments (e.g., HESS, MAGIC, and VERITAS) operate below tens of TeV [5,7,8].Several gamma-ray sources with a photon energy larger than 100 TeV have been revealed by the Large High-Altitude Air Shower Observatory (LHAASO) in recent years [13]; the angular resolution of the LHAASO is around 0.3 • at 100 TeV [14].A gamma-ray detector with an angular resolution of less than 0.1 • operating beyond 100 TeV is needed to study the detailed morphology of those ultra-high-energy gamma-ray sources further [15].The cost-effectiveness is crucial for such large-area detectors.
The design of the detectors is mainly characterised by the configuration of each telescope and by the number and arrangement of these telescopes.The configuration of a single telescope is governed by several key factors.The diameter of the telescope (referred to as aperture below) is critical for obtaining sufficient Cherenkov statistics from air showers.Additionally, in order to reject the night-sky background, the camera needs to capture the Cherenkov photon pulse in a very short time (on the scale of nanoseconds).A good optical point spread function and finely pixelated image over a large field of view are critical for achieving excellent angular resolution of gamma-ray direction and good γ/proton discrimination abilities.Currently, most Cherenkov telescopes are based on single-dish optical systems, with mirror facets attached to either a spherical dish (e.g., Davies-Cotton) or a parabolic dish.The parabolic dish can reduce the time spread of the Cherenkov signal but suffers from significant off-axis aberrations, such as coma.The Davies-Cotton design provides compensation against spherical aberrations and coma.However, global coma still dominates for off-axis images and has significant consequences for the design of a wide field-of-view telescope.In general, off-axis distortions can be reduced by increasing the f-number (f/D, where f is the focal length and D is the aperture of the telescope), as coma scales use 1/f 2 .These single-dish designs are appealing due to their relatively low cost.In recent years, Schwarzschild-Couder telescopes with two mirror surfaces have been developed.The dual-mirror setup corrects spherical and coma aberrations, allowing for finer shower image pixelation and enhancing the optical point spread function and off-axis performance over a large field of view.Significant computing resources and time are required to optimise numerous parameters.
The impacts of the parameters of the telescope, including the pixel size of the camera, aperture, and trigger threshold, on the performances of the whole array were studied mainly for the optimisation of the Cherenkov Telescope Array [16][17][18][19][20].The muon detector of the LHAASO provides excellent γ/proton discrimination capabilities, with a high detection efficiency above ∼20 TeV [14].By operating in synergy with the LHAASO, an imaging Cherenkov telescope array at the LHAASO site combines the γ/proton discrimination capability of the LHAASO's muon detector with the excellent angular resolution of the IACT.This is very important for identifying processes responsible for gamma-ray production.Given that the muon detector has a fixed detection area, its γ/proton discrimination capability is primarily provided above ∼20 TeV.This paper investigates the impact of telescope parameters on the angular resolution of the IACT at the LHAASO site.However, the detection area and γ/proton discrimination ability of the IACT are also important.These results can provide essential elements for the design of a telescope and help narrow down the parameter space in a relatively simple manner without detailed simulations of specific detector parameters and the layout of the array.It is important to note that the detailed optimisation of the layout of telescope arrays is beyond the scope of this study.
Several methods have been developed to reconstruct the direction of gamma rays for IACT experiments.They can be divided into two classes: one is the transitional stereo-reconstruction method, which is based on the second moment parameterisation of the Cherenkov images [21].The shower direction is determined by a weighted mean of all pairwise intersections of the major axes of two suitable images mapped into a common co-ordinate system [22], which is easy to implement and is computationally inexpensive.The other one is based on the maximum likelihood method, which predicts the Cherenkov images of all telescopes based on 3D air shower models or MC templates; the shower direction is derived by performing a maximum likelihood fitting, which is more complicated and consumes more computing sources [23,24].In this study, the transitional stereo-reconstruction method will be used to reconstruct the shower direction.This paper is organised as follows: the simulation procedures, reconstruction method, and image cleaning are introduced in Section 2. The results regarding the impact of several parameters on the angular resolution are presented in Section 3. The differences in the results for the vertical incidence and 45 • zenith angle events, as well as the possible application of the results, are also discussed in Section 3, and, finally, Section 4 presents the summary.

Approach Used
There are numerous parameters to be determined for the design of a telescope; therefore, it is highly beneficial to narrow down the parameter space of the array in a manner that does not consume significant computing resources.In this study, our main focus is on the impact of the parameters of a single telescope on angular resolution.The angular resolution of gamma rays heavily depends on the number of triggered telescopes, the coordinates of the telescopes relative to the shower core position , and the array layouts.In reality, the array layout is fixed, and the core position of the shower is uniformly distributed on the detection plane.The positions of telescopes with respect to the shower core change for each event.As a result, the angular resolution of a sample is an average of many different arrangements of the telescopes with respect to the shower core.It is important to note that the shower core position in this study is always fixed at the origin, and the positions of the telescopes are located at various impact distances and different directions with respect to the shower core.The detector response of a single telescope is modelled using several parameters, which is similar to the method used in reference [19].The main idea is to avoid reproducing the data for each layout and each parameter of a single telescope.

Simulation
The simulation of extensive air showers (EASs) was performed using CORSIKA software (version v7.7410) [25].Only gamma rays were simulated; the EGS4 model was used for the electromagnetic interactions.The energy of the primary particles was fixed at log 10 (E/GeV) = 3.1, 3.7, 4.1, 4.5, and 5.1, representing an energy range from ∼1 TeV to ∼100 TeV.The detection plane was set at an altitude of 4400 m above sea level (the LHAASO's altitude), and the horizontal and vertical components of the Earth's magnetic field at the observation site were set to be 34.618µT and 36.13 µT, respectively.The US standard model was applied for the atmospheric profile.The simulation was performed at two fixed zenith angles (referred to as θ hereafter), which were θ = 0 • and θ = 45 • , corresponding to an atmosphere depth of approximately 600 g/cm 2 and 850 g/cm 2 , respectively, for the altitude chosen in this study.The azimuth angle of the incident gamma ray was fixed at 0 • to fix the perpendicular distance (referred to as Rp) between the telescope and the shower axis; the telescope's pointing direction remained parallel to the direction of the incoming primary particles.For both zenith angles, the number of simulated gamma-ray events was 2 × 10 4 , 2 × 10 4 , 2 × 10 4 , 1 × 10 4 , and 5 × 10 3 , and the energy log 10 (E/GeV) = 3.1, 3.7, 4.1, 4.5, and 5.1, respectively.
In extensive air showers, various secondary particles and Cherenkov photons are produced, including electrons, gamma rays, muons, neutrons, and more [26,27].However, in this study, our focus is solely on Cherenkov photons, specifically in the wavelength band of 200-1000 nm.In order to simulate reality as closely as possible, the atmospheric absorption and the quantum efficiency of the SiPM were applied; the atmospheric extinction coefficient was calculated using MODTRAN [28], and the result in reference [29] was used for quantum efficiency.The number of photons was converted into the number of photoelectrons after applying the quantum efficiency.Meanwhile, seven Rp values (Rp = 20, 50, 100, 150, 200, 300, and 400 m) and 12 positions (uniformly distributed, Rp = 20 m only has 4 positions to avoid shadowing between adjacent telescopes) were selected to sample the Cherenkov photons.The position of the telescopes for θ = 0 • is shown in Figure 1; the configuration "A/B/C/D", defined in Figure 1, was selected as the typical configuration so as to study the impact of telescope parameters on the angular resolution of gamma rays.
A circular region around each telescope (corresponding to the aperture of the telescope), was defined, and only photons within these specific regions at the observation level were recorded.The diameter of the circular region was chosen to be 20 m for log 10 (E/GeV) = 3.1, 12 m for log 10 (E/GeV) = 3.7 and 4.1, and 6 m for log 10 (E/GeV) = 4.5 and 5.1, due to limited disk space.Only a portion of the Cherenkov photons enter the detector response to simulate the different aperture and detection efficiency of the telescope.In order to study the telescope parameter space efficiently, a simplified telescope simulation tool similar to the one presented in reference [19] was developed.This tool mainly describes the telescope characteristics by three groups of parameters, without involving ray tracing or camera technology.The first group is related to design parameters, including the effective aperture, field of view, and the pixel size of the camera.The second group is related to the optical properties of the telescope, specifically the optical point spread function of the telescope.The last group is related to the analysis method, specifically the signal integration time window of the Cherenkov pulse.
The variable of the effective aperture (referred to as aperture hereafter) was used to represent the effects, including the physical aperture of the telescope, the reflecting efficiency of the mirrors, the collecting efficiency of the Winston cone, and any other efficiencies, such as the efficiency of the filter, if applicable.In order to capture the imaging response of the telescope optical system, a straightforward conversion of the Cherenkov photon direction to co-ordinates in the camera frame was employed.Each Cherenkov photon was given a random offset drawn from the optical point spread function, which was parameterised as a 2D Gaussian function with a 68.3% containment radius (referred to as the spot size below).The pixel shape was square, and there were no gaps between the pixels.The telescope's field of view is determined by removing the pixels with an angular distance to the origin of the camera frame that is less than the radius (referred to as the radius of the field of view).The night-sky background (referred to as the NSB) generated for each pixel follows a Poisson distribution, with an average most probable value of around 0.2 pe/m 2 /200 ns/(0.1 • × 0.1 • ) [30,31].This is proportional to the signal integration time window.The result of the simulated events in the camera plane is shown in Figure 2. The default parameters in this study were 6 m for the aperture, 0.1 • for the pixel size, 0.06 • for the spot size, 5 • for the radius of the field of view, and 200 ns for the signal integration time window.
It is important to note that the simulation did not account for the timing of signals and assumed an ideal data acquisition system limited only by irreducible noise.Therefore, the crosstalk, after-pulsing, timing jitter, nonlinearity, and saturation effects of SiPMs were also not simulated; only photo-electrons were simulated by applying quantum efficiency for the camera technology.The reliability of this simple model was verified with the well-tested sim_telarray package [19].

Image Analysis
The NSB was added to the Cherenkov images observed by the telescopes as noise; this needed to be cleaned to obtain a maximum signal from the shower itself and reject the maximum NSB contamination in each image.Traditionally, the tail-cut method is the standard image-cleaning method for the IACT [22,32]; in the tail-cut method, two thresholds (called thr1 and thr2) are applied to the image in the following way: all pixels with a photo-electron (referred to as p.e. below) signal above the higher threshold (thr1) are kept; in addition, pixels with signals above the lower threshold (thr2) are also kept if they are adjacent to at least one pixel with content above thr1.In recent years, other alternative cleaning methods have been developed [33], such as the single-threshold method [34].In the single-threshold method, only pixels with signals above the single threshold are kept.The advantage of the single-threshold method is that it is easy and quick to implement and provides reasonable cleaning regarding the image.The Cherenkov image of one telescope with a different cleaning method is shown in Figure 2. The images were cleaned and then subsequently analysed to reconstruct the shower direction of the gamma ray.The methods mentioned in this paper [22] that are currently in common use are referred to as traditional methods in this work.The direction of the gamma ray was determined by using a traditional method which involves calculating a weighted mean of all the pairwise intersections of the major axes of two suitable images mapped into a common co-ordinate system."Suitable", in this context, refers to images with at least three pixels that survive the cleaning procedures.This traditional method is based on the Hillas parameters of the Cherenkov image, derived from zero-order (amplitude or size, referred to as A), first-order (centre of gravity position), and second-order (width, length, and orientations, referred to as w, l, and ϕ, respectively) moments of the images.The weights for the intersection pairs were calculated using w ij = A 2 red sin 2 (ϕ i − ϕ j )δ 2 i δ 2 j , according to reference [22], where w ij is the weight for the pairwise intersection of i th , j th is the telescope, A red =A i A j /(A i + A j ) is the reduced amplitude of the paired telescope, and δ i = 1 − w i /l i is the simplified variant of the disp parameter for the i th telescope.
The distribution of ∆θ 2 for gamma-ray events with different energies detected by configuration A (see Figure 1) is shown in Figure 3.The blue vertical lines correspond to the region that contains 68.3% of the events, and the associated angle is defined as angular resolution.The shaded red area that contains the blue vertical line represents the region with event fractions of 68.3%−δ and 68.3%+δ; δ is the statistical error associated with the event fraction of 68.3%.The lower and upper errors of angular resolution can be derived from the shaded red region. .The distribution of ∆θ 2 for the gamma-ray events detected by configuration A. The energies of the samples are 10 3.1 GeV (black) and 10 5.1 GeV (purple), and the zenith angle is 45 degrees.∆θ represents the angle between the reconstructed and true gamma-ray directions.The blue vertical lines at ∼0.009 and ∼ 3 × 10 −4 correspond to the region that contains 68.3% of the events for E = 10 3. 1 GeV and E = 10 5.1 GeV, respectively.The shaded red area defines the error band of the region containing 68.3% of the events.Further details can be found in the main text.
The comparisons between the derived angular resolution vs. the threshold and between the tail-cut method and the single-threshold method are shown in Figure 4.The layout of the telescopes is configuration A and θ = 0 • .In the tail-cut method, the x-axis represents the threshold thr1, while the threshold thr2 is fixed at half of the high threshold.As demonstrated, both methods indicate that angular resolution initially worsens with the increasing threshold and then improves with further threshold increases before stabilising; this agrees with the results shown in reference [34].As explained in this reference, the single-threshold method is slightly better than the tail-cut method with a threshold of 15 p.e.; this is due to the improved exclusion of clumps of pixels at the periphery of the image caused by multiple scattered particles in the EAS.The threshold should be optimised to achieve better performance.A threshold that is too high results in a smaller collection area, especially for low-energy events, while a threshold that is too low results in poor angular resolution.In this study, the threshold was selected when the angular resolution stabilised, achieving a relatively good angular resolution, and the fraction of events that did not trigger the arrays A/B/C/D was negligible.The optimised threshold of 35 p.e. was found to be similar for all energies and zenith angles.The associated error of angular resolution (due to the difference between the optimal threshold and 35 p.e.) is indicated with the filled horizontal band in Figure 4.This error is around 10% at E = 10 3.1 GeV and 4% at E = 10 5.1 GeV; it contributes to the main error of angular resolution.It is important to note that for pixel sizes other than 0.

Results
In addition to the telescope parameters, the number of telescopes triggered and the distance between the telescopes and the shower core are also important for angular resolution.In this section, first, the impact of the number of telescopes triggered and the shower impact distance on angular resolution is studied.Then, the study of the impact of the aperture, field of view, pixel size, spot size, and the signal integration time window of the telescope on angular resolution is presented.

Number of Telescopes and Rp
In addition to the parameters of a single telescope, the layout of the telescope also plays a key role in determining angular resolution.Two particularly important factors are the number of telescopes triggered for stereo reconstruction and the distance between the telescopes and the core position.
The left panel of Figure 5 illustrates the angular resolution vs. the number of telescopes triggered for θ = 0 • and θ = 45 • , with the Rp of the telescopes fixed at 100 m and the positions of the telescopes uniformly distributed.Meanwhile, the middle panel of Figure 5 depicts the angular resolution vs. the Rp of the telescope for θ = 0 • and θ = 45 • .There are a total of four telescopes, and the directions of the telescopes with respect to the shower core position are the same as those in configuration A, as shown in Figure 1.As illustrated, the angular resolution improves when a higher number of telescopes are triggered.There is an optimal range of Rp values for angular resolution, which is around 100 m for θ = 0 • and 100-150 m for θ = 45 • .This is due to the telescopes that are far from the core being affected by statistical fluctuations.Conversely, in the region with small Rp, it is more challenging to reconstruct the shower detection plane (or SDP) accurately.It is also observed that at larger zenith angles, the angular resolutions of telescopes at large Rp values are significantly improved.This can be understood based on the lateral distribution (density of photoelectrons vs. Rp) comparison between θ = 0 • and θ = 45 • shown in the right panel of Figure 5.As seen, the density of photo-electrons is higher when the telescope is close to the shower axis for θ = 0 • compared to θ = 45 • , and it is lower when the telescope is further away from the shower axis.This is explained in reference [35].This is due to the fact that the observation plane is far away from the shower maximum for θ = 45 • , and the "Cherenkov pool" effect is more obvious.A configuration for the telescopes, namely configuration A with θ = 45 • , was selected to study the impact of singular telescope parameters on angular resolution.Three additional configurations (Rp = 20-200 m and the number of telescopes being four and six), as shown in Figure 1, were selected to study the dependency on the configuration.The results with different zenith angles and the possible application of the results are also discussed.

Signal Integration Time Window
We integrated the signal from extensive air showers and the night-sky background within a time window to mimic reality.This time window should capture the EAS signals while being as short as possible to minimise the NSB.The duration time determines the noise level of the telescope, which was around 16 ns for the HESS experiment [36].In order to simulate the high NSB level, a time window ranging from 0 ns to 200 ns was simulated, and the resulting angular resolution vs. the signal integration time window is depicted in Figure 6.As illustrated, the angular resolution does not appear to be affected much by the signal integration time window for all the energies and configurations studied in this work.

Field of View
The field of view of a telescope is typically limited due to various technical limitations; it was ±2.5 • for the HESS experiment [4].Increasing the field of view of telescopes is costly.For high-energy events and those with large Rp values, it is highly probable that the Cherenkov image may not be fully contained by the camera.In such instances, the event's reconstruction is likely to be affected.On the other hand, the shower-to-shower fluctuation of the Cherenkov image tail also contributes to worsening angular resolution.Figure 7 illustrates the angular resolution vs. the radius of the field of view for different energies and configurations.As observed, the angular resolution significantly improves with an increasing field of view, particularly at small field-of-view angles, and the angular resolution stabilised when the radius of the field of view was larger than 3 • for all the energies and configurations studied in this work.

Effective Aperture
As described in the simulation Section 2.1, the aperture in this study refers to the effective aperture, encompassing the physical aperture and all the associated efficiencies.A smaller aperture results in the collection of fewer Cherenkov photons, which affects the threshold energy and the detection area and increases the likelihood of the Cherenkov image being affected by statistical fluctuations.This can result in the deterioration of the angular resolution.Additionally, it is anticipated that the angular resolution will stabilise after reaching a certain aperture size when statistical fluctuation is not dominated.Given the significant cost associated with increasing the aperture of telescopes, it is crucial to carefully select an appropriate aperture for the telescope.The aperture for current Cherenkov telescopes ranges from less than 5 m to larger than 20 m.
Figure 8 depicts the angular resolution vs. the effective aperture of the telescope for different energies and configurations.It is evident that the angular resolution improves with the increasing aperture and has a tendency to stabilise after reaching a certain aperture size.The dotted lines in the left and right panels represent an exponential function, corresponding to a 5% improvement in the angular resolution for every 1 m increase in the aperture; it is used to guide the eye, although there are no specific physical reasons for its application.It is worth noting that at approximately 1 TeV, there is an improvement of about 20% in the angular resolution when transitioning from an aperture of 6 m to 10 m, after which the angular resolution stabilises., The dotted lines in both panels represent an exponential function, corresponding to a 5% improvement in angular resolution for every 1 m increase in the aperture; they are used to guide the eye, and there are no physical reasons for their application.

Pixel Size
The pixel size represents the field of view for each PMT or SiPM.The number of SiPMs (or PMTs) is inversely proportional to the pixel size for a given field of view.A larger pixel size will lead to images with worse resolution, thereby worsening angular resolution.Conversely, a smaller pixel size results in more SiPMs, but decreasing the pixel size of telescopes is also costly.The pixel size depends on the design of telescope; the typical pixel size for current Cherenkov telescopes is in the range of ∼0.07 • to ∼0.2 • [37].
Figure 9 illustrates angular resolution vs. the pixel size of the telescope for different energies and configurations.It is observed that the difference in angular resolution between a pixel size of 0.1 • and 0.15 • is minimal, but it worsens the angular resolution with further increases in pixel size.The dotted lines in Figure 9 represent an exponential function, indicating a 35% worsening in angular resolution for every 0.1 • increase in pixel size.Meanwhile, a solid green line is included for comparison in the right panel of Figure 9, which corresponds to a 15% decline in angular resolution for every 0.1 • increase in pixel size.As observed, the worsening of the angular resolution with increasing pixel size is around 15-35% for every 0.1 • increase in pixel size, and this worsening depends on energy and configuration.

Spot Size
The spot size, representing the radius that contains 68.3% of the optical point spread function, is one of the most critical parameters for telescopic performance.A larger spot size can lead to a worsening in the telescope's angular resolution, as it causes the light from a point source to be dispersed over a larger area on the camera, resulting in reduced image sharpness and detail.However, reducing the spot size can be very costly, making it crucial to thoroughly study the impact of spot size on angular resolution.The spot size depends on the design of the telescope; the typical spot size for current Cherenkov telescopes is less than 1 mrad.
Figure 10 depicts angular resolution vs. the spot size of the telescope for different energies and configurations.Since the pixel size may be correlated with spot size, the results for two pixel sizes (pixel size = 0.1 • and pixel size = 0.2 • ) are shown.It is evident that the angular resolution worsens with increasing spot size, particularly when the spot size is very small (less than 0.05 • ).To guide the eye, the black dotted lines shown in Figure 10 represent an exponential function, indicating a 30% worsening in the angular resolution for every 0.05 • increase in spot size.Meanwhile, the dotted orange line in the middle panel and the green lines in the right panel are included for comparison, which represent a 20% and 25% worsening in the angular resolution for every 0.05 • increase in spot size, respectively.As observed, the worsening of the angular resolution with increasing spot size is around 20-30% for every 0.05 • increase in spot size, and this worsening slightly depends on the energy.

Zenith Angle Dependence
The impact of the signal integration time window, field of view, aperture, pixel size, and spot size on angular resolution for θ = 45 • was studied.It was found that angular resolution is highly sensitive to the aperture, pixel size, and spot size of the telescope.For different gamma-ray sources located at varying zenith angles, it is meaningful to study those effects at different zenith angles.Comparisons of angular resolution vs. aperture, pixel size, and spot size between θ = 0 • and θ = 45 • are shown in Figures 11-13, respectively.GeV.The upper panel corresponds to configurations A and B , while the lower panel corresponds to configurations C and D (cfg in the legend is abbreviated from configurations).The magenta line represents an exponential function, corresponding to a 5% improvement in angular resolution for every 1 m increase in aperture; it is the same as the one shown in the left panel of Figure 8, and it is used to guide the eye.For the aperture, the dependencies of angular resolution on the aperture are similar between θ = 0 • (solid lines) and θ = 45 • (dotted lines), with the exception of configuration B at θ = 0 • .The greater dependency on the aperture for configuration B may be attributed to the fact that all telescopes have an Rp of 150 m; this is located outside the Cherenkov pool region for θ = 0 • and inside the Cherenkov pool for θ = 45 • with an energy of 10 3.1 GeV, as shown in the right panel of Figure 5.For pixel size and spot size, the angular resolution is less dependent on pixel size (or spot size) for θ = 0 • compared to θ = 45 • , especially at low energies.

Application
The worsening rate of angular resolution due to aperture, pixel size, and spot size can be utilised during the design of telescopes, with input from the cost model.A large field of view is important for observing extended sources, and it is also crucial for containing the Cherenkov image for high-energy and large-Rp events.However, as described in the Introduction section, single-dish telescopes have the advantage of lower cost, but the pixel size and the spot size for off-axis observations are limited.To achieve better off-axis performance, the f/D parameter needs to be optimised.The typical pixel size is around or less than 0.2 degrees.Alternatively, Schwarzschild-Couder telescopes with two mirror surfaces allow for finer shower image pixelation and a smaller spot size for off-axis instances over a large field of view, for example, a pixel size of 0.07 degrees, but they come with a much higher cost.An increase in pixel size leads to poorer angular resolution but requires fewer pixels (SiPMs or PMTs) and reduces costs, allowing for the possibility of more telescopes.Consequently, a larger detection area or enhanced γ/proton discrimination ability can be achieved for a fixed total cost of the array.This also applies to spot size and aperture.By incorporating a cost model as a function of aperture, pixel size, and spot size for single-dish and two-dish telescopes, the parameter space can be narrowed down for improved overall performance.

Summary
Several key telescope array parameters were investigated, and their influence on angular resolution was studied in detail.The parameters include the aperture, field of view, pixel size, spot size, signal integration time window of the telescope, number of telescopes used, and the distance between telescopes.It was found that the angular resolution is almost independent of the signal integration time window and field of view (except for a very small field of view).The angular resolution is very sensitive to the array layout, aperture, pixel size, and spot size.The angular resolution performs better when the Rp of the telescope is around 100 m for θ = 0 • and 100-150 m for θ = 45 • .The angular resolution improves with increasing aperture and has a tendency to stabilise after reaching a certain aperture size; at approximately 1 TeV, the improvement in the angular resolution is about 20% when transitioning from an aperture of 6 m to 10 m.The worsening rate of the angular resolution with pixel size and spot size is around 15-35% for every 0.1 • increase in pixel size and around 20-30% for every 0.05 • increase in spot size.These results provide important information for the design of the telescope array.
Author Contributions: Generation of Monte Carlo simulation data, F.Z.; the effect of the number and distance of telescopes on angular resolution, H.W.; the effect of angular resolution on the shape of the telescope placement, Q.L.; the effect of angular resolution on the telescope field of view, R.X.; the effect of angular resolution on the telescope pixel size, J.L.; data check, Y.J.; draft version and verification of results, H.L. All authors have read and agreed to the published version of the manuscript.
Funding: This work was funded by the National Natural Science Foundation of China (Grant number 12205244) and "the Fundamental Research Funds for the Central Universities" (Project Number 2023160).
Data Availability Statement: Data will be made available upon request.

Figure 1 .
Figure 1.The location of telescopes on the observational plane (indicated by the black dots), along with the four configurations of telescopes (labelled A, B, C, and D in the legend).These four configurations are used as typical configurations for the study in this work.The red cross at the centre indicates the position of the shower core.

Figure 2 .
Figure 2.The simulated Cherenkov image is depicted with no image cleaning (left), cleaned using the single-threshold cleaning method (middle, with a threshold of 35 p.e.), and cleaned using the tail-cut clean method (right, with thresholds of 35 p.e. and 17.5 p.e.).The image was observed by using a telescope located at Rp = 100 m, and the primary energy of this gamma-ray event is 10 4.1 GeV.The color bar represents the number of photo-electrons detected per pixel (referred to as Npe).

1 E=10Figure 3
Figure3.The distribution of ∆θ 2 for the gamma-ray events detected by configuration A. The energies of the samples are 10 3.1 GeV (black) and 10 5.1 GeV (purple), and the zenith angle is 45 degrees.∆θ represents the angle between the reconstructed and true gamma-ray directions.The blue vertical lines at ∼0.009 and ∼ 3 × 10 −4 correspond to the region that contains 68.3% of the events for E = 10 3.1 GeV and E = 10 5.1 GeV, respectively.The shaded red area defines the error band of the region containing 68.3% of the events.Further details can be found in the main text.

1 •Figure 4 .
Figure 4.The angular resolution against the threshold used during the image-cleaning procedures for configuration A and θ = 0 • .The red lines indicate an energy of 10 3.1 GeV, while the blue lines indicate an energy of 10 5.1 GeV.The solid lines represent the application of the tail-cut clean method, and the dashed lines represent the application of the single-threshold cleaning method.Details can be found in the main text.

Figure 5 .
Figure 5.The angular resolution vs. the number of telescopes triggered (left), and the angular resolution vs. the Rp (perpendicular distance between the telescope and shower axis) of the telescope (middle).The different colors represent different energy levels, with the solid lines corresponding to θ = 0 • and the dashed lines corresponding to θ = 45 • .Detailed information about the configuration can be seen in the main text.The right panel is the mean density of the photo-electrons (total photo-electrons divided by the effective area of the telescope) vs. the Rp of the telescope.

Figure 6 .
Figure 6.The angular resolution vs. signal integration time window for different energy levels (left) and for different configurations of telescopes (right).

Figure 7 .
Figure 7.The angular resolution vs. the radius of the field of view of the telescope for different energy levels (left) and different telescope configurations (right).

Figure 8 .
Figure8.Angular resolution vs. the aperture of the telescope for different energy levels (left) and different telescope configurations (right).The dotted lines in both panels represent an exponential function, corresponding to a 5% improvement in angular resolution for every 1 m increase in the aperture; they are used to guide the eye, and there are no physical reasons for their application.

Figure 9 .
Figure 9. Angular resolution vs. telescope pixel size for different energy levels (left) and different telescope configurations (right).All the dotted lines represent an exponential function; the black dotted lines in both panels indicate a 35% degradation in angular resolution for every 0.1 • increase in pixel size.The dotted green lines in the right panel represent a 15% degradation in angular resolution for every 0.1 • increase in pixel size; they are used to guide the eye.

Figure 10 .
Figure10.Angular resolution vs. telescope spot size for different energy levels with a pixel size of 0.1 • (left), different energy levels with a pixel size of 0.2 • (middle), and different telescope configurations (right).All the dotted lines represent an exponential function; the black dotted lines in all three panels indicate a 30% degradation in angular resolution for every 0.05 • increase in spot size.The dotted orange line in the middle panel represents a 20% degradation in angular resolution for every 0.05 • increase in spot size.The dotted green lines in the right panel represent a 25% degradation in angular resolution for every 0.05 • increase in spot size; they are used to guide the eye.

Figure 11 .
Figure 11.The comparison of angular resolution vs. telescope aperture between θ = 0 • (solid lines) and θ = 45 • (dotted lines).The left panel corresponds to an energy of 10 3.1 GeV, while the right panel corresponds to an energy of 10 5.1 GeV.The upper panel corresponds to configurations A and B , while the lower panel corresponds to configurations C and D (cfg in the legend is abbreviated from configurations).The magenta line represents an exponential function, corresponding to a 5% improvement in angular resolution for every 1 m increase in aperture; it is the same as the one shown in the left panel of Figure8, and it is used to guide the eye.

Figure 12 .Figure 13 .
Figure 12.The comparison of angular resolution vs. telescope pixel size between θ = 0 • (solid lines) and θ = 45 • (dotted lines).The left panel corresponds to an energy of 10 3.1 GeV, while the right panel corresponds to an energy of 10 5.1 GeV.The upper panel corresponds to configurations A and B, while the lower panel corresponds to configurations C and D (cfg in the legend is abbreviated from configurations).The magenta line represents an exponential function, indicating a 35% worsening in angular resolution for every 0.1 • increase in pixel size; it is the same as the one shown in the left panel of Figure9, and it is used to guide the eye.