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.