pyptdi.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pyptdi.f

 
 
C*********************************************************************
 
C...PYPTDI
C...Generates transverse momentum according to a Gaussian.
 
      SUBROUTINE pyptdi(KFL,PX,PY)
 
C...Double precision and integer declarations.
      IMPLICIT DOUBLE PRECISION(a-h, o-z)
      IMPLICIT INTEGER(i-n)
      INTEGER pyk,pychge,pycomp
C...Commonblocks.
      common/pydat1/mstu(200),paru(200),mstj(200),parj(200)
      SAVE /pydat1/
 
C...Generate p_T and azimuthal angle, gives p_x and p_y.
      kfla=iabs(kfl)
      pt=parj(21)*sqrt(-log(max(1d-10,pyr(0))))
      IF(parj(23).GT.pyr(0)) pt=parj(24)*pt
      IF(mstj(91).EQ.1) pt=parj(22)*pt
      IF(kfla.EQ.0.AND.mstj(13).LE.0) pt=0d0
      phi=paru(2)*pyr(0)
      px=pt*cos(phi)
      py=pt*sin(phi)
 
      RETURN
      END