hijing/src/fjet.f

0001 C
0002 C
0003 C
0004         FUNCTION FJET(X,WGT)
0005         IMPLICIT REAL*8(A-H,O-Z)
0006         REAL HIPR1(100),HINT1(100)
0007         COMMON/HIPARNT/HIPR1,IHPR2(50),HINT1,IHNT2(50)
0008         SAVE  /HIPARNT/
0009         DIMENSION X(10)
0010         PT2=(HINT1(1)**2/4.0-HIPR1(8)**2)*X(1)+HIPR1(8)**2
0011         XT=2.0*DSQRT(PT2)/HINT1(1)
0012         YMX1=DLOG(1.0/XT+DSQRT(1.0/XT**2-1.0))
0013         Y1=2.0*YMX1*X(2)-YMX1
0014         YMX2=DLOG(2.0/XT-DEXP(Y1))
0015         YMN2=DLOG(2.0/XT-DEXP(-Y1))
0016         Y2=(YMX2+YMN2)*X(3)-YMN2
0017
0018         FJET=2.0*YMX1*(YMX2+YMN2)*(HINT1(1)**2/4.0-HIPR1(8)**2)
0019      &          *G(Y1,Y2,PT2)/2.0
0020
0021         RETURN
0022         END