Factor out rho limits, increase errorbar zorder
Factor out sep calcs into get_seps(), use full superset matrices more
Add psthcorrdiff