anchor.diff:  !++  Include collimator tilt
   ! this is propably not correct
   !
  -!             xp =  (xp_pencil0*cos_mb(c_rotation)+                         &
  -!     &            sin_mb(c_rotation)*yp_pencil0)
  +!             xp =  (xp_pencil0*cRot+                         &
  +!     &            sRot*yp_pencil0)
   !             if (tiltangle.gt.0.) then
   !                xp = xp - tiltangle
   !!             endif
  
positive.code: !     &                                   .and. mirror.lt.0.) then
  ! changed xp_pencil0(icoll) to xp_pencil0 due to definition mismatch
  ! this has to be solved if necassary and understood
  !                 c_tilt(2) = -1.*(xp_pencil0(icoll)*cos_mb(c_rotation)+    &
  !     &                       sin_mb(c_rotation)*yp_pencil0(icoll))
  !                 c_tilt(2) = -1.*(xp_pencil0*cos_mb(c_rotation)+           &
  !     &                       sin_mb(c_rotation)*yp_pencil0)
  !                 write(*,*) "INFO> Changed tilt2  ICOLL  to  ANGLE  ",  &
  !     &                   icoll, c_tilt(2), j
  !                 changed_tilt2(icoll) = .true.
  
negative.code: 
    real(kind=fPrec) x_flk,xp_flk,y_flk,yp_flk,zpj
    real(kind=fPrec) x_Dump,xpDump,y_Dump,ypDump,s_Dump
    real(kind=fPrec) cRot,sRot
  
    real(kind=fPrec) s_impact
    integer flagsec(npart)
  