Making a scatter plot for a population of events

With PESummary you can scatter plot displaying the distributions of different parameters estimated from the gravitational waves reeleased by the LIGO Collaboration. Here we show an example for the redshift vs total mass distribution for all the GWTC-1 data. This is done including errorbars on the plot. Below the image example and the script used to produce it.

# Here we will show you how to create a population scatter plot using the
# summarypublication executable for the GWTC1 data release.

# 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 make the population plot
summarypublication --plot population_scatter_error \
                   --webdir ./GWTC-1_sample_release \
                   --samples ${FILES[@]} \
                   --labels ${LABELS[@]} \
                   --parameters total_mass redshift \