"no dynamic range"
When running pesummary
on a simple example where only 4 parameters are sampled over, I get the classic no dynamic range error. I believe this is because corner is being called where some parameters are just a single repeated value.
I see that in some places you are jittering to avoid this. An alternative is to add
default_kwargs['range'] = default_kwargs.get('range', [1] * len(corner_parameters))
Which ensures the corner plots generate and just show a single spike.
Of course its up to you how to solve it, but currently I can't get pesummary to run on non-full-PE examples without that.