The impact of electroweak Sudakov logarithms on the endpoint of the photon spectrum for wino dark matter annihilation was studied intensively over the last several years. In this work, we extend these results to Higgsino dark matter $\chi_1^0$. We achieve NLL' resummation accuracy for narrow and intermediate spectral energy resolutions, of order $m_W^2 / m_\chi$ and $m_W$, respectively. This is the most accurate prediction to date for the yield of high-energy $\gamma$-rays from $\chi_1^0 \chi_1^0 \to \gamma + X$ annihilation for the energy resolutions realized by current and next-generation telescopes. We also discuss for the first time the effect of power corrections in $m_W/m_\chi$ in this context and argue why they are not sizeable.

The annihilation cross section of TeV scale dark matter particles $\chi^0$ with electroweak charges into photons is affected by large quantum corrections due to Sudakov logarithms and the Sommerfeld effect. We calculate the semi-inclusive photon energy spectrum in $\chi^0\chi^0\to \gamma+X$ in the vicinity of the maximal photon energy $E_\gamma = m_\chi$ with NLL' accuracy in an all-order summation of the electroweak perturbative expansion adopting the pure wino model. This results in the most precise theoretical prediction of the annihilation rate for $\gamma$-ray telescopes with photon energy resolution of parametric order $m_W^2/m_\chi$ for photons with TeV energies.

Bayesian inference is often used in cosmology and astrophysics to derive constraints on model parameters from observations. This approach relies on the ability to compute the likelihood of the data given a choice of model parameters. In many practical situations, the likelihood function may however be unavailable or intractable due to non-gaussian errors, non-linear measurements processes, or complex data formats such as catalogs and maps. In these cases, the simulation of mock data sets can often be made through forward modeling. We discuss how Approximate Bayesian Computation (ABC) can be used in these cases to derive an approximation to the posterior constraints using simulated data sets. This technique relies on the sampling of the parameter set, a distance metric to quantify the difference between the observation and the simulations and summary statistics to compress the information in the data. We first review the principles of ABC and discuss its implementation using a Population Monte-Carlo (PMC) algorithm and the Mahalanobis distance metric. We test the performance of the implementation using a Gaussian toy model. We then apply the ABC technique to the practical case of the calibration of image simulations for wide field cosmological surveys. We find that the ABC analysis is able to provide reliable parameter constraints for this problem and is therefore a promising technique for other applications in cosmology and astrophysics. Our implementation of the ABC PMC method is made available via a public code release.