hijing/src/hiptdi.f

0001 C
0002 C
0003 C
0004         SUBROUTINE HIPTDI(PT,PTMAX,IOPT)
0005         COMMON/HIPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)
0006         SAVE  /HIPARNT/
0007         COMMON/RANSEED/NSEED
0008         SAVE  /RANSEED/
0009         IF(IOPT.EQ.2) THEN
0010                 PT=HIRND2(7,0.0,PTMAX)
0011                 IF(PT.GT.HIPR1(8))
0012      &          PT=HIPR1(2)*SQRT(-ALOG(EXP(-HIPR1(8)**2/HIPR1(2)**2)
0013      &                  -ATL_RAN(NSEED)*(EXP(-HIPR1(8)**2/HIPR1(2)**2)-
0014      &                  EXP(-PTMAX**2/HIPR1(2)**2))))
0015
0016         ELSE
0017                 PT=HIPR1(2)*SQRT(-ALOG(1.0-ATL_RAN(NSEED)*
0018      &                  (1.0-EXP(-PTMAX**2/HIPR1(2)**2))))
0019         ENDIF
0020         PTMAX0=MAX(PTMAX,0.01)
0021         PT=MIN(PTMAX0-0.01,PT)
0022         RETURN
0023         END