Repository navigation
Fixes a bug pointed out by Isobel Romero-Shaw. The issue was confusion betweeen classic det-frame/source-frame masses conversion. Basically, the redshift shouldn't have been used anywhere except while extracting the eccentricity at a specified frequency.
To reiterate, the general algorithm is:
- Use the source frame quantities (initial separation, initial eccentricity, masses) to get the e vs f curve. Interpolate this curve.
- If you want to extract e at fixed (det frame) f22, evaluate the interpolant at
f22 * (1 + redshift) - f instead you want e at fixed
M * f22, then evaluate the interpolant atM_times_f22 / Mtot_src
Thanks to Isobel for pointing this out!
Full Changelog: v0.1...v0.2