1 #include "isajet/pilot.h"
4 C SET UP W+ W- -> Z0 Z0 AMPLITUDES AS RATIONAL FUNCTIONS OF Z
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
12 C *NOTE* A FACTOR OF SIN(THETA)/SQRT(2) IS REMOVED FROM F01
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
23 C USE UNITS OF WM TO AVOID LARGE NUMBERS - NOTE ANWWWW/ADWWWW
24 C AND AIWWWW ARE DIMENSIONLESS
42 PROPH=(S-HM2)**2+(HM*HG)**2
43 C CORRECT SIGN OF HIGGS AMPLITUDE.
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
64 ANWWWW(1,1,2) = 4.00E+00 * CW2 * S * ZM2
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 $ * ZM3) / (CW * PROPH)) + 8.00E+00 * ((HM2 * S * ZM3) / (CW
68 $ * PROPH)) - 3.20E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * 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) = 6.40E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * ZM2
75 $ - 4.00E+00 * CW2 * S * ZM4 - 2.00E+00 * CW2 * S2 * ZM2
76 ANWWWW(1,2,3) = -6.40E+01 * CW2 * ZM2 + 2.00E+01 * CW2 * S * ZM2
77 $ + 4.00E+00 * CW2 * S * ZM4 - 6.00E+00 * CW2 * S2 * ZM2
78 ANWWWW(1,2,4) = 1.92E+02 * CW2 * QQ0 * QQF * QQI * ZM
79 $ - 3.20E+01 * CW2 * QQ0 * QQF * QQI * ZM3 + 3.20E+01 * CW2
80 $ * QQ0 * QQF * QQI * ZM5 + 1.60E+01 * CW2 * QQ0 * QQF * QQI * S
82 ANWWWW(1,3,1) = -4.00E+00 * CW2 * S2 + 3.00E+00 * CW2 * S3
83 $ - 9.60E+01 * CW2 * ZM4 + 3.20E+01 * CW2 * ZM6 + 8.00E+00 * CW2
84 $ * S * ZM2 + 1.60E+01 * CW2 * S * ZM4 - 8.00E+00 * CW2 * S
85 $ * ZM6 - 1.00E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
86 ANWWWW(1,3,2) = 6.40E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * ZM2
87 $ - 4.00E+00 * CW2 * S * ZM4 - 2.00E+00 * CW2 * S2 * ZM2
88 ANWWWW(1,3,3) = -6.40E+01 * CW2 * ZM2 + 2.00E+01 * CW2 * S * ZM2
89 $ + 4.00E+00 * CW2 * S * ZM4 - 6.00E+00 * CW2 * S2 * ZM2
90 ANWWWW(1,3,4) = -1.92E+02 * CW2 * QQ0 * QQF * QQI * ZM
91 $ + 3.20E+01 * CW2 * QQ0 * QQF * QQI * ZM3 - 3.20E+01 * CW2
92 $ * QQ0 * QQF * QQI * ZM5 - 1.60E+01 * CW2 * QQ0 * QQF * QQI * S
94 ANWWWW(1,4,1) = 0.00E+00
95 ANWWWW(1,4,2) = 0.00E+00
96 ANWWWW(1,4,3) = 0.00E+00
97 ANWWWW(1,4,4) = 0.00E+00
98 ANWWWW(2,1,1) = 0.00E+00
99 ANWWWW(2,1,2) = 0.00E+00
100 ANWWWW(2,1,3) = 0.00E+00
101 ANWWWW(2,1,4) = 8.00E+00 * CW2 * QQ0 * S * ZM
102 ANWWWW(2,2,1) = 4.80E+01 * CW2 * QQF * QQI * S - 2.00E+01 * CW2
103 $ * QQF * QQI * S2 + 6.40E+01 * CW2 * QQF * QQI * ZM2 + 9.60E+01
104 $ * CW2 * QQF * QQI * S * ZM2 - 1.60E+01 * CW2 * QQF * QQI * S
106 ANWWWW(2,2,2) = 8.00E+00 * CW2 * QQF * QQI * S * ZM2
107 ANWWWW(2,2,3) = 5.60E+01 * CW2 * QQF * QQI * S * ZM2
108 ANWWWW(2,2,4) = 1.28E+02 * CW2 * QQ0 * ZM + 1.92E+02 * CW2 * QQ0
109 $ * ZM3 - 3.20E+01 * CW2 * QQ0 * S * ZM + 2.40E+01 * CW2 * QQ0
110 $ * S * ZM3 - 8.00E+00 * CW2 * QQ0 * S * ZM5 - 8.00E+00 * CW2
112 ANWWWW(2,3,1) = -4.80E+01 * CW2 * QQF * QQI * S + 2.00E+01 * CW2
113 $ * QQF * QQI * S2 - 6.40E+01 * CW2 * QQF * QQI * ZM2
114 $ - 9.60E+01 * CW2 * QQF * QQI * S * ZM2 + 1.60E+01 * CW2 * QQF
116 ANWWWW(2,3,2) = -8.00E+00 * CW2 * QQF * QQI * S * ZM2
117 ANWWWW(2,3,3) = -5.60E+01 * CW2 * QQF * QQI * S * ZM2
118 ANWWWW(2,3,4) = 1.28E+02 * CW2 * QQ0 * ZM + 1.92E+02 * CW2 * QQ0
119 $ * ZM3 - 3.20E+01 * CW2 * QQ0 * S * ZM + 2.40E+01 * CW2 * QQ0
120 $ * S * ZM3 - 8.00E+00 * CW2 * QQ0 * S * ZM5 - 8.00E+00 * CW2
122 ANWWWW(2,4,1) = 0.00E+00
123 ANWWWW(2,4,2) = 0.00E+00
124 ANWWWW(2,4,3) = 0.00E+00
125 ANWWWW(2,4,4) = 0.00E+00
126 ANWWWW(3,1,1) = 2.00E+00 * CW2 * S2
127 ANWWWW(3,1,2) = -4.00E+00 * CW2 * S * ZM2
128 ANWWWW(3,1,3) = 4.00E+00 * CW2 * S * ZM2
129 ANWWWW(3,1,4) = 0.00E+00
130 ANWWWW(3,2,1) = 3.20E+01 * CW2 * S - 1.60E+01 * CW2 * S2 + CW2
131 $ * S3 + 9.60E+01 * CW2 * S * ZM2 + 3.20E+01 * CW2 * S * ZM4
132 $ - 2.20E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
133 ANWWWW(3,2,2) = -6.40E+01 * CW2 * ZM2 - 1.20E+01 * CW2 * S
134 $ * ZM2 + 4.00E+00 * CW2 * S * ZM4 + 2.00E+00 * CW2 * S2 * ZM2
135 ANWWWW(3,2,3) = 6.40E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * ZM2
136 $ + 2.80E+01 * CW2 * S * ZM4 - 1.00E+01 * CW2 * S2 * ZM2
137 ANWWWW(3,2,4) = 1.60E+01 * CW2 * QQ0 * QQF * QQI * S * ZM
138 ANWWWW(3,3,1) = 3.20E+01 * CW2 * S - 1.60E+01 * CW2 * S2 + CW2
139 $ * S3 + 9.60E+01 * CW2 * S * ZM2 + 3.20E+01 * CW2 * S * ZM4
140 $ - 2.20E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
141 ANWWWW(3,3,2) = -6.40E+01 * CW2 * ZM2 - 1.20E+01 * CW2 * S
142 $ * ZM2 + 4.00E+00 * CW2 * S * ZM4 + 2.00E+00 * CW2 * S2 * ZM2
143 ANWWWW(3,3,3) = 6.40E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * ZM2
144 $ + 2.80E+01 * CW2 * S * ZM4 - 1.00E+01 * CW2 * S2 * ZM2
145 ANWWWW(3,3,4) = -1.60E+01 * CW2 * QQ0 * QQF * QQI * S * ZM
146 ANWWWW(3,4,1) = 0.00E+00
147 ANWWWW(3,4,2) = 0.00E+00
148 ANWWWW(3,4,3) = 0.00E+00
149 ANWWWW(3,4,4) = 0.00E+00
150 ANWWWW(4,1,1) = 0.00E+00
151 ANWWWW(4,1,2) = 0.00E+00
152 ANWWWW(4,1,3) = 0.00E+00
153 ANWWWW(4,1,4) = 0.00E+00
154 ANWWWW(4,2,1) = 4.00E+00 * CW2 * QQF * QQI * S2
155 ANWWWW(4,2,2) = -8.00E+00 * CW2 * QQF * QQI * S * ZM2
156 ANWWWW(4,2,3) = 8.00E+00 * CW2 * QQF * QQI * S * ZM2
157 ANWWWW(4,2,4) = 0.00E+00
158 ANWWWW(4,3,1) = -4.00E+00 * CW2 * QQF * QQI * S2
159 ANWWWW(4,3,2) = 8.00E+00 * CW2 * QQF * QQI * S * ZM2
160 ANWWWW(4,3,3) = -8.00E+00 * CW2 * QQF * QQI * S * ZM2
161 ANWWWW(4,3,4) = 0.00E+00
162 ANWWWW(4,4,1) = 0.00E+00
163 ANWWWW(4,4,2) = 0.00E+00
164 ANWWWW(4,4,3) = 0.00E+00
165 ANWWWW(4,4,4) = 0.00E+00
167 ADWWWW(1,1) = 1.00E+00
168 ADWWWW(1,2) = S - 2.00E+00 * ZM2
169 ADWWWW(1,3) = S - 2.00E+00 * ZM2
170 ADWWWW(1,4) = 1.00E+00
171 ADWWWW(2,1) = 0.00E+00
172 ADWWWW(2,2) = -1.00E+00 * S + 4.00E+00 * (S / (S + 4.00E+00
173 $ * QQF * QQI)) - 1.60E+01 * (ZM2 / (S + 4.00E+00 * QQF * QQI))
174 $ + 4.00E+00 * ((S * ZM2) / (S + 4.00E+00 * QQF * QQI))
175 ADWWWW(2,3) = S - 4.00E+00 * (S / (S + 4.00E+00 * QQF * QQI))
176 $ + 1.60E+01 * (ZM2 / (S + 4.00E+00 * QQF * QQI)) - 4.00E+00
177 $ * ((S * ZM2) / (S + 4.00E+00 * QQF * QQI))
178 ADWWWW(2,4) = 0.00E+00
180 AIWWWW(1) = 1.60E+01 * ((HG * HM * ZM3) / (CW * PROPH))
181 $ - 8.00E+00 * ((HG * HM * S * ZM) / (CW * PROPH)) - 8.00E+00
182 $ * ((HG * HM * S * ZM3) / (CW * PROPH)) + 4.00E+00 * ((HG * HM
183 $ * S2 * ZM) / (CW * PROPH))
185 AIWWWW(3) = 1.60E+01 * ((HG * HM * ZM3) / (CW * PROPH))
186 $ - 8.00E+00 * ((HG * HM * S * ZM3) / (CW * PROPH))
189 C RESTORE MISSING FACTORS
191 AIWWWW(J)=AIWWWW(J)*GSQ/(16.*ZM2)
194 110 ANWWWW(K,I,J)=ANWWWW(K,I,J)*GSQ/(16.*ZM2)