1 #include "isajet/pilot.h"
3 C-----------------------------------------------------------------------
5 C Calculate H -> gl gl decays including both SM particles and
6 C SUSY particles in loop.
8 C This subroutine uses the tau variable of the Higgs Hunters'
9 C Guide. Many other authors, including the paper cited in
10 C Higgs Hunters' Guide (PR. D. 38(11): 3481) and Collider Physics
11 C by Barger and Phillips use the variable lambda
12 C LAMBDA = ( MASS OF PARTICLE IN LOOP / MASS OF HIGGS )**2
16 C-----------------------------------------------------------------------
17 #if defined(CERNLIB_IMPNONE)
20 #include "isajet/sssm.inc"
21 #include "isajet/sspar.inc"
22 #include "isajet/sstype.inc"
24 DOUBLE PRECISION ETAH,IITOT,RITOT,TAU,IFFF,RFFF
25 $,IFHALF,RFHALF,IF1,RF1,IF0,RF0,TW2,RHF,RHSF,RHSFL,RHSFR
26 $,IIHF,RIHF,IIHSFL,RIHSFL,IIHSFR,RIHSFR,AS,SUMISQ,DW
27 $,RHSF1,RHSF2,IIHSF1,IIHSF2,RIHSF1,RIHSF2
28 DOUBLE PRECISION PI,SR2,XM,THETX,YM,THETY,SGL,CGL,SGR,CGR
29 $,MW1,MW2,THETM,THETP,G2,BETA,ALPHA,SW2,CW2,MH,AMSQ
30 DOUBLE PRECISION MFL(3),MFD(3),MFU(3)
31 DOUBLE PRECISION SSALFS
33 REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS
34 DOUBLE PRECISION SSMQCD
35 INTEGER IJ,II,NUMOUT,NUMH,IDHHA
37 C Mass matrix parameters
45 SGL=1/(DSQRT(1+XM**2))
47 SGR=1/(DSQRT(1+YM**2))
53 G2=4.0*PI*ALFAEM/SN2THW
59 C Loop over neutral Higgs bosons
65 ELSEIF(NUMH.EQ.2) THEN
76 ASMB=SUALFS(AMBT**2,.36,AMTP,3)
77 MBMB=AMBT*(1.-4*ASMB/3./PI)
78 MBQ=SSMQCD(DBLE(MBMB),DBLE(MH))
79 ASMT=SUALFS(AMTP**2,.36,AMTP,3)
80 MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
82 MTQ=SSMQCD(DBLE(MTMT),DBLE(MH))
96 C Down-type quark loops
99 TAU=4.0*MFD(II)**2/MH**2
100 CALL SSHGM1(TAU,IFFF,RFFF)
101 IFHALF=0.5*TAU*(1.0-TAU*ETAH)*IFFF
102 RFHALF=0.5*TAU*(ETAH+(1.0-TAU*ETAH)*RFFF)
104 RHF=SIN(ALPHA)/COS(BETA)
105 ELSEIF(NUMH.EQ.2) THEN
106 RHF=COS(ALPHA)/COS(BETA)
116 C Up-type quark loops
119 TAU=4.0*MFU(II)**2/MH**2
120 CALL SSHGM1(TAU,IFFF,RFFF)
121 IFHALF=0.5*TAU*(1.0-TAU*ETAH)*IFFF
122 RFHALF=0.5*TAU*(ETAH+(1.0-TAU*ETAH)*RFFF)
124 RHF=COS(ALPHA)/SIN(BETA)
125 ELSEIF(NUMH.EQ.2) THEN
126 RHF=-SIN(ALPHA)/SIN(BETA)
137 CALL SSHGM1(TAU,IFFF,RFFF)
138 IFHALF=0.5*TAU*(1.0-TAU*ETAH)*IFFF
139 RFHALF=0.5*TAU*(ETAH+(1.0-TAU*ETAH)*RFFF)
141 RHF=COS(ALPHA)/SIN(BETA)
142 ELSEIF(NUMH.EQ.2) THEN
143 RHF=-SIN(ALPHA)/SIN(BETA)
152 C Down-type squark loops
153 C Mixing between the sbottom squarks is included, so
154 C masses used here are the mixed masses (AMB1SS & AMB2SS)
155 C First do d_L and s_L squarks
159 RHSF=2.0*(MFD(II)/AMW)**2*SIN(ALPHA)/COS(BETA)
160 RHSFL=(-1.0-TW2/3.0)*SIN(BETA-ALPHA)-RHSF
161 ELSEIF(NUMH.EQ.2) THEN
162 RHSF=2.0*(MFD(II)/AMW)**2*COS(ALPHA)/COS(BETA)
163 RHSFL=(-1.0-TW2/3.0)*(-1.0)*COS(BETA-ALPHA)-RHSF
168 IF (II.EQ.1) AMSQ=AMDLSS
169 IF (II.EQ.2) AMSQ=AMSLSS
170 TAU=4.0*AMSQ**2/MH**2
171 CALL SSHGM1(TAU,IFFF,RFFF)
173 RF0=TAU*(1.0-TAU*RFFF)
174 IIHSFL=RHSFL*IF0*(AMW/AMSQ)**2/8.0
175 RIHSFL=RHSFL*RF0*(AMW/AMSQ)**2/8.0
182 RHSF=2.0*(MFD(II)/AMW)**2*SIN(ALPHA)/COS(BETA)
183 RHSFR=(-2.0*TW2/3.0)*SIN(BETA-ALPHA)-RHSF
184 ELSEIF(NUMH.EQ.2) THEN
185 RHSF=2.0*(MFD(II)/AMW)**2*COS(ALPHA)/COS(BETA)
186 RHSFR=(-2.0*TW2/3.0)*(-1.0)*COS(BETA-ALPHA)-RHSF
191 IF (II.EQ.1) AMSQ=AMDRSS
192 IF (II.EQ.2) AMSQ=AMSRSS
193 TAU=4.0*AMSQ**2/MH**2
194 CALL SSHGM1(TAU,IFFF,RFFF)
196 RF0=TAU*(1.0-TAU*RFFF)
197 IIHSFR=RHSFR*IF0*(AMW/AMSQ)**2/8.0
198 RIHSFR=RHSFR*RF0*(AMW/AMSQ)**2/8.0
203 RHSF=2.0*(MBQ/AMW)**2*SIN(ALPHA)/COS(BETA)
204 RHSFL=(-1.0-TW2/3.0)*SIN(BETA-ALPHA)-RHSF
205 RHSFR=(-2.0*TW2/3.0)*SIN(BETA-ALPHA)-RHSF
206 ELSEIF(NUMH.EQ.2) THEN
207 RHSF=2.0*(MBQ/AMW)**2*COS(ALPHA)/COS(BETA)
208 RHSFL=(-1.0-TW2/3.0)*(-1.0)*COS(BETA-ALPHA)-RHSF
209 RHSFR=(-2.0*TW2/3.0)*(-1.0)*COS(BETA-ALPHA)-RHSF
215 RHSF1=RHSFL*COS(THETAB)-RHSFR*SIN(THETAB)
216 RHSF2=RHSFL*SIN(THETAB)+RHSFR*COS(THETAB)
217 TAU=4.0*AMB1SS**2/MH**2
218 CALL SSHGM1(TAU,IFFF,RFFF)
220 RF0=TAU*(1.0-TAU*RFFF)
221 IIHSF1=RHSF1*IF0*(AMW/AMB1SS)**2/8.0
222 RIHSF1=RHSF1*RF0*(AMW/AMB1SS)**2/8.0
225 TAU=4.0*AMB2SS**2/MH**2
226 CALL SSHGM1(TAU,IFFF,RFFF)
228 RF0=TAU*(1.0-TAU*RFFF)
229 IIHSF2=RHSF2*IF0*(AMW/AMB2SS)**2/8.0
230 RIHSF2=RHSF2*RF0*(AMW/AMB2SS)**2/8.0
234 C Up-type squark loops
235 C Mixing between the stop squarks is included, so
236 C masses used here are the mixed masses (AMT1SS & AMT2SS)
237 C First do u_L and c_L
240 RHSF=2.0*(MFU(II)/AMW)**2*COS(ALPHA)/SIN(BETA)
241 RHSFL=(1.0-TW2/3.0)*SIN(BETA-ALPHA)-RHSF
242 ELSEIF(NUMH.EQ.2) THEN
243 RHSF=2.0*(MFU(II)/AMW)**2
244 RHSF=RHSF*(-1.0)*SIN(ALPHA)/SIN(BETA)
245 RHSFL=(1.0-TW2/3.0)*(-1.0)*COS(BETA-ALPHA)-RHSF
250 IF (II.EQ.1) AMSQ=AMULSS
251 IF (II.EQ.2) AMSQ=AMCLSS
252 TAU=4.0*(AMSQ)**2/MH**2
253 CALL SSHGM1(TAU,IFFF,RFFF)
255 RF0=TAU*(1.0-TAU*RFFF)
256 IIHSFL=RHSFL*IF0*(AMW/AMSQ)**2/8.0
257 RIHSFL=RHSFL*RF0*(AMW/AMSQ)**2/8.0
261 C Next, do u_R and c_R
264 RHSF=2.0*(MFU(II)/AMW)**2*COS(ALPHA)/SIN(BETA)
265 RHSFR=(4.0*TW2/3.0)*SIN(BETA-ALPHA)-RHSF
266 ELSEIF(NUMH.EQ.2) THEN
267 RHSF=2.0*(MFU(II)/AMW)**2
268 RHSF=RHSF*(-1.0)*SIN(ALPHA)/SIN(BETA)
269 RHSFR=(4.0*TW2/3.0)*(-1.0)*COS(BETA-ALPHA)-RHSF
274 IF (II.EQ.1) AMSQ=AMURSS
275 IF (II.EQ.2) AMSQ=AMCRSS
276 TAU=4.0*(AMSQ)**2/MH**2
277 CALL SSHGM1(TAU,IFFF,RFFF)
279 RF0=TAU*(1.0-TAU*RFFF)
280 IIHSFR=RHSFR*IF0*(AMW/AMSQ)**2/8.0
281 RIHSFR=RHSFR*RF0*(AMW/AMSQ)**2/8.0
287 RHSF=2.0*(MTQ/AMW)**2*COS(ALPHA)/SIN(BETA)
288 RHSFL=(1.0-TW2/3.0)*SIN(BETA-ALPHA)-RHSF
289 RHSFR=(4.0*TW2/3.0)*SIN(BETA-ALPHA)-RHSF
290 ELSEIF(NUMH.EQ.2) THEN
291 RHSF=2.0*(MTQ/AMW)**2
292 RHSF=RHSF*(-1.0)*SIN(ALPHA)/SIN(BETA)
293 RHSFL=(1.0-TW2/3.0)*(-1.0)*COS(BETA-ALPHA)-RHSF
294 RHSFR=(4.0*TW2/3.0)*(-1.0)*COS(BETA-ALPHA)-RHSF
300 RHSF1=RHSFL*COS(THETAT)-RHSFR*SIN(THETAT)
301 RHSF2=RHSFL*SIN(THETAT)+RHSFR*COS(THETAT)
302 TAU=4.0*AMT1SS**2/MH**2
303 CALL SSHGM1(TAU,IFFF,RFFF)
305 RF0=TAU*(1.0-TAU*RFFF)
306 IIHSF1=RHSF1*IF0*(AMW/AMT1SS)**2/8.0
307 RIHSF1=RHSF1*RF0*(AMW/AMT1SS)**2/8.0
310 TAU=4.0*AMT2SS**2/MH**2
311 CALL SSHGM1(TAU,IFFF,RFFF)
313 RF0=TAU*(1.0-TAU*RFFF)
314 IIHSF2=RHSF2*IF0*(AMW/AMT2SS)**2/8.0
315 RIHSF2=RHSF2*RF0*(AMW/AMT2SS)**2/8.0
319 C IITOT and RITOT now contain the total imaginary and
320 C real parts of the I function
322 SUMISQ=IITOT**2+RITOT**2
324 DW=AS**2*G2*MH**3/(32.0*(PI**3)*AMW**2)
326 CALL SSSAVE(IDHHA,WID,IDGL,IDGL,0,0,0)