GWTC-4.0: Population Properties of Merging Compact Binaries

GWTC-4.0: Population Properties of Merging Compact Binaries

Welcome to the data release associated with GWTC-4.0: Population Properties of Merging Compact Binaries. This is a large data release; if you run into problems please alert us.

Description of files in this data release

analyses_AllCBC.tar: Contains popsummary format files for all analyses that are part of the section “Binary Merger Population Across All Masses”

analyses_NS.tar: Contains popsummary format files for all analyses that are part of the section “Population Properties of Mergers Containing Neutron Stars”

analyses_BBH.tar: Contains popsummary format files for all analyses that are part of the section “Binary Black Hole Population”, except those that use the Binned Gaussian Processes (BGP) model.

analyses_BGP.tar: Contains popsummary format files for all analyses that are part of the section “Binary Black Hole Population” and use the Binned Gaussian Processes (BGP) model.

figure_scripts.tar: Contains all .py scripts needed to produce figures in the paper.

figures.tar: Contains .pdf files of all figures in the paper

o4a_event_list.tar: Contains list of events used in the paper, along with the label used in the PE samples file.

download_gwtc3_data.py: Download necessary GWTC-3 data to reproduce figures.

popsummary_tutorial.ipynb: A short tutorial that shows how to deal with popsummary files.

Setup instructions

Creating a conda environment

Use the following commands to set up a conda environment called o4a-astro and install required packages:

conda create --name o4a-astro python=3.11 -y
conda activate o4a-astro
pip install numpy==2.3.2 scipy==1.16.1 matplotlib==3.10.5 seaborn==0.13.2 bilby==2.6.0 popsummary==0.0.1 zenodo-get==2.0.0

Setup directory structure

Run the following on the bash terminal

mkdir o4a-astro
cd o4a-astro
mkdir data_release figures figure_scripts

# If you want to get the figure .py scripts
zenodo_get 16911563 -g figure_scripts.tar -o figure_scripts
tar -xvf figure_scripts/figure_scripts.tar -C figure_scripts

# If you want to get the figure pdf files
zenodo_get 16911563 -g figures.tar -o figures
tar -xvf figures/figures.tar -C figures

Optionally, if you want to download GWTC-3 data, run the following:

zenodo_get 16911563 -g download_gwtc3_data.py
python download_gwtc3_data.py

Running download_gwtc3_data.py will create a directory called gwtc3_data in the o4a-astro directory.

The final ingredient to run the figure scripts is downloading the popsummary themselves. Below is a table that says which TAR file each script uses, and which input h5 files from inside the TAR archive it uses. Download the requisite TAR files, and extract them inside the o4a-astro/data_release. Thereafter, you can enter the o4a-astro/figure_scripts directory and run scripts using python figure_<figureID>.py. The figure file will be saved in the o4a-astro/figures directory (any existing files will be overwritten).

Figure Script
TAR File
Input Files

figure_1.py
analyses_AllCBC.tar
AllCBC_FullPop.h5AllCBC_FullPopBGP.h5

figure_2.py
analyses_AllCBC.tar
AllCBC_FullPop.h5AllCBC_FullPopBGP.h5

figure_3.py
analyses_BBH.tar, gwtc3_data
BBHMassSpinRedshift_BrokenPowerLawTwoPeaks_GaussianComponentSpins_PowerLawRedshift.h5BBHMassSpinRedshift_BSplineIID.h5gwtc3_data/PowerLawPeak

figure_4.py
analyses_BBH.tar
BBHMassSpinRedshift_BrokenPowerLawTwoPeaks_GaussianComponentSpins_PowerLawRedshift.h5BBHMassSpinRedshift_BrokenPowerLawTwoPeaks_GaussianComponentSpins_PowerLawRedshift_GWTC3.h5

figure_5.py
analyses_BBH.tar, gwtc3_data
BBHMassSpinRedshift_BrokenPowerLawTwoPeaks_GaussianComponentSpins_PowerLawRedshift.h5BBHMassSpinRedshift_BSplineIID.h5gwtc3_data/PowerLawPeak

figure_6.py
analyses_BBH.tar
BBHMassSpinRedshift_BSplineIsopeakIID.h5BBHMass_VaryingBetaQs_DominantMode.h5

