]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | C | |
6 | FUNCTION G(Y1,Y2,PT2) | |
7 | IMPLICIT REAL*8 (A-H,O-Z) | |
bc676b8e | 8 | #define BLANKET_SAVE |
e74335a4 | 9 | #include "hiparnt.inc" |
10 | DIMENSION F(2,7) | |
11 | SAVE | |
12 | XT=2.0*DSQRT(PT2)/HINT1(1) | |
13 | X1=0.50*XT*(DEXP(Y1)+DEXP(Y2)) | |
14 | X2=0.50*XT*(DEXP(-Y1)+DEXP(-Y2)) | |
9f320b1f | 15 | DUMMY = 1.D0-XT**2/X1/X2 |
16 | if (DUMMY.LT.0.D0) DUMMY = 0.D0 | |
17 | Z=DSQRT(DUMMY) | |
e74335a4 | 18 | SS=X1*X2*HINT1(1)**2 |
19 | T=-(1.0-Z)/2.0 | |
20 | U=-(1.0+Z)/2.0 | |
21 | AF=3.0 | |
22 | DLAM=HIPR1(15) | |
23 | APH=12.0*3.1415926/(33.0-2.0*AF)/DLOG(PT2/DLAM**2) | |
24 | C | |
25 | CALL PARTON(F,X1,X2,PT2) | |
26 | C | |
27 | G11=( (F(1,1)+F(1,2))*(F(2,3)+F(2,4)+F(2,5)+F(2,6)) | |
28 | & +(F(1,3)+F(1,4))*(F(2,5)+F(2,6)) )*SUBCRS1(T,U) | |
29 | C | |
30 | G12=( (F(2,1)+F(2,2))*(F(1,3)+F(1,4)+F(1,5)+F(1,6)) | |
31 | & +(F(2,3)+F(2,4))*(F(1,5)+F(1,6)) )*SUBCRS1(U,T) | |
32 | C | |
33 | G13=(F(1,1)*F(2,1)+F(1,2)*F(2,2)+F(1,3)*F(2,3)+F(1,4)*F(2,4) | |
34 | & +F(1,5)*F(2,5)+F(1,6)*F(2,6))*(SUBCRS1(U,T) | |
35 | & +SUBCRS1(T,U)-8.D0/T/U/27.D0) | |
36 | C | |
37 | G2=(AF-1)*(F(1,1)*F(2,2)+F(2,1)*F(1,2)+F(1,3)*F(2,4) | |
38 | & +F(2,3)*F(1,4)+F(1,5)*F(2,6)+F(2,5)*F(1,6))*SUBCRS2(T,U) | |
39 | C | |
40 | G31=(F(1,1)*F(2,2)+F(1,3)*F(2,4)+F(1,5)*F(2,6))*SUBCRS3(T,U) | |
41 | G32=(F(2,1)*F(1,2)+F(2,3)*F(1,4)+F(2,5)*F(1,6))*SUBCRS3(U,T) | |
42 | C | |
43 | G4=(F(1,1)*F(2,2)+F(2,1)*F(1,2)+F(1,3)*F(2,4)+F(2,3)*F(1,4)+ | |
44 | 1 F(1,5)*F(2,6)+F(2,5)*F(1,6))*SUBCRS4(T,U) | |
45 | C | |
46 | G5=AF*F(1,7)*F(2,7)*SUBCRS5(T,U) | |
47 | C | |
48 | G61=F(1,7)*(F(2,1)+F(2,2)+F(2,3)+F(2,4)+F(2,5) | |
49 | & +F(2,6))*SUBCRS6(T,U) | |
50 | G62=F(2,7)*(F(1,1)+F(1,2)+F(1,3)+F(1,4)+F(1,5) | |
51 | & +F(1,6))*SUBCRS6(U,T) | |
52 | C | |
53 | G7=F(1,7)*F(2,7)*SUBCRS7(T,U) | |
54 | C | |
55 | G=(G11+G12+G13+G2+G31+G32+G4+G5+G61+G62+G7)*HIPR1(17)* | |
56 | 1 3.14159D0*APH**2/SS**2 | |
57 | RETURN | |
58 | END |