Comparison of the fitted SED and profiles with data for the models with p-value >0.003. The best fit is obtained for n = 4.5, η = 0.1, α = 1.8, δ = 1.0, B = 1 µG and is highlighted compared to other models, with fitted parameters of cm² s⁻¹ and TeV. The uncertainty bands for H.E.S.S. indicates the systematic uncertainty with statistical uncertainties indicated by the error bars, while for HAWC the uncertainty band represents the systematic and statistical uncertainty.