]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ISAJET/isasusy/sshgl.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 / sshgl.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE SSHGL
3C-----------------------------------------------------------------------
4C
5C Calculate H -> gl gl decays including both SM particles and
6C SUSY particles in loop.
7C
8C This subroutine uses the tau variable of the Higgs Hunters'
9C Guide. Many other authors, including the paper cited in
10C Higgs Hunters' Guide (PR. D. 38(11): 3481) and Collider Physics
11C by Barger and Phillips use the variable lambda
12C LAMBDA = ( MASS OF PARTICLE IN LOOP / MASS OF HIGGS )**2
13C TAU = 4.0 * LAMBDA
14C
15C Bisset's HGLGL
16C-----------------------------------------------------------------------
17#if defined(CERNLIB_IMPNONE)
18 IMPLICIT NONE
19#endif
20#include "isajet/sssm.inc"
21#include "isajet/sspar.inc"
22#include "isajet/sstype.inc"
23C
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
32 REAL WID
33 REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS
34 DOUBLE PRECISION SSMQCD
35 INTEGER IJ,II,NUMOUT,NUMH,IDHHA
36C
37C Mass matrix parameters
38C
39 PI=4.*ATAN(1.D0)
40 SR2=SQRT(2.D0)
41 XM=1./TAN(GAMMAL)
42 THETX=SIGN(1.D0,XM)
43 YM=1./TAN(GAMMAR)
44 THETY=SIGN(1.D0,YM)
45 SGL=1/(DSQRT(1+XM**2))
46 CGL=SGL*XM
47 SGR=1/(DSQRT(1+YM**2))
48 CGR=SGR*YM
49 MW1=DBLE(ABS(AMW1SS))
50 MW2=DBLE(ABS(AMW2SS))
51 THETM=SIGN(1.,AMW1SS)
52 THETP=SIGN(1.,AMW2SS)
53 G2=4.0*PI*ALFAEM/SN2THW
54 BETA=ATAN(1.0/RV2V1)
55 ALPHA=ALFAH
56 SW2=SN2THW
57 CW2=1.-SN2THW
58C
59C Loop over neutral Higgs bosons
60C
61 DO 100 NUMH=1,3
62 IF(NUMH.EQ.1) THEN
63 MH=AMHL
64 IDHHA=ISHL
65 ELSEIF(NUMH.EQ.2) THEN
66 MH=AMHH
67 IDHHA=ISHH
68 ELSE
69 MH=AMHA
70 IDHHA=ISHA
71 ENDIF
72 ETAH=1.0
73 IITOT=0.0
74 RITOT=0.0
75C
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))*
81 $(ASMT/PI)**2)
82 MTQ=SSMQCD(DBLE(MTMT),DBLE(MH))
83
84C
85 MFL(1)=DBLE(AME)
86 MFL(2)=DBLE(AMMU)
87 MFL(3)=DBLE(AMTAU)
88 MFD(1)=DBLE(AMDN)
89 MFD(2)=DBLE(AMST)
90 MFD(3)=DBLE(MBQ)
91 MFU(1)=DBLE(AMUP)
92 MFU(2)=DBLE(AMCH)
93 MFU(3)=DBLE(MTQ)
94C
95C
96C Down-type quark loops
97C
98 DO 20 II=1,3
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)
103 IF(NUMH.EQ.1) THEN
104 RHF=SIN(ALPHA)/COS(BETA)
105 ELSEIF(NUMH.EQ.2) THEN
106 RHF=COS(ALPHA)/COS(BETA)
107 ELSE
108 RHF=TAN(BETA)
109 ENDIF
110 IIHF=RHF*IFHALF
111 RIHF=RHF*RFHALF
112 IITOT=IITOT+IIHF
113 RITOT=RITOT+RIHF
11420 CONTINUE
115C
116C Up-type quark loops
117C
118 DO 30 II=1,2
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)
123 IF(NUMH.EQ.1) THEN
124 RHF=COS(ALPHA)/SIN(BETA)
125 ELSEIF(NUMH.EQ.2) THEN
126 RHF=-SIN(ALPHA)/SIN(BETA)
127 ELSE
128 RHF=TAN(BETA)
129 ENDIF
130 IIHF=RHF*IFHALF
131 RIHF=RHF*RFHALF
132 IITOT=IITOT+IIHF
133 RITOT=RITOT+RIHF
13430 CONTINUE
135C
136 TAU=4.0*MTQ**2/MH**2
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)
140 IF(NUMH.EQ.1) THEN
141 RHF=COS(ALPHA)/SIN(BETA)
142 ELSEIF(NUMH.EQ.2) THEN
143 RHF=-SIN(ALPHA)/SIN(BETA)
144 ELSE
145 RHF=1.0/TAN(BETA)
146 ENDIF
147 IIHF=RHF*IFHALF
148 RIHF=RHF*RFHALF
149 IITOT=IITOT+IIHF
150 RITOT=RITOT+RIHF
151C
152C Down-type squark loops
153C Mixing between the sbottom squarks is included, so
154C masses used here are the mixed masses (AMB1SS & AMB2SS)
155C First do d_L and s_L squarks
156 TW2=SW2/CW2
157 DO 50 II=1,2
158 IF(NUMH.EQ.1) THEN
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
164 ELSE
165 RHSF=0
166 RHSFL=0
167 ENDIF
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)
172 IF0=-TAU*TAU*IFFF
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
176 IITOT=IITOT+IIHSFL
177 RITOT=RITOT+RIHSFL
17850 CONTINUE
179c Next, do R squarks
180 DO 51 II=1,2
181 IF(NUMH.EQ.1) THEN
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
187 ELSE
188 RHSF=0
189 RHSFR=0
190 ENDIF
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)
195 IF0=-TAU*TAU*IFFF
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
199 IITOT=IITOT+IIHSFR
200 RITOT=RITOT+RIHSFR
20151 CONTINUE
202 IF(NUMH.EQ.1) THEN
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
210 ELSE
211 RHSF=0
212 RHSFL=0
213 RHSFR=0
214 ENDIF
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)
219 IF0=-TAU*TAU*IFFF
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
223 IITOT=IITOT+IIHSF1
224 RITOT=RITOT+RIHSF1
225 TAU=4.0*AMB2SS**2/MH**2
226 CALL SSHGM1(TAU,IFFF,RFFF)
227 IF0=-TAU*TAU*IFFF
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
231 IITOT=IITOT+IIHSF2
232 RITOT=RITOT+RIHSF2
233C
234C Up-type squark loops
235C Mixing between the stop squarks is included, so
236C masses used here are the mixed masses (AMT1SS & AMT2SS)
237C First do u_L and c_L
238 DO 60 II=1,2
239 IF(NUMH.EQ.1) THEN
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
246 ELSE
247 RHSF=0
248 RHSFL=0
249 ENDIF
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)
254 IF0=-TAU*TAU*IFFF
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
258 IITOT=IITOT+IIHSFL
259 RITOT=RITOT+RIHSFL
26060 CONTINUE
261C Next, do u_R and c_R
262 DO 61 II=1,2
263 IF(NUMH.EQ.1) THEN
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
270 ELSE
271 RHSF=0
272 RHSFR=0
273 ENDIF
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)
278 IF0=-TAU*TAU*IFFF
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
282 IITOT=IITOT+IIHSFR
283 RITOT=RITOT+RIHSFR
28461 CONTINUE
285C
286 IF(NUMH.EQ.1) THEN
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
295 ELSE
296 RHSF=0
297 RHSFL=0
298 RHSFR=0
299 ENDIF
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)
304 IF0=-TAU*TAU*IFFF
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
308 IITOT=IITOT+IIHSF1
309 RITOT=RITOT+RIHSF1
310 TAU=4.0*AMT2SS**2/MH**2
311 CALL SSHGM1(TAU,IFFF,RFFF)
312 IF0=-TAU*TAU*IFFF
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
316 IITOT=IITOT+IIHSF2
317 RITOT=RITOT+RIHSF2
318C
319C IITOT and RITOT now contain the total imaginary and
320C real parts of the I function
321C
322 SUMISQ=IITOT**2+RITOT**2
323 AS=SSALFS(MH**2)
324 DW=AS**2*G2*MH**3/(32.0*(PI**3)*AMW**2)
325 WID=DW*SUMISQ
326 CALL SSSAVE(IDHHA,WID,IDGL,IDGL,0,0,0)
327100 CONTINUE
328C
329 RETURN
330 END