Additional protection
[u/mrichter/AliRoot.git] / 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