def METHOD_NAME(
    t: npt.ArrayLike,
    std: float = 1.0,
) -> Tuple[npt.ArrayLike, npt.ArrayLike, int]:
    r"""Gaussian wavelet
    Create a Gaussian wavelet given time axis ``t``
    and standard deviation ``std`` using
    :py:func:`scipy.signal.windows.gaussian`.
    Parameters
    ----------
    t : :obj:`numpy.ndarray`
        Time axis (positive part including zero sample)
    std : :obj:`float`, optional
        Standard deviation of gaussian
    Returns
    -------
    w : :obj:`numpy.ndarray`
        Wavelet
    t : :obj:`numpy.ndarray`
        Symmetric time axis
    wcenter : :obj:`int`
        Index of center of wavelet
    """
    t = _tcrop(t)
    t = np.concatenate((np.flipud(-t[1:]), t), axis=0)
    w = spgauss(len(t), std=std)
    wcenter = np.argmax(np.abs(w))
    return w, t, wcenter