]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | FUNCTION FJET(X,WGT) | |
6 | IMPLICIT REAL*8(A-H,O-Z) | |
bc676b8e | 7 | #define BLANKET_SAVE |
e74335a4 | 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) | |
c87d167b | 13 | IF (XT.GT.1.0) XT=1.0 |
e74335a4 | 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 |