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