def METHOD_NAME(pam,
                    fname_shm,
                    fname_dirs,
                    fname_values,
                    fname_indices,
                    fname_gfa,
                    reshape_dirs=False):
        """ Save SH, directions, indices and values of peaks to Nifti.
        """
        save_nifti(fname_shm, pam.shm_coeff.astype(np.float32), pam.affine)
        if reshape_dirs:
            pam_dirs = reshape_peaks_for_visualization(pam)
        else:
            pam_dirs = pam.peak_dirs.astype(np.float32)
        save_nifti(fname_dirs, pam_dirs, pam.affine)
        save_nifti(fname_values, pam.peak_values.astype(np.float32),
                   pam.affine)
        save_nifti(fname_indices, pam.peak_indices, pam.affine)
        save_nifti(fname_gfa, pam.gfa, pam.affine)