In galaxy clustering measurements, summary statistics, e.g., 2-pt correlation functions and power spectra, have been widely used to extract cosmological information. These statistics summarise the information of the observed galaxy number density field and theoretical models such as perturbation theory can accurately predict the statistics up to mildly non-linear regime. On the other hand, the...