GW data from Abbott et al. (2021b). Top panel: mass distribution of the individual BBHs (represented by a histogram of the mean values, with a bin-width of 2 M_⊙). We can also represent the datapoints as asymmetric Gaussians, as described in Appendix A. The sum of these Gaussians, ∑_idN_i/dm, is given as well (solid line), together with a cumulative distribution function (CDF, dotted line on a different axis). Bottom panel: chirp mass distribution (grey lines), which is made to resemble the distribution from Abbott et al. (2023). In order to do this, we use asymmetric Gaussians as well, taking into account the 90% credible interval. We also show an adjustable-width kernel density estimation (AWKDE, dashed line) from Shimazaki & Shinomoto (2010), Shimazaki et al. (2018). The Gaussians are plotted on a logarithmic axis, while the sum and the AWKDE use a linear axis. Similarly to Abbott et al. (2023), we omitted GW190814 from our analysis.