]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | FUNCTION FJETRIG(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 | PTMIN=ABS(HIPR1(10))-0.25 | |
70090d58 | 12 | PTMIN=MAX(PTMIN,DBLE(HIPR1(8))) |
e74335a4 | 13 | AM2=0.D0 |
14 | IF(IHPR2(3).EQ.3) THEN | |
15 | AM2=HIPR1(7)**2 | |
16 | PTMIN=MAX(0.0,HIPR1(10)) | |
17 | ENDIF | |
18 | PTMAX=ABS(HIPR1(10))+0.25 | |
19 | IF(HIPR1(10).LE.0.0) PTMAX=HINT1(1)/2.0-AM2 | |
20 | IF(PTMAX.LE.PTMIN) PTMAX=PTMIN+0.25 | |
21 | PT2=(PTMAX**2-PTMIN**2)*X(1)+PTMIN**2 | |
22 | AMT2=PT2+AM2 | |
23 | XT=2.0*DSQRT(AMT2)/HINT1(1) | |
24 | YMX1=DLOG(1.0/XT+DSQRT(1.0/XT**2-1.0)) | |
25 | Y1=2.0*YMX1*X(2)-YMX1 | |
26 | YMX2=DLOG(2.0/XT-DEXP(Y1)) | |
27 | YMN2=DLOG(2.0/XT-DEXP(-Y1)) | |
28 | Y2=(YMX2+YMN2)*X(3)-YMN2 | |
29 | IF(IHPR2(3).EQ.3) THEN | |
30 | GTRIG=2.0*GHVQ(Y1,Y2,AMT2) | |
31 | ELSE IF(IHPR2(3).EQ.2) THEN | |
32 | GTRIG=2.0*GPHOTON(Y1,Y2,PT2) | |
33 | ELSE | |
34 | GTRIG=G(Y1,Y2,PT2) | |
35 | ENDIF | |
36 | FJETRIG=2.0*YMX1*(YMX2+YMN2)*(PTMAX**2-PTMIN**2) | |
37 | & *GTRIG/2.0 | |
38 | RETURN | |
39 | END |