Use distance grid instead of marginal_cdf, marginal_ppf

This speeds up find_injection_moc when it is called with a large number of true 3D sky positions because it avoids an expensive n^2 loop.