Faster PET reconstruction with a stochastic primal-dual hybrid gradient method