gwastro / 3-ogc Goto Github PK

Hi again!

I see that you have a description in text on how to create the PSDs used for analysis:

We use 512 seconds of data around each candidate (with the exception of GW190725 174728) to
estimate the local power spectral density (PSD) using
a version of Welch’s method (Brown 2004). The data is
divided into overlapping segments (8s duration for BBH;
16s for BNS) which are Hann-windowed. The final PSD
estimate is the result of taking the mean of the median
average of the odd and even segments’ power spectrum.

Im happy to follow this and re-create the PSDs. However, I was wondering -- would it be possible to upload a script to create a PSD from the ini for an event?

Is it something like the following:

from pycbc.psd import welch, interpolate
ts = Merger("EVENT_NAME").strain(ifo)
psd = interpolate(welch(ts), 1.0 / ts.duration)

Alternatively, would it be possible to share the PSD files that were used for GW151226's analysis?

Thanks! :)

I think there might be a typo in the readme. Instead of the following,

import h5py

catalog_high = h5py.File('./3-OGC_low.hdf', 'r')

# Selecting parts of the catalog
region = candidates['mass1'][:] + candidates['mass2'][:] < 4
lowmass_snr = all_candidates['H1_snr'][:][region]

I think that it might be:

import h5py

candidates = h5py.File('./3-OGC_low.hdf', 'r')

# Selecting parts of the catalog
region = candidates['mass1'][:] + candidates['mass2'][:] < 4
lowmass_snr = candidates['H1_snr'][:][region]