Gravitational waves offer a unique probe into the properties of merging black holes (BHs) and neutron stars.
Since the first such detection in 2015, the LIGO-Virgo network [1]}, [2]}, [3]} has reported the detection of \(\sim 90\)  binary black hole (BBH) mergers, with each gravitational-wave (GW) signal encoding physical information about the BHs involved, such as their masses and angular momenta (spins) [4]}, [5]}, [5]}. Extracting this information has has enabled the study of properties of BBH systems on both an individual and population-level basis. From an astrophysical perspective, combining GW detections to infer the mass, spin, and redshift distributions of BBH systems can help answer questions ranging from binary formation and stellar evolution [7]}, [8]} to the expansion rate of the Universe and possible deviations from General Relativity [9]}, [10]}.
