]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | REAL*8 FUNCTION ZJJ7(P1, P2, P3, P4, P5,NHEL,IM1) | |
3 | C | |
4 | C FUNCTION GENERATED BY MADGRAPH | |
5 | C RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS | |
6 | C FOR THE POINT IN PHASE SPACE P1,P2,P3,P4,... | |
7 | C AND HELICITY NHEL(1),NHEL(2),.... | |
8 | C | |
9 | C FOR PROCESS : g u -> z g u | |
10 | C | |
11 | IMPLICIT NONE | |
12 | C | |
13 | C CONSTANTS | |
14 | C | |
15 | INTEGER NGRAPHS, NEIGEN, NEXTERNAL | |
16 | PARAMETER (NGRAPHS= 8,NEIGEN= 2,NEXTERNAL=5) | |
17 | REAL*8 ZERO | |
18 | PARAMETER (ZERO=0D0) | |
19 | C | |
20 | C ARGUMENTS | |
21 | C | |
22 | REAL*8 P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3) | |
23 | INTEGER NHEL(NEXTERNAL) | |
24 | INTEGER IM1 | |
25 | C | |
26 | C LOCAL VARIABLES | |
27 | C | |
28 | INTEGER I,J | |
29 | REAL*8 EIGEN_VAL(NEIGEN), EIGEN_VEC(NGRAPHS,NEIGEN) | |
30 | COMPLEX*16 ZTEMP | |
31 | COMPLEX*16 AMP(NGRAPHS) | |
32 | COMPLEX*16 W1(6) , W2(6) , W3(6) , W4(6) , W5(6) | |
33 | COMPLEX*16 W6(6) , W7(6) , W8(6) , W9(6) , W10(6) | |
34 | COMPLEX*16 W11(6) , W12(6) , W13(6) , W14(6) , W15(6) | |
35 | REAL*8 GZUD(2) | |
36 | C | |
37 | C GLOBAL VARIABLES | |
38 | C | |
39 | #include "isajet/mgcoms.inc" | |
40 | C | |
41 | C COLOR DATA | |
42 | C | |
43 | DATA EIGEN_VAL(1 )/ 5.8333333333333270D-01 / | |
44 | DATA EIGEN_VEC(1 ,1 )/ 4.0824829046386285D-01 / | |
45 | DATA EIGEN_VEC(2 ,1 )/ 0.0000000000000000D+00 / | |
46 | DATA EIGEN_VEC(3 ,1 )/ 4.0824829046386318D-01 / | |
47 | DATA EIGEN_VEC(4 ,1 )/ 0.0000000000000000D+00 / | |
48 | DATA EIGEN_VEC(5 ,1 )/ 4.0824829046386285D-01 / | |
49 | DATA EIGEN_VEC(6 ,1 )/ 4.0824829046386285D-01 / | |
50 | DATA EIGEN_VEC(7 ,1 )/ 4.0824829046386318D-01 / | |
51 | DATA EIGEN_VEC(8 ,1 )/ 4.0824829046386318D-01 / | |
52 | DATA EIGEN_VAL(2 )/ 1.7499999999999991D+00 / | |
53 | DATA EIGEN_VEC(1 ,2 )/ 2.6726124191242445D-01 / | |
54 | DATA EIGEN_VEC(2 ,2 )/ -5.3452248382484879D-01 / | |
55 | DATA EIGEN_VEC(3 ,2 )/ -2.6726124191242434D-01 / | |
56 | DATA EIGEN_VEC(4 ,2 )/ -5.3452248382484879D-01 / | |
57 | DATA EIGEN_VEC(5 ,2 )/ 2.6726124191242445D-01 / | |
58 | DATA EIGEN_VEC(6 ,2 )/ 2.6726124191242445D-01 / | |
59 | DATA EIGEN_VEC(7 ,2 )/ -2.6726124191242434D-01 / | |
60 | DATA EIGEN_VEC(8 ,2 )/ -2.6726124191242434D-01 / | |
61 | C ---------- | |
62 | C BEGIN CODE | |
63 | C ---------- | |
64 | IF(IM1.EQ.3.OR.IM1.EQ.7.OR.IM1.EQ.11) THEN | |
65 | GZUD(1)=GZD(1) | |
66 | GZUD(2)=GZD(2) | |
67 | ELSE | |
68 | GZUD(1)=GZU(1) | |
69 | GZUD(2)=GZU(2) | |
70 | ENDIF | |
71 | C | |
72 | CALL VXXXXX(P1 , ZERO,NHEL(1 ),-1,W1 ) | |
73 | CALL IXXXXX(P2 ,FMASS(IM1),NHEL(2 ), 1,W2 ) | |
74 | CALL VXXXXX(P3 ,ZMASS,NHEL(3 ), 1,W3 ) | |
75 | CALL VXXXXX(P4 , ZERO,NHEL(4 ), 1,W4 ) | |
76 | CALL OXXXXX(P5 ,FMASS(IM1),NHEL(5 ), 1,W5 ) | |
77 | CALL FVIXXX(W2 ,W3 ,GZUD,FMASS(IM1),FWIDTH(IM1),W6 ) | |
78 | CALL FVIXXX(W6 ,W4 ,GG,FMASS(IM1),FWIDTH(IM1),W7 ) | |
79 | CALL IOVXXX(W7 ,W5 ,W1 ,GG,AMP(1 )) | |
80 | CALL JGGXXX(W4 ,W1 ,G,W8 ) | |
81 | CALL FVOXXX(W5 ,W3 ,GZUD,FMASS(IM1),FWIDTH(IM1),W9 ) | |
82 | CALL IOVXXX(W2 ,W9 ,W8 ,GG,AMP(2 )) | |
83 | CALL FVIXXX(W6 ,W1 ,GG,FMASS(IM1),FWIDTH(IM1),W10 ) | |
84 | CALL IOVXXX(W10 ,W5 ,W4 ,GG,AMP(3 )) | |
85 | CALL IOVXXX(W6 ,W5 ,W8 ,GG,AMP(4 )) | |
86 | CALL FVIXXX(W2 ,W4 ,GG,FMASS(IM1),FWIDTH(IM1),W11 ) | |
87 | CALL FVOXXX(W5 ,W1 ,GG,FMASS(IM1),FWIDTH(IM1),W12 ) | |
88 | CALL IOVXXX(W11 ,W12 ,W3 ,GZUD,AMP(5 )) | |
89 | CALL IOVXXX(W11 ,W9 ,W1 ,GG,AMP(6 )) | |
90 | CALL FVIXXX(W2 ,W1 ,GG,FMASS(IM1),FWIDTH(IM1),W13 ) | |
91 | CALL FVIXXX(W13 ,W4 ,GG,FMASS(IM1),FWIDTH(IM1),W14 ) | |
92 | CALL IOVXXX(W14 ,W5 ,W3 ,GZUD,AMP(7 )) | |
93 | CALL FVIXXX(W13 ,W3 ,GZUD,FMASS(IM1),FWIDTH(IM1),W15 ) | |
94 | CALL IOVXXX(W15 ,W5 ,W4 ,GG,AMP(8 )) | |
95 | ZJJ7 = 0.D0 | |
96 | DO I = 1, NEIGEN | |
97 | ZTEMP = (0.D0,0.D0) | |
98 | DO J = 1, NGRAPHS | |
99 | ZTEMP = ZTEMP + EIGEN_VEC(J,I)*AMP(J) | |
100 | ENDDO | |
101 | ZJJ7 =ZJJ7+ZTEMP*EIGEN_VAL(I)*CONJG(ZTEMP) | |
102 | ENDDO | |
103 | C CALL GAUGECHECK(AMP,ZTEMP,EIGEN_VEC,EIGEN_VAL,NGRAPHS,NEIGEN) | |
104 | END |