Profiled Likelihoods from A Search for Dark Matter Lines at the Galactic Center with 14 Years of Fermi Data

This repository contains the profiled likelihood ratios calculated in "A Search for Dark Matter Lines at the Galactic Center with 14 Years of Fermi Data" [arXiv:2212.07435] for the unmasked analysis. Profiled likelihood ratios are provided for each of the 30 annuli considered in the work joined over the top three quartiles of data in energy dispersion.

The data are provided in an hdf5 archive containing 532 datasets. Line rest energies in units of GeV are provided in the `Line_Energy` dataset, which is a length 531 1-D array. Each of the remaining 531 datasets contains a 3-D array of shape (30, 2, 1001) corresponding to the profiled likelihood ratio as a function of line intensity, in units of counts/cm²/s/sr. The log of the profiled likelihood ratio for the i^th line rest energy in the j^th annulus may be accessed by:

archive = h5py.File('./Archive.h5', 'r')

line_energy = archive['Line_Energy'][i]
data = archive['Line_' + str(i)][j]

line_intensity = data[0]
log_profiled_likelihood_ratio = data[1]

archive.close()

For details, see https://arxiv.org/abs/2212.07435.