The calibrations of DAMPE $γ$-ray effective area

Zhao-Qiang Shen,Wen-Hao Li,Kai-Kai Duan,Wei Jiang,Zun-Lei Xu,Chuan Yue,Xiang Li
2024-10-03
Abstract:The DArk Matter Particle Explorer (DAMPE) is a cosmic-ray detector as well as a pair-converting $\gamma$-ray telescope. The effective area, reflecting the geometrical cross-section area, the $\gamma$-ray conversion probability and the photon selection efficiency, is important in the $\gamma$-ray analyses. In the work, we find a significant time variation in the effective area, as large as $\sim -4\%/{\rm yr}$ at 2 GeV for the high-energy trigger. We derive the data-based correction factors to the effective areas and apply corrections to both the effective areas and the exposure maps. The calibrated exposure can be $\sim 12\%$ smaller than the Monte Carlo one on average at 2 GeV. The calibration is further verified using the observation of the Vela pulsar, showing the spectral parameters with the correction are more consistent with those in the Fermi-LAT catalog than the ones without correction. All the corrections are now implemented in the latest version of the DAMPE $\gamma$-ray analysis toolkit DmpST.
Instrumentation and Methods for Astrophysics,High Energy Astrophysical Phenomena,Instrumentation and Detectors
What problem does this paper attempt to address?
This paper attempts to solve the problem of the time - varying effective area of DAMPE (Dark Matter Particle Explorer) γ - rays. Specifically, the main objectives of the paper include: 1. **Identify and correct the time - dependent changes in the effective area**: It has been found that there are significant time - varying changes in the effective area of DAMPE's High - Energy Trigger (HET) and Low - Energy Trigger (LET). For example, at an energy of 2 GeV, the HET effective area decreases by about 4% per year. Such changes will affect the results of γ - ray spectral analysis of astrophysical sources. 2. **Propose data - driven correction factors**: In order to correct these time - dependent changes, the authors propose correction factors based on observational data and apply them to the effective area and exposure maps. The corrected exposure is on average about 12% smaller than the results of Monte Carlo simulations at 2 GeV. 3. **Verify the effectiveness of the correction**: Verification is carried out by using the observational data of the Vela pulsar. The results show that the corrected spectral parameters are more consistent with those in the Fermi - LAT catalogue, proving the effectiveness of the correction method. ### Specific problem analysis #### 1. Identification of time - dependent changes The paper first identifies the time - dependent changes in the effective area by comparing the γ - ray photon counts observed in different time periods with the expected counts from Monte Carlo simulations. A detailed analysis was carried out for the HET and LET trigger logics respectively. In particular, the rate of change of the HET effective area is about - 4%/year at 2 GeV, and less than 1%/year above 5 GeV. #### 2. Data - driven correction factors To correct these changes, the authors introduce two correction factors: - **Primary correction factor \( f(E) \)**: Defined according to flight data, it reflects the change in the effective area with energy. - **Secondary correction factor \( g(E,\theta) \)**: Based on MC simulations, it takes into account the effect of changes in the trigger threshold. The final correction formula is: \[ A_{\text{cal}}^{\text{eff}, H}(E, \theta, \phi, t) = A_{\text{MC}}^{\text{eff}, H}(E, \theta, \phi) \times [1 - f(E)g(E, \theta)\Delta t] \] where \( E \) is the incident energy, \( \theta \) is the incident inclination angle, \( \phi \) is the incident azimuth angle, \( \Delta t=t - t_{\text{ref}} \), and \( t_{\text{ref}} \) is the reference time (selected as July 1, 2016). #### 3. Verification of the effectiveness of the correction Verification is carried out by using the data of the Vela pulsar. The results show that the corrected spectral parameters are more consistent with those in the Fermi - LAT catalogue. Specifically, the corrected maximum likelihood value \( \Delta \ln L \) increases significantly, indicating that the correction method effectively solves the problem of mismatch between the observed counts and the expected counts. ### Summary This paper improves the accuracy of γ - ray spectral analysis of astrophysical sources by identifying and correcting the time - dependent changes in the DAMPE γ - ray effective area. This work not only improves the data processing flow of DAMPE, but also provides valuable experience and methods for future similar high - energy astronomical observations.