swh:1:snp:3a699297f000109a1bc833f294a54171df990207
Tip revision: 5d35c811cec3cdd131433fafe86880715993a372 authored by Ian Harry on 01 March 2019, 17:17:05 UTC
Prepping for release (#2535)
Prepping for release (#2535)
Tip revision: 5d35c81
data.py
import pycbc.catalog, pylab
m = pycbc.catalog.Merger("GW170817")
# Retreive data from the Livingston observatory around the BNS merger
ts = m.strain("L1").time_slice(m.time - 30, m.time + 6)
# Whiten the data with a 4s filter
white = ts.whiten(4, 4)
times, freqs, power = white.qtransform(.01, logfsteps=200,
qrange=(110, 110),
frange=(20, 512))
pylab.pcolormesh(times, freqs, power**0.5)
pylab.yscale('log')
pylab.ylabel("Frequency (Hz)")
pylab.xlabel("Time (s)")
pylab.show()