CLUMPY: Jeans analysis, γ-ray and neutrino fluxes from dark matter (sub-)structures

We present an update of the CLUMPY code for the calculation of the astrophysical J-factors (from dark matter annihilation/decay) for any Galactic or extragalactic dark matter halo including substructures: halo-to-halo concentration scatter may now be enabled, boost factors can include several levels of substructures, and triaxiality is a new option for dark matter haloes. This new version takes advantage of the CFITSIO and HEALPix libraries to propose FITS output maps using the HEALPix pixelisation scheme. Skymaps for $\gamma$-ray and $\nu$ signals from generic annihilation/decay spectra are now direct outputs of CLUMPY. Making use of HEALPix routines, smoothing by a user-defined instrumental Gaussian beam and computing the angular power spectrum of the maps is now possible. In addition to these improvements, the main novelty is the implementation of a Jeans analysis module, to obtain dark matter density profiles from kinematic data in relaxed spherical systems (e.g., dwarf spheroidal galaxies). The code is also interfaced with the GreAT toolkit designed for Markov Chain Monte Carlo analyses, from which probability density functions and credible intervals can be obtained for velocity dispersions, dark matter profiles, and J-factors.