1 #include "isajet/pilot.h"
3 C-----------------------------------------------------------------------
4 C Calculates the decays Hi -> Hj + X.
6 C Includes vertex corrections for triple Higgs couplings due
7 C to top and stop quarks effects.
8 C See Kunszt and Zwirner CERN-TH.6150/91 for all but hh-hc-hc
9 C correction which is in our Higgs-->SUSY paper:
10 C Baer et. al. FSU-HEP-920630 or UH-511-749-92.
12 C The hh-hl-hl vertex correction now includes both
13 C top & bottom and stop and sbottom squark
14 C (non-degenerate with mixing) effects.
15 C A-terms and mu=-2m1 are also included.
19 C-----------------------------------------------------------------------
20 #if defined(CERNLIB_IMPNONE)
23 #include "isajet/sspar.inc"
24 #include "isajet/sssm.inc"
25 #include "isajet/sstype.inc"
27 DOUBLE PRECISION PI,SR2,G2,GP2,BETA,ALPHA,SW2,CW2,LGTST,CBMA
28 $,SBMA,LAMB1,DWID,DELLPP,MH,M1,M2,LAMB,TEMP,DTEMPL,DTEMPR
29 $,DELHLL,DELHPP,DELHCC,CAB2,SAB2
30 DOUBLE PRECISION SSDLAM
31 REAL WID,ASMT,MTMT,MTQ,SUALFS,HIGFRZ
32 DOUBLE PRECISION SSMQCD
36 G2=4.0*PI*ALFAEM/SN2THW
37 GP2=4*PI*ALFAEM/(1-SN2THW)
38 HIGFRZ=SQRT(AMTLSS*AMTRSS)
39 ASMT=SUALFS(AMTP**2,.36,AMTP,3)
40 MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
42 MTQ=SSMQCD(DBLE(MTMT),DBLE(HIGFRZ))
49 LGTST=(1+(AMTLSS/MTQ)**2)*(1+(AMTRSS/MTQ)**2)
53 CAB2=(DCOS(ALPHA+BETA))**2
58 IF(AMHL.GT.2*AMHA) THEN
59 LAMB1=AMHL**2-4.0*AMHA**2
60 DWID=SBMA*COS(2.0*BETA)
61 C Now add hl-hp-hp vertex correction
62 DELLPP=3.0*G2*CW2*(MTQ**4)*COS(ALPHA)
63 DELLPP=DELLPP*(COS(BETA)**2)/(16.0*(PI**2))
64 DELLPP=DELLPP/((AMW**4)*(SIN(BETA))**3)
67 DWID=DWID*G2*(AMZ**2)/(128.0*PI*CW2*(AMHL**2))
70 CALL SSSAVE(ISHL,WID,ISHA,ISHA,0,0,0)
75 IF(AMHH.GT.AMHA+AMZ) THEN
79 LAMB=SSDLAM(MH**2,M1**2,M2**2)
80 DWID=SQRT(G2*CW2)+SQRT(GP2*SW2)
81 DWID=DWID**2*SAB2*SQRT(LAMB)
82 DWID=DWID/(64.0*PI*(AMZ**2)*(AMHH**3))
85 CALL SSSAVE(ISHH,WID,ISHA,IDZ,0,0,0)
90 IF(AMHH.GT.2*AMHL) THEN
91 LAMB1=AMHH**2-4.0*AMHL**2
92 TEMP=CBMA*COS(2.0*ALPHA)
93 TEMP=TEMP+2.0*SBMA*SIN(2.0*ALPHA)
95 C Now add hh-hl-hl vertex correction
97 C The following 8 lines calculate the radiative
98 C hh-hl-hl vertex correction including only
99 C effects from tops and stop squarks.
101 C DTEMPL=3.0*LOG(1.0+(AMTLSS/MTQ)**2)
102 C DTEMPL=DTEMPL-2.0*AMTLSS**2/(AMTLSS**2+MTQ**2)
103 C DTEMPR=3.0*LOG(1.0+(AMTRSS/MTQ)**2)
104 C DTEMPR=DTEMPR-2.0*AMTRSS**2/(AMTRSS**2+MTQ**2)
105 C DELHLL=3.0*G2*CW2*(MTQ**4)*SIN(ALPHA)
106 C DELHLL=DELHLL*(COS(ALPHA)**2)/(PI**2)
107 C DELHLL=DELHLL/(16.0*(AMW**4)*(SIN(BETA))**3)
108 C DELHLL=DELHLL*(DTEMPL+DTEMPR)
110 C The subroutine SSHL calculates the radiative
111 C hh-hl-hl vertex correction including both
112 C top & bottom and stop and sbottom squark
113 C (non-degenerate with mixing) effects.
114 C A-terms and mu=-2m1 are also included.
118 C Note: the variable TEMP in the line below
119 C this is the Lagrangian term (as noted on
120 C page 27 of Prof. Tata's personal Lagrangian
121 C term notes. Thus DELHLL must also be the
122 C Lagrangian entry - not the potential entry.
123 C The subroutine SSHLL IS set up to yield the
124 C the Lagrangian entry. (We must be very careful
125 C about the relative sign between TEMP and DELHLL.)
127 DWID=G2*(AMZ**2)*(TEMP+DELHLL)**2
128 DWID=DWID/(128.0*PI*CW2*(AMHH**2))
129 DWID=DWID*SQRT(LAMB1)
131 CALL SSSAVE(ISHH,WID,ISHL,ISHL,0,0,0)
136 IF(AMHH.GT.2*AMHA) THEN
137 LAMB1=AMHH**2-4.0*AMHA**2
138 DWID=CBMA*COS(2*BETA)
139 C Now add hh-hp-hp vertex correction
140 DELHPP=3.0*G2*CW2*(MTQ**4)*SIN(ALPHA)
141 DELHPP=DELHPP*(COS(BETA)**2)/(16.0*(PI**2))
142 DELHPP=DELHPP/((AMW**4)*(SIN(BETA))**3)
144 DWID=G2*(AMZ**2)*(DWID+DELHPP)**2
145 DWID=DWID/(128.0*PI*CW2*(AMHH**2))
146 DWID=DWID*SQRT(LAMB1)
148 CALL SSSAVE(ISHH,WID,ISHA,ISHA,0,0,0)
153 IF(AMHH.GT.2*AMHC) THEN
154 LAMB1=1.0-4.0*(AMHC**2)/(AMHH**2)
155 DWID=CBMA*COS(2.0*BETA)/(2.0*CW2)
156 DWID=COS(BETA+ALPHA)-DWID
157 C Now add hh-hc-hc vertex correction
158 DELHCC=3.0*G2*MTQ**4*SIN(ALPHA)
159 DELHCC=DELHCC/( SIN(BETA)*(DTAN(BETA))**2 )
160 DELHCC=DELHCC/(32.0*PI**2*AMW**4)
162 DWID=G2*AMW**2*(-DWID+DELHCC)**2
163 DWID=DWID*SQRT(LAMB1)/(16.0*PI*AMHH)
165 CALL SSSAVE(ISHH,WID,ISHC,-ISHC,0,0,0)
170 IF(AMHA.GT.AMHL+AMZ) THEN
174 LAMB=SSDLAM(MH**2,M1**2,M2**2)
175 DWID=SQRT(G2*CW2)+SQRT(GP2*SW2)
176 DWID=DWID**2*CAB2*SQRT(LAMB)
177 DWID=DWID/(64.0*PI*(AMZ**2)*(AMHA**3))
180 CALL SSSAVE(ISHA,WID,ISHL,IDZ,0,0,0)
185 IF(AMHC.GT.AMW+AMHL) THEN
189 LAMB=SSDLAM(MH**2,M1**2,M2**2)
190 DWID=G2*CAB2*SQRT(LAMB)
191 DWID=DWID/( 64.0*PI*(AMW**2)*(AMHC**3) )
194 CALL SSSAVE(ISHC,WID,ISHL,IDW,0,0,0)