scri.SpEC.com_motion

scri.SpEC.com_motion(path_to_horizons_h5, path_to_matter_h5=None, m_A=None, m_B=None)

Calculate the center-of-mass motion from Horizons.h5 and/or Matter.h5

Note that no input checking is done here. If the data sets are missing, incomplete, or corrupted, you should expect to get exceptions or silently incorrect data.

Parameters:

path_to_horizons_h5string, optional [default is ‘Horizons.h5’]

Absolute or relative path to ‘Horizons.h5’

path_to_matter_h5: None or string [default is None]

Absolute or relative path to ‘Matter.h5’ for neutron star trajectories. If None, we assume that only black holes are involved, so that both A and B trajectories will be in Horizons.h5.

m_A: None or float [default is None]

If None, the mass will be read from Horizons.h5 (for BHs) or Matter.h5 (for NSs).

m_B: None or float [default is None]

If None, the mass will be read from Horizons.h5 (for BHs) or Matter.h5 (for NSs).

Returns:

tfloat array

These are simply the coordinate times

com3xN float array

Position of CoM at each instant of time