PYTHIA/jetset/luptdi.F

   1
   2 C*********************************************************************
   3
   4       SUBROUTINE LUPTDI(KFL,PX,PY)
   5
   6 C...Purpose: to generate transverse momentum according to a Gaussian.
   7       COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
   8       SAVE /LUDAT1/
   9
  10 C...Generate p_T and azimuthal angle, gives p_x and p_y.
  11       KFLA=IABS(KFL)
  12       PT=PARJ(21)*SQRT(-LOG(MAX(1E-10,RLU(0))))
  13       IF(PARJ(23).GT.RLU(0)) PT=PARJ(24)*PT
  14       IF(MSTJ(91).EQ.1) PT=PARJ(22)*PT
  15       IF(KFLA.EQ.0.AND.MSTJ(13).LE.0) PT=0.
  16       PHI=PARU(2)*RLU(0)
  17       PX=PT*COS(PHI)
  18       PY=PT*SIN(PHI)
  19
  20       RETURN
  21       END