Remaking plots that appear in GWTC1

Here we will go through step by step how to make all plots that appear in the GWTC1 paper produced by the LIGO Virgo collaboration. This involves first downloading the data, and then using the summarypublication executable provided by PESummary.

# Here I will show you how to recreate the plots in GWTC1 using the
# summarypublication executable. Note that we cannot reproduce all of the
# plots in GWTC1 because the data release does not hold all of the required
# parameters

# First, we need to download the data from the DCC
curl -O https://dcc.ligo.org/public/0157/P1800370/004/GWTC-1_sample_release.tar.gz
tar -xf GWTC-1_sample_release.tar.gz

# Lets remove the priorChoice files as they do not contain posterior samples
FILES_TO_REMOVE=($(ls ./GWTC-1_sample_release/*priorChoices*))
for i in ${FILES_TO_REMOVE[@]}; do
    rm ${i}
done

# Now lets get a list of all of the files
FILES=($(ls ./GWTC-1_sample_release/*_GWTC-1.hdf5))

# Now lets get a list of labels
LABELS=()
for i in ${FILES[@]}; do
    label=`python -c "print('${i}'.split('_GWTC-1.hdf5')[0])"`
    LABELS+=(`python -c "print('${label}'.split('./GWTC-1_sample_release/')[1])"`)
done

# Now lets assign the same colors and linestyles that were used in GWTC1
COLORS=('#00a7f0' '#9da500' '#c59700' '#55b300' '#f77b00' '#ea65ff' '#00b1a4' '#00b47d' '#ff6c91' '#00aec2' '#9f8eff')
LINESTYLES=(solid dashed solid solid dashed solid solid dashed solid dashed dashed)

Now that we have all of the data and all of the variables setup, we can now run the summarypublication executable and make a 2d bounded contour of the mass_1 and mass_2 parameter space

summarypublication --plot 2d_contour \
                   --webdir ./GWTC-1_sample_release \
                   --samples ${FILES[@]} \
                   --labels ${LABELS[@]} \
                   --parameters mass_1 mass_2 \
                   --colors ${COLORS[@]} \
                   --linestyles ${LINESTYLES[@]} \

Now we can produced a violin plot showing the variation in mass_ratio

# Top left panel of figure 5
summarypublication --plot violin \
                   --webdir ./GWTC-1_sample_release \
                   --samples ${FILES[@]} \
                   --labels ${LABELS[@]} \
                   --parameters mass_ratio \

Now we can produce a 2d bounded contour of the theta_jn and luminosity_distance parameter space

# Left panel of Figure 7
summarypublication --plot 2d_contour \
                   --webdir ./GWTC-1_sample_release \
                   --samples ${FILES[@]} \
                   --labels ${LABELS[@]} \
                   --parameters theta_jn luminosity_distance \
                   --colors ${COLORS[@]} \

Now we can produce a 2d bounded contour of the luminosity_distance chirp_mass parameter space

# Right panel of Figure 8
summarypublication --plot 2d_contour \
                   --webdir ./GWTC-1_sample_release \
                   --samples ${FILES[@]} \
                   --labels ${LABELS[@]} \
                   --parameters luminosity_distance chirp_mass \
                   --colors ${COLORS[@]} \