IMPLICIT DOUBLE PRECISION(A-H, O-Z)
IMPLICIT INTEGER(I-N)
INTEGER PYK,PYCHGE,PYCOMP
- external plthik, pln, plt, pls, gauss, gluang
+ external plthik, pln, plt, pls, pygauss, gluang
external pyp,pyr,pyk
common /plpar1/ tau0,T0,TC,sigqq,AW,RA,mvisc,nf
common /plpar2/ pln0,taupl,tauh,sigpl,sigh,sigplh,sigqqh,rg,rgn
return
end
-* function to generate gauss distribution
- double precision function gauss(x0,sig)
- IMPLICIT DOUBLE PRECISION(A-H, O-Z)
- IMPLICIT INTEGER(I-N)
- 41 u1=pyr(0)
- u2=pyr(0)
- v1=2.d0*u1-1.d0
- v2=2.d0*u2-1.d0
- s=v1**2+v2**2
- if(s.gt.1) go to 41
- gauss=v1*dsqrt(-2.d0*dlog(s)/s)*sig+x0
- return
- end
-
* function to calculate angular distribution of emitted gluons
double precision function gluang(x)
IMPLICIT DOUBLE PRECISION(A-H, O-Z)