def final_remnant_properties_from_NRSurrogate(*args, f_low=20.0, f_ref=20.0,
    model='NRSur7dq4Remnant', return_fits_used=False, properties=[
    'final_mass', 'final_spin', 'final_kick'], approximant='SEOBNRv4PHM'):
    """Return the properties of the final remnant resulting from a BBH merger using
    NRSurrogate fits

    Parameters
    ---------
    f_low: float/np.ndarray
        The low frequency cut-off used in the analysis. Default is 20Hz
    f_ref: float/np.ndarray
        The reference frequency used in the analysis. Default is 20Hz
    model: str, optional
        The name of the NRSurrogate model you wish to use
    return_fits_used: Bool, optional
        if True, return the approximant that was used.
    properties: list, optional
        The list of properties you wish to calculate
    approximant: str, optional
        The approximant that was used to generate the posterior samples
    """
    from pesummary.gw.file.nrutils import NRSur_fit
    fit = NRSur_fit(*args, f_low=f_low, f_ref=f_ref, model=model, fits=
        properties, approximant=approximant)
    if return_fits_used:
        return fit, [model]
    return fit
