hijing/src/fjetrig.f

0001 C
0002 C
0003 C
0004         FUNCTION FJETRIG(X,WGT)
0005         IMPLICIT REAL*8(A-H,O-Z)
0006         REAL HIPR1(100),HINT1(100),PTMAX,PTMIN
0007         COMMON/HIPARNT/HIPR1,IHPR2(50),HINT1,IHNT2(50)
0008         SAVE  /HIPARNT/
0009         DIMENSION X(10)
0010         PTMIN=ABS(HIPR1(10))-0.25
0011         PTMIN=MAX(PTMIN,HIPR1(8))
0012         AM2=0.D0
0013         IF(IHPR2(3).EQ.3) THEN
0014            AM2=HIPR1(7)**2
0015            PTMIN=MAX(0.0,HIPR1(10))
0016         ENDIF
0017         PTMAX=ABS(HIPR1(10))+0.25
0018         IF(HIPR1(10).LE.0.0) PTMAX=HINT1(1)/2.0-AM2
0019         IF(PTMAX.LE.PTMIN) PTMAX=PTMIN+0.25
0020         PT2=(PTMAX**2-PTMIN**2)*X(1)+PTMIN**2
0021         AMT2=PT2+AM2
0022         XT=2.0*DSQRT(AMT2)/HINT1(1)
0023         YMX1=DLOG(1.0/XT+DSQRT(1.0/XT**2-1.0))
0024         Y1=2.0*YMX1*X(2)-YMX1
0025         YMX2=DLOG(2.0/XT-DEXP(Y1))
0026         YMN2=DLOG(2.0/XT-DEXP(-Y1))
0027         Y2=(YMX2+YMN2)*X(3)-YMN2
0028         IF(IHPR2(3).EQ.3) THEN
0029            GTRIG=2.0*GHVQ(Y1,Y2,AMT2)
0030         ELSE IF(IHPR2(3).EQ.2) THEN
0031            GTRIG=2.0*GPHOTON(Y1,Y2,PT2)
0032         ELSE
0033            GTRIG=G(Y1,Y2,PT2)
0034         ENDIF
0035         FJETRIG=2.0*YMX1*(YMX2+YMN2)*(PTMAX**2-PTMIN**2)
0036      &          *GTRIG/2.0
0037         RETURN
0038         END