We test the Higgs dilaton inflation model (HDM) using the latest cosmological datasets, including the cosmic microwave background temperature, polarization and lensing data from the Planck satellite (2015), the BICEP and Keck Array experiments, the type Ia supernovae from the JLA catalogue, the baryon acoustic oscillations from CMASS, LOWZ and 6dF, the weak lensing data from the CFHTLenS survey and the matter power spectrum measurements from the latest SDSS data release. We find that the values of all cosmological parameters allowed by the HDM are well within the Planck satellite (2015) constraints. In particular, we determine w(0) = -1.0001(-0.0074)(+0.0072), w(a) = 0.00(-0.16)(+0.15), n(s) = 0.9693(-0.0082)(+0.0083), a(s) = -0.001(-0.014)(+0.013) and r(0.05) = 0.0025(-0.0016)(+0.0017) (at 95.5% c.l.). We also place new stringent constraints on the couplings of the HDM, xi(chi) < 0.00328 and xi(h)/root lambda = 59200(-20000)(+30000) (at 95.5% c.l.). We find that the HDM is only slightly better than the w(0)w(a)CDM model, with Delta chi(2) = chi(w0waCDM)-C-2 - chi(2)(HDM) = 0.18. Given that the HDM has two fewer parameters, we find Bayesian evidence favouring the HDM over the w(0)w(a)CDM model. We also study the critical Higgs inflation model, taking into account the running of both the self-coupling lambda(mu) and the non-minimal coupling to gravity xi(mu). We find peaks in the curvature power spectrum at scales corresponding to the critical value mu that re-enter during the radiation era and collapse to form a broad distribution of clustered primordial black holes, which could constitute today the main component of dark matter. This article is part of the Theo Murphy meeting issue 'Higgs cosmology'.