anchor.diff:  !   but it might be then that only one jaw is hit on the first turn, thus only by half of the particles
   !   the particle generated on the other side will then hit the same jaw several turns later, possibly smearing the impact parameter
   !   This could possibly be improved in the future.
  -    call makedis_coll(napx,myalphax,myalphay, mybetax, mybetay, myemitx0_collgap, myemity0_collgap, &
  +    call makedis_coll(myalphax,myalphay, mybetax, mybetay, myemitx0_collgap, myemity0_collgap, &
    &                    myenom, mynex2, mdex, myney2, mdey, myx, myxp, myy, myyp, myp, mys)
   
       do j = 1, napx
  
positive.code: !!     YIL EDIT 2010: particle 0 is always on orbit...
  !!
  !! @author Javier Barranco <jbarranc@cern.ch>
  !! @param mynp
  !! @param myalphax
  !! @param myalphay
  !! @param mybetax
  
negative.code: !     centred in the aperture centre are generated. (SR, 08-05-2005)
  
    use crcoall
    use mod_common, only : napx
    implicit none
  
    integer :: j
    real(kind=fPrec), allocatable :: myx(:) !(npart)
    real(kind=fPrec), allocatable :: myxp(:) !(npart)
    real(kind=fPrec), allocatable :: myy(:) !(npart)
    real(kind=fPrec), allocatable :: myyp(:) !(npart)
    real(kind=fPrec), allocatable :: myp(:) !(npart)
    real(kind=fPrec), allocatable :: mys(:) !(npart)
  
    real(kind=fPrec) myalphax,mybetax,myemitx0,myemitx,mynex,mdex, &
    &mygammax,myalphay,mybetay,myemity0,myemity,myney,mdey,mygammay,   &
  