This is the data release associated with https://arxiv.org/abs/2209.06978 Contains all of the hyper posterior samples for the runs listed in Tables G.1. The files are in json format. Bilby offers a dedicated routine to read them in import bilby data= bilby.core.result.read_in_result(path_to_json) See the Bilby documentation for what is contained in the result object. For each run, we report the posterior hyper samples for the mass model, reshift model, spin magnitude model, spin tilt model and merger rate [Gpc^-3 yr^-1] Here the name used to store and a short description of each parameter: Primary mass model (Power Law + Peak for all runs) power_law_slope_m1, slope of the primary mass power law component minmass_m1, minimum BH mass maxmass_m1, maximum BH mass low_end_smoothing_m1, smoothing at the low-mass end peak_branchingratio_m1, branching ratio between Gaussian peak and power law (1= 100% peak) peak_mean_m1, mean of the Gaussian peak peak_sigma_m1, sigma of the Gaussian peak Mass ratio model (power law for all runs) power_law_slope_mass_ratio, slope of the mass ratio Redshift (power law for all runs) power_law_slope_redshift, slope of the redshift Spin magnitude (IID beta distributions for all runs) alpha_chi, first argument of beta distribution beta_chi, second argument of beta distribution Cosine of tilt angle Gaussian models mu_0_costilt, for Gaussian models w/o correlation, the mean of the left (or only) Gaussian sigma_0_costilt, for Gaussian models w/o correlation, the sigma of the left (or only) Gaussian mu_1_costilt, for Gaussian models w/o correlation, the mean of the right Gaussian sigma_1_costilt, for Gaussian models w/o correlation, the sigma of the right Gaussian mu_a_costilt, for Gaussian model with correlation, the constant part of the Gaussian mean mu_b_costilt, for Gaussian model with correlation, the coefficient of the linearly evolving part of the Gaussian mean sigma_a_costilt, for Gaussian model with correlation, the constant part of the Gaussian sigma sigma_b_costilt, for Gaussian model with correlation, the coefficient of the linearly evolving part of the Gaussian sigma Beta models alpha_a_costilt, for all Beta models, the constant part of the first parameter of the Beta distribution alpha_b_costilt, for all Beta models, the coefficient of the linearly evolving part of the first parameter of the Beta distribution beta_a_costilt, for all Beta models, the constant part of the second parameter of the Beta distribution beta_b_costilt, for all Beta models, the coefficient of the linearly evolving part of the second parameter of the Beta distribution Tukey models: tukey_x0, the center of the Tukey as defined in appendix E of the paper tukey_k, Tk as defined in appendix E of the paper tukey_r, Tk as defined in appendix E of the paper Branching ratios: spin_mixture_0, for 2-component models, this is the branching ratio of the non-isotropic component spin_mixture_1, for Isotropic + Gaussian + Tukey and Isotropic + Gaussian + Beta this is the branching ratio of the Gaussian component; for Isotropic + 2 Gaussian this is the branching ratio of the Gaussian on the right. Merger rate rates, merger rate per unit Gpc cubed per unit year Note that some of the parameters for the tilt models might not be used, but still stored (and fixed to - usually - zero). This can be checked by verifying what priors were used for each parameter. For example the Isotropic run was obtained from the Isotropic + Gaussian model by setting the branching ratio of the Gaussian component to zero (at which point the values of mu and sigma costitl are irrelevant) > data['prior'] {'alpha_chi': Uniform(minimum=1, maximum=5, name='alphachi', latex_label='$\\alpha_\\chi$', unit=None, boundary=None), 'beta_chi': Uniform(minimum=1, maximum=5, name='betachi', latex_label='$\\beta_\\chi$', unit=None, boundary=None), 'spin_mixture_0': DeltaFunction(peak=0, name=None, latex_label=None, unit=None), 'power_law_slope_mass_ratio': Uniform(minimum=-2, maximum=7, name=None, latex_label='$\\beta$', unit=None, boundary=None), 'power_law_slope_m1': Uniform(minimum=1, maximum=6, name=None, latex_label='$\\alpha$', unit=None, boundary=None), 'mu_0_costilt': DeltaFunction(peak=0, name=None, latex_label=None, unit=None), 'mu_1_costilt': DeltaFunction(peak=0.0, name=None, latex_label=None, unit=None), 'sigma_0_costilt': DeltaFunction(peak=1, name=None, latex_label=None, unit=None), 'sigma_1_costilt': DeltaFunction(peak=0.0, name=None, latex_label=None, unit=None), 'peak_branchingratio_m1': Uniform(minimum=0, maximum=0.25, name=None, latex_label='$\\lambda$', unit=None, boundary=None), 'power_law_slope_redshift': Uniform(minimum=-2, maximum=10, name='lamb', latex_label='$\\lambda_z$', unit=None, boundary=None), 'peak_mean_m1': Uniform(minimum=20, maximum=45, name=None, latex_label='$\\mu_{m}$', unit=None, boundary=None), 'peak_sigma_m1': Uniform(minimum=1, maximum=10, name=None, latex_label='$\\sigma_{m}$', unit=None, boundary=None), 'minmass_m1': Uniform(minimum=2, maximum=10, name=None, latex_label='$m_{\\min}$', unit=None, boundary=None), 'maxmass_m1': Uniform(minimum=60, maximum=100, name=None, latex_label='$m_{\\max}$', unit=None, boundary=None), 'low_end_smoothing_m1': Uniform(minimum=0, maximum=10, name=None, latex_label='$\\delta_{m}$', unit=None, boundary=None)} Drop me (Salvatore Vitale) an email if anything doesn't work or if you spot issues. Thanks!

{"references": ["https://arxiv.org/abs/2209.06978"]}