Skip to content

Draft: Amanda's GWTC distribution

Michael William Coughlin requested to merge michael-coughlin/ligo.skymap:massdist into main

Potential PR to add Amanda's GWTC3 distribution to ligo.skymap.

Can read and plot her file as:

filename = "GWTC3_powerlaw_redshift_maxP_pdf.h5"
f = h5py.File(filename, 'r')
mass = f['mass']
p_mass_1 = f['p_mass_1']
plt.figure()
plt.semilogy(mass, p_mass_1)
plt.xlabel('primary mass [Msun]')
plt.ylabel('differential merger rate [Msun^-1 Gpc^-3 yr^-1]')
plt.savefig('p_m1_maxP.png')
plt.close()
Edited by Michael William Coughlin

Merge request reports