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)

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.RectBivariateSpline that computes g(n,e) from Peters (1964). 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)