anchor.diff:        n_chan  = 0          !valentina :initialize to zero the counters for crystal effects
         n_VR    = 0          !
         n_amorphous = 0      !
  +
  +      cRot  = cos(c_rotation)
  +      sRot  = sin(c_rotation)
  +      cMRot = cos(-c_rotation)
  +      sMRot = sin(-c_rotation)
   
         do j = 1, nev
   !
  
positive.code: !
  !++  Last do rotation into collimator frame
  !
           X_IN(J)  = X  *COS(-1.*C_ROTATION) + Z  *SIN(-1.*C_ROTATION)
           Y_IN(J)  = Z  *COS(-1.*C_ROTATION) - X  *SIN(-1.*C_ROTATION)
           XP_IN(J) = XP *COS(-1.*C_ROTATION) + ZP *SIN(-1.*C_ROTATION)
           YP_IN(J) = ZP *COS(-1.*C_ROTATION) - XP *SIN(-1.*C_ROTATION)
  
  !----- other pencil beam stuff-------
           IF ( ICOLL.EQ.IPENCIL) then
             X00  = MIRROR * X00
             X_IN(J)  = X00  *COS(-1.*C_ROTATION) + Z00  *SIN(-1.*C_ROTATION)
             Y_IN(J)  = Z00  *COS(-1.*C_ROTATION) - X00  *SIN(-1.*C_ROTATION)
  !
             XP_IN(J) = XP_IN(J) + MIRROR*XP_PENCIL0
             YP_IN(J) = YP_IN(J) + MIRROR*YP_PENCIL0
  
negative.code: 
             IF (.NOT. CHANGED_TILT1(ICOLL) .AND. MIRROR.GT.0.) THEN
                     WRITE (*,*) 'NEVER!!!'
                     C_TILT(1) = XP_PENCIL0*cRot+sRot*YP_PENCIL0
                     WRITE(*,*) 'INFO> Changed tilt1  ICOLL  to  ANGLE  ',ICOLL, C_TILT(1)
  !
                     CHANGED_TILT1(ICOLL) = .true.
             ELSEIF (.NOT. CHANGED_TILT2(ICOLL) .AND. MIRROR.LT.0.) THEN
                     C_TILT(2) = -1.*(XP_PENCIL0*cRot+sRot*YP_PENCIL0)
                     WRITE(*,*) 'INFO> Changed tilt2  ICOLL  to  ANGLE  ',ICOLL, C_TILT(2)
  !
                     CHANGED_TILT2(ICOLL) = .true.
  