Evaluate "damping" of higher-order transfer function - something similar to settling time without the discontinuities

I have an 11th order transfer function in the continuous domain (simplified linearized model of a grid-forming converter) and I’m looking for a way to characterize it with a single number for evaluation of a parameter sweep. I’ve so far been using settling time because that’s a typical requirement in grid codes, but there are discontinuities in a parameter sweep (when a peak moves past the threshold the function jumps). That for example makes the use of search algorithms difficult and also makes it harder to define an optimum over all operating points (the location of the discontinuities of course varies depending on the operating point). There also wouldn’t really be a point in linearizing if I end up evaluating something that I can extract from a nonlinear model 🙂

So I’m looking for something that is similar to settling time, but without the discontinuities. I have several pole pairs in close proximity to each other and to zeros, so simply evaluating the damping ratio of the dominant pole pair won’t do. I know that I can weigh pole pairs by frequency, but that would still neglect the zeros. I need to somehow incorporate the zeros as well.

I hope my problem is comprehensible. Any hint in the right direction is appreciated!