In the advent of large scale surveys, individually modelling strong-gravitational lenses and their counterpart time-delays in order to precisely measure H 0 will become computationally expensive, and highly complex. A complimentary approach is to study the cumulative distribution function (CDF) of time-delays where the global population of lenses is modelled along with H 0 . In this paper we use a suite of hydro-dynamical simulations to estimate the CDF of time-delays from doubly-imaged quasars for a realistic distribution of lenses. We find that the CDFs exhibit large amounts of halohalo variance, regulated by the density profile inner slope and the total mass within 5kpc. With the objective of fitting to data, we compress the CDFs using Principal Component Analysis and fit a Gaussian Processes Regressor consisting of three physical features: the redshift of the lens, z L ; the power law index of the halo, α, and the mass within 5kpc, plus four cosmological features. Assuming a flat Universe, we fit our model to 27 doubly-imaged quasars finding Hand Ω Λ = 0.7 +0.04 −0.04 . We compare our estimates of z L and log(M (< 5kpc)/M ) to the data and find that within the sensitivity of the data, they are not systematically biased. We generate mock CDFs and find with that the Vera Rubin Observatory (VRO) could measure σ/H 0 to < 3%, limited by the precision of the model. If we are to exploit fully VRO, we require simulations that sample a larger proportion of the lens population, with a variety of feedback models, exploring all possible systematics.