Tracking Spin in PTC

For the full code, please see: Sec. (10.9).

The first thing to make sure is that you have the right kink of particle. The fibre parameter fibre%AG ( which is $g-2$) must have the value associated with the kind of particle you are using. This is automatic for electron, protons and muons but one can always print fibre%AG of any fibre to double check. PTC has the following $g-2$ constants:

  real(dp),parameter::A_ELECTRON=1.15965218111e-3_dp  !frs NIST CODATA 2006
  real(dp),parameter::A_MUON=1.16592069e-3_dp               !frs NIST CODATA 2006
  real(dp),parameter::A_PROTON=1.79284735e-0_dp          !frs (approx) NIST CODATA 2006