def final_spin_of_merger_from_NR(*args, NRfit='average', return_fits_used=False
    ):
    """Return the final spin resulting from a BBH merger using NR fits

    Parameters
    ----------
    NRfit: str
        Name of the fit you wish to use. If you wish to use a precessing fit
        please use the syntax 'precessing_{}'.format(fit_name). If you wish
        to have an average NR fit, then pass 'average'
    return_fits_used: Bool, optional
        if True, return the fits that were used. Only used when NRfit='average'
    """
    from pesummary.gw.file import nrutils
    if NRfit.lower() == 'average':
        func = getattr(nrutils, 'bbh_final_spin_average_precessing')
    elif 'pan' in NRfit.lower():
        func = getattr(nrutils, 'bbh_final_spin_non_spinning_Panetal')
    elif 'precessing' in NRfit.lower():
        func = getattr(nrutils, 'bbh_final_spin_precessing_{}'.format(NRfit
            .split('precessing_')[1]))
    else:
        func = getattr(nrutils, 'bbh_final_spin_non_precessing_{}'.format(
            NRfit))
    if NRfit.lower() == 'average':
        return func(*args, return_fits_used=return_fits_used)
    return func(*args)
