anchor.diff:          TLdech1 = (const_dech/c2e2)*PC*(one-one/ratio)**2  ! Updated dechanneling length (m)
         end if
   
  -      Dechan = -log(rndm4())  ! probability of dechanneling
  -      Ldech  = TLdech1*Dechan ! actual dechan. length
  +      Dechan = -log_mb(rndm4()) ! Probability of dechanneling
  +      Ldech  = TLdech1*Dechan   ! Actual dechan. length
   
         ! careful: the dechanneling lentgh is along the trajectory
         ! of the particle -not along the longitudinal coordinate...
         if(ldech < l_chan) then
           iProc = proc_DC
  -        Dxp   = Ldech/Rcurv ! change angle from channeling [mrad]
  -        Sdech = Ldech*cos(miscut + half*Dxp)
  -        x     = x + Ldech*(sin(half*Dxp+miscut)) ! trajectory at channeling exit
  +        Dxp   = Ldech/Rcurv ! Change angle from channeling [mrad]
  +        Sdech = Ldech*cos_mb(miscut + half*Dxp)
  +        x     = x + Ldech*(sin_mb(half*Dxp+miscut)) ! Trajectory at channeling exit
           xp    = xp + Dxp + two*(rndm4()-half)*xpcrit
           y     = y + yp * Sdech
   
           call calc_ion_loss_cry(is,pc,ldech,dest)
  -        pc = pc - half*dest*Ldech ! energy loss to ionization while in CH [GeV]
  +        pc = pc - half*dest*Ldech ! Energy loss to ionization while in CH [GeV]
           x  = x  + half*(s_length-Sdech)*xp
           y  = y  + half*(s_length-Sdech)*yp
   
  
positive.code:   use floatPrecision
    use coll_crystal
    use coll_materials, only : zatom, exenergy, rho, anuc
  
    implicit none
  
  
negative.code:   ! New treatment of scattering routine based on standard sixtrack routine
    ! useful calculations for cross-section and event topology calculation
    ecmsq  = two*pmap*c1m3*pc
    xln15s = log_mb(0.15_fPrec*ecmsq)
  
    ! New models, see Claudia's thesis
    pptot = 0.041084_fPrec - 0.0023302_fPrec*log_mb(ecmsq) + 0.00031514_fPrec*log_mb(ecmsq)**2
    ppel  = (11.7_fPrec - 1.59_fPrec*log_mb(ecmsq) + 0.134_fPrec*log_mb(ecmsq)**2)/c1e3
    ppsd  = (4.3_fPrec + 0.3_fPrec*log_mb(ecmsq))/c1e3
    bpp   = 7.156_fPrec + 1.439_fPrec*log_mb(sqrt(ecmsq))
  
    ! Distribution for Ruth. scatt.
    tlow      = tlcut_cry
  