(also `events[np.in1d(events[:, 2], [1, 3])]` is a quite complicated line of numpy ;))