3 C*********************************************************************
5 SUBROUTINE LUPTDI_HIJING(KFL,PX,PY)
7 C...Purpose: to generate transverse momentum according to a Gaussian.
8 #include "ludat1_hijing.inc"
10 C...Generate p_T and azimuthal angle, gives p_x and p_y.
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)