In [1]:
using GR
inline()
Out[1]:
"svg"
In [2]:
tevent = 1126259462.422

d = readdlm("strain.dat")
time, strain_H1_whitenbp, strain_L1_shift = d[:,1], d[:,2], d[:,3]

d = readdlm("nr_waveform.dat")
NRtime, NR_H1_whitenbp = d[:,1], d[:,2]
Out[2]:
([-0.620404,-0.62016,-0.619916,-0.619672,-0.619428,-0.619183,-0.618939,-0.618695,-0.618451,-0.618207  …  0.05318,0.053424,0.053668,0.053912,0.054156,0.054401,0.054645,0.054889,0.055133,0.055377],[-0.017659,-0.047861,-0.074998,-0.096486,-0.110634,-0.116817,-0.115441,-0.10772,-0.095359,-0.080203  …  0.432571,0.42737,0.411599,0.385199,0.348498,0.30223,0.247548,0.18603,0.119649,0.050696])
In [3]:
xlim([-0.1,0.05])
ylim([-4,4])

legend("H1 strain", "L1 strain", "matched NR waveform")
xlabel(@sprintf("time (s) since %.3f", tevent))
ylabel("whitented strain")
title("Advanced LIGO WHITENED strain data near GW150914")

plot(time-tevent, strain_H1_whitenbp, "r",
     time-tevent, strain_L1_shift, "g",
     NRtime+0.002, NR_H1_whitenbp, "k")
Out[3]:
- 4 - 2 0 2 4 - 0.10 - 0.05 0 0.05 Advanced LIGO WHITENED strain data near GW150914 time (s) since 1126259462.422 whitented strain H1 strain L1 strain matched NR waveform