def remove_straddlers(events, time, s_freq, toler=0.1):
    """Reject an event if it straddles a cision point, by comparing its
    duration to its timespan

    Parameters
    ----------
    events : ndarray (dtype='int')
        N x 3 matrix with start, peak, end samples
    time : ndarray (dtype='float')
        vector with time points
    s_freq : float
        sampling frequency
    toler : float, def=0.1
        maximum tolerated difference between event duration and timespan

    Returns
    -------
    ndarray (dtype='int')
        N x 3 matrix with start , peak, end samples
    """
    duration = (events[:, 2] - 1 - events[:, 0]) / s_freq
    continuous = time[events[:, 2] - 1] - time[events[:, 0]] - duration < toler

    return events[continuous, :]