]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | SUBROUTINE XWWWW | |
3 | C | |
4 | C SET UP W+ W- -> 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 | |
44 | CW=1./ZM | |
45 | CW2=CW**2 | |
46 | SW2=1.-CW2 | |
47 | SW=SQRT(SW2) | |
48 | QQ0=.5*RTS | |
49 | QQI=.5*SQRT(S-4.) | |
50 | QQF=.5*SQRT(S-4.) | |
51 | GSQ=4.*PI*ALFA/SW2 | |
52 | C | |
53 | C FROM WWWW3.EX | |
54 | ANWWWW(1,1,1) = 8.00E+00 * S - 3.00E+00 * S2 - 1.60E+01 | |
55 | $ * ( HM2 / PROPH) + 1.60E+01 * (S / PROPH) - 1.60E+01 * (S2 | |
56 | $ / PROPH) + 4.00E+00 * (S3 / PROPH) + 1.60E+01 * ((HM2 * S) | |
57 | $ / PROPH) - 4.00E+00 * ((HM2 * S2) / PROPH) | |
58 | ANWWWW(1,1,2) = 2.00E+00 * S | |
59 | ANWWWW(1,1,3) = -1.60E+01 + 6.00E+00 * S - 1.60E+01 * (HM2 | |
60 | $ / PROPH) + 1.60E+01 * (S / PROPH) - 8.00E+00 * (S2 / PROPH) | |
61 | $ + 8.00E+00 * ((HM2 * S) / PROPH) | |
62 | ANWWWW(1,1,4) = -2.40E+01 * RTS + 6.40E+01 * (RTS / (S | |
63 | $ - 1.00E+00 * ZM2)) + 1.60E+01 * ((RTS * S) / (S - 1.00E+00 | |
64 | $ * ZM2)) - 8.00E+00 * ((RTS * S2) / (S - 1.00E+00 * ZM2)) | |
65 | $ + 6.40E+01 * ((RTS * SW2) / S) - 6.40E+01 * ((RTS * SW2) / (S | |
66 | $ - 1.00E+00 * ZM2)) - 1.60E+01 * ((RTS * S * SW2) / (S | |
67 | $ - 1.00E+00 * ZM2)) + 8.00E+00 * ((RTS * S2 * SW2) / (S | |
68 | $ - 1.00E+00 * ZM2)) + 6.00E+00 * RTS * S + 1.60E+01 * RTS | |
69 | $ * SW2 - 8.00E+00 * RTS * S * SW2 | |
70 | ANWWWW(1,2,1) = -6.40E+01 + 1.60E+01 * S - 1.20E+01 * S2 | |
71 | $ + 3.00E+00 * S3 + 6.40E+01 * SW2 - 1.60E+01 * S * SW2 | |
72 | $ + 1.20E+01 * S2 * SW2 - 3.00E+00 * S3 * SW2 | |
73 | ANWWWW(1,2,2) = 6.40E+01 + 8.00E+00 * S - 2.00E+00 * S2 | |
74 | $ - 6.40E+01 * SW2 - 8.00E+00 * S * SW2 + 2.00E+00 * S2 * SW2 | |
75 | ANWWWW(1,2,3) = -6.40E+01 + 2.40E+01 * S - 6.00E+00 * S2 | |
76 | $ + 6.40E+01 * SW2 - 2.40E+01 * S * SW2 + 6.00E+00 * S2 * SW2 | |
77 | ANWWWW(1,2,4) = -9.60E+01 * RTS + 1.60E+01 * RTS * S + 2.00E+00 | |
78 | $ * RTS * S2 + 9.60E+01 * RTS * SW2 - 1.60E+01 * RTS * S * SW2 | |
79 | $ - 2.00E+00 * RTS * S2 * SW2 | |
80 | ANWWWW(1,3,1) = -6.40E+01 * SW2 + 1.60E+01 * S * SW2 - 1.20E+01 | |
81 | $ * S2 * SW2 + 3.00E+00 * S3 * SW2 | |
82 | ANWWWW(1,3,2) = 6.40E+01 * SW2 + 8.00E+00 * S * SW2 - 2.00E+00 | |
83 | $ * S2 * SW2 | |
84 | ANWWWW(1,3,3) = -6.40E+01 * SW2 + 2.40E+01 * S * SW2 - 6.00E+00 | |
85 | $ * S2 * SW2 | |
86 | ANWWWW(1,3,4) = -9.60E+01 * RTS * SW2 + 1.60E+01 * RTS * S * SW2 | |
87 | $ + 2.00E+00 * RTS * S2 * SW2 | |
88 | ANWWWW(1,4,1) = -3.20E+01 + 1.60E+01 * S - 2.00E+00 * S2 | |
89 | ANWWWW(1,4,2) = -4.00E+00 * S | |
90 | ANWWWW(1,4,3) = 4.00E+00 * S | |
91 | ANWWWW(1,4,4) = -1.60E+01 * RTS + 4.00E+00 * RTS * S | |
92 | ANWWWW(2,1,1) = -2.40E+01 * S + 6.00E+00 * S2 + 4.80E+01 * SW2 | |
93 | $ + 6.40E+01 * (1.00E+00 / (S - 1.00E+00 * ZM2)) + 4.80E+01 * (S | |
94 | $ / (S - 1.00E+00 * ZM2)) - 4.00E+00 * (S3 / (S - 1.00E+00 | |
95 | $ * ZM2)) + 6.40E+01 * (SW2 / S) - 6.40E+01 * (SW2 / (S | |
96 | $ - 1.00E+00 * ZM2)) - 4.80E+01 * ((S * SW2) / (S - 1.00E+00 | |
97 | $ * ZM2)) + 4.00E+00 * ((S3 * SW2) / (S - 1.00E+00 * ZM2)) | |
98 | $ - 4.00E+00 * S2 * SW2 | |
99 | ANWWWW(2,1,2) = 0.00E+00 | |
100 | ANWWWW(2,1,3) = 1.60E+01 * SW2 + 6.40E+01 * (1.00E+00 / (S | |
101 | $ - 1.00E+00 * ZM2)) + 1.60E+01 * (S / (S - 1.00E+00 * ZM2)) | |
102 | $ - 8.00E+00 * (S2 / (S - 1.00E+00 * ZM2)) + 6.40E+01 * (SW2 | |
103 | $ / S) - 6.40E+01 * (SW2 / (S - 1.00E+00 * ZM2)) - 1.60E+01 | |
104 | $ * ((S * SW2) / (S - 1.00E+00 * ZM2)) + 8.00E+00 * ((S2 * SW2) | |
105 | $ / (S - 1.00E+00 * ZM2)) - 8.00E+00 * S * SW2 | |
106 | ANWWWW(2,1,4) = 2.00E+00 * RTS * S | |
107 | ANWWWW(2,2,1) = -6.40E+01 - 1.12E+02 * S + 5.20E+01 * S2 | |
108 | $ - 5.00E+00 * S3 + 6.40E+01 * SW2 + 1.12E+02 * S * SW2 | |
109 | $ - 5.20E+01 * S2 * SW2 + 5.00E+00 * S3 * SW2 | |
110 | ANWWWW(2,2,2) = -8.00E+00 * S + 2.00E+00 * S2 + 8.00E+00 * S | |
111 | $ * SW2 - 2.00E+00 * S2 * SW2 | |
112 | ANWWWW(2,2,3) = -5.60E+01 * S + 1.40E+01 * S2 + 5.60E+01 * S | |
113 | $ * SW2 - 1.40E+01 * S2 * SW2 | |
114 | ANWWWW(2,2,4) = 1.60E+02 * RTS - 8.00E+00 * RTS * S - 4.00E+00 | |
115 | $ * RTS * S2 - 1.60E+02 * RTS * SW2 + 8.00E+00 * RTS * S * SW2 | |
116 | $ + 4.00E+00 * RTS * S2 * SW2 | |
117 | ANWWWW(2,3,1) = -6.40E+01 * SW2 - 1.12E+02 * S * SW2 + 5.20E+01 | |
118 | $ * S2 * SW2 - 5.00E+00 * S3 * SW2 | |
119 | ANWWWW(2,3,2) = -8.00E+00 * S * SW2 + 2.00E+00 * S2 * SW2 | |
120 | ANWWWW(2,3,3) = -5.60E+01 * S * SW2 + 1.40E+01 * S2 * SW2 | |
121 | ANWWWW(2,3,4) = 1.60E+02 * RTS * SW2 - 8.00E+00 * RTS * S * SW2 | |
122 | $ - 4.00E+00 * RTS * S2 * SW2 | |
123 | ANWWWW(2,4,1) = -1.60E+01 * S + 4.00E+00 * S2 | |
124 | ANWWWW(2,4,2) = 0.00E+00 | |
125 | ANWWWW(2,4,3) = 0.00E+00 | |
126 | ANWWWW(2,4,4) = -4.00E+00 * RTS * S | |
127 | ANWWWW(3,1,1) = S2 | |
128 | ANWWWW(3,1,2) = -2.00E+00 * S | |
129 | ANWWWW(3,1,3) = 2.00E+00 * S | |
130 | ANWWWW(3,1,4) = 0.00E+00 | |
131 | ANWWWW(3,2,1) = 1.60E+02 * S - 3.60E+01 * S2 + S3 - 1.60E+02 | |
132 | $ * S * SW2 + 3.60E+01 * S2 * SW2 - 1.00E+00 * S3 * SW2 | |
133 | ANWWWW(3,2,2) = -6.40E+01 - 8.00E+00 * S + 2.00E+00 * S2 | |
134 | $ + 6.40E+01 * SW2 + 8.00E+00 * S * SW2 - 2.00E+00 * S2 * SW2 | |
135 | ANWWWW(3,2,3) = 6.40E+01 + 4.00E+01 * S - 1.00E+01 * S2 | |
136 | $ - 6.40E+01 * SW2 - 4.00E+01 * S * SW2 + 1.00E+01 * S2 * SW2 | |
137 | ANWWWW(3,2,4) = -8.00E+00 * RTS * S + 2.00E+00 * RTS * S2 | |
138 | $ + 8.00E+00 * RTS * S * SW2 - 2.00E+00 * RTS * S2 * SW2 | |
139 | ANWWWW(3,3,1) = 1.60E+02 * S * SW2 - 3.60E+01 * S2 * SW2 + S3 | |
140 | $ * SW2 | |
141 | ANWWWW(3,3,2) = -6.40E+01 * SW2 - 8.00E+00 * S * SW2 + 2.00E+00 | |
142 | $ * S2 * SW2 | |
143 | ANWWWW(3,3,3) = 6.40E+01 * SW2 + 4.00E+01 * S * SW2 - 1.00E+01 | |
144 | $ * S2 * SW2 | |
145 | ANWWWW(3,3,4) = -8.00E+00 * RTS * S * SW2 + 2.00E+00 * RTS * S2 | |
146 | $ * SW2 | |
147 | ANWWWW(3,4,1) = -2.00E+00 * S2 | |
148 | ANWWWW(3,4,2) = 4.00E+00 * S | |
149 | ANWWWW(3,4,3) = -4.00E+00 * S | |
150 | ANWWWW(3,4,4) = 0.00E+00 | |
151 | ANWWWW(4,1,1) = 0.00E+00 | |
152 | ANWWWW(4,1,2) = 0.00E+00 | |
153 | ANWWWW(4,1,3) = 0.00E+00 | |
154 | ANWWWW(4,1,4) = 0.00E+00 | |
155 | ANWWWW(4,2,1) = -4.00E+00 * S2 + S3 + 4.00E+00 * S2 * SW2 | |
156 | $ - 1.00E+00 * S3 * SW2 | |
157 | ANWWWW(4,2,2) = 8.00E+00 * S - 2.00E+00 * S2 - 8.00E+00 * S | |
158 | $ * SW2 + 2.00E+00 * S2 * SW2 | |
159 | ANWWWW(4,2,3) = -8.00E+00 * S + 2.00E+00 * S2 + 8.00E+00 * S | |
160 | $ * SW2 - 2.00E+00 * S2 * SW2 | |
161 | ANWWWW(4,2,4) = 0.00E+00 | |
162 | ANWWWW(4,3,1) = -4.00E+00 * S2 * SW2 + S3 * SW2 | |
163 | ANWWWW(4,3,2) = 8.00E+00 * S * SW2 - 2.00E+00 * S2 * SW2 | |
164 | ANWWWW(4,3,3) = -8.00E+00 * S * SW2 + 2.00E+00 * S2 * SW2 | |
165 | ANWWWW(4,3,4) = 0.00E+00 | |
166 | ANWWWW(4,4,1) = 0.00E+00 | |
167 | ANWWWW(4,4,2) = 0.00E+00 | |
168 | ANWWWW(4,4,3) = 0.00E+00 | |
169 | ANWWWW(4,4,4) = 0.00E+00 | |
170 | C | |
171 | ADWWWW(1,1) = 1.00E+00 | |
172 | ADWWWW(1,2) = -4.00E+00 + S + 2.00E+00 * ZM2 | |
173 | ADWWWW(1,3) = -4.00E+00 + S | |
174 | ADWWWW(1,4) = -4.00E+00 + 2.00E+00 * HM2 + S | |
175 | ADWWWW(2,1) = 0.00E+00 | |
176 | ADWWWW(2,2) = 4.00E+00 - 1.00E+00 * S | |
177 | ADWWWW(2,3) = 4.00E+00 - 1.00E+00 * S | |
178 | ADWWWW(2,4) = 4.00E+00 - 1.00E+00 * S | |
179 | C | |
180 | AIWWWW(1) = 1.60E+01 * ((HG * HM) / PROPH) - 1.60E+01 * ((HG | |
181 | $ * HM * S) / PROPH) + 4.00E+00 * ((HG * HM * S2) / PROPH) | |
182 | AIWWWW(2) = 0.00E+00 | |
183 | AIWWWW(3) = 1.60E+01 * ((HG * HM) / PROPH) - 8.00E+00 * ((HG | |
184 | $ * HM * S) / PROPH) | |
185 | AIWWWW(4) = 0.00E+00 | |
186 | C | |
187 | C RESTORE MISSING FACTORS | |
188 | DO 100 J=1,4 | |
189 | AIWWWW(J)=AIWWWW(J)*GSQ/(16.) | |
190 | DO 100 I=1,4 | |
191 | DO 110 K=1,4 | |
192 | 110 ANWWWW(K,I,J)=ANWWWW(K,I,J)*GSQ/(16.) | |
193 | 100 CONTINUE | |
194 | C | |
195 | RETURN | |
196 | END |