figure_7.py
analyses_BBH.tar, gwtc3_data
BBHMassSpinRedshift_BSplineIID.h5BBHSpin_MagTruncnormIidTiltIsotropicTruncnormNid.h5gwtc3_data/PowerLawPeak

figure_8.py
analyses_BBH.tar
BBHSpin_MagTruncnormIidTiltIsotropicTruncnormSpinSorting.h5

figure_9.py
analyses_BBH.tar, gwtc3_data
BBHSpin_EpsSkewNormalChiEff.h5gwtc3_data/GaussianSpin

figure_10.py
analyses_BBH.tar, gwtc3_data
BBHMassSpinRedshift_BrokenPowerLawTwoPeaks_GaussianComponentSpins_PowerLawRedshift.h5BBHMassSpinRedshift_BSplineIID.h5gwtc3_data/PowerLawPeak

figure_11_and_14.py
analyses_BBH.tar
BBHCorr_qchieffLinearCorrelationModel.h5BBHCorr_zchieffLinearCorrelationModel.h5BBHCorr_qchieffSplineCorrelationModel.h5BBHCorr_zchieffSplineCorrelationModel.h5

figure_12.py
analyses_BBH.tar
BBHCorr_qchieffCopulaCorrelationModel.h5BBHCorr_qchieffLinearCorrelationModel.h5

figure_13.py
analyses_BGP.tar
BBHCorr_MassSpinCorrelatedBGPModel.h5

figure_15.py
analyses_BBH.tar, gwtc3_data
BBHMassSpinRedshift_BrokenPowerLawTwoPeaks_GaussianComponentSpins_PowerLawRedshift.h5BBHMassSpinRedshift_BrokenPowerLawOnePeak_GaussianComponentSpins_PowerLawRedshift.h5BBHMassSpinRedshift_BrokenPowerLawThreePeaks_GaussianComponentSpins_PowerLawRedshift.h5BBHMassSpinRedshift_PowerLawPeak_GaussianComponentSpins_PowerLawRedshift.h5BBHMassSpinRedshift_BSplineIID.h5gwtc3_data/PowerLawPeak

figure_16.py
analyses_BBH.tar
BBHMassSpinRedshift_BrokenPowerLawTwoPeaks_GaussianComponentSpins_PowerLawRedshift_DominantMode.h5BBHMassSpinRedshift_BrokenPowerLawTwoPeaks_GaussianComponentSpins_PowerLawRedshift_SubdominantMode.h5BBHMassSpinRedshift_BrokenPowerLawOnePeak_GaussianComponentSpins_PowerLawRedshift.h5

figure_17.py
analyses_BBH.tar
BBHSpin_MagTruncnormIndTiltIsotropicTruncnormNnd.h5BBHSpin_MagTruncnormIidTiltIsotropicTruncnormNid.h5

figure_18_and_19.py
analyses_BBH.tar, gwtc3_data
BBHSpin_GaussianChiEffChiP.h5BBHSpin_GaussianChiEffChiP_NeffCut.h5gwtc3_data/GaussianSpin

figure_20.py
analyses_BBH.tar, analyses_BGP.tar
BBHMass_Autoregressive.h5BBHCorr_MassRedshiftCorrelatedBGPModel.h5BBHMassSpinRedshift_BSplineIID.h5BBHMass_FlexibleMixtures.h5BBHMassSpinRedshift_BrokenPowerLawTwoPeaks_GaussianComponentSpins_PowerLawRedshift.h5BBHMassSpinRedshift_BSplineIID_NoMinM.h5

figure_21.py
analyses_BBH.tar, analyses_BGP.tar, gwtc3_data
BBHCorr_MassSpinCorrelatedBGPModel.h5BBHSpin_EpsSkewNormalChiEff.h5BBHSpin_BSplineChiEff.h5BBHSpin_GaussianChiEffChiP.h5BBHCorr_qchieffSplineCorrelationModel.h5gwtc3_data/GaussianSpin

figure_22.py
analyses_BBH.tar, analyses_BGP.tar
BBHCorr_zchieffCopulaCorrelationModel.h5BBHCorr_MassRedshiftCorrelatedBGPModel.h5BBHMass_FlexibleMixtures.h5