anchor.diff:    real(kind=fPrec) xIn,xpIn,yIn,ypIn,xOut,xpOut,yOut,ypOut,sImp,sOut
     real(kind=fPrec) x_in0,xp_in0
   
  +  ! ien0,ien1: particle energy entering/leaving the collimator
  +  ! energy in MeV
  +  real(kind=fPrec)    :: ien0, ien1
  +  integer(kind=int16) :: nnuc0,nnuc1
  +
     ! Initilaisation
     mat    = cdb_cMaterialID(icoll)
     length = c_length
  
positive.code:   sRot   = sin_mb(c_rotation)
    cRRot  = cos_mb(-c_rotation)
    sRRot  = sin_mb(-c_rotation)
  
    do j=1,napx
  
  
negative.code:     end if
  
    end do ! End of loop over all particles
  
  ! write out energy change over this collimator
    if((ien0-ien1) > one) then
  #ifdef ROOT
      if(root_flag .and. root_Collimation == 1) then
        call root_EnergyDeposition(icoll, nnuc0-nnuc1,c1m3*(ien0-ien1))
      end if
  #endif
      write(unit208,"(2(i6,1x),e24.16)") icoll, (nnuc0-nnuc1), c1m3*(ien0-ien1)
      flush(unit208)
    end if
  
  end subroutine k2coll_collimate
  
  