CosmoStatGW / gwfast

A Fisher information matrix python package for GW detector networks.
GNU General Public License v3.0
38 stars 13 forks source link

In Signals dictionary, if useEarthMotion=True, then astropy fails to compute LSMT from binary coalescence time if it is too long, i.e ~> Myr #4

Open nsteinle opened 5 months ago

nsteinle commented 5 months ago

I'm using gwfast to compute the SNR of binary neutron stars that merge within the age of the Universe. These BNSs are produced from a population synthesis code and thus each comes with a GW coalescence time, t_coal.

For a given detector, ie the CE with 40km baseline sensitivity, in the Signals dictionary if useEarthMotion=True, then astropy fails to compute LSMT from t_coal if it is too long, i.e ~> Myr, and gives an error because astropy does not know the motion of the Earth far in the future. See the astropy issue discussion here: https://github.com/astropy/astropy/issues/16294

Am I misusing this flag/feature? For now I can just set useEarthMotion=False and pass t_coal to the waveform to avoid the error from astropy when computing LMST from t_GPS via gwfastUtils.GPSt_to_LMST().

Is this feature not appropriate for astrophysically derived coalescence timescales, or am I missing something here?

Thanks!