anchor.diff:    pc  = pc - dei*dz ! Energy lost because of ionization process[GeV]
   
     dya   = (13.6_fPrec/pc)*sqrt(dz/dlr) ! RMS of coloumb scattering MCS (mrad)
  -  kxmcs = dya*ran_gauss(one)
  -  kymcs = dya*ran_gauss(one)
  +  kxmcs = dya*ran_gauss(zero)
  +  kymcs = dya*ran_gauss(zero)
   
     xp = xp+kxmcs
     yp = yp+kymcs
  
positive.code:         if(Ldech < Length-Lrefl) then
  
            iProc = proc_DC
            Dxp   = Ldech/c_rcurv + (half*ran_gauss(one))*xpcrit
            x     = x + Ldech*(sin_mb(half*Dxp+xp)) ! Trajectory at channeling exit
            y     = y + Sdech*yp
            xp    =  Dxp
  
negative.code:       teta = sqrt(t)/pc
      end if
  
      tx = (teta*ran_gauss(zero))*c1e3
      tz = (teta*ran_gauss(zero))*c1e3
  
      ! Change p angle
      xp = xp + tx
  