1 #include "isajet/pilot.h"
3 C-----------------------------------------------------------------------
4 C This subroutine calculates the chargino (wi) subset of
5 C SSWZBF, which was too long.
6 C Valid for all scalar masses (functions in double precision)
7 C Includes Higgs sector radiative corrections (Aug. 31)
9 C Auxiliary functions are called SSWxyi, SSZxyi, where normally
10 C x indicates the SUSY particle, y the SM particle(s), and i is
13 C Part of Baer's GAUGBF
15 C-----------------------------------------------------------------------
16 #if defined(CERNLIB_IMPNONE)
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"
27 EXTERNAL SSZWF1,SSZZF1,SSZZF2,SSZZF3,SSWZF1,SSWZF2,SSWZF3
28 $,SSWZF4,SSWZF5,SSWZF6,SSWZF7,SSWWF1,SSZZF4,SSZZF5,SSGX1
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)
38 REAL MWIW,SL,PP,SP,PL,PH,SH,COSB,SINB,COSA,SINA
39 $,FACTOR,MZJZ,ULIM,XINTGL,COSBE,SINBE
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
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
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/
66 C Constants from neutralino mass matrix
71 G=SQRT(4*PI*ALFAEM/SN2THW)
72 GP=G*SQRT(SN2THW/(1.-SN2THW))
75 TANW=SQRT(SN2THW/(1.-SN2THW))
83 APU=-5*TANW/12.+COTW/4.
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))*
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)
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)
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))
133 C Constants from Higgs mass matrix
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)
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.
162 AUWI(2)=G*THX*COS(GAMMAL)
163 ADWI(2)=SNW2*G*THY*COS(GAMMAR)
166 AUWI(1)=G*SIN(GAMMAL)
167 ADWI(1)=SNW1*G*SIN(GAMMAR)
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)))
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.
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))
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))
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))
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)
243 C-----------------------------------------------------------------------
244 C Generate Chargino Branching Fractions
245 C-----------------------------------------------------------------------
250 C 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)
275 SNJZ=SIGN(1.,AMZISS(JZ))
277 IF (AMZISS(JZ).LT.0.) THJZ=1.
278 IF(MWIW.LE.FUDGE*MZJZ) GOTO 310
288 C wi --> f + fbar + zj
290 IF (MWIW.GT.(MZJZ+AMUP+AMDN)) THEN
291 IF (IW.EQ.1.AND.JZ.EQ.1.AND.(MW1.LT.(MZ1+1.))) THEN
292 C EXIT HADRONIC MODE IF EXCLUSIVE PION DECAY ACTIVE
295 IF (MWIW.LT.(AMW+MZJZ)) THEN
296 TMP(1)=XIPM**2+YIPM**2
297 TMP(2)=XIPM**2-YIPM**2
300 TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
301 $ 2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
305 ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
308 IF (MWIW.LT.AMULSS) THEN
310 PSIINT=SSXINT(0.,SSWZF2,ULIM)
311 T1=AUI(JZ)**2*ADWI(IW)**2*PSIINT
315 IF (MWIW.LT.AMDLSS) THEN
317 PSIINT=SSXINT(0.,SSWZF2,ULIM)
318 T2=ADI(JZ)**2*AUWI(IW)**2*PSIINT
322 IF (MWIW.LT.AMDLSS.AND.MWIW.LT.AMULSS) THEN
326 PHIINT=SSXINT(0.,SSGX2,ULIM)
327 CC=2*SIGN(1.,AMZISS(JZ))*AUWI(IW)*ADWI(IW)*AUI(JZ)*ADI(JZ)
332 TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
333 FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
336 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMULSS) THEN
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))
345 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMDLSS) THEN
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
354 WID=3*(TERM1+TERM2+TERM3+TERM4)
355 CALL SSSAVE(ISWIW,WID,ISZJZ,IDUP,-IDDN,0,0)
356 C Enter information for decay matrix element
357 Z1(1)=ZI**THJZ*G*XIPM
358 Z1(2)=ZI**THJZ*G*YIPM
361 CALL SSME3(1,AMW,Z1,Z2)
364 Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ))/2.
366 CALL SSME3(2,AMDLSS,Z1,Z2)
367 Z1(1)=CONJG(ZI*ADWI(IW))/2.
369 Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.
371 CALL SSME3(3,AMULSS,Z1,Z2)
373 C wi --> c + sbar + zj
374 200 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
380 TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
381 $ 2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
385 ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
388 IF (MWIW.LT.AMCLSS) THEN
390 PSIINT=SSXINT(0.,SSWZF2,ULIM)
391 T1=AUI(JZ)**2*ADWI(IW)**2*PSIINT
395 IF (MWIW.LT.AMSLSS) THEN
397 PSIINT=SSXINT(0.,SSWZF2,ULIM)
398 T2=ADI(JZ)**2*AUWI(IW)**2*PSIINT
402 IF (MWIW.LT.AMSLSS.AND.MWIW.LT.AMCLSS) THEN
406 PHIINT=SSXINT(0.,SSGX2,ULIM)
407 CC=2*SIGN(1.,AMZISS(JZ))*AUWI(IW)*ADWI(IW)*AUI(JZ)*ADI(JZ)
412 TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
413 FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
416 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMCLSS) THEN
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))
425 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMSLSS) THEN
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
434 WID=3*(TERM1+TERM2+TERM3+TERM4)
435 CALL SSSAVE(ISWIW,WID,ISZJZ,IDCH,-IDST,0,0)
436 C Enter information for decay matrix element
437 Z1(1)=ZI**THJZ*G*XIPM
438 Z1(2)=ZI**THJZ*G*YIPM
441 CALL SSME3(1,AMW,Z1,Z2)
444 Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ))/2.
446 CALL SSME3(2,AMSLSS,Z1,Z2)
447 Z1(1)=CONJG(ZI*ADWI(IW))/2.
449 Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.
451 CALL SSME3(3,AMCLSS,Z1,Z2)
453 C wi -> t + bbar + zj neglected since 2-body modes should dominate
454 C 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
461 TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
462 $ 2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
466 ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
469 IF (MWIW.LT.AMN1SS) THEN
471 T1=ANI(JZ)**2*ALWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
475 IF (MWIW.LT.AMELSS) THEN
477 T2=ALI(JZ)**2*ANWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
481 IF (MWIW.LT.AMELSS.AND.MWIW.LT.AMN1SS) THEN
485 T3=-2*SIGN(1.,AMZISS(JZ))*ANWI(IW)*ALWI(IW)*ANI(JZ)*
486 $ ALI(JZ)*SSXINT(0.,SSGX2,ULIM)
490 TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
491 FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
494 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMN1SS) THEN
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))
503 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMELSS) THEN
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
512 WID=TERM1+TERM2+TERM3+TERM4
513 CALL SSSAVE(ISWIW,WID,ISZJZ,-IDE,IDNE,0,0)
514 C Enter information for decay matrix element
515 Z1(1)=ZI**THJZ*G*XIPM
516 Z1(2)=ZI**THJZ*G*YIPM
519 CALL SSME3(1,AMW,Z1,Z2)
522 Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ))/2.
524 CALL SSME3(2,AMELSS,Z1,Z2)
525 Z1(1)=CONJG(ZI*ALWI(IW))/2.
527 Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.
529 CALL SSME3(3,AMN1SS,Z1,Z2)
531 C 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
538 TERM1=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)/
539 $ 2./MWIW/(2*PI)**5*2*G**4*PI**2/3.
543 ULIM=MWIW/2.*(1.-MZJZ**2/MWIW**2)
546 IF (MWIW.LT.AMN2SS) THEN
548 T1=ANI(JZ)**2*ALWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
552 IF (MWIW.LT.AMMLSS) THEN
554 T2=ALI(JZ)**2*ANWI(IW)**2*SSXINT(0.,SSWZF2,ULIM)
558 IF (MWIW.LT.AMMLSS.AND.MWIW.LT.AMN2SS) THEN
562 T3=-2*SIGN(1.,AMZISS(JZ))*ANWI(IW)*ALWI(IW)*ANI(JZ)*
563 $ ALI(JZ)*SSXINT(0.,SSGX2,ULIM)
567 TERM2=(T1+T2+T3)/2./MWIW/(2*PI)**5
568 FACTOR=1./2./MWIW/(2*PI)**5*2*SR2*G**2
571 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMN2SS) THEN
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))
580 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMMLSS) THEN
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
589 WID=TERM1+TERM2+TERM3+TERM4
590 CALL SSSAVE(ISWIW,WID,ISZJZ,-IDMU,IDNM,0,0)
591 C Enter information for decay matrix element
592 Z1(1)=ZI**THJZ*G*XIPM
593 Z1(2)=ZI**THJZ*G*YIPM
596 CALL SSME3(1,AMW,Z1,Z2)
599 Z2(1)=-CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ))/2.
601 CALL SSME3(2,AMMLSS,Z1,Z2)
602 Z1(1)=CONJG(ZI*ALWI(IW))/2.
604 Z2(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.
606 CALL SSME3(3,AMN2SS,Z1,Z2)
608 C 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))
615 C Change ALWI def'ns in accord with Drees note
616 ALWI(1)=-G*SIN(GAMMAR)
617 ALWI(2)=-G*THY*COS(GAMMAR)
618 C Polarization for stau_1 -> z1ss+tau.
619 C 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)
625 ALIW1=-G*SIN(GAMMAL)*COSL+FL*COS(GAMMAL)*SINL
626 ALIW2=-G*SIN(GAMMAL)*SINL-FL*COS(GAMMAL)*COSL
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
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
641 TERMW=SSXINT(MZJZ,SSWZF1,(MWIW**2+MZJZ**2)/2./MWIW)*
642 $ FACT*2*G**4*PI**2/3.
646 ULIM=(MWIW**2-MZJZ**2)/2./MWIW
649 IF (MWIW.LT.AMN3SS) THEN
651 POLNL=FACT*ANI(JZ)**2*ALWI(IW)**2*
652 $ SSXINT(0.,SSWZF2,ULIM)
653 POLNR=POLNL*BPLWI(IW)**2/ALWI(IW)**2
660 IF (MWIW.LT.AML1SS) THEN
662 POL1L=FACT*ALJZ1**2*ALIW1**2*SSXINT(0.,SSWZF2,ULIM)
663 POL1R=POL1L*BEJZ1**2/ALJZ1**2
670 IF (MWIW.LT.AML2SS) THEN
672 POL2L=FACT*ALJZ2**2*ALIW2**2*SSXINT(0.,SSWZF2,ULIM)
673 POL2R=POL2L*BEJZ2**2/ALJZ2**2
680 IF (MWIW.LT.AML1SS) THEN
685 POL12L=FACT*2*ALIW1*ALIW2*ALJZ1*ALJZ2*
686 $ SSXINT(0.,SSGX1,ULIM)
687 POL12R=POL12L*BEJZ1*BEJZ2/ALJZ1/ALJZ2
694 IF (MWIW.LT.(AMHC+MZJZ)) THEN
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)
705 IF (MWIW.LT.AML1SS.AND.MWIW.LT.AMN3SS) THEN
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)
720 IF (MWIW.LT.AML2SS.AND.MWIW.LT.AMN3SS) THEN
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)
737 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AMN3SS) THEN
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))
746 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AML1SS) THEN
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)
755 IF (MWIW.LT.(MZJZ+AMW).AND.MWIW.LT.AML2SS) THEN
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)
769 IF (MWIW.LT.(AMHC+MZJZ).AND.MWIW.LT.AML1SS) THEN
771 TERMH1=PI**2/2./MWIW*FACT*SR2*ALIW1*BEJZ1*G*MTAMZ*
772 $ TANB/AMW*SSXINT(0.,SSWZF7,(MWIW-MZJZ)**2)
776 IF (MWIW.LT.(AMHC+MZJZ).AND.MWIW.LT.AML2SS) THEN
778 TERMH2=PI**2/2./MWIW*FACT*SR2*ALIW2*BEJZ2*G*MTAMZ*
779 $ TANB/AMW*SSXINT(0.,SSWZF7,(MWIW-MZJZ)**2)
783 IF (MWIW.LT.(AMHC+MZJZ).AND.MWIW.LT.AMN3SS) THEN
785 TERMHN=PI**2/2./MWIW*FACT*SR2*ANI(JZ)*BPLWI(IW)*G*
786 $ MTAMZ*TANB/AMW*SSXINT(0.,SSWZF7,(MWIW-MZJZ)**2)
790 WID=TERMW+TERMN+TERM1+TERM2+TERMH+TERMWN+TERMW1+
791 $ TERMW2+TERM12+TERMN1+TERMN2+TERMH1+TERMH2+
793 C 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
800 CALL SSSAVE(ISWIW,WID,ISZJZ,-IDTAU,IDNT,0,0)
801 Z1(1)=ZI**THJZ*G*XIPM
802 Z1(2)=ZI**THJZ*G*YIPM
805 CALL SSME3(1,AMW,Z1,Z2)
806 Z1(1)=-(ZI)**THJZ*ALIW1/2.
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.
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.
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.
823 Z2(1)=(SNIW*AHCJZ-SNJZ*BHCJZ)/2.
824 Z2(2)=(SNIW*AHCJZ+SNJZ*BHCJZ)/2.
825 CALL SSME3(4,AMHC,Z1,Z2)
830 IF (MWIW.GT.(MZJZ+AMW)) THEN
831 EF=MWIW**2+MZJZ**2-AMW**2+((MWIW**2-MZJZ**2)**2
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*
836 CALL SSSAVE(ISWIW,WID,ISZJZ,IDW,0,0,0)
841 IF (MWIW.GT.(MZJZ+AMHC)) 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.
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.
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)
861 C wi --> quark + squark
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)
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)
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)
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)
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)
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)
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)
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)
912 C wi --> lepton + slepton
914 IF(MWIW.GT.(AMN1SS+AME)) THEN
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)
920 IF(MWIW.GT.(AMN2SS+AMMU)) THEN
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)
926 IF(MWIW.GT.(AMN3SS+AMTAU)) THEN
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)
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)
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)
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)
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)
954 C w2 --> w1 + f + fbar
956 IF (MW2.GT.(MW1+AMZ)) THEN
957 EF=MW2**2+MW1**2-AMZ**2+((MW2**2-MW1**2)**2-
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)
971 CONST=E**4*(COTW+TANW)**2/96./PI**3/MW2
972 UPPER=(MW2**2+MW1**2)/2./MW2
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
985 C do w2 ->w1+q+qbar via sq'
988 UPPER=(MW2**2-MW1**2)/2./MW2
989 IF (MW2.LT.AMULSS) THEN
991 XINTGL=SSXINT(0.,SSWZF2,UPPER)
992 W21D=3*(ADWI(2)*ADWI(1))**2*XINTGL/2./MW2/(2*PI)**5
996 IF (MW2.LT.AMDLSS) THEN
998 XINTGL=SSXINT(0.,SSWZF2,UPPER)
999 W21U=3*(AUWI(2)*AUWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1003 IF (MW2.LT.AMCLSS) THEN
1005 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1006 W21S=3*(ADWI(2)*ADWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1010 IF (MW2.LT.AMSLSS) THEN
1012 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1013 W21C=3*(AUWI(2)*AUWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1017 IF (MW2.LT.AMN1SS) THEN
1019 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1020 W21E=(ALWI(2)*ALWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1024 IF (MW2.LT.AMN2SS) THEN
1026 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1027 W21M=(ALWI(2)*ALWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1031 IF (MW2.LT.AMN3SS) THEN
1033 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1034 W21L=(ALWI(2)*ALWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1038 IF (MW2.LT.AMELSS) THEN
1040 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1041 W21N1=(ANWI(2)*ANWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1045 IF (MW2.LT.AMMLSS) THEN
1047 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1048 W21N2=(ANWI(2)*ANWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1052 C !!! W2->W1+NU_TAU+NU_TAUBAR NEEDS UPDATING FOR STAU MIXING
1053 IF (MW2.LT.AML1SS) THEN
1055 XINTGL=SSXINT(0.,SSWZF2,UPPER)
1056 W21N3=(ANWI(2)*ANWI(1))**2*XINTGL/2./MW2/(2*PI)**5
1060 C-----WINO-2 ->WINO-1 +BBBAR NEEDS UPDATE -------------------
1061 C-----WINO-2 ->WINO-1 +TTBAR NEEDS UPDATE -------------------
1062 C-----THESE ALL LACK INTERFERENCE TERMS AS WELL
1074 CALL SSSAVE(ISW2,W21D,ISW1,IDDN,-IDDN,0,0)
1076 IF(MW2.GT.(MW1+2*AMST+1.)) THEN
1077 CALL SSSAVE(ISW2,W21S,ISW1,IDST,-IDST,0,0)
1079 C IF (MW2.GT.(MW1+2*AMBT+2.)) THEN
1080 C CALL SSSAVE(ISW2,W21D,ISW1,IDBT,-IDBT,0,0)
1083 CALL SSSAVE(ISW2,W21U,ISW1,IDUP,-IDUP,0,0)
1085 IF (MW2.GT.(MW1+2*AMCH+1.)) THEN
1086 CALL SSSAVE(ISW2,W21C,ISW1,IDCH,-IDCH,0,0)
1088 C IF (MW2.GT.(MW1+2*AMTP+2.)) THEN
1089 C CALL SSSAVE(ISW2,W21U,ISW1,IDTP,-IDTP,0,0)
1091 IF(W21N1.GT.0.) THEN
1092 CALL SSSAVE(ISW2,W21N1,ISW1,IDNE,-IDNE,0,0)
1094 IF(W21N2.GT.0.) THEN
1095 CALL SSSAVE(ISW2,W21N2,ISW1,IDNM,-IDNM,0,0)
1097 IF(W21N3.GT.0.) THEN
1098 CALL SSSAVE(ISW2,W21N3,ISW1,IDNT,-IDNT,0,0)
1101 CALL SSSAVE(ISW2,W21E,ISW1,IDE,-IDE,0,0)
1103 IF(MW2.GT.(MW1+2*AMMU+1.)) THEN
1104 CALL SSSAVE(ISW2,W21M,ISW1,IDMU,-IDMU,0,0)
1106 IF(MW2.GT.(MW1+2*AMTAU+1.)) THEN
1107 CALL SSSAVE(ISW2,W21L,ISW1,IDTAU,-IDTAU,0,0)
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)
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)
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)
1134 C Normalize wi branching ratios