Gaussian noise is wrong!!
The implementation of Gaussian noise is not working. Minimal example:
%pylab inline
import peyote
time_duration = 1.
sampling_frequency = 4096.
H1 = peyote.detector.H1
H1_hf_noise, frequencies = H1.power_spectral_density.get_noise_realisation(
sampling_frequency, time_duration)
H1.set_data(sampling_frequency, time_duration,
frequency_domain_strain=H1_hf_noise)
plt.loglog(H1.frequency_array, np.abs(H1.data))
plt.loglog(H1.frequency_array, np.abs(H1.amplitude_spectral_density_array))
You will see that the noise realisation this creates does not overlap with the ASD.
@colm.talbot or @nikhil.sarin can you please take a look at this? @gregory.ashton and I are sitting in Seattle trying to work this out. We think there is some double counting of terms between create_white_noise (in utils.py) and get_noise_realisation (in detector.py).