Add 90% contours to residual plots
At the moment PESummary can generate residual plots where the maxL waveform is plotted on top of the real gravitational wave strain (see here as an example). The function that generates this plot is here. However it would be better if it plotted the maxL waveform, with the 90% confidence region highlighted. This is already done in bilby
. See https://git.ligo.org/lscsoft/bilby/blob/master/bilby/gw/result.py#L243 for details.