Bug in constant_solenoid beam with angular momentum
By email, from Ao, confirmed by Rogers...
As I was going through a way to do tracking in G4BL, I also read your
I noticed in __birth_transverse_ellipse, for the constant_solenoid mode,
beta_4d = (1.+beam_def['normalised_angular_momentum'])**0.5/k_s
while from Kenn's note 71, it should be
beta_4d = (1.+beam_def['normalised_angular_momentum']***2*)**0.5/k_s
Was it because when you define the data card, you want people to include