This commit was generated by cvs2svn to compensate for changes in r1073,
[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / hiptdi.F
CommitLineData
e74335a4 1* $Id$
2C
3C
4C
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