]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ISAJET/isasusy/sswibf.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
EffC++ warnings corrected.
[u/mrichter/AliRoot.git] / ISAJET / isasusy / sswibf.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE SSWIBF
3C-----------------------------------------------------------------------
4C This subroutine calculates the chargino (wi) subset of
5C SSWZBF, which was too long.
6C Valid for all scalar masses (functions in double precision)
7C Includes Higgs sector radiative corrections (Aug. 31)
8C
9C Auxiliary functions are called SSWxyi, SSZxyi, where normally
10C x indicates the SUSY particle, y the SM particle(s), and i is
11C a counter.
12C
13C Part of Baer's GAUGBF
14C
15C-----------------------------------------------------------------------
16#if defined(CERNLIB_IMPNONE)
17 IMPLICIT NONE
18#endif
19#include "isajet/sslun.inc"
20#include "isajet/ssmode.inc"
21#include "isajet/sssm.inc"
22#include "isajet/sspar.inc"
23#include "isajet/sstype.inc"
24#include "isajet/sstmp.inc"
25#include "isajet/sspols.inc"
26C
27 EXTERNAL SSZWF1,SSZZF1,SSZZF2,SSZZF3,SSWZF1,SSWZF2,SSWZF3
28 $,SSWZF4,SSWZF5,SSWZF6,SSWZF7,SSWWF1,SSZZF4,SSZZF5,SSGX1
29 $,SSGX2,SSGX8
30C
31 REAL AUI(4),BUI(4),ADI(4),BDI(4),ALI(4),BLI(4),ANI(4),BNI(4)
32 $,WIJ(4,4),AUWI(4),ADWI(4),ANWI(4),ALWI(4),XIM(4),YIM(4)
33 $,XIP(4),YIP(4),SNIJ(4,4),XLIJ(4,4),HIJ(4,4)
34 $,V1I(4),V2I(4),V3I(4),V4I(4),XHIJ(4,4),XPIJ(4,4),AMWISS(2)
35C
36 INTEGER ISZ(4),THJZ
37C
38 REAL MWIW,SL,PP,SP,PL,PH,SH,COSB,SINB,COSA,SINA
39 $,FACTOR,MZJZ,ULIM,XINTGL,COSBE,SINBE
40 $,UPPER,CONST
41 REAL W21ZU,W21ZN,W21ZL,W21ZD,W21U,W21D,W21S,W21C,W21N1
42 $,W21N2,W21N3,W21E,W21M,W21L,STHW,CTHW
43 REAL T3,XI2,CC,PSIINT,T2,T1,PHIINT,XI1,EF,A,Z,B
44 $,TANB,FB,FT,SR2,G,PI,GP,FL
45 $,MW2,SNW1,MW1,YM,BE,SNW2,XM,THX,THY
46 $,BTN,APD,APL,APU,BTD,BTL,APN,BTU,Y,MZ1,FPI
47 REAL TANW,COTW,XWINO,YWINO,SNIW,SNJZ
48 REAL SSXINT,SSXLAM
49 REAL WID,TERM1,TERM2,TERM3,TERM4,E,TERMH
50 REAL FACT,ALJZ1,ALJZ2,BEJZ1,BEJZ2
51 REAL ALIW1,ALIW2,AHCJZ,BHCJZ,TERMW,TERMN,TERM12,TERMN1,
52 $TERMN2,TERMWN,TERMW1,TERMW2,TERMH1,TERMH2,TERMHN
53 REAL XIPM,YIPM,COSL,SINL,BPWI(2),BPLWI(2)
54 REAL BWI(2),AS,COST,SINT
55 REAL POLNL,POLNR,POL1L,POL1R,POL2L,POL2R,POL12L,POL12R,
56 $POLN1L,POLN1R,POLN2L,POLN2R
57 REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS
58 REAL SUALFE,MTAMTA,MTAMB,MTAMZ,AMPL,AMPI
59 REAL FUDGE
60 DOUBLE PRECISION SSMQCD
61 COMPLEX ZI,ZONE,Z1(2),Z2(2)
62 INTEGER IW,JZ,IZ,ISZJZ,ISWIW
63 DATA FUDGE/1.0/,AMPI/.140/
64 DATA ZONE/(1.,0.)/,ZI/(0.,1.)/,FPI/.1315/
65C
66C Constants from neutralino mass matrix
67C
68 AMPL=2.4E18
69 PI=4.*ATAN(1.)
70 SR2=SQRT(2.)
71 G=SQRT(4*PI*ALFAEM/SN2THW)
72 GP=G*SQRT(SN2THW/(1.-SN2THW))
73 E=SQRT(4*PI/128.)
74C
75 TANW=SQRT(SN2THW/(1.-SN2THW))
76 COTW=1./TANW
77 STHW=SQRT(SN2THW)
78 CTHW=SQRT(1.-SN2THW)
79 APL=.25*(3*TANW-COTW)
80 BTL=.25*(COTW+TANW)
81 APN=.25*(TANW+COTW)
82 BTN=-.25*(COTW+TANW)
83 APU=-5*TANW/12.+COTW/4.
84 BTU=-.25*(COTW+TANW)
85 APD=-COTW/4.+TANW/12.
86 BTD=.25*(COTW+TANW)
87C
88 TANB=1./RV2V1
89 BE=ATAN(1./RV2V1)
90 SINBE=SIN(BE)
91 COSBE=COS(BE)
92 XM=1./TAN(GAMMAL)
93 YM=1./TAN(GAMMAR)
94 THX=SIGN(1.,XM)
95 THY=SIGN(1.,YM)
96 ASMB=SUALFS(AMBT**2,.36,AMTP,3)
97 MBMB=AMBT*(1.-4*ASMB/3./PI)
98 MBQ=SSMQCD(DBLE(MBMB),DBLE(AMTP))
99 ASMT=SUALFS(AMTP**2,.36,AMTP,3)
100 MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
101 $(ASMT/PI)**2)
102 MTQ=SSMQCD(DBLE(MTMT),DBLE(AMTP))
103 FB=G*MBQ/SR2/AMW/COS(BE)
104 FT=G*MTQ/SR2/AMW/SIN(BE)
105 MTAMTA=AMTAU*(1.-SUALFE(AMTAU**2)/PI)
106 MTAMB=MTAMTA*(SUALFE(AMBT**2)/SUALFE(AMTAU**2))**(-27./76.)
107 MTAMZ=MTAMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.)
108 FL=G*MTAMZ/SR2/AMW/COS(BE)
109 SNW1=SIGN(1.,AMW1SS)
110 SNW2=SIGN(1.,AMW2SS)
111 AMWISS(1)=AMW1SS
112 AMWISS(2)=AMW2SS
113 BWI(1)=-FT*SNW1*COS(GAMMAR)
114 BWI(2)=FT*SNW2*THY*SIN(GAMMAR)
115 BPWI(1)=-FB*COS(GAMMAL)
116 BPWI(2)=FB*THX*SIN(GAMMAL)
117 BPLWI(1)=-FL*COS(GAMMAL)
118 BPLWI(2)=FL*THX*SIN(GAMMAL)
119 MW1=ABS(AMW1SS)
120 MW2=ABS(AMW2SS)
121 MZ1=ABS(AMZ1SS)
122 XWINO=.5*(THX*SIN(GAMMAL)*COS(GAMMAL)
123 $-THY*SIN(GAMMAR)*COS(GAMMAR))
124 YWINO=.5*(THX*SIN(GAMMAL)*COS(GAMMAL)
125 $+THY*SIN(GAMMAR)*COS(GAMMAR))
126 COST=COS(THETAT)
127 SINT=SIN(THETAT)
128 COSB=COS(THETAB)
129 SINB=SIN(THETAB)
130 COSL=COS(THETAL)
131 SINL=SIN(THETAL)
132C
133C Constants from Higgs mass matrix
134C
135 SINA=SIN(ALFAH)
136 COSA=COS(ALFAH)
137C
138C Gaugino couplings
139C
140 DO 100 IZ=1,4
141 AUI(IZ)=G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
142 BUI(IZ)=4.*GP/3./SR2*ZMIXSS(4,IZ)
143 ADI(IZ)=-G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ)
144 BDI(IZ)=-2.*GP/3./SR2*ZMIXSS(4,IZ)
145 ALI(IZ)=G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ)
146 BLI(IZ)=-SR2*GP*ZMIXSS(4,IZ)
147 ANI(IZ)=G/SR2*ZMIXSS(3,IZ)-GP/SR2*ZMIXSS(4,IZ)
148 BNI(IZ)=0.0
149100 CONTINUE
150C
151 DO 110 IZ=1,4
152 DO 110 JZ=1,4
153 IF(IZ.LT.JZ) THEN
154 WIJ(IZ,JZ)=SQRT(G**2+GP**2)
155 $ *(ZMIXSS(1,IZ)*ZMIXSS(1,JZ)-ZMIXSS(2,IZ)*ZMIXSS(2,JZ))/4.
156 ELSEIF(IZ.GT.JZ) THEN
157 WIJ(IZ,JZ)=-SQRT(G**2+GP**2)
158 $ *(ZMIXSS(1,IZ)*ZMIXSS(1,JZ)-ZMIXSS(2,IZ)*ZMIXSS(2,JZ))/4.
159 ENDIF
160110 CONTINUE
161C
162 AUWI(2)=G*THX*COS(GAMMAL)
163 ADWI(2)=SNW2*G*THY*COS(GAMMAR)
164 ALWI(2)=ADWI(2)
165 ANWI(2)=AUWI(2)
166 AUWI(1)=G*SIN(GAMMAL)
167 ADWI(1)=SNW1*G*SIN(GAMMAR)
168 ALWI(1)=ADWI(1)
169 ANWI(1)=AUWI(1)
170C
171 DO 120 IZ=1,4
172 XIM(IZ)=.5*(SNW1*SIGN(1.,AMZISS(IZ))*(COS(GAMMAR)
173 $ *ZMIXSS(1,IZ)/SR2+SIN(GAMMAR)*ZMIXSS(3,IZ))-COS(GAMMAL)
174 $ *ZMIXSS(2,IZ)/SR2+SIN(GAMMAL)*ZMIXSS(3,IZ))
175 YIM(IZ)=.5*(-SNW1*SIGN(1.,AMZISS(IZ))*(COS(GAMMAR)
176 $ *ZMIXSS(1,IZ)/SR2+SIN(GAMMAR)*ZMIXSS(3,IZ))-COS(GAMMAL)
177 $ *ZMIXSS(2,IZ)/SR2+SIN(GAMMAL)*ZMIXSS(3,IZ))
178 XIP(IZ)=.5*(SNW2*SIGN(1.,AMZISS(IZ))*THY*(-SIN(GAMMAR)
179 $ *ZMIXSS(1,IZ)/SR2+COS(GAMMAR)*ZMIXSS(3,IZ))+THX*(SIN(GAMMAL)
180 $ *ZMIXSS(2,IZ)/SR2+COS(GAMMAL)*ZMIXSS(3,IZ)))
181 YIP(IZ)=.5*(-SNW2*SIGN(1.,AMZISS(IZ))*THY*(-SIN(GAMMAR)
182 $ *ZMIXSS(1,IZ)/SR2+COS(GAMMAR)*ZMIXSS(3,IZ))+THX*(SIN(GAMMAL)
183 $ *ZMIXSS(2,IZ)/SR2+COS(GAMMAL)*ZMIXSS(3,IZ)))
184120 CONTINUE
185C
186 DO 130 IZ=1,4
187 DO 130 JZ=1,4
188 IF(IZ.NE.JZ) THEN
189 SNIJ(IZ,JZ)=-1.*SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
190 XLIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
191 $ *(ZMIXSS(2,IZ)*SINA-ZMIXSS(1,IZ)*COSA)
192 $ *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
193 XHIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
194 $ *(ZMIXSS(2,IZ)*COSA+ZMIXSS(1,IZ)*SINA)
195 $ *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
196 XPIJ(IZ,JZ)=SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
197 $ *(ZMIXSS(2,IZ)*SINBE-ZMIXSS(1,IZ)*COSBE)
198 $ *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
199 HIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ))
200 $ *(ZMIXSS(2,IZ)*COSA+ZMIXSS(1,IZ)*SINA)
201 $ *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2.
202 ENDIF
203130 CONTINUE
204C
205 SP=-.5*(-THY*SNW2*COSBE*SIN(GAMMAL)*SIN(GAMMAR)+
206 $THY*SNW2*SINBE*COS(GAMMAL)*COS(GAMMAR)-
207 $THX*SNW1*COSBE*COS(GAMMAL)*COS(GAMMAR)+
208 $THX*SNW1*SINBE*SIN(GAMMAL)*SIN(GAMMAR))
209 PP=-.5*(-THY*SNW2*COSBE*SIN(GAMMAL)*SIN(GAMMAR)+
210 $THY*SNW2*SINBE*COS(GAMMAL)*COS(GAMMAR)+
211 $THX*SNW1*COSBE*COS(GAMMAL)*COS(GAMMAR)-
212 $THX*SNW1*SINBE*SIN(GAMMAL)*SIN(GAMMAR))
213C
214 SL=.5*(THY*SNW2*SINA*COS(GAMMAL)*COS(GAMMAR)-
215 $THY*SNW2*COSA*SIN(GAMMAL)*SIN(GAMMAR)+
216 $THX*SNW1*COSA*COS(GAMMAL)*COS(GAMMAR)-
217 $THX*SNW1*SINA*SIN(GAMMAL)*SIN(GAMMAR))
218 PL=.5*(THY*SNW2*SINA*COS(GAMMAL)*COS(GAMMAR)-
219 $THY*SNW2*COSA*SIN(GAMMAL)*SIN(GAMMAR)-
220 $THX*SNW1*COSA*COS(GAMMAL)*COS(GAMMAR)+
221 $THX*SNW1*SINA*SIN(GAMMAL)*SIN(GAMMAR))
222C
223 SH=.5*(THY*SNW2*COSA*COS(GAMMAL)*COS(GAMMAR)+
224 $THY*SNW2*SINA*SIN(GAMMAL)*SIN(GAMMAR)-
225 $THX*SNW1*SINA*COS(GAMMAL)*COS(GAMMAR)-
226 $THX*SNW1*COSA*SIN(GAMMAL)*SIN(GAMMAR))
227 PH=.5*(THY*SNW2*COSA*COS(GAMMAL)*COS(GAMMAR)+
228 $THY*SNW2*SINA*SIN(GAMMAL)*SIN(GAMMAR)+
229 $THX*SNW1*SINA*COS(GAMMAL)*COS(GAMMAR)+
230 $THX*SNW1*COSA*SIN(GAMMAL)*SIN(GAMMAR))
231C
232 DO 140 IZ=1,4
233 V1I(IZ)=-SIN(GAMMAR)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
234 $ -G*COS(GAMMAR)*ZMIXSS(1,IZ)
235 V2I(IZ)=COS(GAMMAR)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
236 $ -G*SIN(GAMMAR)*ZMIXSS(1,IZ)
237 V3I(IZ)=-SIN(GAMMAL)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
238 $ +G*COS(GAMMAL)*ZMIXSS(2,IZ)
239 V4I(IZ)=COS(GAMMAL)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))
240 $ +G*SIN(GAMMAL)*ZMIXSS(2,IZ)
241140 CONTINUE
242C
243C-----------------------------------------------------------------------
244C Generate Chargino Branching Fractions
245C-----------------------------------------------------------------------
246 ISZ(1)=ISZ1
247 ISZ(2)=ISZ2
248 ISZ(3)=ISZ3
249 ISZ(4)=ISZ4
250C FIRST TRY EXCLUSIVE DECAY TO SINGLE PION
251 IF (MW1.GT.(MZ1+AMPI).AND.MW1.LT.(MZ1+1.)) THEN
252 WID=G**4*FPI**2*SQRT(SSXLAM(MW1**2,MZ1**2,AMPI**2))/
253 $ 128./MW1**3/PI/AMW**4*((XIM(1)**2+YIM(1)**2)*
254 $ (MW1-MZ1)**2*(MW1+MZ1)**2-AMPI**2*(XIM(1)**2*(MW1-MZ1)**2+
255 $ YIM(1)**2*(MW1+MZ1)**2))
256 CALL SSSAVE(ISW1,WID,ISZ(1),120,0,0,0)
257 END IF
258 DO 300 IW=1,2
259C Loop over w1, w2
260 IF(IW.EQ.1) THEN
261 MWIW=MW1
262 SNIW=SNW1
263 ISWIW=ISW1
264 ELSE
265 MWIW=MW2
266 SNIW=SNW2
267 ISWIW=ISW2
268 ENDIF
269C
270C Decays to zj
271C
272 DO 310 JZ=1,4
273 MZJZ=ABS(AMZISS(JZ))
274 ISZJZ=ISZ(JZ)
275 SNJZ=SIGN(1.,AMZISS(JZ))
276 THJZ=0.
277 IF (AMZISS(JZ).LT.0.) THJZ=1.
278 IF(MWIW.LE.FUDGE*MZJZ) GOTO 310
279C Couplings
280 IF(IW.EQ.1) THEN
281 XIPM=XIM(JZ)
282 YIPM=YIM(JZ)
283 ELSE
284 XIPM=XIP(JZ)
285 YIPM=YIP(JZ)
286 ENDIF
287C
288C wi --> f + fbar + zj
289C
290 IF (MWIW.GT.(MZJZ+AMUP+AMDN)) THEN
291 IF (IW.EQ.1.AND.JZ.EQ.1.AND.(MW1.LT.(MZ1+1.))) THEN
292C EXIT HADRONIC MODE IF EXCLUSIVE PION DECAY ACTIVE
293 GO TO 200
294 END IF
295 IF (MWIW.LT.(AMW+MZJZ)) THEN
296 TMP(1)=XIPM**2+YIPM**2
297 TMP(2)=XIPM**2-YIPM**2
298 TMP(3)=MWIW
299 TMP(4)=MZJZ
300 TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
301 $ 2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
302 ELSE
303 TERM1=0.
304 END IF
305 ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
306 TMP(1)=MWIW
307 TMP(3)=MZJZ
308 IF (MWIW.LT.AMULSS) THEN
309 TMP(2)=AMULSS
310 PSIINT=SSXINT(0.,SSWZF2,ULIM)
311 T1=AUI(JZ)**2*ADWI(IW)**2*PSIINT
312 ELSE
313 T1=0.
314 END IF
315 IF (MWIW.LT.AMDLSS) THEN
316 TMP(2)=AMDLSS
317 PSIINT=SSXINT(0.,SSWZF2,ULIM)
318 T2=ADI(JZ)**2*AUWI(IW)**2*PSIINT
319 ELSE
320 T2=0.
321 END IF
322 IF (MWIW.LT.AMDLSS.AND.MWIW.LT.AMULSS) THEN
323 TMP(2)=0.
324 TMP(4)=AMDLSS
325 TMP(5)=AMULSS
326 PHIINT=SSXINT(0.,SSGX2,ULIM)
327 CC=2*SIGN(1.,AMZISS(JZ))*AUWI(IW)*ADWI(IW)*AUI(JZ)*ADI(JZ)
328 T3=CC*PHIINT
329 ELSE
330 T3=0.
331 END IF
332 TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
333 FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
334 TMP(1)=MWIW
335 TMP(3)=MZJZ
336 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMULSS) THEN
337 TMP(2)=AMULSS
338 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
339 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
340 TERM3=FACTOR*ADWI(IW)*AUI(JZ)*((XIPM-YIPM)*XI1
341 $ -(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
342 ELSE
343 TERM3=0.
344 END IF
345 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMDLSS) THEN
346 TMP(2)=AMDLSS
347 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
348 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
349 TERM4=-FACTOR*AUWI(IW)*ADI(JZ)*((XIPM+YIPM)*XI1
350 $ -(XIPM-YIPM)*XI2)
351 ELSE
352 TERM4=0.
353 END IF
354 WID=3*(TERM1+TERM2+TERM3+TERM4)
355 CALL SSSAVE(ISWIW,WID,ISZJZ,IDUP,-IDDN,0,0)
356C Enter information for decay matrix element
357 Z1(1)=ZI**THJZ*G*XIPM
358 Z1(2)=ZI**THJZ*G*YIPM
359 Z2(1)=G/2./SR2
360 Z2(2)=-G/2./SR2
361 CALL SSME3(1,AMW,Z1,Z2)
362 Z1(1)=ZI*AUWI(IW)/2.
363 Z1(2)=Z1(1)
364 Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ))/2.
365 Z2(2)=-Z2(1)
366 CALL SSME3(2,AMDLSS,Z1,Z2)
367 Z1(1)=CONJG(ZI*ADWI(IW))/2.
368 Z1(2)=-Z1(1)
369 Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.
370 Z2(2)=Z2(1)
371 CALL SSME3(3,AMULSS,Z1,Z2)
372 END IF
373C wi --> c + sbar + zj
374200 IF (MWIW.GT.(MZJZ+AMCH+AMST)) THEN
375 IF (MWIW.LT.(AMW+MZJZ)) THEN
376 TMP(1)=XIPM**2+YIPM**2
377 TMP(2)=XIPM**2-YIPM**2
378 TMP(3)=MWIW
379 TMP(4)=MZJZ
380 TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
381 $ 2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
382 ELSE
383 TERM1=0.
384 END IF
385 ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
386 TMP(1)=MWIW
387 TMP(3)=MZJZ
388 IF (MWIW.LT.AMCLSS) THEN
389 TMP(2)=AMCLSS
390 PSIINT=SSXINT(0.,SSWZF2,ULIM)
391 T1=AUI(JZ)**2*ADWI(IW)**2*PSIINT
392 ELSE
393 T1=0.
394 END IF
395 IF (MWIW.LT.AMSLSS) THEN
396 TMP(2)=AMSLSS
397 PSIINT=SSXINT(0.,SSWZF2,ULIM)
398 T2=ADI(JZ)**2*AUWI(IW)**2*PSIINT
399 ELSE
400 T2=0.
401 END IF
402 IF (MWIW.LT.AMSLSS.AND.MWIW.LT.AMCLSS) THEN
403 TMP(2)=0.
404 TMP(4)=AMSLSS
405 TMP(5)=AMCLSS
406 PHIINT=SSXINT(0.,SSGX2,ULIM)
407 CC=2*SIGN(1.,AMZISS(JZ))*AUWI(IW)*ADWI(IW)*AUI(JZ)*ADI(JZ)
408 T3=CC*PHIINT
409 ELSE
410 T3=0.
411 END IF
412 TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
413 FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
414 TMP(1)=MWIW
415 TMP(3)=MZJZ
416 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMCLSS) THEN
417 TMP(2)=AMCLSS
418 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
419 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
420 TERM3=FACTOR*ADWI(IW)*AUI(JZ)*((XIPM-YIPM)*XI1
421 $ -(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
422 ELSE
423 TERM3=0.
424 END IF
425 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMSLSS) THEN
426 TMP(2)=AMSLSS
427 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
428 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
429 TERM4=-FACTOR*AUWI(IW)*ADI(JZ)*((XIPM+YIPM)*XI1
430 $ -(XIPM-YIPM)*XI2)
431 ELSE
432 TERM4=0.
433 END IF
434 WID=3*(TERM1+TERM2+TERM3+TERM4)
435 CALL SSSAVE(ISWIW,WID,ISZJZ,IDCH,-IDST,0,0)
436C Enter information for decay matrix element
437 Z1(1)=ZI**THJZ*G*XIPM
438 Z1(2)=ZI**THJZ*G*YIPM
439 Z2(1)=G/2./SR2
440 Z2(2)=-G/2./SR2
441 CALL SSME3(1,AMW,Z1,Z2)
442 Z1(1)=ZI*AUWI(IW)/2.
443 Z1(2)=Z1(1)
444 Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ))/2.
445 Z2(2)=-Z2(1)
446 CALL SSME3(2,AMSLSS,Z1,Z2)
447 Z1(1)=CONJG(ZI*ADWI(IW))/2.
448 Z1(2)=-Z1(1)
449 Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.
450 Z2(2)=Z2(1)
451 CALL SSME3(3,AMCLSS,Z1,Z2)
452 END IF
453C wi -> t + bbar + zj neglected since 2-body modes should dominate
454C wi --> nu_e + e + zj
455 IF (MWIW.GT.(MZJZ+AME)) THEN
456 IF (MWIW.LT.(AMW+MZJZ)) THEN
457 TMP(1)=XIPM**2+YIPM**2
458 TMP(2)=XIPM**2-YIPM**2
459 TMP(3)=MWIW
460 TMP(4)=MZJZ
461 TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
462 $ 2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
463 ELSE
464 TERM1=0.
465 END IF
466 ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
467 TMP(1)=MWIW
468 TMP(3)=MZJZ
469 IF (MWIW.LT.AMN1SS) THEN
470 TMP(2)=AMN1SS
471 T1=ANI(JZ)**2*ALWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
472 ELSE
473 T1=0.
474 END IF
475 IF (MWIW.LT.AMELSS) THEN
476 TMP(2)=AMELSS
477 T2=ALI(JZ)**2*ANWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
478 ELSE
479 T2=0.
480 END IF
481 IF (MWIW.LT.AMELSS.AND.MWIW.LT.AMN1SS) THEN
482 TMP(2)=0.
483 TMP(4)=AMELSS
484 TMP(5)=AMN1SS
485 T3=-2*SIGN(1.,AMZISS(JZ))*ANWI(IW)*ALWI(IW)*ANI(JZ)*
486 $ ALI(JZ)*SSXINT(0.,SSGX2,ULIM)
487 ELSE
488 T3=0.
489 END IF
490 TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
491 FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
492 TMP(1)=MWIW
493 TMP(3)=MZJZ
494 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMN1SS) THEN
495 TMP(2)=AMN1SS
496 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
497 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
498 TERM3=FACTOR*ALWI(IW)*ANI(JZ)*((XIPM-YIPM)*XI1
499 $ -(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
500 ELSE
501 TERM3=0.
502 END IF
503 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMELSS) THEN
504 TMP(2)=AMELSS
505 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
506 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
507 TERM4=FACTOR*ANWI(IW)*ALI(JZ)*((XIPM+YIPM)*XI1
508 $ -(XIPM-YIPM)*XI2)
509 ELSE
510 TERM4=0.
511 END IF
512 WID=TERM1+TERM2+TERM3+TERM4
513 CALL SSSAVE(ISWIW,WID,ISZJZ,-IDE,IDNE,0,0)
514C Enter information for decay matrix element
515 Z1(1)=ZI**THJZ*G*XIPM
516 Z1(2)=ZI**THJZ*G*YIPM
517 Z2(1)=G/2./SR2
518 Z2(2)=-G/2./SR2
519 CALL SSME3(1,AMW,Z1,Z2)
520 Z1(1)=ZI*ANWI(IW)/2.
521 Z1(2)=Z1(1)
522 Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ))/2.
523 Z2(2)=-Z2(1)
524 CALL SSME3(2,AMELSS,Z1,Z2)
525 Z1(1)=CONJG(ZI*ALWI(IW))/2.
526 Z1(2)=-Z1(1)
527 Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.
528 Z2(2)=Z2(1)
529 CALL SSME3(3,AMN1SS,Z1,Z2)
530 END IF
531C wi --> nu_mu + mu + zj
532 IF (MWIW.GT.(MZJZ+AMMU)) THEN
533 IF (MWIW.LT.(AMW+MZJZ)) THEN
534 TMP(1)=XIPM**2+YIPM**2
535 TMP(2)=XIPM**2-YIPM**2
536 TMP(3)=MWIW
537 TMP(4)=MZJZ
538 TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
539 $ 2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
540 ELSE
541 TERM1=0.
542 END IF
543 ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
544 TMP(1)=MWIW
545 TMP(3)=MZJZ
546 IF (MWIW.LT.AMN2SS) THEN
547 TMP(2)=AMN2SS
548 T1=ANI(JZ)**2*ALWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
549 ELSE
550 T1=0.
551 END IF
552 IF (MWIW.LT.AMMLSS) THEN
553 TMP(2)=AMMLSS
554 T2=ALI(JZ)**2*ANWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
555 ELSE
556 T2=0.
557 END IF
558 IF (MWIW.LT.AMMLSS.AND.MWIW.LT.AMN2SS) THEN
559 TMP(2)=0.
560 TMP(4)=AMMLSS
561 TMP(5)=AMN2SS
562 T3=-2*SIGN(1.,AMZISS(JZ))*ANWI(IW)*ALWI(IW)*ANI(JZ)*
563 $ ALI(JZ)*SSXINT(0.,SSGX2,ULIM)
564 ELSE
565 T3=0.
566 END IF
567 TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
568 FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
569 TMP(1)=MWIW
570 TMP(3)=MZJZ
571 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMN2SS) THEN
572 TMP(2)=AMN2SS
573 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
574 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
575 TERM3=FACTOR*ALWI(IW)*ANI(JZ)*((XIPM-YIPM)*XI1
576 $ -(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
577 ELSE
578 TERM3=0.
579 END IF
580 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMMLSS) THEN
581 TMP(2)=AMMLSS
582 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
583 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
584 TERM4=FACTOR*ANWI(IW)*ALI(JZ)*((XIPM+YIPM)*XI1
585 $ -(XIPM-YIPM)*XI2)
586 ELSE
587 TERM4=0.
588 END IF
589 WID=TERM1+TERM2+TERM3+TERM4
590 CALL SSSAVE(ISWIW,WID,ISZJZ,-IDMU,IDNM,0,0)
591C Enter information for decay matrix element
592 Z1(1)=ZI**THJZ*G*XIPM
593 Z1(2)=ZI**THJZ*G*YIPM
594 Z2(1)=G/2./SR2
595 Z2(2)=-G/2./SR2
596 CALL SSME3(1,AMW,Z1,Z2)
597 Z1(1)=ZI*ANWI(IW)/2.
598 Z1(2)=Z1(1)
599 Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ))/2.
600 Z2(2)=-Z2(1)
601 CALL SSME3(2,AMMLSS,Z1,Z2)
602 Z1(1)=CONJG(ZI*ALWI(IW))/2.
603 Z1(2)=-Z1(1)
604 Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.
605 Z2(2)=Z2(1)
606 CALL SSME3(3,AMN2SS,Z1,Z2)
607 END IF
608C wi --> nu_tau + tau + zj ; includes mixing and Yukawas
609 FACT=1./2./MWIW/(2*PI)**5
610 ALJZ1=-ALI(JZ)*COSL-FL*ZMIXSS(2,JZ)*SINL
611 ALJZ2=-ALI(JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
612 BEJZ1=BLI(JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
613 BEJZ2=-BLI(JZ)*COSL+FL*ZMIXSS(2,JZ)*SINL
614 SNJZ=SIGN(1.,AMZISS(JZ))
615C Change ALWI def'ns in accord with Drees note
616 ALWI(1)=-G*SIN(GAMMAR)
617 ALWI(2)=-G*THY*COS(GAMMAR)
618C Polarization for stau_1 -> z1ss+tau.
619C See above for other cases.
620 IF(IW.EQ.1.AND.JZ.EQ.1) THEN
621 PTAU1(JZ)=(BEJZ1**2-ALJZ1**2)/(BEJZ1**2+ALJZ1**2)
622 PTAU2(JZ)=(BEJZ2**2-ALJZ2**2)/(BEJZ2**2+ALJZ2**2)
623 ENDIF
624 IF (IW.EQ.1) THEN
625 ALIW1=-G*SIN(GAMMAL)*COSL+FL*COS(GAMMAL)*SINL
626 ALIW2=-G*SIN(GAMMAL)*SINL-FL*COS(GAMMAL)*COSL
627 AHCJZ=COSBE*V2I(JZ)
628 BHCJZ=-SINBE*V4I(JZ)
629 ELSE IF (IW.EQ.2) THEN
630 ALIW1=(-G*COS(GAMMAL)*COSL-FL*SIN(GAMMAL)*SINL)*THX
631 ALIW2=(-G*COS(GAMMAL)*SINL+FL*SIN(GAMMAL)*COSL)*THX
632 AHCJZ=COSBE*V1I(JZ)*THY
633 BHCJZ=-SINBE*V3I(JZ)*THX
634 END IF
635 IF (MWIW.GT.(MZJZ+AMTAU)) THEN
636 IF (MWIW.LT.(AMW+MZJZ)) THEN
637 TMP(1)=XIPM**2+YIPM**2
638 TMP(2)=XIPM**2-YIPM**2
639 TMP(3)=MWIW
640 TMP(4)=MZJZ
641 TERMW=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)*
642 $ FACT*2*G**4*PI**2/3.
643 ELSE
644 TERMW=0.
645 END IF
646 ULIM=(MWIW**2-MZJZ**2)/2./MWIW
647 TMP(1)=MWIW
648 TMP(3)=MZJZ
649 IF (MWIW.LT.AMN3SS) THEN
650 TMP(2)=AMN3SS
651 POLNL=FACT*ANI(JZ)**2*ALWI(IW)**2*
652 $ SSXINT(0.,SSWZF2,ULIM)
653 POLNR=POLNL*BPLWI(IW)**2/ALWI(IW)**2
654 TERMN=POLNL+POLNR
655 ELSE
656 POLNL=0.
657 POLNR=0.
658 TERMN=0.
659 END IF
660 IF (MWIW.LT.AML1SS) THEN
661 TMP(2)=AML1SS
662 POL1L=FACT*ALJZ1**2*ALIW1**2*SSXINT(0.,SSWZF2,ULIM)
663 POL1R=POL1L*BEJZ1**2/ALJZ1**2
664 TERM1=POL1L+POL1R
665 ELSE
666 POL1L=0.
667 POL1R=0.
668 TERM1=0.
669 END IF
670 IF (MWIW.LT.AML2SS) THEN
671 TMP(2)=AML2SS
672 POL2L=FACT*ALJZ2**2*ALIW2**2*SSXINT(0.,SSWZF2,ULIM)
673 POL2R=POL2L*BEJZ2**2/ALJZ2**2
674 TERM2=POL2L+POL2R
675 ELSE
676 POL2L=0.
677 POL2R=0.
678 TERM2=0.
679 END IF
680 IF (MWIW.LT.AML1SS) THEN
681 TMP(2)=0.
682 TMP(3)=MZJZ
683 TMP(4)=AML1SS
684 TMP(5)=AML2SS
685 POL12L=FACT*2*ALIW1*ALIW2*ALJZ1*ALJZ2*
686 $ SSXINT(0.,SSGX1,ULIM)
687 POL12R=POL12L*BEJZ1*BEJZ2/ALJZ1/ALJZ2
688 TERM12=POL12L+POL12R
689 ELSE
690 POL12L=0.
691 POL12R=0.
692 TERM12=0.
693 END IF
694 IF (MWIW.LT.(AMHC+MZJZ)) THEN
695 TMP(2)=AMHC
696 TMP(3)=MZJZ
697 TMP(4)=AHCJZ
698 TMP(5)=BHCJZ
699 TMP(6)=SIGN(1.,AMZISS(JZ))*SIGN(1.,AMWISS(IW))
700 TERMH=FACT*PI**2*MWIW*(G*MTAMZ*TANB/AMW)**2/2.*
701 $ SSXINT(MZJZ,SSWZF6,(MWIW**2+MZJZ**2)/2./MWIW)
702 ELSE
703 TERMH=0.
704 END IF
705 IF (MWIW.LT.AML1SS.AND.MWIW.LT.AMN3SS) THEN
706 TMP(2)=0.
707 TMP(3)=MZJZ
708 TMP(4)=AMN3SS
709 TMP(5)=AML1SS
710 POLN1L=+2*FACT*ANI(JZ)*ALIW1*SNJZ*SNIW*ALWI(IW)*
711 $ ALJZ1*SSXINT(0.,SSGX2,ULIM)
712 POLN1R=-2*FACT*ANI(JZ)*ALIW1*BPLWI(IW)*BEJZ1*
713 $ SSXINT(0.,SSGX8,ULIM)
714 TERMN1=POLN1L+POLN1R
715 ELSE
716 POLN1L=0.
717 POLN1R=0.
718 TERMN1=0.
719 END IF
720 IF (MWIW.LT.AML2SS.AND.MWIW.LT.AMN3SS) THEN
721 TMP(2)=0.
722 TMP(3)=MZJZ
723 TMP(4)=AMN3SS
724 TMP(5)=AML2SS
725 POLN2L=+2*FACT*ANI(JZ)*ALIW2*SNJZ*SNIW*ALWI(IW)*ALJZ2*
726 $ SSXINT(0.,SSGX2,ULIM)
727 POLN2R=-2*FACT*ANI(JZ)*ALIW2*BPLWI(IW)*BEJZ2*
728 $ SSXINT(0.,SSGX8,ULIM)
729 TERMN2=POLN2L+POLN2R
730 ELSE
731 POLN2L=0.
732 POLN2R=0.
733 TERMN2=0.
734 END IF
735 TMP(1)=MWIW
736 TMP(3)=MZJZ
737 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMN3SS) THEN
738 TMP(2)=AMN3SS
739 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
740 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
741 TERMWN=2*SR2*G**2*FACT*ALWI(IW)*ANI(JZ)*((XIPM-
742 $ YIPM)*XI1-(XIPM+YIPM)*XI2)*SIGN(1.,AMZISS(JZ))
743 ELSE
744 TERMWN=0.
745 END IF
746 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AML1SS) THEN
747 TMP(2)=AML1SS
748 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
749 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
750 TERMW1=2*SR2*G**2*FACT*ALIW1*ALJZ1*
751 $ ((XIPM+YIPM)*XI1-(XIPM-YIPM)*XI2)
752 ELSE
753 TERMW1=0.
754 END IF
755 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AML2SS) THEN
756 TMP(2)=AML2SS
757 XI1=SSXINT(0.,SSWZF4,(MWIW-MZJZ)**2)
758 XI2=SSXINT(0.,SSWZF5,(MWIW-MZJZ)**2)
759 TERMW2=2*SR2*G**2*FACT*ALIW2*ALJZ2*
760 $ ((XIPM+YIPM)*XI1-(XIPM-YIPM)*XI2)
761 ELSE
762 TERMW2=0.
763 END IF
764 TMP(2)=MZJZ
765 TMP(3)=AMHC
766 TMP(5)=AHCJZ
767 TMP(6)=BHCJZ
768 TMP(7)=SNJZ*SNIW
769 IF (MWIW.LT.(AMHC+MZJZ).AND.MWIW.LT.AML1SS) THEN
770 TMP(4)=AML1SS
771 TERMH1=PI**2/2./MWIW*FACT*SR2*ALIW1*BEJZ1*G*MTAMZ*
772 $ TANB/AMW*SSXINT(0.,SSWZF7,(MWIW-MZJZ)**2)
773 ELSE
774 TERMH1=0.
775 END IF
776 IF (MWIW.LT.(AMHC+MZJZ).AND.MWIW.LT.AML2SS) THEN
777 TMP(4)=AML2SS
778 TERMH2=PI**2/2./MWIW*FACT*SR2*ALIW2*BEJZ2*G*MTAMZ*
779 $ TANB/AMW*SSXINT(0.,SSWZF7,(MWIW-MZJZ)**2)
780 ELSE
781 TERMH2=0.
782 END IF
783 IF (MWIW.LT.(AMHC+MZJZ).AND.MWIW.LT.AMN3SS) THEN
784 TMP(4)=AMN3SS
785 TERMHN=PI**2/2./MWIW*FACT*SR2*ANI(JZ)*BPLWI(IW)*G*
786 $ MTAMZ*TANB/AMW*SSXINT(0.,SSWZF7,(MWIW-MZJZ)**2)
787 ELSE
788 TERMHN=0.
789 END IF
790 WID=TERMW+TERMN+TERM1+TERM2+TERMH+TERMWN+TERMW1+
791 $ TERMW2+TERM12+TERMN1+TERMN2+TERMH1+TERMH2+
792 $ TERMHN
793C tau polarization for 3-body w1 -> z1 tau nu
794 IF (IW.EQ.1.AND.JZ.EQ.1.AND.WID.GT.0.) THEN
795 PTAUWZ=(POLNR+POL1R+POL2R+POL12R+TERMH+POLN1R+
796 $ POLN2R+TERMHN+TERMH1+TERMH2-(TERMW+POLNL+
797 $ POL1L+POL2L+POL12L+POLN1L+POLN2L+TERMWN+
798 $ TERMW1+TERMW2))/WID
799 END IF
800 CALL SSSAVE(ISWIW,WID,ISZJZ,-IDTAU,IDNT,0,0)
801 Z1(1)=ZI**THJZ*G*XIPM
802 Z1(2)=ZI**THJZ*G*YIPM
803 Z2(1)=-G/2./SR2
804 Z2(2)=-Z2(1)
805 CALL SSME3(1,AMW,Z1,Z2)
806 Z1(1)=-(ZI)**THJZ*ALIW1/2.
807 Z1(2)=-Z1(1)
808 Z2(1)=(ALJZ1+(-1.)**THJZ*BEJZ1)/2.
809 Z2(2)=(ALJZ1-(-1.)**THJZ*BEJZ1)/2.
810 CALL SSME3(2,AML1SS,Z1,Z2)
811 Z1(1)=-(ZI)**THJZ*ALIW2/2.
812 Z1(2)=-Z1(1)
813 Z2(1)=(ALJZ2+(-1.)**THJZ*BEJZ2)/2.
814 Z2(2)=(ALJZ2-(-1.)**THJZ*BEJZ2)/2.
815 CALL SSME3(2,AML2SS,Z1,Z2)
816 Z1(1)=(-1.)**THJZ*ANI(JZ)/2.
817 Z1(2)=-Z1(1)
818 Z2(1)=(SNIW*ALWI(IW)+BPLWI(IW))/2.
819 Z2(2)=(SNIW*ALWI(IW)-BPLWI(IW))/2.
820 CALL SSME3(3,AMN3SS,Z1,Z2)
821 Z1(1)=ZI**THJZ*G*MTAMZ*TANB/SR2/AMW/2.
822 Z1(2)=-Z1(1)
823 Z2(1)=(SNIW*AHCJZ-SNJZ*BHCJZ)/2.
824 Z2(2)=(SNIW*AHCJZ+SNJZ*BHCJZ)/2.
825 CALL SSME3(4,AMHC,Z1,Z2)
826 END IF
827C
828C wi --> w + zj
829C
830 IF (MWIW.GT.(MZJZ+AMW)) THEN
831 EF=MWIW**2+MZJZ**2-AMW**2+((MWIW**2-MZJZ**2)**2
832 $ -AMW**4)/AMW/AMW
833 WID=G*G*SQRT(SSXLAM(MWIW**2,MZJZ**2,AMW**2))/32./PI/
834 $ MWIW**3*(2.*EF*(XIPM**2+YIPM**2)-12*MWIW*MZJZ*
835 $ (XIPM**2-YIPM**2))
836 CALL SSSAVE(ISWIW,WID,ISZJZ,IDW,0,0,0)
837 END IF
838C
839C wi --> h+ + zj
840C
841 IF (MWIW.GT.(MZJZ+AMHC)) THEN
842 IF(IW.EQ.1) THEN
843 A=(SNW1*COSBE*V2I(JZ)
844 $ -SIGN(1.,AMZISS(JZ))*SINBE*V4I(JZ))/2.
845 B=(SNW1*COSBE*V2I(JZ)
846 $ +SIGN(1.,AMZISS(JZ))*SINBE*V4I(JZ))/2.
847 ELSE
848 A=(THY*SNW2*COSBE*V1I(JZ)
849 $ -THX*SIGN(1.,AMZISS(JZ))*SINBE*V3I(JZ))/2.
850 B=(THY*SNW2*COSBE*V1I(JZ)
851 $ +THX*SIGN(1.,AMZISS(JZ))*SINBE*V3I(JZ))/2.
852 ENDIF
853 WID=SQRT(MWIW**4+MZJZ**4+AMHC**4-2.*(MWIW*MZJZ)**2-
854 $ 2.*(MWIW*AMHC)**2-2.*(MZJZ*AMHC)**2)/8./PI/MWIW**3*
855 $ ((A*A+B*B)*(MWIW*MWIW+MZJZ*MZJZ-AMHC*AMHC)/2.
856 $ +(A*A-B*B)*MWIW*MZJZ)
857 CALL SSSAVE(ISWIW,WID,ISZJZ,ISHC,0,0,0)
858 ENDIF
859310 CONTINUE
860C
861C wi --> quark + squark
862C
863 IF(MWIW.GT.(AMULSS+AMDN)) THEN
864 WID=3.*ADWI(IW)**2/32./PI/MWIW*(1.+AMDN**2/MWIW**2-
865 $ AMULSS**2/MWIW**2)*SQRT(SSXLAM(MWIW**2,AMDN**2,AMULSS**2))
866 CALL SSSAVE(ISWIW,WID,+ISUPL,-IDDN,0,0,0)
867 END IF
868 IF(MWIW.GT.(AMDLSS+AMUP)) THEN
869 WID=3.*AUWI(IW)**2/32./PI/MWIW*(1.+AMUP**2/MWIW**2-
870 $ AMDLSS**2/MWIW**2)*SQRT(SSXLAM(MWIW**2,AMUP**2,AMDLSS**2))
871 CALL SSSAVE(ISWIW,WID,-ISDNL,+IDUP,0,0,0)
872 END IF
873 IF(MWIW.GT.(AMCLSS+AMST)) THEN
874 WID=3.*ADWI(IW)**2/32./PI/MWIW*(1.+AMST**2/MWIW**2-
875 $ AMCLSS**2/MWIW**2)*SQRT(SSXLAM(MWIW**2,AMST**2,AMCLSS**2))
876 CALL SSSAVE(ISWIW,WID,+ISCHL,-IDST,0,0,0)
877 END IF
878 IF(MWIW.GT.(AMSLSS+AMCH)) THEN
879 WID=3.*AUWI(IW)**2/32./PI/MWIW*(1.+AMCH**2/MWIW**2-
880 $ AMCLSS**2/MWIW**2)*SQRT(SSXLAM(MWIW**2,AMCH**2,AMCLSS**2))
881 CALL SSSAVE(ISWIW,WID,-ISSTL,+IDCH,0,0,0)
882 ENDIF
883 IF(MWIW.GT.(AMT1SS+AMBT)) THEN
884 AS=(-ADWI(IW)*COST-BWI(IW)*SINT)**2
885 WID=3*((AS+BPWI(IW)**2*COST**2)*(MWIW**2+AMBT**2-AMT1SS**2)
886 $ +4*SQRT(AS)*BPWI(IW)*COST*AMBT*MWIW)/32./PI/MWIW**3*
887 $ SQRT(SSXLAM(MWIW**2,AMBT**2,AMT1SS**2))
888 CALL SSSAVE(ISWIW,WID,ISTP1,-IDBT,0,0,0)
889 ENDIF
890 IF(MWIW.GT.(AMT2SS+AMBT)) THEN
891 AS=(-ADWI(IW)*SINT+BWI(IW)*COST)**2
892 WID=3*((AS+BPWI(IW)**2*SINT**2)*(MWIW**2+AMBT**2-AMT2SS**2)
893 $ +4*SQRT(AS)*BPWI(IW)*SINT*AMBT*MWIW)/32./PI/MWIW**3*
894 $ SQRT(SSXLAM(MWIW**2,AMBT**2,AMT2SS**2))
895 CALL SSSAVE(ISWIW,WID,ISTP2,-IDBT,0,0,0)
896 ENDIF
897 IF(MWIW.GT.(AMB1SS+AMTP)) THEN
898 AS=(-AUWI(IW)*COSB-BPWI(IW)*SINB)**2
899 WID=3*((AS+BWI(IW)**2*COSB**2)*(MWIW**2+AMTP**2-AMB1SS**2)
900 $ +4*SQRT(AS)*BWI(IW)*COSB*AMTP*MWIW)/32./PI/MWIW**3*
901 $ SQRT(SSXLAM(MWIW**2,AMTP**2,AMB1SS**2))
902 CALL SSSAVE(ISWIW,WID,-ISBT1,IDTP,0,0,0)
903 ENDIF
904 IF(MWIW.GT.(AMB2SS+AMTP)) THEN
905 AS=(-AUWI(IW)*SINB+BPWI(IW)*COSB)**2
906 WID=3*((AS+BWI(IW)**2*SINB**2)*(MWIW**2+AMTP**2-AMB2SS**2)
907 $ +4*SQRT(AS)*BWI(IW)*SINB*AMTP*MWIW)/32./PI/MWIW**3*
908 $ SQRT(SSXLAM(MWIW**2,AMTP**2,AMB2SS**2))
909 CALL SSSAVE(ISWIW,WID,-ISBT2,IDTP,0,0,0)
910 ENDIF
911C
912C wi --> lepton + slepton
913C
914 IF(MWIW.GT.(AMN1SS+AME)) THEN
915 AS=(-ALWI(IW))**2
916 WID=AS*(MWIW**2+AME**2-AMN1SS**2)/32./PI/MWIW**3*
917 $ SQRT(SSXLAM(MWIW**2,AME**2,AMN1SS**2))
918 CALL SSSAVE(ISWIW,WID,ISNEL,-IDE,0,0,0)
919 END IF
920 IF(MWIW.GT.(AMN2SS+AMMU)) THEN
921 AS=(-ALWI(IW))**2
922 WID=AS*(MWIW**2+AMMU**2-AMN2SS**2)/32./PI/MWIW**3*
923 $ SQRT(SSXLAM(MWIW**2,AMMU**2,AMN2SS**2))
924 CALL SSSAVE(ISWIW,WID,ISNML,-IDMU,0,0,0)
925 END IF
926 IF(MWIW.GT.(AMN3SS+AMTAU)) THEN
927 AS=(-ALWI(IW))**2
928 WID=((AS+BPLWI(IW)**2)*(MWIW**2+AMTAU**2-AMN3SS**2)+
929 $ 4*SQRT(AS)*BPLWI(IW)*AMTAU*MWIW)/32./PI/MWIW**3*
930 $ SQRT(SSXLAM(MWIW**2,AMTAU**2,AMN3SS**2))
931 CALL SSSAVE(ISWIW,WID,ISNTL,-IDTAU,0,0,0)
932 ENDIF
933 IF(MWIW.GT.AMELSS) THEN
934 WID=ANWI(IW)**2*(MWIW**2-AMELSS**2)**2/32./PI/MWIW**3
935 CALL SSSAVE(ISWIW,WID,-ISEL,IDNE,0,0,0)
936 ENDIF
937 IF(MWIW.GT.AMMLSS) THEN
938 WID=ANWI(IW)**2*(MWIW**2-AMMLSS**2)**2/32./PI/MWIW**3
939 CALL SSSAVE(ISWIW,WID,-ISMUL,IDNM,0,0,0)
940 END IF
941 IF(MWIW.GT.AML1SS) THEN
942 AS=(-ANWI(IW)*COSL-BPLWI(IW)*SINL)**2
943 WID=AS*(MWIW**2-AML1SS**2)**2/32./PI/MWIW**3
944 CALL SSSAVE(ISWIW,WID,-ISTAU1,IDNT,0,0,0)
945 END IF
946 IF(MWIW.GT.AML2SS) THEN
947 AS=(-ANWI(IW)*SINL+BPLWI(IW)*COSL)**2
948 WID=AS*(MWIW**2-AML2SS**2)**2/32./PI/MWIW**3
949 CALL SSSAVE(ISWIW,WID,-ISTAU2,IDNT,0,0,0)
950 END IF
951300 CONTINUE
952C
953C w2 --> w1 + z
954C w2 --> w1 + f + fbar
955C
956 IF (MW2.GT.(MW1+AMZ)) THEN
957 EF=MW2**2+MW1**2-AMZ**2+((MW2**2-MW1**2)**2-
958 $ AMZ**4)/AMZ/AMZ
959 Y=(THX*SIN(GAMMAL)*COS(GAMMAL)-THY*SIN(GAMMAR)*COS(GAMMAR))/2.
960 Z=(THX*SIN(GAMMAL)*COS(GAMMAL)+THY*SIN(GAMMAR)*COS(GAMMAR))/2.
961 WID=(COTW+TANW)**2*SQRT(SSXLAM(MW2**2,MW1**2,AMZ**2))/
962 $ 128./32./MW2**3*(2*EF*(Y*Y+Z*Z)+
963 $ 12*MW2*MW1*(Y*Y-Z*Z)*SNW2*SNW1)
964 CALL SSSAVE(ISW2,WID,ISW1,IDZ,0,0,0)
965 W21ZL=0.
966 W21ZN=0.
967 W21ZU=0.
968 W21ZD=0.
969C ...w1 + f + fbar
970 ELSE
971 CONST=E**4*(COTW+TANW)**2/96./PI**3/MW2
972 UPPER=(MW2**2+MW1**2)/2./MW2
973 TMP(1)=MW2
974 TMP(2)=MW1
975 TMP(3)=AMZ
976 TMP(4)=SNW1*SNW2
977 TMP(5)=XWINO
978 TMP(6)=YWINO
979 XINTGL=SSXINT(MW1,SSWWF1,UPPER)
980 W21ZL=(APL**2+BTL**2)*CONST*XINTGL
981 W21ZN=(APN**2+BTN**2)*CONST*XINTGL
982 W21ZU=3*(APU**2+BTU**2)*CONST*XINTGL
983 W21ZD=3*(APD**2+BTD**2)*CONST*XINTGL
984 END IF
985C do w2 ->w1+q+qbar via sq'
986 TMP(1)=MW2
987 TMP(3)=MW1
988 UPPER=(MW2**2-MW1**2)/2./MW2
989 IF (MW2.LT.AMULSS) THEN
990 TMP(2)=AMULSS
991 XINTGL=SSXINT(0.,SSWZF2,UPPER)
992 W21D=3*(ADWI(2)*ADWI(1))**2*XINTGL/2./MW2/(2*PI)**5
993 ELSE
994 W21D=0.
995 END IF
996 IF (MW2.LT.AMDLSS) THEN
997 TMP(2)=AMDLSS
998 XINTGL=SSXINT(0.,SSWZF2,UPPER)
999 W21U=3*(AUWI(2)*AUWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1000 ELSE
1001 W21U=0.
1002 END IF
1003 IF (MW2.LT.AMCLSS) THEN
1004 TMP(2)=AMCLSS
1005 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1006 W21S=3*(ADWI(2)*ADWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1007 ELSE
1008 W21S=0.
1009 END IF
1010 IF (MW2.LT.AMSLSS) THEN
1011 TMP(2)=AMSLSS
1012 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1013 W21C=3*(AUWI(2)*AUWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1014 ELSE
1015 W21C=0.
1016 END IF
1017 IF (MW2.LT.AMN1SS) THEN
1018 TMP(2)=AMN1SS
1019 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1020 W21E=(ALWI(2)*ALWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1021 ELSE
1022 W21E=0.
1023 END IF
1024 IF (MW2.LT.AMN2SS) THEN
1025 TMP(2)=AMN2SS
1026 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1027 W21M=(ALWI(2)*ALWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1028 ELSE
1029 W21M=0.
1030 END IF
1031 IF (MW2.LT.AMN3SS) THEN
1032 TMP(2)=AMN3SS
1033 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1034 W21L=(ALWI(2)*ALWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1035 ELSE
1036 W21L=0.
1037 END IF
1038 IF (MW2.LT.AMELSS) THEN
1039 TMP(2)=AMELSS
1040 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1041 W21N1=(ANWI(2)*ANWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1042 ELSE
1043 W21N1=0.
1044 END IF
1045 IF (MW2.LT.AMMLSS) THEN
1046 TMP(2)=AMMLSS
1047 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1048 W21N2=(ANWI(2)*ANWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1049 ELSE
1050 W21N2=0.
1051 END IF
1052C !!! W2->W1+NU_TAU+NU_TAUBAR NEEDS UPDATING FOR STAU MIXING
1053 IF (MW2.LT.AML1SS) THEN
1054 TMP(2)=AML1SS
1055 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1056 W21N3=(ANWI(2)*ANWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1057 ELSE
1058 W21N3=0.
1059 END IF
1060C-----WINO-2 ->WINO-1 +BBBAR NEEDS UPDATE -------------------
1061C-----WINO-2 ->WINO-1 +TTBAR NEEDS UPDATE -------------------
1062C-----THESE ALL LACK INTERFERENCE TERMS AS WELL
1063 W21D=W21D+W21ZD
1064 W21U=W21U+W21ZU
1065 W21S=W21S+W21ZD
1066 W21C=W21C+W21ZU
1067 W21N1=W21N1+W21ZN
1068 W21N2=W21N2+W21ZN
1069 W21N3=W21N3+W21ZN
1070 W21E=W21E+W21ZL
1071 W21M=W21M+W21ZL
1072 W21L=W21L+W21ZL
1073 IF(W21D.GT.0.) THEN
1074 CALL SSSAVE(ISW2,W21D,ISW1,IDDN,-IDDN,0,0)
1075 END IF
1076 IF(MW2.GT.(MW1+2*AMST+1.)) THEN
1077 CALL SSSAVE(ISW2,W21S,ISW1,IDST,-IDST,0,0)
1078 END IF
1079C IF (MW2.GT.(MW1+2*AMBT+2.)) THEN
1080C CALL SSSAVE(ISW2,W21D,ISW1,IDBT,-IDBT,0,0)
1081C END IF
1082 IF(W21U.GT.0.) THEN
1083 CALL SSSAVE(ISW2,W21U,ISW1,IDUP,-IDUP,0,0)
1084 END IF
1085 IF (MW2.GT.(MW1+2*AMCH+1.)) THEN
1086 CALL SSSAVE(ISW2,W21C,ISW1,IDCH,-IDCH,0,0)
1087 ENDIF
1088C IF (MW2.GT.(MW1+2*AMTP+2.)) THEN
1089C CALL SSSAVE(ISW2,W21U,ISW1,IDTP,-IDTP,0,0)
1090C END IF
1091 IF(W21N1.GT.0.) THEN
1092 CALL SSSAVE(ISW2,W21N1,ISW1,IDNE,-IDNE,0,0)
1093 ENDIF
1094 IF(W21N2.GT.0.) THEN
1095 CALL SSSAVE(ISW2,W21N2,ISW1,IDNM,-IDNM,0,0)
1096 ENDIF
1097 IF(W21N3.GT.0.) THEN
1098 CALL SSSAVE(ISW2,W21N3,ISW1,IDNT,-IDNT,0,0)
1099 ENDIF
1100 IF(W21E.GT.0.) THEN
1101 CALL SSSAVE(ISW2,W21E,ISW1,IDE,-IDE,0,0)
1102 ENDIF
1103 IF(MW2.GT.(MW1+2*AMMU+1.)) THEN
1104 CALL SSSAVE(ISW2,W21M,ISW1,IDMU,-IDMU,0,0)
1105 ENDIF
1106 IF(MW2.GT.(MW1+2*AMTAU+1.)) THEN
1107 CALL SSSAVE(ISW2,W21L,ISW1,IDTAU,-IDTAU,0,0)
1108 ENDIF
1109C
1110C w2 --> w1 + higgs
1111C
1112C w2 --> w1 + hl
1113 IF(MW2.GT.(MW1+AMHL)) THEN
1114 WID=G*G/16./PI/MW2**3*SQRT(SSXLAM(MW2**2,MW1**2,AMHL**2))*
1115 $ ((SL*SL+PL*PL)*(MW2*MW2+MW1*MW1-AMHL*AMHL)/2.+
1116 $ (SL*SL-PL*PL)*MW2*MW1)
1117 CALL SSSAVE(ISW2,WID,ISW1,ISHL,0,0,0)
1118 ENDIF
1119C w2 --> w1 + hh
1120 IF(MW2.GT.(MW1+AMHH)) THEN
1121 WID=G*G/16./PI/MW2**3*SQRT(SSXLAM(MW2**2,MW1**2,AMHH**2))*
1122 $ ((SH*SH+PH*PH)*(MW2*MW2+MW1*MW1-AMHH*AMHH)/2.+
1123 $ (SH*SH-PH*PH)*MW2*MW1)
1124 CALL SSSAVE(ISW2,WID,ISW1,ISHH,0,0,0)
1125 ENDIF
1126C w2 --> w1 + ha
1127 IF(MW2.GT.(MW1+AMHA)) THEN
1128 WID=G*G/16./PI/MW2**3*SQRT(SSXLAM(MW2**2,MW1**2,AMHA**2))*
1129 $ ((SP*SP+PP*PP)*(MW2*MW2+MW1*MW1-AMHA*AMHA)/2.+
1130 $ (SP*SP-PP*PP)*MW2*MW1)
1131 CALL SSSAVE(ISW2,WID,ISW1,ISHA,0,0,0)
1132 END IF
1133C
1134C Normalize wi branching ratios
1135C
1136 CALL SSNORM(ISW1)
1137 CALL SSNORM(ISW2)
1138C
1139 RETURN
1140 END