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 as t_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)