anchor.diff:    integer, parameter :: dist_fmtYP          = 14 ! Vertical angle
     integer, parameter :: dist_fmtPX          = 15 ! Horizontal momentum
     integer, parameter :: dist_fmtPY          = 16 ! Vertical momentum
  -  integer, parameter :: dist_fmtSIGMA       = 17 ! Longitudinal relative position
  -  integer, parameter :: dist_fmtDT          = 18 ! Time delay
  -  integer, parameter :: dist_fmtE           = 19 ! Particle energy
  -  integer, parameter :: dist_fmtP           = 20 ! Particle momentum
  -  integer, parameter :: dist_fmtDEE0        = 21 ! Relative particle energy (to reference particle)
  -  integer, parameter :: dist_fmtDPP0        = 22 ! Relative particle momentum (to reference particle)
  +  integer, parameter :: dist_fmtPXP0        = 17 ! Relative horizontal momentum
  +  integer, parameter :: dist_fmtPYP0        = 18 ! Relative vertical momentum
  +  integer, parameter :: dist_fmtSIGMA       = 19 ! Longitudinal relative position
  +  integer, parameter :: dist_fmtDT          = 20 ! Time delay
  +  integer, parameter :: dist_fmtE           = 21 ! Particle energy
  +  integer, parameter :: dist_fmtP           = 22 ! Particle momentum
  +  integer, parameter :: dist_fmtDEE0        = 23 ! Relative particle energy (to reference particle)
  +  integer, parameter :: dist_fmtDPP0        = 24 ! Relative particle momentum (to reference particle)
   
     ! Normalised Coordinates
     integer, parameter :: dist_fmtX_NORM      = 31 ! Normalised horizontal position
  
positive.code:   mtc(1:napx)   = (nqq(1:napx)*nucm0)/(qq0*nucm(1:napx))
  
    ! If we have no energy arrays, we set all energies from deltaP = 0, that is reference momentum/energy
    call part_updatePartEnergy(dist_updtEFrom)
  
    if(dist_readIonA .and. dist_readIonZ .and. .not. dist_readPDGID) then
      do j=1,napx
  
negative.code:     call dist_appendFormat(dist_fmtDEE0, one)
    case("DP/P0","DPP0") ! Relative particle momentum (to reference particle)
      call dist_appendFormat(dist_fmtDPP0, one)
  
    case("MASS","M") ! Particle mass
      call dist_unitScale(fmtName, fmtUnit, 3, uFac, fErr)
  