Gammapy: present status and future roadmap

TL;DR

Gammapy is an open-source Python library for high-energy astrophysics data analysis, now at version 1.0, supporting multi-instrument data reduction, modeling, and future enhancements like unbinned analysis and transient detection.

Contribution

This paper presents the current features and future roadmap of Gammapy, a widely-used analysis tool in high-energy astrophysics, highlighting new capabilities and planned developments.

Findings

Released Gammapy v1.0 with new analysis features

Demonstrated multi-instrument data reduction and modeling

Outlined future plans for unbinned analysis and transient detection

Abstract

Since its start in 2014, the lightweight open source Python library Gammapy has come a long way to become a popular data analysis package for high-energy astrophysics. Selected as the official CTAO Science Analysis tool, it is also an approved analysis software within the H.E.S.S. and MAGIC collaborations. The first long-term version, Gammapy v1.0 was released on late 2022. It is compliant with several well-established data conventions in high-energy astrophysics, and provides serialised data products that are interoperable with other software. Event lists and instrument response functions curated within the same format from various instruments can be reduced to data binned in energy, time or spatial coordinates. Thereafter, the flux and morphology of one or more gamma-ray sources can be estimated using Poisson maximum likelihood fitting and assuming a variety of spectral, temporal and spatial models. Flux points, likelihood profiles and light curves extractions are supported. Complex user defined likelihoods and models can also be implemented. In this contribution, we will highlight the main features of Gammapy v1.0, including data reduction and analysis examples from different space and ground-based instruments, applications of various background rejection techniques, and a simultaneous fitting across multiple instruments with astrophysical models. We will also present our plans for the future, showcasing new features such as the support of different event types, unbinned likelihood analysis, spectral unfolding and transient source detections. In addition to an improved API with distributed computing for scalable analysis, enhanced support for all-sky instruments like Fermi-LAT and HAWC is foreseen.