1 #include "isajet/pilot.h"
3 C--------------------------------------------------------
5 C This program gives slepton branching fractions to gauginos
6 C according to Baer,Bartl,Karatas,Majerotto,Tata
7 C (Int. J. Mod. Phys. A4,4111 (1989); updated 10/21/96
10 C-------------------------------------------------------
11 #if defined(CERNLIB_IMPNONE)
14 #include "isajet/sslun.inc"
15 #include "isajet/ssmode.inc"
16 #include "isajet/sssm.inc"
17 #include "isajet/sspar.inc"
18 #include "isajet/sstype.inc"
19 #include "isajet/sstmp.inc"
21 EXTERNAL SSLRT1,SSL1ST,SSN1ST,SSSNWS
22 REAL SSXLAM,SSXINT,SSLRT1,SSN1ST,SSL1ST,SSSNWS
23 REAL WID,SNZI,THIZ,XM,YM,THX,THY,BPLWI(2),CS2THW
24 REAL VS,PI,SR2,G,GP,VP,V,TANB,ANWI(2),ALWI(2)
25 REAL ANIZ,AEIZ,BEIZ,MZIZ,SINL,COSL,BETA,FL,AMPL
26 REAL MTAMTA,MTAMB,MTAMZ,SUALFE
27 REAL MW1,MW2,TN2THW,SNW1,SNW2,AS,BS,BH,COSA,SINA,A
28 REAL XLO,SUM1,SUM2,AEJZ,TERM,WID1,WID2,SINB,COSB,COS2B,ANJZ,EMAX
29 REAL MEME,MEMB,MEMZ,MMMM,MMUMB,MMUMZ,TANTHE,TANTHM,THETAE,THETAM
30 REAL AAE,FE,AAM,FM,AME1,AML1,APE1,APL1,TM1,TM2,TM3,AMM1,APM1
34 COMPLEX ZI,ZONE,ZA,ZB,ZALIZ,ZBLIZ,ZPP,ZPM
35 DATA ZONE/(1.,0.)/,ZI/(0.,1.)/
37 C Partly duplicated from SSMASS
42 G=SQRT(4*PI*ALFAEM/SN2THW)
43 GP=G*SQRT(SN2THW/(1.-SN2THW))
46 VS=2*AMW**2/G/G/(1.+RV2V1**2)
60 MEME=AME*(1.-SUALFE(AME**2)/PI)
61 MEMB=MEME*(SUALFE(AMBT**2)/SUALFE(AME**2))**(-27./76.)
62 MEMZ=MEMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.)
63 FE=G*MEMZ/SR2/AMW/COS(BETA)
64 C SINCE A_e not defined in ISAJET, use A_tau as approximation
66 TANTHE=(AMERSS**2-MEMZ**2+AMZ**2*COS2B*(.5-SN2THW)-
67 $AMELSS**2)/MEMZ/(TWOM1*SINB/COSB+AAE)
69 MMMM=AMMU*(1.-SUALFE(AMMU**2)/PI)
70 MMUMB=MMMM*(SUALFE(AMBT**2)/SUALFE(AMMU**2))**(-27./76.)
71 MMUMZ=MMUMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.)
72 FM=G*MMUMZ/SR2/AMW/COS(BETA)
74 TANTHM=(AMMRSS**2-MMUMZ**2+AMZ**2*COS2B*(.5-SN2THW)-
75 $AMMLSS**2)/MMUMZ/(TWOM1*SINB/COSB+AAM)
77 MTAMTA=AMTAU*(1.-SUALFE(AMTAU**2)/PI)
78 MTAMB=MTAMTA*(SUALFE(AMBT**2)/SUALFE(AMTAU**2))**(-27./76.)
79 MTAMZ=MTAMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.)
80 FL=G*MTAMZ/SR2/AMW/COS(BETA)
85 BPLWI(1)=-FL*COS(GAMMAL)
86 BPLWI(2)=FL*THX*SIN(GAMMAL)
88 ALWI(1)=SNW1*G*SIN(GAMMAR)
89 ANWI(2)=G*THX*COS(GAMMAL)
90 ALWI(2)=SNW2*G*THY*COS(GAMMAR)
91 C Reconstruct masses from SSMASS
95 C Compute slepton branching fractions to zi
102 SNZI=SIGN(1.,AMZISS(IZ))
109 ANIZ=G/SR2*ZMIXSS(3,IZ)-GP/SR2*ZMIXSS(4,IZ)
110 AEIZ=G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ)
111 BEIZ=SR2*GP*ZMIXSS(4,IZ)
112 ZALIZ=ZI**(THIZ-1.)*SNZI
113 $ *(G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ))
114 ZBLIZ=-1*ZI**(THIZ-1.)*SR2*GP*ZMIXSS(4,IZ)
117 C sLEPTON --> LEPTON + zi
118 IF (AMELSS.GT.(MZIZ+AME)) THEN
119 WID=AEIZ**2*(AMELSS**2-MZIZ**2-AME**2)/AMELSS**3
120 $ /16./PI*SQRT(SSXLAM(AMELSS**2,MZIZ**2,AME**2))
121 CALL SSSAVE(ISEL,WID,ISZIZ(IZ),IDE,0,0,0)
123 IF (AMMLSS.GT.(MZIZ+AMMU)) THEN
124 WID=AEIZ**2*(AMMLSS**2-MZIZ**2-AMMU**2)/AMMLSS**3
125 $ /16./PI*SQRT(SSXLAM(AMMLSS**2,MZIZ**2,AMMU**2))
126 CALL SSSAVE(ISMUL,WID,ISZIZ(IZ),IDMU,0,0,0)
128 IF (AMN1SS.GT.MZIZ) THEN
129 WID=ANIZ**2*AMN1SS*(1.-MZIZ**2/AMN1SS**2)**2/16./PI
130 CALL SSSAVE(ISNEL,WID,ISZIZ(IZ),IDNE,0,0,0)
132 IF (AMN2SS.GT.MZIZ) THEN
133 WID=ANIZ**2*AMN2SS*(1.-MZIZ**2/AMN2SS**2)**2/16./PI
134 CALL SSSAVE(ISNML,WID,ISZIZ(IZ),IDNM,0,0,0)
136 IF (AMN3SS.GT.MZIZ) THEN
137 WID=ANIZ**2*AMN3SS*(1.-MZIZ**2/AMN3SS**2)**2/16./PI
138 CALL SSSAVE(ISNTL,WID,ISZIZ(IZ),IDNT,0,0,0)
140 IF (AMERSS.GT.(MZIZ+AME)) THEN
141 WID=BEIZ**2*(AMERSS**2-MZIZ**2-AME**2)/AMERSS**3
142 $ /16./PI*SQRT(SSXLAM(AMERSS**2,MZIZ**2,AME**2))
143 CALL SSSAVE(ISER,WID,ISZIZ(IZ),IDE,0,0,0)
145 IF (AMMRSS.GT.(MZIZ+AMMU)) THEN
146 WID=BEIZ**2*(AMMRSS**2-MZIZ**2-AMMU**2)/AMMRSS**3
147 $ /16./PI*SQRT(SSXLAM(AMMRSS**2,MZIZ**2,AMMU**2))
148 CALL SSSAVE(ISMUR,WID,ISZIZ(IZ),IDMU,0,0,0)
150 C sTAU_1 --> TAU + zi
151 IF (AML1SS.GT.(MZIZ+AMTAU)) THEN
152 ZA=((ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*COSL
153 $ -(ZI*ZBLIZ-ZPM*FL*ZMIXSS(2,IZ))*SINL)/2.
154 ZB=((-ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*COSL
155 $ -(ZI*ZBLIZ+ZPM*FL*ZMIXSS(2,IZ))*SINL)/2.
158 WID=(AS*(AML1SS**2-(AMTAU+MZIZ)**2)+BS*(AML1SS**2-
159 $ (MZIZ-AMTAU)**2))/8./PI/AML1SS**3*
160 $ SQRT(SSXLAM(AML1SS**2,MZIZ**2,AMTAU**2))
161 CALL SSSAVE(ISTAU1,WID,ISZIZ(IZ),IDTAU,0,0,0)
163 IF (AML2SS.GT.(MZIZ+AMTAU)) THEN
164 ZA=((ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*SINL
165 $ +(ZI*ZBLIZ-ZPM*FL*ZMIXSS(2,IZ))*COSL)/2.
166 ZB=((-ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*SINL
167 $ +(ZI*ZBLIZ+ZPM*FL*ZMIXSS(2,IZ))*COSL)/2.
170 WID=(AS*(AML2SS**2-(AMTAU+MZIZ)**2)+BS*(AML2SS**2-
171 $ (MZIZ-AMTAU)**2))/8./PI/AML2SS**3*
172 $ SQRT(SSXLAM(AML2SS**2,MZIZ**2,AMTAU**2))
173 CALL SSSAVE(ISTAU2,WID,ISZIZ(IZ),IDTAU,0,0,0)
177 C Compute branching fractions to wi --- theta-C = 0
179 IF (AMELSS.GT.MW1) THEN
180 WID=ANWI(1)**2*AMELSS*(1.-MW1**2/AMELSS**2)**2/16./PI
181 CALL SSSAVE(ISEL,WID,-ISW1,IDNE,0,0,0)
183 IF (AMMLSS.GT.MW1) THEN
184 WID=ANWI(1)**2*AMMLSS*(1.-MW1**2/AMMLSS**2)**2/16./PI
185 CALL SSSAVE(ISMUL,WID,-ISW1,IDNM,0,0,0)
187 IF (AML1SS.GT.MW1) THEN
188 AS=(-ANWI(1)*COSL-BPLWI(1)*SINL)**2
189 WID=AS*AML1SS*(1.-MW1**2/AML1SS**2)**2/16./PI
190 CALL SSSAVE(ISTAU1,WID,-ISW1,IDNT,0,0,0)
192 IF (AML2SS.GT.MW1) THEN
193 AS=(-ANWI(1)*SINL+BPLWI(1)*COSL)**2
194 WID=AS*AML2SS*(1.-MW1**2/AML2SS**2)**2/16./PI
195 CALL SSSAVE(ISTAU2,WID,-ISW1,IDNT,0,0,0)
198 IF (AMN1SS.GT.(MW1+AME)) THEN
199 WID=ALWI(1)**2*(AMN1SS**2-MW1**2-AME**2)*
200 $ SQRT(SSXLAM(AMN1SS**2,MW1**2,AME**2))/16./PI/AMN1SS**3
201 CALL SSSAVE(ISNEL,WID,ISW1,IDE,0,0,0)
203 IF (AMN2SS.GT.(MW1+AMMU)) THEN
204 WID=ALWI(1)**2*(AMN2SS**2-MW1**2-AMMU**2)*
205 $ SQRT(SSXLAM(AMN2SS**2,MW1**2,AMMU**2))/16./PI/AMN2SS**3
206 CALL SSSAVE(ISNML,WID,ISW1,IDMU,0,0,0)
209 IF (AMN3SS.GT.(MW1+AMTAU)) THEN
210 WID=((ALWI(1)**2+BPLWI(1)**2)*(AMN3SS**2-MW1**2-AMTAU**2)+
211 $ 4*MW1*AMTAU*BPLWI(1)*ALWI(1))*
212 $ SQRT(SSXLAM(AMN3SS**2,MW1**2,AMTAU**2))/16./PI/AMN3SS**3
213 CALL SSSAVE(ISNTL,WID,ISW1,IDTAU,0,0,0)
217 IF (AMELSS.GT.MW2) THEN
218 WID=ANWI(2)**2*AMELSS*(1.-MW2**2/AMELSS**2)**2/16./PI
219 CALL SSSAVE(ISEL,WID,-ISW2,IDNE,0,0,0)
221 IF (AMMLSS.GT.MW2) THEN
222 WID=ANWI(2)**2*AMMLSS*(1.-MW2**2/AMMLSS**2)**2/16./PI
223 CALL SSSAVE(ISMUL,WID,-ISW2,IDNM,0,0,0)
225 IF (AML1SS.GT.MW2) THEN
226 AS=(-ANWI(2)*COSL-BPLWI(2)*SINL)**2
227 WID=AS*AML1SS*(1.-MW2**2/AML1SS**2)**2/16./PI
228 CALL SSSAVE(ISTAU1,WID,-ISW2,IDNT,0,0,0)
230 IF (AML2SS.GT.MW2) THEN
231 AS=(-ANWI(2)*SINL+BPLWI(2)*COSL)**2
232 WID=AS*AML2SS*(1.-MW2**2/AML2SS**2)**2/16./PI
233 CALL SSSAVE(ISTAU2,WID,-ISW2,IDNT,0,0,0)
236 IF (AMN1SS.GT.(MW2+AME)) THEN
237 WID=ALWI(2)**2*(AMN1SS**2-MW2**2-AME**2)*
238 $ SQRT(SSXLAM(AMN1SS**2,MW2**2,AME**2))/16./PI/AMN1SS**3
239 CALL SSSAVE(ISNEL,WID,ISW2,IDE,0,0,0)
241 IF (AMN2SS.GT.(MW2+AMMU)) THEN
242 WID=ALWI(2)**2*(AMN2SS**2-MW2**2-AMMU**2)*
243 $ SQRT(SSXLAM(AMN2SS**2,MW2**2,AMMU**2))/16./PI/AMN2SS**3
244 CALL SSSAVE(ISNML,WID,ISW2,IDMU,0,0,0)
247 IF (AMN3SS.GT.(MW2+AMTAU)) THEN
248 WID=((ALWI(2)**2+BPLWI(2)**2)*(AMN3SS**2-MW2**2-AMTAU**2)+
249 $ 4*MW2*AMTAU*BPLWI(2)*ALWI(2))*
250 $ SQRT(SSXLAM(AMN3SS**2,MW2**2,AMTAU**2))/16./PI/AMN3SS**3
251 CALL SSSAVE(ISNTL,WID,ISW2,IDTAU,0,0,0)
254 C DECAYS TO VECTOR BOSONS
256 IF (AMELSS.GT.(AMN1SS+AMW)) THEN
257 WID=G*G*(SSXLAM(AMELSS**2,AMN1SS**2,AMW**2))**1.5/
258 $ 32./PI/AMELSS**3/AMW**2
259 CALL SSSAVE(ISEL,WID,-IDW,ISNEL,0,0,0)
262 IF (AMMLSS.GT.(AMN2SS+AMW)) THEN
263 WID=G*G*(SSXLAM(AMMLSS**2,AMN2SS**2,AMW**2))**1.5/
264 $ 32./PI/AMMLSS**3/AMW**2
265 CALL SSSAVE(ISMUL,WID,-IDW,ISNML,0,0,0)
268 IF (AML1SS.GT.(AMN3SS+AMW)) THEN
269 WID=G*G*(SSXLAM(AML1SS**2,AMN3SS**2,AMW**2))**1.5/
270 $ 32./PI/AML1SS**3/AMW**2*COSL**2
271 CALL SSSAVE(ISTAU1,WID,-IDW,ISNTL,0,0,0)
274 IF (AML2SS.GT.(AMN3SS+AMW)) THEN
275 WID=G*G*(SSXLAM(AML2SS**2,AMN3SS**2,AMW**2))**1.5/
276 $ 32./PI/AML2SS**3/AMW**2*SINL**2
277 CALL SSSAVE(ISTAU2,WID,-IDW,ISNTL,0,0,0)
280 IF (AMN1SS.GT.(AMELSS+AMW)) THEN
281 WID=G*G*(SSXLAM(AMN1SS**2,AMELSS**2,AMW**2))**1.5/
282 $ 32./PI/AMN1SS**3/AMW**2
283 CALL SSSAVE(ISNEL,WID,IDW,ISEL,0,0,0)
286 IF (AMN2SS.GT.(AMMLSS+AMW)) THEN
287 WID=G*G*(SSXLAM(AMN2SS**2,AMMLSS**2,AMW**2))**1.5/
288 $ 32./PI/AMN2SS**3/AMW**2
289 CALL SSSAVE(ISNML,WID,IDW,ISMUL,0,0,0)
292 IF (AMN3SS.GT.(AML1SS+AMW)) THEN
293 WID=G*G*(SSXLAM(AMN3SS**2,AML1SS**2,AMW**2))**1.5/
294 $ 32./PI/AMN3SS**3/AMW**2*COSL**2
295 CALL SSSAVE(ISNTL,WID,IDW,ISTAU1,0,0,0)
298 IF (AMN3SS.GT.(AML2SS+AMW)) THEN
299 WID=G*G*(SSXLAM(AMN3SS**2,AML2SS**2,AMW**2))**1.5/
300 $ 32./PI/AMN3SS**3/AMW**2*SINL**2
301 CALL SSSAVE(ISNTL,WID,IDW,ISTAU2,0,0,0)
304 IF (AML2SS.GT.(AML1SS+AMZ)) THEN
305 WID=G*G*(SSXLAM(AML2SS**2,AML1SS**2,AMZ**2))**1.5/
306 $ 64./PI/AML2SS**3/CS2THW/AMZ**2*SINL**2*COSL**2
307 CALL SSSAVE(ISTAU2,WID,IDZ,ISTAU1,0,0,0)
310 C 3-body decay of l_R -> l+tau+stau_1
312 XLO=(AML1SS+AMTAU)**2
315 IF (AMERSS.GT.(AML1SS+AMTAU+AME)) THEN
319 TMP(2)=-SR2*GP*ZMIXSS(4,IZ)
320 TMP(3)=-SR2*GP*ZMIXSS(4,JZ)
321 AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
322 AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
323 TMP(4)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
324 TMP(5)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
325 TMP(6)=TMP(2)*SINL+FL*ZMIXSS(2,IZ)*COSL
326 TMP(7)=TMP(3)*SINL+FL*ZMIXSS(2,JZ)*COSL
329 IF (AMERSS.LT.ABS(AMZISS(IZ)).AND.
330 , AMERSS.LT.ABS(AMZISS(JZ))) THEN
331 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMERSS**2)/AMERSS
335 IF (IZ.EQ.JZ) TERM=TERM/2.
339 TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
340 TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
341 IF (AMERSS.LT.ABS(AMZISS(IZ)).AND.
342 , AMERSS.LT.ABS(AMZISS(JZ))) THEN
343 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMERSS**2)/AMERSS
347 IF (IZ.EQ.JZ) TERM=TERM/2.
351 WID1=SUM1/2./AMERSS/(2*PI)**5
352 WID2=SUM2/2./AMERSS/(2*PI)**5
353 CALL SSSAVE(ISER,WID1,ISTAU1,IDE,-IDTAU,0,0)
354 CALL SSSAVE(ISER,WID2,-ISTAU1,IDE,IDTAU,0,0)
358 IF (AMMRSS.GT.(AML1SS+AMTAU+AMMU)) THEN
362 TMP(2)=-SR2*GP*ZMIXSS(4,IZ)
363 TMP(3)=-SR2*GP*ZMIXSS(4,JZ)
364 AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
365 AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
366 TMP(4)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
367 TMP(5)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
368 TMP(6)=TMP(2)*SINL+FL*ZMIXSS(2,IZ)*COSL
369 TMP(7)=TMP(3)*SINL+FL*ZMIXSS(2,JZ)*COSL
372 IF (AMMRSS.LT.ABS(AMZISS(IZ)).AND.
373 , AMMRSS.LT.ABS(AMZISS(JZ))) THEN
374 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMMRSS**2)/AMMRSS
378 IF (IZ.EQ.JZ) TERM=TERM/2.
382 TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
383 TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
384 IF (AMMRSS.LT.ABS(AMZISS(IZ)).AND.
385 , AMMRSS.LT.ABS(AMZISS(JZ))) THEN
386 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMMRSS**2)/AMMRSS
390 IF (IZ.EQ.JZ) TERM=TERM/2.
394 WID1=SUM1/2./AMMRSS/(2*PI)**5
395 WID2=SUM2/2./AMMRSS/(2*PI)**5
396 CALL SSSAVE(ISMUR,WID1,ISTAU1,IDMU,-IDTAU,0,0)
397 CALL SSSAVE(ISMUR,WID2,-ISTAU1,IDMU,IDTAU,0,0)
400 C 3-body decay of l_L -> l+tau+stau_1
404 IF (AMELSS.GT.(AML1SS+AMTAU+AME)) THEN
408 AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
409 AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
412 TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
413 TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
414 TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
415 TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
418 IF (AMELSS.LT.ABS(AMZISS(IZ)).AND.
419 , AMELSS.LT.ABS(AMZISS(JZ))) THEN
420 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMELSS**2)/AMELSS
424 IF (IZ.EQ.JZ) TERM=TERM/2.
428 TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
429 TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
430 IF (AMELSS.LT.ABS(AMZISS(IZ)).AND.
431 , AMELSS.LT.ABS(AMZISS(JZ))) THEN
432 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMELSS**2)/AMELSS
436 IF (IZ.EQ.JZ) TERM=TERM/2.
440 WID1=SUM1/2./AMELSS/(2*PI)**5
441 WID2=SUM2/2./AMELSS/(2*PI)**5
442 CALL SSSAVE(ISEL,WID1,ISTAU1,IDE,-IDTAU,0,0)
443 CALL SSSAVE(ISEL,WID2,-ISTAU1,IDE,IDTAU,0,0)
447 IF (AMMLSS.GT.(AML1SS+AMTAU+AMMU)) THEN
451 AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
452 AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
455 TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
456 TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
457 TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
458 TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
461 IF (AMMLSS.LT.ABS(AMZISS(IZ)).AND.
462 , AMMLSS.LT.ABS(AMZISS(JZ))) THEN
463 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMMLSS**2)/AMMLSS
467 IF (IZ.EQ.JZ) TERM=TERM/2.
471 TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
472 TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
473 IF (AMMLSS.LT.ABS(AMZISS(IZ)).AND.
474 , AMMLSS.LT.ABS(AMZISS(JZ))) THEN
475 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMMLSS**2)/AMMLSS
479 IF (IZ.EQ.JZ) TERM=TERM/2.
483 WID1=SUM1/2./AMMLSS/(2*PI)**5
484 WID2=SUM2/2./AMMLSS/(2*PI)**5
485 CALL SSSAVE(ISMUL,WID1,ISTAU1,IDMU,-IDTAU,0,0)
486 CALL SSSAVE(ISMUL,WID2,-ISTAU1,IDMU,IDTAU,0,0)
489 C 3-body decay of nu_eL -> nu_e+tau+stau_1
493 IF (AMN1SS.GT.(AML1SS+AMTAU)) THEN
497 ANIZ=(G*ZMIXSS(3,IZ)-GP*ZMIXSS(4,IZ))/SR2
498 ANJZ=(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/SR2
499 AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
500 AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
503 TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
504 TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
505 TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
506 TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
509 IF (AMN1SS.LT.ABS(AMZISS(IZ)).AND.
510 , AMN1SS.LT.ABS(AMZISS(JZ))) THEN
511 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN1SS**2)/AMN1SS
515 IF (IZ.EQ.JZ) TERM=TERM/2.
519 TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
520 TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
521 IF (AMN1SS.LT.ABS(AMZISS(IZ)).AND.
522 , AMN1SS.LT.ABS(AMZISS(JZ))) THEN
523 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN1SS**2)/AMN1SS
527 IF (IZ.EQ.JZ) TERM=TERM/2.
531 WID1=SUM1/2./AMN1SS/(2*PI)**5
532 WID2=SUM2/2./AMN1SS/(2*PI)**5
533 CALL SSSAVE(ISNEL,WID1,ISTAU1,IDNE,-IDTAU,0,0)
534 CALL SSSAVE(ISNEL,WID2,-ISTAU1,IDNE,IDTAU,0,0)
537 C 3-body decay of nu_muL -> nu_mu+tau+stau_1
541 IF (AMN2SS.GT.(AML1SS+AMTAU)) THEN
545 ANIZ=(G*ZMIXSS(3,IZ)-GP*ZMIXSS(4,IZ))/SR2
546 ANJZ=(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/SR2
547 AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
548 AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
551 TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
552 TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
553 TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
554 TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
557 IF (AMN2SS.LT.ABS(AMZISS(IZ)).AND.
558 , AMN2SS.LT.ABS(AMZISS(JZ))) THEN
559 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN2SS**2)/AMN2SS
563 IF (IZ.EQ.JZ) TERM=TERM/2.
567 TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
568 TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
569 IF (AMN2SS.LT.ABS(AMZISS(IZ)).AND.
570 , AMN2SS.LT.ABS(AMZISS(JZ))) THEN
571 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN2SS**2)/AMN2SS
575 IF (IZ.EQ.JZ) TERM=TERM/2.
579 WID1=SUM1/2./AMN2SS/(2*PI)**5
580 WID2=SUM2/2./AMN2SS/(2*PI)**5
581 CALL SSSAVE(ISNML,WID1,ISTAU1,IDNM,-IDTAU,0,0)
582 CALL SSSAVE(ISNML,WID2,-ISTAU1,IDNM,IDTAU,0,0)
585 C 3-body decay of nu_tauL -> nu_tau+tau+stau_1
586 C Here, we include decay via Z_i, W_i and W*, but neglect
587 C some interference terms (Z_i-W* and W_i-W* and Z_i-W_i)
588 C until a future date. HB 8/24/98
590 AME1=-G*SIN(GAMMAL)*COS(THETAE)+FE*COS(GAMMAL)*SIN(THETAE)
591 AML1=-G*SIN(GAMMAL)*COS(THETAL)+FL*COS(GAMMAL)*SIN(THETAL)
592 APE1=-G*THX*COS(GAMMAL)*COS(THETAE)-FE*THX*SIN(GAMMAL)*
594 APL1=-G*THX*COS(GAMMAL)*COS(THETAL)-FL*THX*SIN(GAMMAL)*
598 IF (AMN3SS.GT.(AML1SS+AMTAU)) THEN
602 ANIZ=(G*ZMIXSS(3,IZ)-GP*ZMIXSS(4,IZ))/SR2
603 ANJZ=(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/SR2
604 AEIZ=-(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2
605 AEJZ=-(G*ZMIXSS(3,JZ)+GP*ZMIXSS(4,JZ))/SR2
608 TMP(6)=AEIZ*COSL-FL*ZMIXSS(2,IZ)*SINL
609 TMP(7)=AEJZ*COSL-FL*ZMIXSS(2,JZ)*SINL
610 TMP(4)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
611 TMP(5)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
614 IF (AMN3SS.LT.ABS(AMZISS(IZ)).AND.
615 , AMN3SS.LT.ABS(AMZISS(JZ))) THEN
616 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN3SS**2)/AMN3SS
620 IF (IZ.EQ.JZ) TERM=TERM/2.
624 TMP(6)=-SR2*GP*ZMIXSS(4,IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL
625 TMP(7)=-SR2*GP*ZMIXSS(4,JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL
626 IF (AMN3SS.LT.ABS(AMZISS(IZ)).AND.
627 , AMN3SS.LT.ABS(AMZISS(JZ))) THEN
628 TERM=2*PI**2*SSXINT(XLO,SSLRT1,AMN3SS**2)/AMN3SS
632 IF (IZ.EQ.JZ) TERM=TERM/2.
636 WID1=SUM1/2./AMN3SS/(2*PI)**5
637 WID2=SUM2/2./AMN3SS/(2*PI)**5
638 C-----COMPUTE SNU_TAU --> TAU +STAU_1BAR +NU_TAU via W_i DECAYS --------------
640 BPE1=FL*THX*SIN(GAMMAL)
645 IF (AMN3SS.LT.MW1) THEN
646 TM1=AML1**2*(AME1**2*MW1**2*SSXINT(AML1SS**2,SSN1ST,AMN3SS**2)
647 $+BME1**2*SSXINT(AML1SS**2,SSL1ST,AMN3SS**2))
649 TM2=2*AML1*APL1*(SNW1*SNW2*AME1*APE1*
650 $ SSXINT(AML1SS**2,SSN1ST,AMN3SS**2)+BME1*BPE1*
651 $ SSXINT(AML1SS**2,SSL1ST,AMN3SS**2))
658 IF (AMN3SS.LT.MW2) THEN
659 TM3=APL1**2*(APE1**2*MW2**2*SSXINT(AML1SS**2,SSN1ST,AMN3SS**2)
660 $+BPE1**2*SSXINT(AML1SS**2,SSL1ST,AMN3SS**2))
664 WID=PI**2*(TM1+TM2+TM3)/8./2./AMN3SS/(2.*PI)**5
665 CALL SSSAVE(ISNTL,WID1+WID,ISTAU1,IDNT,-IDTAU,0,0)
666 C-----COMPUTE SNU_TAU --> TAUBAR +STAU_1 +NU_TAU via W* DECAYS --
667 EMAX=(AML1SS**2+AMN3SS**2)/2./AMN3SS
672 IF (AMN3SS.LT.(AMW+AML1SS)) THEN
673 TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
677 WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
678 CALL SSSAVE(ISNTL,WID2+WID,-ISTAU1,IDNT,IDTAU,0,0)
680 C-----Now impose all sneutrino 3 body decays via W* ------------
685 EMAX=(AML1SS**2+AMN3SS**2)/2./AMN3SS
686 IF (AMN3SS.GT.(AML1SS+AME).AND.AMN3SS.LT.(AMW+AML1SS)) THEN
687 TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
688 WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
689 CALL SSSAVE(ISNTL,WID,ISTAU1,IDNE,-IDE,0,0)
691 IF (AMN3SS.GT.(AML1SS+AMMU).AND.AMN3SS.LT.(AMW+AML1SS)) THEN
692 TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
693 WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
694 CALL SSSAVE(ISNTL,WID,ISTAU1,IDNM,-IDMU,0,0)
696 IF (AMN3SS.GT.(AML1SS+AMUP+AMDN).AND.AMN3SS.LT.(AMW+AML1SS)) THEN
697 TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
698 WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
699 CALL SSSAVE(ISNTL,3*WID,ISTAU1,IDUP,-IDDN,0,0)
701 IF (AMN3SS.GT.(AML1SS+AMCH+AMST).AND.AMN3SS.LT.(AMW+AML1SS)) THEN
702 TM1=G**4*COS(THETAL)**2*SSXINT(AML1SS,SSSNWS,EMAX)
703 WID=2*PI**2*TM1/3./2./AMN3SS/(2.*PI)**5
704 CALL SSSAVE(ISNTL,3*WID,ISTAU1,IDCH,-IDST,0,0)
707 EMAX=(AMMLSS**2+AMN2SS**2)/2./AMN2SS
708 IF (AMN2SS.GT.(AMMLSS+AME).AND.AMN2SS.LT.(AMW+AML1SS)) THEN
709 TM1=G**4*COS(THETAL)**2*SSXINT(AMMLSS,SSSNWS,EMAX)
710 WID=2*PI**2*TM1/3./2./AMN2SS/(2.*PI)**5
711 CALL SSSAVE(ISNML,WID,ISMUL,IDNE,-IDE,0,0)
713 IF (AMN2SS.GT.(AMMLSS+AMMU).AND.AMN2SS.LT.(AMW+AML1SS)) THEN
714 TM1=G**4*COS(THETAL)**2*SSXINT(AMMLSS,SSSNWS,EMAX)
715 WID=2*PI**2*TM1/3./2./AMN2SS/(2.*PI)**5
716 CALL SSSAVE(ISNML,WID,ISMUL,IDNM,-IDMU,0,0)
718 IF (AMN2SS.GT.(AMMLSS+AMUP+AMDN).AND.AMN2SS.LT.(AMW+AML1SS)) THEN
719 TM1=G**4*COS(THETAL)**2*SSXINT(AMMLSS,SSSNWS,EMAX)
720 WID=2*PI**2*TM1/3./2./AMN2SS/(2.*PI)**5
721 CALL SSSAVE(ISNML,3*WID,ISMUL,IDUP,-IDDN,0,0)
723 IF (AMN2SS.GT.(AMMLSS+AMCH+AMST).AND.AMN2SS.LT.(AMW+AML1SS)) THEN
724 TM1=G**4*COS(THETAL)**2*SSXINT(AMMLSS,SSSNWS,EMAX)
725 WID=2*PI**2*TM1/3./2./AMN2SS/(2.*PI)**5
726 CALL SSSAVE(ISNML,3*WID,ISMUL,IDCH,-IDST,0,0)
729 EMAX=(AMELSS**2+AMN1SS**2)/2./AMN1SS
730 IF (AMN1SS.GT.(AMELSS+AME).AND.AMN1SS.LT.(AMW+AML1SS)) THEN
731 TM1=G**4*COS(THETAL)**2*SSXINT(AMELSS,SSSNWS,EMAX)
732 WID=2*PI**2*TM1/3./2./AMN1SS/(2.*PI)**5
733 CALL SSSAVE(ISNEL,WID,ISEL,IDNE,-IDE,0,0)
735 IF (AMN1SS.GT.(AMELSS+AMMU).AND.AMN1SS.LT.(AMW+AML1SS)) THEN
736 TM1=G**4*COS(THETAL)**2*SSXINT(AMELSS,SSSNWS,EMAX)
737 WID=2*PI**2*TM1/3./2./AMN1SS/(2.*PI)**5
738 CALL SSSAVE(ISNEL,WID,ISEL,IDNM,-IDMU,0,0)
740 IF (AMN1SS.GT.(AMELSS+AMUP+AMDN).AND.AMN1SS.LT.(AMW+AML1SS)) THEN
741 TM1=G**4*COS(THETAL)**2*SSXINT(AMELSS,SSSNWS,EMAX)
742 WID=2*PI**2*TM1/3./2./AMN1SS/(2.*PI)**5
743 CALL SSSAVE(ISNEL,3*WID,ISEL,IDUP,-IDDN,0,0)
745 IF (AMN1SS.GT.(AMELSS+AMCH+AMST).AND.AMN1SS.LT.(AMW+AML1SS)) THEN
746 TM1=G**4*COS(THETAL)**2*SSXINT(AMELSS,SSSNWS,EMAX)
747 WID=2*PI**2*TM1/3./2./AMN1SS/(2.*PI)**5
748 CALL SSSAVE(ISNEL,3*WID,ISEL,IDCH,-IDST,0,0)
750 C-----COMPUTE ER AND MUR--> STAU_1+NEUTRINOS DECAYS --------------
751 IF (AMERSS.LT.AMELSS.AND.AMERSS.GT.AML1SS) THEN
756 IF (AMERSS.LT.MW1) THEN
757 TM1=AML1**2*AME1**2*SSXINT(AML1SS**2,SSL1ST,AMERSS**2)
759 TM2=2*AML1*AME1*APL1*APE1*SSXINT(AML1SS**2,SSL1ST,AMERSS**2)
764 IF (AMERSS.LT.MW2) THEN
767 TM3=APL1**2*APE1**2*SSXINT(AML1SS**2,SSL1ST,AMERSS**2)
771 WID=PI**2*(TM1+TM2+TM3)/8./2./AMERSS/(2.*PI)**5
772 CALL SSSAVE(ISER,WID,ISTAU1,IDNE,-IDNT,0,0)
774 IF (AMMRSS.LT.AMMLSS.AND.AMMRSS.GT.AML1SS) THEN
775 AMM1=-G*SIN(GAMMAL)*COS(THETAM)+FM*COS(GAMMAL)*SIN(THETAM)
776 AML1=-G*SIN(GAMMAL)*COS(THETAL)+FL*COS(GAMMAL)*SIN(THETAL)
777 APM1=-G*THX*COS(GAMMAL)*COS(THETAM)-FM*THX*SIN(GAMMAL)*
779 APL1=-G*THX*COS(GAMMAL)*COS(THETAL)-FL*THX*SIN(GAMMAL)*
785 IF (AMMRSS.LT.MW1) THEN
786 TM1=AML1**2*AMM1**2*SSXINT(AML1SS**2,SSL1ST,AMMRSS**2)
788 TM2=2*AML1*AMM1*APL1*APM1*SSXINT(AML1SS**2,SSL1ST,AMMRSS**2)
793 IF (AMMRSS.LT.MW2) THEN
796 TM3=APL1**2*APM1**2*SSXINT(AML1SS**2,SSL1ST,AMMRSS**2)
800 WID=PI**2*(TM1+TM2+TM3)/8./2./AMMRSS/(2.*PI)**5
801 CALL SSSAVE(ISMUR,WID,ISTAU1,IDNM,-IDNT,0,0)
803 C-----COMPUTE SNU_E --> E +STAU_1BAR +NU_TAU DECAYS --------------
805 AML1=-G*SIN(GAMMAL)*COS(THETAL)+FL*COS(GAMMAL)*SIN(THETAL)
806 APE1=-G*THX*COS(GAMMAL)
807 APL1=-G*THX*COS(GAMMAL)*COS(THETAL)-FL*THX*SIN(GAMMAL)*
809 IF (AMN1SS.GT.(AME+AML1SS)) THEN
811 BPE1=FE*THX*SIN(GAMMAL)
816 IF (AMN1SS.LT.MW1) THEN
817 TM1=AML1**2*(AME1**2*MW1**2*SSXINT(AML1SS**2,SSN1ST,AMN1SS**2)
818 $+BME1**2*SSXINT(AML1SS**2,SSL1ST,AMN1SS**2))
820 TM2=2*AML1*APL1*(SNW1*SNW2*AME1*APE1*
821 $ SSXINT(AML1SS**2,SSN1ST,AMN1SS**2)+BME1*BPE1*
822 $ SSXINT(AML1SS**2,SSL1ST,AMN1SS**2))
829 IF (AMN1SS.LT.MW2) THEN
830 TM3=APL1**2*(APE1**2*MW2**2*SSXINT(AML1SS**2,SSN1ST,AMN1SS**2)
831 $+BPE1**2*SSXINT(AML1SS**2,SSL1ST,AMN1SS**2))
835 WID=PI**2*(TM1+TM2+TM3)/8./2./AMN1SS/(2.*PI)**5
836 CALL SSSAVE(ISNEL,WID,-ISTAU1,IDE,IDNT,0,0)
838 C-----COMPUTE SNU_M --> MU +STAU_1BAR +NU_TAU DECAYS --------------
839 IF (AMN2SS.GT.(AMMU+AML1SS)) THEN
841 BPE1=FM*THX*SIN(GAMMAL)
846 IF (AMN2SS.LT.MW1) THEN
847 TM1=AML1**2*(AME1**2*MW1**2*SSXINT(AML1SS**2,SSN1ST,AMN2SS**2)
848 $+BME1**2*SSXINT(AML1SS**2,SSL1ST,AMN2SS**2))
850 TM2=2*AML1*APL1*(SNW1*SNW2*AME1*APE1*
851 $ SSXINT(AML1SS**2,SSN1ST,AMN2SS**2)+BME1*BPE1*
852 $ SSXINT(AML1SS**2,SSL1ST,AMN2SS**2))
859 IF (AMN2SS.LT.MW2) THEN
860 TM3=APL1**2*(APE1**2*MW2**2*SSXINT(AML1SS**2,SSN1ST,AMN2SS**2)
861 $+BPE1**2*SSXINT(AML1SS**2,SSL1ST,AMN2SS**2))
865 WID=PI**2*(TM1+TM2+TM3)/8./2./AMN2SS/(2.*PI)**5
866 CALL SSSAVE(ISNML,WID,-ISTAU1,IDMU,IDNT,0,0)
870 C decay to neutral higgs bosons
872 IF (AML2SS.GT.(AMHL+AML1SS)) THEN
873 BH=G*AMW*SIN(BETA-ALFAH)*(-1.+3*TN2THW)*SINL*COSL/2.+G*
874 $ AMTAU*(TWOM1*COSA+AAL*SINA)*COS(2*THETAL)/2./AMW/COS(BETA)
875 WID=BH**2*SQRT(SSXLAM(AML2SS**2,AMHL**2,AML1SS**2))/
877 CALL SSSAVE(ISTAU2,WID,ISHL,ISTAU1,0,0,0)
880 IF (AML2SS.GT.(AMHA+AML1SS)) THEN
881 BH=G*AMTAU*(TWOM1-AAL*TANB)/2./AMW
882 WID=BH**2*SQRT(SSXLAM(AML2SS**2,AMHA**2,AML1SS**2))/
884 CALL SSSAVE(ISTAU2,WID,ISHA,ISTAU1,0,0,0)
887 IF (AML2SS.GT.(AMHH+AML1SS)) THEN
888 BH=-G*AMW*COS(BETA-ALFAH)*(-1.+3*TN2THW)*SINL*COSL/2.+G*
889 $ AMTAU*(-TWOM1*SINA+AAL*COSA)*COS(2*THETAL)/2./AMW/COS(BETA)
890 WID=BH**2*SQRT(SSXLAM(AML2SS**2,AMHH**2,AML1SS**2))/
892 CALL SSSAVE(ISTAU2,WID,ISHH,ISTAU1,0,0,0)
895 C stau_i -> H^- stau_j
897 IF (AMN3SS.GT.(AML1SS+AMHC)) THEN
898 A=G/SR2/AMW*((AMTAU**2*TANB-AMW**2*SIN(2*BETA))*COSL-
899 $ AMTAU*(TWOM1-AAL*TANB)*SINL)
900 WID=A*A*SQRT(SSXLAM(AMN3SS**2,AML1SS**2,AMHC**2))/
902 CALL SSSAVE(ISNTL,WID,ISHC,ISTAU1,0,0,0)
905 IF (AMN3SS.GT.(AML2SS+AMHC)) THEN
906 A=G/SR2/AMW*((AMTAU**2*TANB-AMW**2*SIN(2*BETA))*SINL+
907 $ AMTAU*(TWOM1-AAL*TANB)*COSL)
908 WID=A*A*SQRT(SSXLAM(AMN3SS**2,AML2SS**2,AMHC**2))/
910 CALL SSSAVE(ISNTL,WID,ISHC,ISTAU2,0,0,0)
913 IF (AML1SS.GT.(AMN3SS+AMHC)) THEN
914 A=G/SR2/AMW*((AMTAU**2*TANB-AMW**2*SIN(2*BETA))*COSL-
915 $ AMTAU*(TWOM1-AAL*TANB)*SINL)
916 WID=A*A*SQRT(SSXLAM(AML1SS**2,AMN3SS**2,AMHC**2))/
918 CALL SSSAVE(ISTAU1,WID,-ISHC,ISNTL,0,0,0)
921 IF (AML2SS.GT.(AMN3SS+AMHC)) THEN
922 A=G/SR2/AMW*((AMTAU**2*TANB-AMW**2*SIN(2*BETA))*SINL+
923 $ AMTAU*(TWOM1-AAL*TANB)*COSL)
924 WID=A*A*SQRT(SSXLAM(AML2SS**2,AMN3SS**2,AMHC**2))/
926 CALL SSSAVE(ISTAU2,WID,-ISHC,ISNTL,0,0,0)
929 C Gauge mediated SUSY breaking model decays to gravitino
930 IF (AMERSS.GT.(AME+AMGVSS)) THEN
931 WID=(AMERSS**2-AME**2)**4/48./PI/AMERSS**3/(AMPL*AMGVSS)**2
932 CALL SSSAVE(ISER,WID,IDE,91,0,0,0)
934 IF (AMMRSS.GT.(AMMU+AMGVSS)) THEN
935 WID=(AMMRSS**2-AMMU**2)**4/48./PI/AMMRSS**3/(AMPL*AMGVSS)**2
936 CALL SSSAVE(ISMUR,WID,IDMU,91,0,0,0)
938 IF (AML1SS.GT.(AMTAU+AMGVSS)) THEN
939 WID=(AML1SS**2-AMTAU**2)**4/48./PI/AML1SS**3/(AMPL*AMGVSS)**2
940 CALL SSSAVE(ISTAU1,WID,IDTAU,91,0,0,0)
942 IF (AMELSS.GT.(AME+AMGVSS)) THEN
943 WID=(AMELSS**2-AME**2)**4/48./PI/AMELSS**3/(AMPL*AMGVSS)**2
944 CALL SSSAVE(ISEL,WID,IDE,91,0,0,0)
946 IF (AMMLSS.GT.(AMMU+AMGVSS)) THEN
947 WID=(AMMLSS**2-AMMU**2)**4/48./PI/AMMLSS**3/(AMPL*AMGVSS)**2
948 CALL SSSAVE(ISMUL,WID,IDMU,91,0,0,0)
950 IF (AMN1SS.GT.AMGVSS) THEN
951 WID=AMN1SS**5/48./PI/(AMPL*AMGVSS)**2
952 CALL SSSAVE(ISNEL,WID,IDNE,91,0,0,0)
954 IF (AMN2SS.GT.AMGVSS) THEN
955 WID=AMN2SS**5/48./PI/(AMPL*AMGVSS)**2
956 CALL SSSAVE(ISNML,WID,IDNM,91,0,0,0)
958 IF (AMN3SS.GT.AMGVSS) THEN
959 WID=AMN3SS**5/48./PI/(AMPL*AMGVSS)**2
960 CALL SSSAVE(ISNTL,WID,IDNT,91,0,0,0)
963 C Normalize branching ratios