snr_ecc_stationary
- legwork.snr.snr_ecc_stationary(m_c, f_orb, ecc, dist, t_obs, harmonics_required, interpolated_g=None, interpolated_sc=None, ret_max_snr_harmonic=False, ret_snr2_by_harmonic=False, **kwargs)[source]
Computes SNR for eccentric and stationary sources
- Parameters
- m_cfloat/array
Chirp mass
- f_orbfloat/array
Orbital frequency
- eccfloat/array
Eccentricity
- distfloat/array
Distance to the source
- t_obsfloat
Total duration of the observation
- harmonics_requiredinteger
Maximum integer harmonic to compute
- interpolated_gfunction
A function returned by
scipy.interpolate.interp2d
that computes g(n,e) from Peters (1964). The code assumes that the function returns the output sorted as with the interp2d returned functions (and thus unsorts). Default is None and uses exact g(n,e) in this case.- interpolated_scfunction
A function returned by
scipy.interpolate.interp1d
that computes the LISA sensitivity curve. Default is None and uses exact values. Note: take care to ensure that your interpolated function has the same LISA observation time ast_obs
and uses the same instrument.- ret_max_snr_harmonicboolean
Whether to return (in addition to the snr), the harmonic with the maximum SNR
- ret_snr2_by_harmonicboolean
Whether to return the SNR^2 in each individual harmonic rather than the total. The total can be retrieving by summing and then taking the square root.
- **kwargsvarious
Keyword args are passed to
legwork.psd.power_spectral_density()
, see those docs for details on possible arguments.
- Returns
- snrfloat/array
SNR for each binary
- max_snr_harmonicint/array
harmonic with maximum SNR for each binary (only returned if
ret_max_snr_harmonic=True
)