Massively parallelizable list-mode reconstruction using a Monte Carlo-based elliptical Gaussian model