Double check the implementation of `calculate.cutoff` - should be fine but it's sort of unexplained in the paper/code