anchor.diff:    theta    = 13.6e-3_fPrec/(p0*(one+dpop))
     rad_len  = radl(mat)
   
  -  x  = (x/theta)/radl(mat)
  -  xp = xp/theta
  -  z  = (z/theta)/radl(mat)
  -  zp = zp/theta
  +  x     = (x/theta)/radl(mat)
  +  xp    = xp/theta
  +  z     = (z/theta)/radl(mat)
  +  zp    = zp/theta
     rlen0 = zlm1/radl(mat)
     rlen  = rlen0
   
  
positive.code:   vb  = rndm4()
    va2 = va**2
    vb2 = vb**2
    r2 = va2 + vb2
    if(r2 > one) goto 10
    tx = (teta*((two*va)*vb))/r2
    tz = (teta*(va2 - vb2))/r2
  
  end subroutine k2coll_tetat
  
  
negative.code: #endif
  
    ! Claudia new fit for the slope parameter with new data at sqrt(s)=7 TeV from TOTEM
    bpp = 7.156_fPrec + 1.439_fPrec*log_mb(sqrt(ecmsq))
  
    ! unmeasured tungsten data,computed with lead data and power laws
    bnref(4) = (bnref(5)*(anuc(4))/anuc(5))**(two/three)
    emr(4)   = (emr(5)  *(anuc(4))/anuc(5))**(one/three)
  
    ! Compute cross-sections (CS) and probabilities + Interaction length
    ! Last two material treated below statement number 100
  