HIJING/hipyset1_35/luptdi_hijing.F

   1 * $Id$
   2
   3 C*********************************************************************
   4
   5       SUBROUTINE LUPTDI_HIJING(KFL,PX,PY)
   6
   7 C...Purpose: to generate transverse momentum according to a Gaussian.
   8 #include "ludat1_hijing.inc"
   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_HIJING(0))))
  13       IF(MSTJ(91).EQ.1) PT=PARJ(22)*PT
  14       IF(KFLA.EQ.0.AND.MSTJ(13).LE.0) PT=0.
  15       PHI=PARU(2)*RLU_HIJING(0)
  16       PX=PT*COS(PHI)
  17       PY=PT*SIN(PHI)
  18
  19       RETURN
  20       END