]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | SUBROUTINE XZZWW | |
3 | C | |
4 | C SET UP Z0 Z0 -> W+ W- AMPLITUDES AS RATIONAL FUNCTIONS OF Z | |
5 | C | |
6 | C RE(F(Z,L)) = SUM(I,J)(ANWWWW(I+1,J,L)*Z**I | |
7 | C /(ADWWWW(1,J)+ADWWWW(2,J)*Z)) | |
8 | C IM(F(Z,L)) = AIWWWW(L) (INDEPENDENT OF Z) | |
9 | C J LABELS PIECES WITH SAME DENOMINATOR. | |
10 | C L=1 FOR 0,0; L=2 FOR 1,-1; L=3 FOR 1,1; L=4 FOR 0,1 | |
11 | C | |
12 | C *NOTE* A FACTOR OF SIN(THETA)/SQRT(2) IS REMOVED FROM F01 | |
13 | C | |
14 | #include "isajet/const.inc" | |
15 | #include "isajet/jetpar.inc" | |
16 | #include "isajet/wcon.inc" | |
17 | #include "isajet/hcon.inc" | |
18 | #if defined(CERNLIB_DOUBLE) | |
19 | DOUBLE PRECISION WM,ZM,ZM2,ZM3,ZM4,ZM5,ZM6,HM,HM2,HM3,HM4,HG,HG2 | |
20 | $,PROPH,RTS,S,S2,S3,SW,QQ0,QQI,QQF | |
21 | #endif | |
22 | C | |
23 | C USE UNITS OF WM TO AVOID LARGE NUMBERS - NOTE ANWWWW/ADWWWW | |
24 | C AND AIWWWW ARE DIMENSIONLESS | |
25 | WM=WMASS(2) | |
26 | ZM=WMASS(4)/WM | |
27 | ZM2=ZM**2 | |
28 | ZM3=ZM**3 | |
29 | ZM4=ZM**4 | |
30 | ZM5=ZM**5 | |
31 | ZM6=ZM**6 | |
32 | HM=HMASS/WM | |
33 | HM2=HM**2 | |
34 | HM3=HM**3 | |
35 | HM4=HM**4 | |
36 | HG=HGAM/WM | |
37 | HG2=HG**2 | |
38 | RTS=QMW/WM | |
39 | S=RTS**2 | |
40 | S2=S**2 | |
41 | S3=S**3 | |
42 | PROPH=(S-HM2)**2+(HM*HG)**2 | |
43 | C CORRECT SIGN OF HIGGS AMPLITUDE. | |
44 | PROPH=-PROPH | |
45 | C | |
46 | CW=1./ZM | |
47 | CW2=CW**2 | |
48 | SW2=1.-CW2 | |
49 | SW=SQRT(SW2) | |
50 | QQ0=.5*RTS | |
51 | QQI=.5*SQRT(S-4.*ZM2) | |
52 | QQF=.5*SQRT(S-4.) | |
53 | GSQ=4.*PI*ALFA/SW2 | |
54 | C | |
55 | C FROM ZZWW3.EX | |
56 | ANWWWW(1,1,1) = -1.60E+01 * ((HM2 * ZM3) / (CW * PROPH)) | |
57 | $ + 1.60E+01 * ((S * ZM3) / (CW * PROPH) ) - 8.00E+00 * ((S2 | |
58 | $ * ZM) / (CW * PROPH)) - 8.00E+00 * ((S2 * ZM3) / (CW * PROPH)) | |
59 | $ + 4.00E+00 * ((S3 * ZM) / (CW * PROPH)) + 8.00E+00 * ((HM2 * S | |
60 | $ * ZM) / (CW * PROPH)) + 8.00E+00 * ((HM2 * S * ZM3) / (CW | |
61 | $ * PROPH)) - 4.00E+00 * ((HM2 * S2 * ZM) / (CW * PROPH)) | |
62 | $ + 8.00E+00 * CW2 * S - 6.00E+00 * CW2 * S2 + 8.00E+00 * CW2 | |
63 | $ * S * ZM2 | |
64 | ANWWWW(1,1,2) = 4.00E+00 * CW2 * S | |
65 | ANWWWW(1,1,3) = -1.60E+01 * ((HM2 * ZM3) / (CW * PROPH)) | |
66 | $ + 1.60E+01 * ((S * ZM3) / (CW * PROPH) ) - 8.00E+00 * ((S2 | |
67 | $ * ZM) / (CW * PROPH)) + 8.00E+00 * ((HM2 * S * ZM) / (CW | |
68 | $ * PROPH)) + 1.20E+01 * CW2 * S - 3.20E+01 * CW2 * ZM2 | |
69 | ANWWWW(1,1,4) = 0.00E+00 | |
70 | ANWWWW(1,2,1) = -4.00E+00 * CW2 * S2 + 3.00E+00 * CW2 * S3 | |
71 | $ - 9.60E+01 * CW2 * ZM4 + 3.20E+01 * CW2 * ZM6 + 8.00E+00 * CW2 | |
72 | $ * S * ZM2 + 1.60E+01 * CW2 * S * ZM4 - 8.00E+00 * CW2 * S | |
73 | $ * ZM6 - 1.00E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4 | |
74 | ANWWWW(1,2,2) = -4.00E+00 * CW2 * S - 2.00E+00 * CW2 * S2 | |
75 | $ + 6.40E+01 * CW2 * ZM4 + 1.20E+01 * CW2 * S * ZM2 | |
76 | ANWWWW(1,2,3) = 4.00E+00 * CW2 * S - 6.00E+00 * CW2 * S2 | |
77 | $ - 6.40E+01 * CW2 * ZM4 + 2.00E+01 * CW2 * S * ZM2 | |
78 | ANWWWW(1,2,4) = 6.40E+01 * CW2 * QQ0 * QQF * QQI + 1.60E+01 | |
79 | $ * CW2 * QQ0 * QQF * QQI * S + 9.60E+01 * CW2 * QQ0 * QQF * QQI | |
80 | $ * ZM2 + 3.20E+01 * CW2 * QQ0 * QQF * QQI * ZM4 | |
81 | ANWWWW(1,3,1) = -4.00E+00 * CW2 * S2 + 3.00E+00 * CW2 * S3 | |
82 | $ - 9.60E+01 * CW2 * ZM4 + 3.20E+01 * CW2 * ZM6 + 8.00E+00 * CW2 | |
83 | $ * S * ZM2 + 1.60E+01 * CW2 * S * ZM4 - 8.00E+00 * CW2 * S | |
84 | $ * ZM6 - 1.00E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4 | |
85 | ANWWWW(1,3,2) = -4.00E+00 * CW2 * S - 2.00E+00 * CW2 * S2 | |
86 | $ + 6.40E+01 * CW2 * ZM4 + 1.20E+01 * CW2 * S * ZM2 | |
87 | ANWWWW(1,3,3) = 4.00E+00 * CW2 * S - 6.00E+00 * CW2 * S2 | |
88 | $ - 6.40E+01 * CW2 * ZM4 + 2.00E+01 * CW2 * S * ZM2 | |
89 | ANWWWW(1,3,4) = -6.40E+01 * CW2 * QQ0 * QQF * QQI - 1.60E+01 | |
90 | $ * CW2 * QQ0 * QQF * QQI * S - 9.60E+01 * CW2 * QQ0 * QQF * QQI | |
91 | $ * ZM2 - 3.20E+01 * CW2 * QQ0 * QQF * QQI * ZM4 | |
92 | ANWWWW(1,4,1) = 0.00E+00 | |
93 | ANWWWW(1,4,2) = 0.00E+00 | |
94 | ANWWWW(1,4,3) = 0.00E+00 | |
95 | ANWWWW(1,4,4) = 0.00E+00 | |
96 | ANWWWW(2,1,1) = 0.00E+00 | |
97 | ANWWWW(2,1,2) = 0.00E+00 | |
98 | ANWWWW(2,1,3) = 0.00E+00 | |
99 | ANWWWW(2,1,4) = 8.00E+00 * CW2 * QQ0 * S | |
100 | ANWWWW(2,2,1) = 4.80E+01 * CW2 * QQF * QQI * S - 2.00E+01 * CW2 | |
101 | $ * QQF * QQI * S2 + 6.40E+01 * CW2 * QQF * QQI * ZM2 + 9.60E+01 | |
102 | $ * CW2 * QQF * QQI * S * ZM2 - 1.60E+01 * CW2 * QQF * QQI * S | |
103 | $ * ZM4 | |
104 | ANWWWW(2,2,2) = 8.00E+00 * CW2 * QQF * QQI * S | |
105 | ANWWWW(2,2,3) = 5.60E+01 * CW2 * QQF * QQI * S | |
106 | ANWWWW(2,2,4) = -8.00E+00 * CW2 * QQ0 * S2 + 1.92E+02 * CW2 | |
107 | $ * QQ0 * ZM2 + 1.28E+02 * CW2 * QQ0 * ZM4 - 8.00E+00 * CW2 | |
108 | $ * QQ0 * S * ZM2 - 8.00E+00 * CW2 * QQ0 * S * ZM4 | |
109 | ANWWWW(2,3,1) = -4.80E+01 * CW2 * QQF * QQI * S + 2.00E+01 * CW2 | |
110 | $ * QQF * QQI * S2 - 6.40E+01 * CW2 * QQF * QQI * ZM2 | |
111 | $ - 9.60E+01 * CW2 * QQF * QQI * S * ZM2 + 1.60E+01 * CW2 * QQF | |
112 | $ * QQI * S * ZM4 | |
113 | ANWWWW(2,3,2) = -8.00E+00 * CW2 * QQF * QQI * S | |
114 | ANWWWW(2,3,3) = -5.60E+01 * CW2 * QQF * QQI * S | |
115 | ANWWWW(2,3,4) = -8.00E+00 * CW2 * QQ0 * S2 + 1.92E+02 * CW2 | |
116 | $ * QQ0 * ZM2 + 1.28E+02 * CW2 * QQ0 * ZM4 - 8.00E+00 * CW2 | |
117 | $ * QQ0 * S * ZM2 - 8.00E+00 * CW2 * QQ0 * S * ZM4 | |
118 | ANWWWW(2,4,1) = 0.00E+00 | |
119 | ANWWWW(2,4,2) = 0.00E+00 | |
120 | ANWWWW(2,4,3) = 0.00E+00 | |
121 | ANWWWW(2,4,4) = 0.00E+00 | |
122 | ANWWWW(3,1,1) = 2.00E+00 * CW2 * S2 | |
123 | ANWWWW(3,1,2) = -4.00E+00 * CW2 * S | |
124 | ANWWWW(3,1,3) = 4.00E+00 * CW2 * S | |
125 | ANWWWW(3,1,4) = 0.00E+00 | |
126 | ANWWWW(3,2,1) = 3.20E+01 * CW2 * S - 1.60E+01 * CW2 * S2 + CW2 | |
127 | $ * S3 + 9.60E+01 * CW2 * S * ZM2 + 3.20E+01 * CW2 * S * ZM4 | |
128 | $ - 2.20E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4 | |
129 | ANWWWW(3,2,2) = 4.00E+00 * CW2 * S + 2.00E+00 * CW2 * S2 | |
130 | $ - 6.40E+01 * CW2 * ZM4 - 1.20E+01 * CW2 * S * ZM2 | |
131 | ANWWWW(3,2,3) = 2.80E+01 * CW2 * S - 1.00E+01 * CW2 * S2 | |
132 | $ + 6.40E+01 * CW2 * ZM4 + 1.20E+01 * CW2 * S * ZM2 | |
133 | ANWWWW(3,2,4) = 1.60E+01 * CW2 * QQ0 * QQF * QQI * S | |
134 | ANWWWW(3,3,1) = 3.20E+01 * CW2 * S - 1.60E+01 * CW2 * S2 + CW2 | |
135 | $ * S3 + 9.60E+01 * CW2 * S * ZM2 + 3.20E+01 * CW2 * S * ZM4 | |
136 | $ - 2.20E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4 | |
137 | ANWWWW(3,3,2) = 4.00E+00 * CW2 * S + 2.00E+00 * CW2 * S2 | |
138 | $ - 6.40E+01 * CW2 * ZM4 - 1.20E+01 * CW2 * S * ZM2 | |
139 | ANWWWW(3,3,3) = 2.80E+01 * CW2 * S - 1.00E+01 * CW2 * S2 | |
140 | $ + 6.40E+01 * CW2 * ZM4 + 1.20E+01 * CW2 * S * ZM2 | |
141 | ANWWWW(3,3,4) = -1.60E+01 * CW2 * QQ0 * QQF * QQI * S | |
142 | ANWWWW(3,4,1) = 0.00E+00 | |
143 | ANWWWW(3,4,2) = 0.00E+00 | |
144 | ANWWWW(3,4,3) = 0.00E+00 | |
145 | ANWWWW(3,4,4) = 0.00E+00 | |
146 | ANWWWW(4,1,1) = 0.00E+00 | |
147 | ANWWWW(4,1,2) = 0.00E+00 | |
148 | ANWWWW(4,1,3) = 0.00E+00 | |
149 | ANWWWW(4,1,4) = 0.00E+00 | |
150 | ANWWWW(4,2,1) = 4.00E+00 * CW2 * QQF * QQI * S2 | |
151 | ANWWWW(4,2,2) = -8.00E+00 * CW2 * QQF * QQI * S | |
152 | ANWWWW(4,2,3) = 8.00E+00 * CW2 * QQF * QQI * S | |
153 | ANWWWW(4,2,4) = 0.00E+00 | |
154 | ANWWWW(4,3,1) = -4.00E+00 * CW2 * QQF * QQI * S2 | |
155 | ANWWWW(4,3,2) = 8.00E+00 * CW2 * QQF * QQI * S | |
156 | ANWWWW(4,3,3) = -8.00E+00 * CW2 * QQF * QQI * S | |
157 | ANWWWW(4,3,4) = 0.00E+00 | |
158 | ANWWWW(4,4,1) = 0.00E+00 | |
159 | ANWWWW(4,4,2) = 0.00E+00 | |
160 | ANWWWW(4,4,3) = 0.00E+00 | |
161 | ANWWWW(4,4,4) = 0.00E+00 | |
162 | C | |
163 | ADWWWW(1,1) = 1.00E+00 | |
164 | ADWWWW(1,2) = S - 2.00E+00 * ZM2 | |
165 | ADWWWW(1,3) = S - 2.00E+00 * ZM2 | |
166 | ADWWWW(1,4) = 1.00E+00 | |
167 | ADWWWW(2,1) = 0.00E+00 | |
168 | ADWWWW(2,2) = -1.00E+00 * S + 4.00E+00 * (S / (S + 4.00E+00 | |
169 | $ * QQF * QQI)) - 1.60E+01 * (ZM2 / (S + 4.00E+00 * QQF * QQI)) | |
170 | $ + 4.00E+00 * ((S * ZM2) / (S + 4.00E+00 * QQF * QQI)) | |
171 | ADWWWW(2,3) = S - 4.00E+00 * (S / (S + 4.00E+00 * QQF * QQI)) | |
172 | $ + 1.60E+01 * (ZM2 / (S + 4.00E+00 * QQF * QQI)) - 4.00E+00 | |
173 | $ * ((S * ZM2) / (S + 4.00E+00 * QQF * QQI)) | |
174 | ADWWWW(2,4) = 0.00E+00 | |
175 | C | |
176 | AIWWWW(1) = 1.60E+01 * ((HG * HM * ZM3) / (CW * PROPH)) | |
177 | $ - 8.00E+00 * ((HG * HM * S * ZM) / (CW * PROPH)) | |
178 | $ - 8.00E+00 * ((HG * HM * S * ZM3) / (CW * PROPH)) + 4.00E+00 | |
179 | $ * ((HG * HM * S2 * ZM) / (CW * PROPH)) | |
180 | AIWWWW(2) = 0.00E+00 | |
181 | AIWWWW(3) = 1.60E+01 * ((HG * HM * ZM3) / (CW * PROPH)) | |
182 | $ - 8.00E+00 * ((HG * HM * S * ZM) / (CW * PROPH)) | |
183 | AIWWWW(4) = 0.00E+00 | |
184 | C | |
185 | C RESTORE MISSING FACTORS | |
186 | DO 100 J=1,4 | |
187 | AIWWWW(J)=AIWWWW(J)*GSQ/(16.*ZM2) | |
188 | DO 100 I=1,4 | |
189 | DO 110 K=1,4 | |
190 | 110 ANWWWW(K,I,J)=ANWWWW(K,I,J)*GSQ/(16.*ZM2) | |
191 | 100 CONTINUE | |
192 | C | |
193 | RETURN | |
194 | END |