anchor.diff:  
         if (cdb_isCrystal(icoll)) then
           p_in(j) = p
  +        s_in(j) = s_in(j) + s
         else
           p_in(j) = (one + dpop) * p0
  -      end if
  -
  -      s_in(j) = sp + (real(j_slices,fPrec)-one) * c_length
  +        s_in(j) = sp + (real(j_slices,fPrec)-one) * c_length
  +      end if
  +
       else
         x_in(j) = x
         y_in(j) = z
  
positive.code:       x  = x + sin_mb(tiltangle) * c_length
        xp = xp - tiltangle
      end if
  
      ! For selected collimator, first turn reset particle distribution to simple pencil beam
      nprim = 3
  
negative.code:     zlm   = -one * length
  
      if (cdb_isCrystal(icoll)) then
        call cry_doCrystal(ie,iturn,j,mat,x,xp,z,zp,s,p,x_in0,xp_in0,zlm,nhit,nabs,lhit_pos,lhit_turn,&
          part_abs_pos_local,part_abs_turn_local,impact,indiv,c_length)
  !      p = sqrt(p**2+pmap**2)
      else
  