For the mirror plot (e.g. here), I've expected a plot like this:
I think Orel has a script to do this.