I'm actually a bit annoyed that there's no way to ask `numpy.random.choice` for multiple samples without replacement (i.e. sample 5 from 10, 5 times). But it doesn't really apply here because the poisson is needed first. But the poisson sampling can be vectorised a bit better.