Currently the maximum value of ASq (used for the accept/reject method of toy generation) for a given DP model has to be determined by trial and error and then set by hand.
It should be possible to determine it via minimising the `-ASq` over the DP coordinates.
If so, add a method to do so.
D85 implements this for `LauIsobarDynamics`
Should next be done for more complex cases such as `LauTimeDepFitModel`.