HIJING/hijing1_36/fjet.F

   1 * $Id$
   2 C
   3 C
   4 C
   5         FUNCTION FJET(X,WGT)
   6         IMPLICIT REAL*8(A-H,O-Z)
   7 #define BLANKET_SAVE
   8 #include "hiparnt.inc"
   9         DIMENSION X(10)
  10         SAVE
  11         PT2=(HINT1(1)**2/4.0-HIPR1(8)**2)*X(1)+HIPR1(8)**2
  12         XT=2.0*DSQRT(PT2)/HINT1(1)
  13         IF (XT.GT.1.0) XT=1.0
  14         YMX1=DLOG(1.0/XT+DSQRT(1.0/XT**2-1.0))
  15         Y1=2.0*YMX1*X(2)-YMX1
  16         YMX2=DLOG(2.0/XT-DEXP(Y1))
  17         YMN2=DLOG(2.0/XT-DEXP(-Y1))
  18         Y2=(YMX2+YMN2)*X(3)-YMN2
  19         FJET=2.0*YMX1*(YMX2+YMN2)*(HINT1(1)**2/4.0-HIPR1(8)**2)
  20      &          *G(Y1,Y2,PT2)/2.0
  21         RETURN
  22         END