The set of log-concave densities on $\mathbb{R}^d$ forms a very attractive infinite-dimensional class. It is large enough to include many standard parametric families and yet it retains many of the desirable properties of Gaussian densities, e.g. closure under marginalisation, conditioning and convolution operations. Importantly, the class is also small enough to allow fully automatic estimation...