Galaxy N-point function statistics provide a powerful tool for extracting non-Gaussian information. To efficiently compute N-point functions, various approaches have been proposed, including radial-angular decomposition in harmonic space and compressed statistics. In particular, the statistics involving four fields can be designed to be an effective tool for extracting parity-odd features from the field, where these odd features may originate from the beyond standard model physics. In this talk, I will discuss the application of the galaxy four-point function, the connections between different algorithms involving four fields, and the potential implications of using the galaxy four-point functions to link late-time and early-time physics.