]>
Commit | Line | Data |
---|---|---|
e74335a4 | 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 |