• Energy Research

Context.Traditionally, TeV-γ-ray astronomy has been conducted by experiments employing proprietary data and analysis software. However, the next generation ofγ-ray instruments, such as the Cherenkov Telescope Array Observatory (CTAO), will be operated as open observatories. Alongside the data, they will also make the associated software tools available to a wider community. This necessity prompted the development of open, high-level, astronomical software customized for high-energy astrophysics.Aims.In this article, we present Gammapy, an open-source Python package for the analysis of astronomicalγ-ray data, and illustrate the functionalities of its first long-term-support release, version 1.0. Built on the modern Python scientific ecosystem, Gammapy provides a uniform platform for reducing and modeling data from differentγ-ray instruments for many analysis scenarios. Gammapy complies with several well-established data conventions in high-energy astrophysics, providing serialized data products that are interoperable with other software packages.Methods.Starting from event lists and instrument response functions, Gammapy provides functionalities to reduce these data by binning them in energy and sky coordinates. Several techniques for background estimation are implemented in the package to handle the residual hadronic background affecting γ-ray instruments. After the data are binned, the flux and morphology of one or moreγ-ray sources can be estimated using Poisson maximum likelihood fitting and assuming a variety of spectral, temporal, and spatial models. Estimation of flux points, likelihood profiles, and light curves is also supported.Results.After describing the structure of the package, we show, using publicly available gamma-ray data, the capabilities of Gammapy in multiple traditional and novelγ-ray analysis scenarios, such as spectral and spectro-morphological modeling and estimations of a spectral energy distribution and a light curve. Its flexibility and its power are displayed in a final multi-instrument example, where datasets from different instruments, at different stages of data reduction, are simultaneously fitted with an astrophysical flux model.

Context.Traditionally, TeV-γ-ray astronomy has been conducted by experiments employing proprietary data and analysis software. However, the next generation ofγ-ray instruments, such as the Cherenkov Telescope Array Observatory (CTAO), will be operated as open observatories. Alongside the data, they will also make the associated software tools available to a wider community. This necessity prompted the development of open, high-level, astronomical software customized for high-energy astrophysics.Aims.In this article, we present Gammapy, an open-source Python package for the analysis of astronomicalγ-ray data, and illustrate the functionalities of its first long-term-support release, version 1.0. Built on the modern Python scientific ecosystem, Gammapy provides a uniform platform for reducing and modeling data from differentγ-ray instruments for many analysis scenarios. Gammapy complies with several well-established data conventions in high-energy astrophysics, providing serialized data products that are interoperable with other software packages.Methods.Starting from event lists and instrument response functions, Gammapy provides functionalities to reduce these data by binning them in energy and sky coordinates. Several techniques for background estimation are implemented in the package to handle the residual hadronic background affecting γ-ray instruments. After the data are binned, the flux and morphology of one or moreγ-ray sources can be estimated using Poisson maximum likelihood fitting and assuming a variety of spectral, temporal, and spatial models. Estimation of flux points, likelihood profiles, and light curves is also supported.Results.After describing the structure of the package, we show, using publicly available gamma-ray data, the capabilities of Gammapy in multiple traditional and novelγ-ray analysis scenarios, such as spectral and spectro-morphological modeling and estimations of a spectral energy distribution and a light curve. Its flexibility and its power are displayed in a final multi-instrument example, where datasets from different instruments, at different stages of data reduction, are simultaneously fitted with an astrophysical flux model.

Abstract The Cherenkov Telescope Array (CTA) is a next-generation ground-based observatory for gamma-ray astronomy at very high energies. The Large-Sized Telescope prototype (LST-1) is located at the CTA-North site, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to ≃20 GeV. LST-1 started performing astronomical observations in 2019 November, during its commissioning phase, and it has been taking data ever since. We present the first LST-1 observations of the Crab Nebula, the standard candle of very-high-energy gamma-ray astronomy, and use them, together with simulations, to assess the performance of the telescope. LST-1 has reached the expected performance during its commissioning period—only a minor adjustment of the preexisting simulations was needed to match the telescope’s behavior. The energy threshold at trigger level is around 20 GeV, rising to ≃30 GeV after data analysis. Performance parameters depend strongly on energy, and on the strength of the gamma-ray selection cuts in the analysis: angular resolution ranges from 0.°12–0.°40, and energy resolution from 15%–50%. Flux sensitivity is around 1.1% of the Crab Nebula flux above 250 GeV for a 50 hr observation (12% for 30 minutes). The spectral energy distribution (in the 0.03–30 TeV range) and the light curve obtained for the Crab Nebula agree with previous measurements, considering statistical and systematic uncertainties. A clear periodic signal is also detected from the pulsar at the center of the Nebula.

Abstract The Cherenkov Telescope Array (CTA) is a next-generation ground-based observatory for gamma-ray astronomy at very high energies. The Large-Sized Telescope prototype (LST-1) is located at the CTA-North site, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to ≃20 GeV. LST-1 started performing astronomical observations in 2019 November, during its commissioning phase, and it has been taking data ever since. We present the first LST-1 observations of the Crab Nebula, the standard candle of very-high-energy gamma-ray astronomy, and use them, together with simulations, to assess the performance of the telescope. LST-1 has reached the expected performance during its commissioning period—only a minor adjustment of the preexisting simulations was needed to match the telescope’s behavior. The energy threshold at trigger level is around 20 GeV, rising to ≃30 GeV after data analysis. Performance parameters depend strongly on energy, and on the strength of the gamma-ray selection cuts in the analysis: angular resolution ranges from 0.°12–0.°40, and energy resolution from 15%–50%. Flux sensitivity is around 1.1% of the Crab Nebula flux above 250 GeV for a 50 hr observation (12% for 30 minutes). The spectral energy distribution (in the 0.03–30 TeV range) and the light curve obtained for the Crab Nebula agree with previous measurements, considering statistical and systematic uncertainties. A clear periodic signal is also detected from the pulsar at the center of the Nebula.