fd_model_difference

GWCorrect.wfu.parameterization.fd_model_difference(hf1,hf2,injection=None,npoints=None,
polarization='plus',psd_data=None,correction_parameter=0.0001,ref_amplitude=None)

Generates frequency domain waveform differences between two models hf1 and hf2.

\[\begin{split}\Delta\mathcal{A}_{\mu}(f;\vartheta)= \begin{cases} \frac{\mathcal{A}_2(f;\vartheta)}{\mathcal{A}_1(f;\vartheta)}-1 & f \leq f_{\mathrm{disc}} \\ \Delta\mathcal{A}_\mu(f_{\mathrm{disc}};\vartheta) & f > f_{\mathrm{disc}} \end{cases}\end{split}\]

\[\begin{split}\Delta\phi_{\mu}(f;\vartheta)= \begin{cases} \phi_2(f;\vartheta)-\phi_1(f;\vartheta)-2\pi ft_0-\phi_0 & f \leq f_{\mathrm{disc}} \\ \Delta\phi_\mu(f_{\mathrm{disc}};\vartheta) & f > f_{\mathrm{disc}} \end{cases}\end{split}\]

Parameters:

hf1: bilby.gw.waveform_generator.WaveformGenerator: frequency domain waveform generator object
hf2: bilby.gw.waveform_generator.WaveformGenerator: frequency domain waveform generator object
injection: dictionary, optional, None: dictionary of injection parameters if waveform generators do not have parameters; if they do not, this argument is not optional
npoints: int, optional, (None): length of the desired frequency grid; if None, this will be set automatically by the frequency grid
polarization: string, optional, (‘plus’): polarization of the strain data {‘plus’,’cross’}
psd_data: numpy.ndarray, optional, (None): array containing the power spectral density data and their corresponding frequencies; if None, psd_data will be set to GW170817 PSD data
correction_parameter: float, optional, (0.0001): fraction of the peak frequency domain amplitudes at which to cut off the amplitudes; part of the f_disc calculation
ref_amplitude: numpy.ndarray, optional, (None): reference amplitude for aligned phase calculation; will be generated automatically if not given

Returns:

frequency_grid: numpy.ndarray: frequency grid that corresponds to the uncertainty arrays
amplitude_difference: numpy.ndarray: array of amplitude difference values
phase_difference: numpy.ndarray: array of phase difference values; if psd data is None, unaligned_phase_difference will be returned, aligned_phase_difference otherwise