ISAJET/isasusy/sshhx.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE SSHHX
   3 C-----------------------------------------------------------------------
   4 C     Calculates the decays Hi -> Hj + X.
   5 C
   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.
  11 C
  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.
  16 C
  17 C
  18 C     Bisset's HIGPRO
  19 C-----------------------------------------------------------------------
  20 #if defined(CERNLIB_IMPNONE)
  21       IMPLICIT NONE
  22 #endif
  23 #include "isajet/sspar.inc"
  24 #include "isajet/sssm.inc"
  25 #include "isajet/sstype.inc"
  26 C
  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
  33 C
  34       PI=4.*ATAN(1.D0)
  35       SR2=SQRT(2.D0)
  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))*
  41      $(ASMT/PI)**2)
  42       MTQ=SSMQCD(DBLE(MTMT),DBLE(HIGFRZ))
  43
  44       BETA=ATAN(1.0/RV2V1)
  45       ALPHA=ALFAH
  46       SW2=SN2THW
  47       CW2=1.-SN2THW
  48 C
  49       LGTST=(1+(AMTLSS/MTQ)**2)*(1+(AMTRSS/MTQ)**2)
  50       LGTST=LOG(LGTST)
  51       SBMA=SIN(BETA-ALPHA)
  52       CBMA=COS(BETA-ALPHA)
  53       CAB2=(DCOS(ALPHA+BETA))**2
  54       SAB2=1.0-CAB2
  55 C
  56 C          hl0 -> ha0 + ha0
  57 C
  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)
  65         DELLPP=DELLPP*LGTST
  66         DWID=(DWID+DELLPP)**2
  67         DWID=DWID*G2*(AMZ**2)/(128.0*PI*CW2*(AMHL**2))
  68         DWID=DWID*SQRT(LAMB1)
  69         WID=DWID
  70         CALL SSSAVE(ISHL,WID,ISHA,ISHA,0,0,0)
  71       ENDIF
  72 C
  73 C          hh -> ha + z
  74 C
  75       IF(AMHH.GT.AMHA+AMZ) THEN
  76         MH=AMHH
  77         M1=AMHA
  78         M2=AMZ
  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))
  83         DWID=DWID*LAMB
  84         WID=DWID
  85         CALL SSSAVE(ISHH,WID,ISHA,IDZ,0,0,0)
  86       ENDIF
  87 C
  88 C          hh -> hl + hl
  89 C
  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)
  94 C
  95 C          Now add hh-hl-hl vertex correction
  96 C
  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.
 100 C
 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)
 109 C
 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.
 115 C
 116         CALL SSDHLL(DELHLL)
 117 C
 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.)
 126 C
 127         DWID=G2*(AMZ**2)*(TEMP+DELHLL)**2
 128         DWID=DWID/(128.0*PI*CW2*(AMHH**2))
 129         DWID=DWID*SQRT(LAMB1)
 130         WID=DWID
 131         CALL SSSAVE(ISHH,WID,ISHL,ISHL,0,0,0)
 132       ENDIF
 133 C
 134 C          hh -> ha + ha
 135 C
 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)
 143         DELHPP=DELHPP*LGTST
 144         DWID=G2*(AMZ**2)*(DWID+DELHPP)**2
 145         DWID=DWID/(128.0*PI*CW2*(AMHH**2))
 146         DWID=DWID*SQRT(LAMB1)
 147         WID=DWID
 148         CALL SSSAVE(ISHH,WID,ISHA,ISHA,0,0,0)
 149       ENDIF
 150 C
 151 C          hh -> hc+ + hc-
 152 C
 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)
 161         DELHCC=DELHCC*LGTST
 162         DWID=G2*AMW**2*(-DWID+DELHCC)**2
 163         DWID=DWID*SQRT(LAMB1)/(16.0*PI*AMHH)
 164         WID=DWID
 165         CALL SSSAVE(ISHH,WID,ISHC,-ISHC,0,0,0)
 166       ENDIF
 167 C
 168 C          ha -> hl + z
 169 C
 170       IF(AMHA.GT.AMHL+AMZ) THEN
 171          MH=AMHA
 172          M1=AMHL
 173          M2=AMZ
 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))
 178          DWID=DWID*LAMB
 179          WID=DWID
 180          CALL SSSAVE(ISHA,WID,ISHL,IDZ,0,0,0)
 181       ENDIF
 182 C
 183 C          hc+ -> w+ + hl
 184 C
 185       IF(AMHC.GT.AMW+AMHL) THEN
 186         MH=AMHC
 187         M1=AMW
 188         M2=AMHL
 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) )
 192         DWID=DWID*LAMB
 193         WID=DWID
 194         CALL SSSAVE(ISHC,WID,ISHL,IDW,0,0,0)
 195       ENDIF
 196 C
 197       RETURN
 198       END