Dodgy Hack Virtual Planes. Have written an algorithm directly into the stepping action Takes either end of step, tests to see if it straddles a virtual plane If yes, performs linear interpolation over the step Rubbish hack, but I need data. Results - 5.5pi emittance.
Turquoise - phase space emittance
Pink - unnormalised trace space emittance (I think)
Results - equilibrium emittance & full magnetic lattice(a) no absorbers or electrostatics (b) electrostatic field but no absorbers (c) absorbers but no electrostatics (d) absorbers & electrostatics