]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | SUBROUTINE HIPTDI(PT,PTMAX,IOPT) | |
6 | #include "hiparnt.inc" | |
7 | SAVE | |
8 | IF(IOPT.EQ.2) THEN | |
9 | PT=HIRND2(7,0.0,PTMAX) | |
10 | IF(PT.GT.HIPR1(8)) | |
11 | & PT=HIPR1(2)*SQRT(-ALOG(EXP(-HIPR1(8)**2/HIPR1(2)**2) | |
12 | & -RLU_HIJING(0)*(EXP(-HIPR1(8)**2/HIPR1(2)**2)- | |
13 | & EXP(-PTMAX**2/HIPR1(2)**2)))) | |
14 | ||
15 | ELSE | |
16 | PT=HIPR1(2)*SQRT(-ALOG(1.0-RLU_HIJING(0)* | |
17 | & (1.0-EXP(-PTMAX**2/HIPR1(2)**2)))) | |
18 | ENDIF | |
19 | PTMAX0=MAX(PTMAX,0.01) | |
20 | PT=MIN(PTMAX0-0.01,PT) | |
21 | RETURN | |
22 | END |