Help with Data Challenge 3

Open Data Workshop
2

Hi! Thank you for your question!

Are you using the q_transform from pycbc or gwpy?

  1. The q-value should correspond to the number of cycles in the LIGO band. You can see some notes in the following places:

quickview app - Click “see notes” under the q-transform plot.

intro tutorials

  1. For pycbc, you could try something like the code below.
    For gwpy, you might look instead at the Quickview notebook 
t0 = 1204
data = pycbc_strain.time_slice(t0-10, t0+10)
times, freqs, power = data.qtransform(.001, logfsteps=100,
                                                qrange=(3, 10),
                                                frange=(10, 1024)
                                         )
plt.figure(figsize=[15, 3])
plt.pcolormesh(times, freqs, power**0.5)
plt.xlim(t0-2, t0+2)
plt.ylim(10,1024)
plt.yscale('log')
plt.clim(0,20)
plt.colorbar()
plt.show()