GammaBayes: a Bayesian pipeline for dark matter detection with CTA

Liam Pinchbeck,Eric Thrane,Csaba Balazs

DOI: https://doi.org/10.1088/1475-7516/2024/05/020

2024-05-03

Journal of Cosmology and Astroparticle Physics

Abstract:We present GammaBayes, a Bayesian Python package for dark matter detection with the Cherenkov Telescope Array (CTA). GammaBayes takes as input the CTA measurements of gamma rays and a user-specified dark-matter particle model. It outputs the posterior distribution for parameters of the dark-matter model including the velocity-averaged cross section for dark-matter self interactions 〈σv〉 and the dark-matter mass m χ . It also outputs the Bayesian evidence, which can be used for model selection. We demonstrate GammaBayes using 525 hours of simulated data, corresponding to 10 8 observed gamma-ray events. The vast majority of this simulated data consists of noise, but 100000 events arise from the annihilation of scalar singlet dark matter with m χ = 1 TeV. We recover the dark matter mass within a 95% credible interval of m χ ∼ 0.96–1.07 TeV. Meanwhile, the velocity averaged cross section is constrained to 〈σv〉 ∼ 1.4–2.1 × 10 -25 cm 3 s -1 (95% credibility). This is equivalent to measuring the number of dark-matter annihilation events to be N S ∼ 1.1 -0.2 +0.2 × 10 5 . The no-signal hypothesis 〈σv〉 = 0 is ruled out with about 5σ credibility. We discuss how GammaBayes can be extended to include more sophisticated signal and background models and the computational challenges that must be addressed to facilitate these upgrades. The source code is publicly available here.

astronomy & astrophysics,physics, particles & fields