Using SomeTF2->SetNpx and SomeTF2->SetNpy solves the problem for these parameters, but the integral takes much longer. And I also find it hard to trust the result. So my other questions remain:
Is there another integrator which should work better? (or/and faster?)
Is there a better way of implementing non-square integration limits?