DC1 PIT

This short notebook aims to show the difference in PIT values computed now and in the past (for DC1 paper).

Three sets of PIT values are comapared:

  1. PITs computed now using qp CDF method
  2. PITs computed now using scipy.integrate.quad
  3. PITs computed in the past at the time of writing DC1 paper (reference sample).

A subset of the data (1%) is used for the purpose of illustration. Computing the PITs using scipy.integrate.quad would take many hours for the ~399k galaxies present in DC1 sample.

Conclusion: using qp causes small discrepancies in the PIT values compared to the reference results (DC1), which would be negligible if using the scipy integration method. However, the latter is not as efficient as qp, and can make the computation unfeasible if the dataset is larger.