]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ISAJET/isasusy/sshsf.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Adding MUON HLT code to the repository.
[u/mrichter/AliRoot.git] / ISAJET / isasusy / sshsf.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE SSHSF
3C-----------------------------------------------------------------------
4C
5C Calculates the partial decay widths of
6C the Higgs bosons into sfermions.
7C calculated by X. Tata
8C program by M. Bisset
9C
10C 10/23/93: modified by H. Baer, 10/8/96
11C Intra-flavor sfermion mixing is neglected
12C for all flavors EXCEPT for stops, sbottoms and staus.
13C
14C
15C 10/23/93
16C It is assumed that the A-terms are real.
17C In addition, all coefficients of the sfermion
18C trilinear terms from the superpotential
19C EXCEPT the stop (AAT), sbottom (AAB) and stau (AAL)
20C coefficients are set to zero.
21C
22C ===> Code for the general case removing all these
23C artificial restrictions is present below.
24C The preceeding restrictions are specified
25C by giving special values to some variables
26C This is discussed in two sections beginning
27C with the symbols (*@&*) in the code below.
28C
29C-----------------------------------------------------------------------
30#if defined(CERNLIB_IMPNONE)
31 IMPLICIT NONE
32#endif
33#include "isajet/sslun.inc"
34#include "isajet/sssm.inc"
35#include "isajet/sspar.inc"
36#include "isajet/sstype.inc"
37C
38C
39 REAL SR2,PI,GG,TW2,BETA,DSA,DCA,DSB,DCB,MH
40 REAL EP,TANB,COTB,ATERM,MSFMIX,THETSF,SIN2B
41 REAL TEMP,TEMP1,TEMP2,YA1,YA2
42 REAL SINA,COSA,SINA2,COSA2,M1,M2,M12,LAMB
43 REAL SINAU,COSAU,SINAD,COSAD
44 REAL A11,A22,A12,B11,B22,B12,C11,C12,C21,C22
45 REAL ASQ,BSQ,CSQ,DWSF
46 REAL DWSFL,DWSFH,DWSFP,DWSFC,SSXLAM
47 REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS
48 DOUBLE PRECISION SSMQCD
49 DIMENSION ATERM(12),MSFMIX(12,2),THETSF(12)
50 DIMENSION ASQ(10,3),BSQ(9),CSQ(6,4)
51 DIMENSION DWSF(12,4),DWSFL(12,4),DWSFH(12,4)
52 DIMENSION DWSFP(12,4),DWSFC(6,4)
53 INTEGER II,IJ,JJ,IC,IJU,IJD,NUMH
54C
55C
56 SR2=SQRT(2.0)
57 PI=4.0*ATAN(1.0)
58 TW2=SN2THW/(1.0-SN2THW)
59 GG=SQRT(4.0*PI*ALFAEM/SN2THW)
60 EP=TWOM1
61C
62 TANB=1.0/RV2V1
63 COTB=RV2V1
64 BETA=ATAN(1.0/RV2V1)
65 DSA=SIN(ALFAH)
66 DCA=COS(ALFAH)
67 DSB=SIN(BETA)
68 DCB=COS(BETA)
69 SIN2B=2.0*DSB*DCB
70C
71C Set A-terms.
72C (all A-terms are assumed to be real)
73C The A-terms are loaded into the array ATERM(12)
74C in the following way:
75C ATERM(1)=selectron A-term
76C ATERM(2)=smuon A-term
77C ATERM(3)=stau A-term
78C ATERM(4)=up squark A-term
79C ATERM(5)=charm squark A-term
80C ATERM(6)=down squark A-term
81C ATERM(7)=strange squark A-term
82C ATERM(8)=sbottom A-term
83C ATERM(9)=stop A-term
84C ATERM(10)=selectronic sneutrino A-term
85C ATERM(11)=smuonic sneutrino A-term
86C ATERM(12)=stauonic sneutrino A-term
87C
88 DO 10 II=1,7
89 ATERM(II)=0.0
9010 CONTINUE
91 ATERM(3)=AAL
92 ATERM(8)=AAB
93 ATERM(9)=AAT
94 DO 20 II=10,12
95 ATERM(II)=0.0
9620 CONTINUE
97C
98C Set mixing parameters.
99C The intra-flavor-mixed sfermion masses are loaded into
100C the array MSFMIX(12,2) where (#,1) is the lighter
101C mixed sfermion mass of a given flavor and (#,2) is the
102C heavier sfermion mass. The sfermionic mixing angles are
103C loaded into the array THETSF(12). The identities of the
104C elements of these arrays are given below:
105C MSFMIX(1,*)=mixed selectron masses
106C THETSF(1)=selectron mixing angle
107C MSFMIX(2,*)=mixed smuon masses
108C THETSF(2)=smuon mixing angle
109C MSFMIX(3,*)=mixed stau masses
110C THETSF(3)=stau mixing angle
111C MSFMIX(4,*)=mixed up squark masses
112C THETSF(4)=up squark mixing angle
113C MSFMIX(5,*)=mixed charm squark masses
114C THETSF(5)=charm squark mixing angle
115C MSFMIX(6,*)=mixed down squark masses
116C THETSF(6)=down squark mixing angle
117C MSFMIX(7,*)=mixed strange squark masses
118C THETSF(7)=strange squark mixing angle
119C MSFMIX(8,*)=mixed sbottom masses
120C THETSF(8)=sbottom mixing angle
121C MSFMIX(9,*)=mixed stop masses
122C THETSF(9)=stop mixing angle
123C For sneuterinos MSFMIX(#,2)=0.0, THETSF(#)=0.0 ; #=10-12
124C Yukawa contributions from D-terms to the sneutrino masses
125C are supposed to be added in here.
126C MSFMIX(10,1)= selectronic sneutrino mass with D-terms
127C MSFMIX(11,1)= smuonic sneutrino mass with D-terms
128C MSFMIX(12,1)= stauonic sneutrino mass with D-terms
129C
130 DO 30 II=10,12
131 MSFMIX(II,2)=0.0
132 THETSF(II)=0.0
13330 CONTINUE
134C
135C
136C (*@&*) 10/24/93 - Special conditions used ---
137C set all mixing angles EXCEPT stop, sbottom, stau to zero.
138C For all EXCEPT st, sb and stau, set mixed sfermion masses
139C to bare sfermion masses:
140C MSFMIX(#,1) = Left sfermion mass
141C MSFMIX(#,2) = Right sfermion mass ; # = 1-8
142C but
143C MSFMIX(9,1) = AMT1SS
144C MSFMIX(9,2) = AMT2SS , etc.
145C
146C (The choice of which to call Left and which to call
147C Right is based on the definition of the sfermion
148C mixing angle theta_sf :
149C sfermion_1 = cos(theta_sf) * sfermion_L
150C - sin(theta_sf) * sfermion_R
151C sfermion_2 = sin(theta_sf) * sfermion_L
152C + cos(theta_sf) * sfermion_R
153C Thus if we set theta_sf = 0, then
154C sfermion_1 = sfermion_L
155C and sfermion_2 = sfermion_R . )
156C
157 DO 40 II=1,7
158 THETSF(II)=0.0
15940 CONTINUE
160 MSFMIX(1,1)=AMELSS
161 MSFMIX(1,2)=AMERSS
162 MSFMIX(2,1)=AMMLSS
163 MSFMIX(2,2)=AMMRSS
164 MSFMIX(3,1)=AML1SS
165 MSFMIX(3,2)=AML2SS
166 THETSF(3)=THETAL
167 MSFMIX(4,1)=AMULSS
168 MSFMIX(4,2)=AMURSS
169 MSFMIX(5,1)=AMCLSS
170 MSFMIX(5,2)=AMCRSS
171 MSFMIX(6,1)=AMDLSS
172 MSFMIX(6,2)=AMDRSS
173 MSFMIX(7,1)=AMSLSS
174 MSFMIX(7,2)=AMSRSS
175 MSFMIX(8,1)=AMB1SS
176 MSFMIX(8,2)=AMB2SS
177 THETSF(8)=THETAB
178 MSFMIX(9,1)=AMT1SS
179 MSFMIX(9,2)=AMT2SS
180 THETSF(9)=THETAT
181 MSFMIX(10,1)=AMN1SS
182 MSFMIX(11,1)=AMN2SS
183 MSFMIX(12,1)=AMN3SS
184C
185 DO 1000 NUMH=1,4
186 IF(NUMH.EQ.1) THEN
187 MH=AMHL
188 ELSE IF(NUMH.EQ.2) THEN
189 MH=AMHH
190 ELSE IF(NUMH.EQ.3) THEN
191 MH=AMHA
192 GO TO 233
193 ELSE IF(NUMH.EQ.4) THEN
194 MH=AMHC
195 GO TO 333
196 ENDIF
197 ASMB=SUALFS(AMBT**2,.36,AMTP,3)
198 MBMB=AMBT*(1.-4*ASMB/3./PI)
199 MBQ=SSMQCD(DBLE(MBMB),DBLE(MH))
200 ASMT=SUALFS(AMTP**2,.36,AMTP,3)
201 MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
202 $(ASMT/PI)**2)
203 MTQ=SSMQCD(DBLE(MTMT),DBLE(MH))
204
205C
206C Scalar neutral Higgses --> sfermions
207C partial decay widths
208C
209 IF(NUMH.EQ.1) THEN
210 TEMP=GG*AMW*SIN(BETA-ALFAH)/2.0
211 YA1=DCA
212 YA2=DSA
213 ELSE IF(NUMH.EQ.2) THEN
214 TEMP=-GG*AMW*COS(BETA-ALFAH)/2.0
215 YA1=-DSA
216 YA2=DCA
217 ENDIF
218C
219 TEMP1=TEMP*(1.0-TW2/3.0)
220 TEMP2=GG*YA1/(AMW*DSB)
221 ASQ(4,1)=TEMP1-TEMP2*AMUP**2
222 ASQ(5,1)=TEMP1-TEMP2*AMCH**2
223 ASQ(9,1)=TEMP1-TEMP2*MTQ**2
224C
225 TEMP1=-TEMP*(1.0+TW2/3.0)
226 TEMP2=GG*YA2/(AMW*DCB)
227 ASQ(6,1)=-TEMP1-TEMP2*AMDN**2
228 ASQ(7,1)=-TEMP1-TEMP2*AMST**2
229 ASQ(8,1)=-TEMP1-TEMP2*MBQ**2
230C
231 ASQ(10,1)=TEMP*(1.0+TW2)
232 TEMP1=TEMP*(TW2-1.0)
233 TEMP2=GG*YA2/(AMW*DCB)
234 ASQ(1,1)=TEMP1-TEMP2*AME**2
235 ASQ(2,1)=TEMP1-TEMP2*AMMU**2
236 ASQ(3,1)=TEMP1-TEMP2*AMTAU**2
237C
238 TEMP1=4.0*TEMP*TW2/3.0
239 TEMP2=GG*YA1/(AMW*DSB)
240 ASQ(4,2)=TEMP1-TEMP2*AMUP**2
241 ASQ(5,2)=TEMP1-TEMP2*AMCH**2
242 ASQ(9,2)=TEMP1-TEMP2*MTQ**2
243C
244 TEMP1=-2.0*TEMP*TW2/3.0
245 TEMP2=GG*YA2/(AMW*DCB)
246 ASQ(6,2)=TEMP1-TEMP2*AMDN**2
247 ASQ(7,2)=TEMP1-TEMP2*AMST**2
248 ASQ(8,2)=TEMP1-TEMP2*MBQ**2
249C
250 ASQ(10,2)=0.0
251 TEMP1=-2.0*TEMP*TW2
252 TEMP2=GG*YA2/(AMW*DCB)
253 ASQ(1,2)=TEMP1-TEMP2*AME**2
254 ASQ(2,2)=TEMP1-TEMP2*AMMU**2
255 ASQ(3,2)=TEMP1-TEMP2*AMTAU**2
256C
257 TEMP1=GG/(2.0*AMW*DSB)
258 ASQ(4,3)=(EP*YA2 + ATERM(4)*YA1)*TEMP1*AMUP
259 ASQ(5,3)=(EP*YA2 + ATERM(5)*YA1)*TEMP1*AMCH
260 ASQ(9,3)=(EP*YA2 + ATERM(9)*YA1)*TEMP1*MTQ
261C
262 TEMP1=GG/(2.0*AMW*DCB)
263 ASQ(6,3)=(ATERM(6)*YA2 + EP*YA1)*TEMP1*AMDN
264 ASQ(7,3)=(ATERM(7)*YA2 + EP*YA1)*TEMP1*AMST
265 ASQ(8,3)=(ATERM(8)*YA2 + EP*YA1)*TEMP1*MBQ
266C
267 ASQ(10,3)=0.0
268 ASQ(1,3)=(ATERM(1)*YA2 + EP*YA1)*TEMP1*AME
269 ASQ(2,3)=(ATERM(2)*YA2 + EP*YA1)*TEMP1*AMMU
270 ASQ(3,3)=(ATERM(3)*YA2 + EP*YA1)*TEMP1*AMTAU
271C
272C
273 DO 150 IJ=1,9
274 IF(IJ.LT.4) THEN
275 TEMP1=1.0/(16.0*PI*MH**3)
276 ELSE
277 TEMP1=3.0/(16.0*PI*MH**3)
278 ENDIF
279 SINA=SIN(THETSF(IJ))
280 COSA=COS(THETSF(IJ))
281 SINA2=SINA**2
282 COSA2=COSA**2
283 M1=MSFMIX(IJ,1)
284 M2=MSFMIX(IJ,1)
285 M12=M1+M2
286 IF(MH.GT.M12) THEN
287 A11=ASQ(IJ,1)*COSA2+ASQ(IJ,2)*SINA2
288 $ -2.0*ASQ(IJ,3)*SINA*COSA
289 LAMB=SSXLAM(MH**2,M1**2,M2**2)
290 DWSF(IJ,1)=TEMP1*SQRT(LAMB)*A11**2
291 ELSE IF(MH.LE.M12) THEN
292 DWSF(IJ,1)=0.0
293 ENDIF
294C
295 M1=MSFMIX(IJ,2)
296 M2=MSFMIX(IJ,2)
297 M12=M1+M2
298 IF(MH.GT.M12) THEN
299 A22=ASQ(IJ,1)*SINA2+ASQ(IJ,2)*COSA2
300 $ +2.0*ASQ(IJ,3)*SINA*COSA
301 LAMB=SSXLAM(MH**2,M1**2,M2**2)
302 DWSF(IJ,2)=TEMP1*SQRT(LAMB)*A22**2
303 ELSE IF(MH.LE.M12) THEN
304 DWSF(IJ,2)=0.0
305 ENDIF
306C
307 M1=MSFMIX(IJ,1)
308 M2=MSFMIX(IJ,2)
309 M12=M1+M2
310 IF(MH.GT.M12) THEN
311 A12=(ASQ(IJ,1)-ASQ(IJ,2))*SINA*COSA
312 $ +ASQ(IJ,3)*(COSA2-SINA2)
313 LAMB=SSXLAM(MH**2,M1**2,M2**2)
314 DWSF(IJ,3)=TEMP1*SQRT(LAMB)*A12**2
315 ELSE IF(MH.LE.M12) THEN
316 DWSF(IJ,3)=0.0
317 ENDIF
318C
319 DWSF(IJ,4)=DWSF(IJ,3)
320C
321 IF(NUMH.EQ.1) THEN
322 DO 121 JJ=1,4
323 DWSFL(IJ,JJ)=DWSF(IJ,JJ)
324121 CONTINUE
325 ELSE IF(NUMH.EQ.2) THEN
326 DO 122 JJ=1,4
327 DWSFH(IJ,JJ)=DWSF(IJ,JJ)
328122 CONTINUE
329 ENDIF
330C
331150 CONTINUE
332C
333C Now take care of sneutrinos.
334C
335 DO 155 IJ=10,12
336 M1=MSFMIX(IJ,1)
337 M2=MSFMIX(IJ,1)
338 M12=M1+M2
339 IF(MH.GT.M12) THEN
340 LAMB=SSXLAM(MH**2,M1**2,M2**2)
341 DWSF(IJ,1)=SQRT(LAMB)*(ASQ(10,1))**2
342 $ /(16.0*PI*MH**3)
343 ELSE IF(MH.LE.M12) THEN
344 DWSF(IJ,1) = 0.0
345 ENDIF
346 DWSF(IJ,2)=0.0
347 DWSF(IJ,3)=0.0
348 DWSF(IJ,4)=0.0
349 IF(NUMH.EQ.1) THEN
350 DO 151 JJ=1,4
351 DWSFL(IJ,JJ)=DWSF(IJ,JJ)
352151 CONTINUE
353 ELSE IF(NUMH.EQ.2) THEN
354 DO 152 JJ=1,4
355 DWSFH(IJ,JJ)=DWSF(IJ,JJ)
356152 CONTINUE
357 ENDIF
358C
359155 CONTINUE
360 GO TO 1000
361C
362C
363C Pseudocalar neutral Higgses --> sfermions
364C partial decay widths
365C
366233 TEMP1=GG/(2.0*AMW)
367 BSQ(1)=TEMP1*AME*(EP-TANB*ATERM(1))
368 BSQ(2)=TEMP1*AMMU*(EP-TANB*ATERM(2))
369 BSQ(3)=TEMP1*AMTAU*(EP-TANB*ATERM(3))
370 BSQ(4)=TEMP1*AMUP*(EP-COTB*ATERM(4))
371 BSQ(5)=TEMP1*AMCH*(EP-COTB*ATERM(5))
372 BSQ(6)=TEMP1*AMDN*(EP-TANB*ATERM(6))
373 BSQ(7)=TEMP1*AMST*(EP-TANB*ATERM(7))
374 BSQ(8)=TEMP1*MBQ*(EP-TANB*ATERM(8))
375 BSQ(9)=TEMP1*MTQ*(EP-COTB*ATERM(9))
376C
377 DO 260 IJ=1,9
378 IF(IJ.LT.4) THEN
379 TEMP1=1.0/(16.0*PI*MH**3)
380 ELSE
381 TEMP1=3.0/(16.0*PI*MH**3)
382 ENDIF
383 SINA=SIN(THETSF(IJ))
384 COSA=COS(THETSF(IJ))
385 SINA2=SINA**2
386 COSA2=COSA**2
387 M1=MSFMIX(IJ,1)
388 M2=MSFMIX(IJ,1)
389 M12=M1+M2
390 IF(MH.GT.M12) THEN
391 B11=-2.0*COSA*SINA*BSQ(IJ)
392 LAMB=SSXLAM(MH**2,M1**2,M2**2)
393 DWSFP(IJ,1)=TEMP1*SQRT(LAMB)*B11**2
394 ELSE IF(MH.LE.M12) THEN
395 DWSFP(IJ,1)=0.0
396 ENDIF
397C
398 M1=MSFMIX(IJ,2)
399 M2=MSFMIX(IJ,2)
400 M12=M1+M2
401 IF(MH.GT.M12) THEN
402 B22=-B11
403 LAMB=SSXLAM(MH**2,M1**2,M2**2)
404 DWSFP(IJ,2)=TEMP1*SQRT(LAMB)*B22**2
405 ELSE IF(MH.LE.M12) THEN
406 DWSFP(IJ,2)=0.0
407 ENDIF
408 M1=MSFMIX(IJ,1)
409 M2=MSFMIX(IJ,2)
410 M12=M1+M2
411 IF(MH.GT.M12) THEN
412 B12=(COSA2-SINA2)*BSQ(IJ)
413 LAMB=SSXLAM(MH**2,M1**2,M2**2)
414 DWSFP(IJ,3)=TEMP1*SQRT(LAMB)*B12**2
415 ELSE IF(MH.LE.M12) THEN
416 DWSFP(IJ,3)=0.0
417 ENDIF
418 DWSFP(IJ,4)=DWSFP(IJ,3)
419260 CONTINUE
420 DO 265 IJ=10,12
421 DO 264 JJ=1,4
422 DWSFP(IJ,JJ)=0.0
423264 CONTINUE
424265 CONTINUE
425 GO TO 1000
426C
427C Charged Higgses --> sfermions
428C partial decay widths
429C
430333 TEMP1=-AMW*SIN2B
431 CSQ(1,1)=GG*(TEMP1+(TANB*AMDN**2 + COTB*AMUP**2)/AMW)/SR2
432 CSQ(2,1)=GG*(TEMP1+(TANB*AMST**2 + COTB*AMCH**2)/AMW)/SR2
433 CSQ(3,1)=GG*(TEMP1+(TANB*MBQ**2 + COTB*MTQ**2)/AMW)/SR2
434 CSQ(4,1)=GG*(TEMP1 + (TANB*AME**2)/AMW)/SR2
435 CSQ(5,1)=GG*(TEMP1 + (TANB*AMMU**2)/AMW)/SR2
436 CSQ(6,1)=GG*(TEMP1 + (TANB*AMTAU**2)/AMW)/SR2
437C
438 TEMP1=GG*(COTB+TANB)/(SR2*AMW)
439 CSQ(1,2)=TEMP1*AMUP*AMDN
440 CSQ(2,2)=TEMP1*AMCH*AMST
441 CSQ(3,2)=TEMP1*MTQ*MBQ
442 CSQ(4,2)=0.0
443 CSQ(5,2)=0.0
444 CSQ(6,2)=0.0
445C
446 TEMP1=GG/(SR2*AMW)
447 CSQ(1,3)=TEMP1*AMUP*(EP-COTB*ATERM(4))
448 CSQ(2,3)=TEMP1*AMCH*(EP-COTB*ATERM(5))
449 CSQ(3,3)=TEMP1*MTQ*(EP-COTB*ATERM(9))
450 CSQ(4,3)=0.0
451 CSQ(5,3)=0.0
452 CSQ(6,3)=0.0
453C
454 CSQ(1,4)=TEMP1* AMDN*(EP-TANB*ATERM(6))
455 CSQ(2,4)=TEMP1* AMST*(EP-TANB*ATERM(7))
456 CSQ(3,4)=TEMP1* MBQ*(EP-TANB*ATERM(8))
457 CSQ(4,4)=TEMP1* AME*(EP-TANB*ATERM(1))
458 CSQ(5,4)=TEMP1* AMMU*(EP-TANB*ATERM(2))
459 CSQ(6,4)=TEMP1* AMTAU*(EP-TANB*ATERM(3))
460C
461 DO 350 IC=1,3
462 TEMP1=3.0/(16.0*PI*MH**3)
463 IF(IC.EQ.1) THEN
464 IJU=4
465 IJD=6
466 ELSE IF(IC.EQ.2) THEN
467 IJU=5
468 IJD=7
469 ELSE IF(IC.EQ.3) THEN
470 IJU=9
471 IJD=8
472 ENDIF
473 SINAU=SIN(THETSF(IJU))
474 COSAU=COS(THETSF(IJU))
475 SINAD=SIN(THETSF(IJD))
476 COSAD=COS(THETSF(IJD))
477C
478 M1=MSFMIX(IJU,1)
479 M2=MSFMIX(IJD,1)
480 M12=M1+M2
481 IF(MH.GT.M12) THEN
482 C11=COSAU*COSAD*CSQ(IC,1)
483 $ + SINAU*SINAD*CSQ(IC,2)
484 $ - SINAU*COSAD*CSQ(IC,3)
485 $ - COSAU*SINAD*CSQ(IC,4)
486 LAMB=SSXLAM(MH**2,M1**2,M2**2)
487 DWSFC(IC,1)=TEMP1*SQRT(LAMB)*C11**2
488 ELSE IF(MH.LE.M12) THEN
489 DWSFC(IC,1) = 0.0
490 ENDIF
491C
492 M1=MSFMIX(IJU,1)
493 M2=MSFMIX(IJD,2)
494 M12=M1+M2
495 IF(MH.GT.M12) THEN
496 C12=COSAU*SINAD*CSQ(IC,1)
497 $ - SINAU*COSAD*CSQ(IC,2)
498 $ - SINAU*SINAD*CSQ(IC,3)
499 $ + COSAU*COSAD*CSQ(IC,4)
500 LAMB=SSXLAM(MH**2,M1**2,M2**2)
501 DWSFC(IC,2)=TEMP1*SQRT(LAMB)*C12**2
502 ELSE IF(MH.LE.M12) THEN
503 DWSFC(IC,2)=0.0
504 ENDIF
505C
506 M1=MSFMIX(IJU,2)
507 M2=MSFMIX(IJD,1)
508 M12=M1+M2
509 IF(MH.GT.M12) THEN
510 C21=SINAU*COSAD*CSQ(IC,1)
511 $ - COSAU*SINAD*CSQ(IC,2)
512 $ + COSAU*COSAD*CSQ(IC,3)
513 $ - SINAU*SINAD*CSQ(IC,4)
514 LAMB=SSXLAM(MH**2,M1**2,M2**2)
515 DWSFC(IC,3)=TEMP1*SQRT(LAMB)*C21**2
516 ELSE IF(MH.LE.M12) THEN
517 DWSFC(IC,3)=0.0
518 ENDIF
519C
520 M1=MSFMIX(IJU,2)
521 M2=MSFMIX(IJD,2)
522 M12=M1+M2
523 IF(MH.GT.M12) THEN
524 C22=SINAU*SINAD*CSQ(IC,1)
525 $ + COSAU*COSAD*CSQ(IC,2)
526 $ + COSAU*SINAD*CSQ(IC,3)
527 $ - SINAU*COSAD*CSQ(IC,4)
528 LAMB=SSXLAM(MH**2,M1**2,M2**2)
529 DWSFC(IC,4)=TEMP1*SQRT(LAMB)*C22**2
530 ELSE IF(MH.LE.M12) THEN
531 DWSFC(IC,4)=0.0
532 ENDIF
533C
534350 CONTINUE
535C
536C
537C Now calculate the sleptonic
538C partial decay widths of the
539C charged Higgs.
540C
541 DO 355 IC = 4,6
542 TEMP1=1.0/(16.0*PI*MH**3)
543 IF(IC.EQ.4) THEN
544 IJU=10
545 IJD=1
546 ELSE IF(IC.EQ.5) THEN
547 IJU=11
548 IJD=2
549 ELSE IF(IC.EQ.6) THEN
550 IJU=12
551 IJD=3
552 ENDIF
553 SINAD=SIN(THETSF(IJD))
554 COSAD=COS(THETSF(IJD))
555C
556 M1=MSFMIX(IJU,1)
557 M2=MSFMIX(IJD,1)
558 M12=M1+M2
559 IF(MH.GT.M12) THEN
560 C11=COSAD*CSQ(IC,1)-SINAD*CSQ(IC,4)
561 LAMB=SSXLAM(MH**2,M1**2,M2**2)
562 DWSFC(IC,1)=TEMP1*SQRT(LAMB)*C11**2
563 ELSE IF(MH.LE.M12) THEN
564 DWSFC(IC,1)=0.0
565 ENDIF
566C
567 M1=MSFMIX(IJU,1)
568 M2=MSFMIX(IJD,2)
569 M12=M1+M2
570 IF(MH.GT.M12) THEN
571 C12=SINAD*CSQ(IC,1)+COSAD*CSQ(IC,4)
572 LAMB=SSXLAM(MH**2,M1**2,M2**2)
573 DWSFC(IC,2)=TEMP1*SQRT(LAMB)*C12**2
574 ELSE IF(MH.LE.M12) THEN
575 DWSFC(IC,2)=0.0
576 ENDIF
577 DWSFC(IC,3)=0.0
578 DWSFC(IC,4)=0.0
579355 CONTINUE
5801000 CONTINUE
581C H_l decays
582 CALL SSSAVE(ISHL,DWSFL(1,1),ISEL,-ISEL,0,0,0)
583 CALL SSSAVE(ISHL,DWSFL(1,2),ISER,-ISER,0,0,0)
584 CALL SSSAVE(ISHL,DWSFL(2,1),ISMUL,-ISMUL,0,0,0)
585 CALL SSSAVE(ISHL,DWSFL(2,2),ISMUR,-ISMUR,0,0,0)
586 CALL SSSAVE(ISHL,DWSFL(3,1),ISTAU1,-ISTAU1,0,0,0)
587 CALL SSSAVE(ISHL,DWSFL(3,2),ISTAU2,-ISTAU2,0,0,0)
588 CALL SSSAVE(ISHL,DWSFL(3,3),ISTAU1,-ISTAU2,0,0,0)
589 CALL SSSAVE(ISHL,DWSFL(3,4),ISTAU2,-ISTAU1,0,0,0)
590 CALL SSSAVE(ISHL,DWSFL(4,1),ISUPL,-ISUPL,0,0,0)
591 CALL SSSAVE(ISHL,DWSFL(4,2),ISUPR,-ISUPR,0,0,0)
592 CALL SSSAVE(ISHL,DWSFL(5,1),ISCHL,-ISCHL,0,0,0)
593 CALL SSSAVE(ISHL,DWSFL(5,2),ISCHR,-ISCHR,0,0,0)
594 CALL SSSAVE(ISHL,DWSFL(6,1),ISDNL,-ISDNL,0,0,0)
595 CALL SSSAVE(ISHL,DWSFL(6,2),ISDNR,-ISDNR,0,0,0)
596 CALL SSSAVE(ISHL,DWSFL(7,1),ISSTL,-ISSTL,0,0,0)
597 CALL SSSAVE(ISHL,DWSFL(7,2),ISSTR,-ISSTR,0,0,0)
598 CALL SSSAVE(ISHL,DWSFL(8,1),ISBT1,-ISBT1,0,0,0)
599 CALL SSSAVE(ISHL,DWSFL(8,2),ISBT2,-ISBT2,0,0,0)
600 CALL SSSAVE(ISHL,DWSFL(8,3),ISBT1,-ISBT2,0,0,0)
601 CALL SSSAVE(ISHL,DWSFL(8,4),ISBT2,-ISBT1,0,0,0)
602 CALL SSSAVE(ISHL,DWSFL(9,1),ISTP1,-ISTP1,0,0,0)
603 CALL SSSAVE(ISHL,DWSFL(9,2),ISTP2,-ISTP2,0,0,0)
604 CALL SSSAVE(ISHL,DWSFL(9,3),ISTP1,-ISTP2,0,0,0)
605 CALL SSSAVE(ISHL,DWSFL(9,4),ISTP2,-ISTP1,0,0,0)
606 CALL SSSAVE(ISHL,DWSFL(10,1),ISNEL,-ISNEL,0,0,0)
607 CALL SSSAVE(ISHL,DWSFL(11,1),ISNML,-ISNML,0,0,0)
608 CALL SSSAVE(ISHL,DWSFL(12,1),ISNTL,-ISNTL,0,0,0)
609C H_h decays
610 CALL SSSAVE(ISHH,DWSFH(1,1),ISEL,-ISEL,0,0,0)
611 CALL SSSAVE(ISHH,DWSFH(1,2),ISER,-ISER,0,0,0)
612 CALL SSSAVE(ISHH,DWSFH(2,1),ISMUL,-ISMUL,0,0,0)
613 CALL SSSAVE(ISHH,DWSFH(2,2),ISMUR,-ISMUR,0,0,0)
614 CALL SSSAVE(ISHH,DWSFH(3,1),ISTAU1,-ISTAU1,0,0,0)
615 CALL SSSAVE(ISHH,DWSFH(3,2),ISTAU2,-ISTAU2,0,0,0)
616 CALL SSSAVE(ISHH,DWSFH(3,3),ISTAU1,-ISTAU2,0,0,0)
617 CALL SSSAVE(ISHH,DWSFH(3,4),ISTAU2,-ISTAU1,0,0,0)
618 CALL SSSAVE(ISHH,DWSFH(4,1),ISUPL,-ISUPL,0,0,0)
619 CALL SSSAVE(ISHH,DWSFH(4,2),ISUPR,-ISUPR,0,0,0)
620 CALL SSSAVE(ISHH,DWSFH(5,1),ISCHL,-ISCHL,0,0,0)
621 CALL SSSAVE(ISHH,DWSFH(5,2),ISCHR,-ISCHR,0,0,0)
622 CALL SSSAVE(ISHH,DWSFH(6,1),ISDNL,-ISDNL,0,0,0)
623 CALL SSSAVE(ISHH,DWSFH(6,2),ISDNR,-ISDNR,0,0,0)
624 CALL SSSAVE(ISHH,DWSFH(7,1),ISSTL,-ISSTL,0,0,0)
625 CALL SSSAVE(ISHH,DWSFH(7,2),ISSTR,-ISSTR,0,0,0)
626 CALL SSSAVE(ISHH,DWSFH(8,1),ISBT1,-ISBT1,0,0,0)
627 CALL SSSAVE(ISHH,DWSFH(8,2),ISBT2,-ISBT2,0,0,0)
628 CALL SSSAVE(ISHH,DWSFH(8,3),ISBT1,-ISBT2,0,0,0)
629 CALL SSSAVE(ISHH,DWSFH(8,4),ISBT2,-ISBT1,0,0,0)
630 CALL SSSAVE(ISHH,DWSFH(9,1),ISTP1,-ISTP1,0,0,0)
631 CALL SSSAVE(ISHH,DWSFH(9,2),ISTP2,-ISTP2,0,0,0)
632 CALL SSSAVE(ISHH,DWSFH(9,3),ISTP1,-ISTP2,0,0,0)
633 CALL SSSAVE(ISHH,DWSFH(9,4),ISTP2,-ISTP1,0,0,0)
634 CALL SSSAVE(ISHH,DWSFH(10,1),ISNEL,-ISNEL,0,0,0)
635 CALL SSSAVE(ISHH,DWSFH(11,1),ISNML,-ISNML,0,0,0)
636 CALL SSSAVE(ISHH,DWSFH(12,1),ISNTL,-ISNTL,0,0,0)
637C Decay of H_p
638 CALL SSSAVE(ISHA,DWSFP(1,3),ISEL,-ISER,0,0,0)
639 CALL SSSAVE(ISHA,DWSFP(1,4),ISER,-ISEL,0,0,0)
640 CALL SSSAVE(ISHA,DWSFP(2,3),ISMUL,-ISMUR,0,0,0)
641 CALL SSSAVE(ISHA,DWSFP(2,4),ISMUR,-ISMUL,0,0,0)
642 CALL SSSAVE(ISHA,DWSFP(3,1),ISTAU1,-ISTAU1,0,0,0)
643 CALL SSSAVE(ISHA,DWSFP(3,2),ISTAU2,-ISTAU2,0,0,0)
644 CALL SSSAVE(ISHA,DWSFP(3,3),ISTAU1,-ISTAU2,0,0,0)
645 CALL SSSAVE(ISHA,DWSFP(3,4),ISTAU2,-ISTAU1,0,0,0)
646 CALL SSSAVE(ISHA,DWSFP(4,3),ISUPL,-ISUPR,0,0,0)
647 CALL SSSAVE(ISHA,DWSFP(4,4),ISUPR,-ISUPL,0,0,0)
648 CALL SSSAVE(ISHA,DWSFP(5,3),ISCHL,-ISCHR,0,0,0)
649 CALL SSSAVE(ISHA,DWSFP(5,4),ISCHR,-ISCHL,0,0,0)
650 CALL SSSAVE(ISHA,DWSFP(6,3),ISDNL,-ISDNR,0,0,0)
651 CALL SSSAVE(ISHA,DWSFP(6,4),ISDNR,-ISDNL,0,0,0)
652 CALL SSSAVE(ISHA,DWSFP(7,3),ISSTL,-ISSTR,0,0,0)
653 CALL SSSAVE(ISHA,DWSFP(7,4),ISSTR,-ISSTL,0,0,0)
654 CALL SSSAVE(ISHA,DWSFP(8,1),ISBT1,-ISBT1,0,0,0)
655 CALL SSSAVE(ISHA,DWSFP(8,2),ISBT2,-ISBT2,0,0,0)
656 CALL SSSAVE(ISHA,DWSFP(8,3),ISBT1,-ISBT2,0,0,0)
657 CALL SSSAVE(ISHA,DWSFP(8,4),ISBT2,-ISBT1,0,0,0)
658 CALL SSSAVE(ISHA,DWSFP(9,1),ISTP1,-ISTP1,0,0,0)
659 CALL SSSAVE(ISHA,DWSFP(9,2),ISTP2,-ISTP2,0,0,0)
660 CALL SSSAVE(ISHA,DWSFP(9,3),ISTP1,-ISTP2,0,0,0)
661 CALL SSSAVE(ISHA,DWSFP(9,4),ISTP2,-ISTP1,0,0,0)
662C Decay of H+
663 CALL SSSAVE(ISHC,DWSFC(1,1),ISUPL,-ISDNL,0,0,0)
664 CALL SSSAVE(ISHC,DWSFC(1,2),ISUPR,-ISDNR,0,0,0)
665 CALL SSSAVE(ISHC,DWSFC(2,1),ISCHL,-ISSTL,0,0,0)
666 CALL SSSAVE(ISHC,DWSFC(2,2),ISCHR,-ISSTR,0,0,0)
667 CALL SSSAVE(ISHC,DWSFC(3,1),ISTP1,-ISBT1,0,0,0)
668 CALL SSSAVE(ISHC,DWSFC(3,2),ISTP1,-ISBT2,0,0,0)
669 CALL SSSAVE(ISHC,DWSFC(3,3),ISTP2,-ISBT1,0,0,0)
670 CALL SSSAVE(ISHC,DWSFC(3,4),ISTP2,-ISBT2,0,0,0)
671 CALL SSSAVE(ISHC,DWSFC(4,1),-ISEL,ISNEL,0,0,0)
672 CALL SSSAVE(ISHC,DWSFC(5,1),-ISMUL,ISNML,0,0,0)
673 CALL SSSAVE(ISHC,DWSFC(6,1),-ISTAU1,ISNTL,0,0,0)
674 CALL SSSAVE(ISHC,DWSFC(6,2),-ISTAU2,ISNTL,0,0,0)
675 RETURN
676 END