ISAJET/isasusy/sshff.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE SSHFF
   3 C-----------------------------------------------------------------------
   4 C
   5 C     Calculate all decays higgs -> f fbar, including QCD radiative
   6 C     corrections for quarks.
   7 C
   8 C     Bisset's SETFAC, WDHFFN, QCDRAD
   9 C
  10 C-----------------------------------------------------------------------
  11 #if defined(CERNLIB_IMPNONE)
  12       IMPLICIT NONE
  13 #endif
  14 #include "isajet/sslun.inc"
  15 #include "isajet/sssm.inc"
  16 #include "isajet/sspar.inc"
  17 #include "isajet/sstype.inc"
  18 C
  19       DOUBLE PRECISION PI,SR2,G2,DWID,MHIH,BETA,BEFAC,ALFAC,MH,MF
  20      $,MFRUN,FACTOR,ALAM,MF1,MF2,SUM,MF1RUN,MF2RUN,COLOR,TEMP1
  21      $,QCDFAC
  22       DOUBLE PRECISION MFIFF(9),MFIF1(6),MFIF2(6)
  23       DOUBLE PRECISION SSDLAM,SSMQCD,SSHFF1
  24       REAL WID
  25       INTEGER IH,IDIH,IFF,IDF,ID1,ID2
  26       INTEGER IDIFF(9),IDIF1(6),IDIF2(6)
  27 C
  28       PI=4.*ATAN(1.)
  29       SR2=SQRT(2.)
  30       BETA=ATAN(1./RV2V1)
  31       G2=4.0*PI*ALFAEM/SN2THW
  32 C
  33 C          Loop over HL, HH, HA and fermions
  34 C
  35       MFIFF(1)=AME
  36       IDIFF(1)=IDE
  37       MFIFF(2)=AMMU
  38       IDIFF(2)=IDMU
  39       MFIFF(3)=AMTAU
  40       IDIFF(3)=IDTAU
  41       MFIFF(4)=AMDN
  42       IDIFF(4)=IDDN
  43       MFIFF(5)=AMST
  44       IDIFF(5)=IDST
  45       MFIFF(6)=AMBT
  46       IDIFF(6)=IDBT
  47       MFIFF(7)=AMUP
  48       IDIFF(7)=IDUP
  49       MFIFF(8)=AMCH
  50       IDIFF(8)=IDCH
  51       MFIFF(9)=AMTP
  52       IDIFF(9)=IDTP
  53 C
  54       DO 100 IH=1,3
  55         IF(IH.EQ.1) THEN
  56           MH=AMHL
  57           IDIH=ISHL
  58           BEFAC=COS(BETA)
  59           ALFAC=SIN(ALFAH)
  60         ELSEIF(IH.EQ.2) THEN
  61           MH=AMHH
  62           IDIH=ISHH
  63           BEFAC=COS(BETA)
  64           ALFAC=COS(ALFAH)
  65         ELSE
  66           MH=AMHA
  67           IDIH=ISHA
  68           BEFAC=1/TAN(BETA)
  69           ALFAC=1.
  70         ENDIF
  71 C
  72 C          Down type fermions
  73 C
  74         DO 110 IFF=1,6
  75           MF=MFIFF(IFF)
  76           IDF=IDIFF(IFF)
  77           FACTOR=1.-4.*MF**2/MH**2
  78           IF(FACTOR.LE.0) GO TO 110
  79           FACTOR=SQRT(FACTOR)
  80           IF(IFF.GE.4) THEN
  81              COLOR=3.
  82              MFRUN=SSMQCD(MF,MH)
  83              QCDFAC=SSHFF1(MH,MF,IH)
  84           ELSE
  85              COLOR=1.
  86              MFRUN=MF
  87              QCDFAC=1.
  88           ENDIF
  89           DWID=G2*MFRUN**2*MH*ALFAC**2/(32.*PI*AMW**2*BEFAC**2)
  90           IF(IH.EQ.1.OR.IH.EQ.2) THEN
  91             DWID=DWID*FACTOR**3
  92           ELSEIF(IH.EQ.3) THEN
  93             DWID=DWID*FACTOR
  94           ENDIF
  95           DWID=DWID*COLOR*QCDFAC
  96           WID=DWID
  97           CALL SSSAVE(IDIH,WID,IDF,-IDF,0,0,0)
  98 110     CONTINUE
  99 C
 100 C          Up type fermions
 101 C
 102         IF(IH.EQ.1) THEN
 103           BEFAC=SIN(BETA)
 104           ALFAC=COS(ALFAH)
 105         ELSEIF(IH.EQ.2) THEN
 106           BEFAC=SIN(BETA)
 107           ALFAC=SIN(ALFAH)
 108         ELSE
 109           BEFAC=TAN(BETA)
 110           ALFAC=1.
 111         ENDIF
 112         DO 120 IFF=7,9
 113           MF=MFIFF(IFF)
 114           IDF=IDIFF(IFF)
 115           FACTOR=1.-4.*MF**2/MH**2
 116           IF(FACTOR.LE.0) GO TO 120
 117           FACTOR=SQRT(FACTOR)
 118           MFRUN=SSMQCD(MF,MH)
 119           QCDFAC=SSHFF1(MH,MF,IH)
 120           DWID=G2*MFRUN**2*MH*ALFAC**2/(32.*PI*AMW**2*BEFAC**2)
 121           IF(IH.EQ.1.OR.IH.EQ.2) THEN
 122             DWID=DWID*FACTOR**3
 123           ELSEIF(IH.EQ.3) THEN
 124             DWID=DWID*FACTOR
 125           ENDIF
 126           DWID=3.*DWID*QCDFAC
 127           WID=DWID
 128           CALL SSSAVE(IDIH,WID,IDF,-IDF,0,0,0)
 129 120     CONTINUE
 130 100   CONTINUE
 131 C
 132 C           HC decays. F1 has Iz=+1/2, F2 has Iz=-1/2
 133 C
 134       MFIF1(1)=0
 135       IDIF1(1)=IDNE
 136       MFIF2(1)=AME
 137       IDIF2(1)=IDE
 138       MFIF1(2)=0
 139       IDIF1(2)=IDNM
 140       MFIF2(2)=AMMU
 141       IDIF2(2)=IDMU
 142       MFIF1(3)=0
 143       IDIF1(3)=IDNT
 144       MFIF2(3)=AMTAU
 145       IDIF2(3)=IDTAU
 146       MFIF1(4)=AMUP
 147       IDIF1(4)=IDUP
 148       MFIF2(4)=AMDN
 149       IDIF2(4)=IDDN
 150       MFIF1(5)=AMCH
 151       IDIF1(5)=IDCH
 152       MFIF2(5)=AMST
 153       IDIF2(5)=IDST
 154       MFIF1(6)=AMTP
 155       IDIF1(6)=IDTP
 156       MFIF2(6)=AMBT
 157       IDIF2(6)=IDBT
 158       MH=AMHC
 159 C
 160       DO 200 IFF=1,6
 161         MF1=MFIF1(IFF)
 162         MF2=MFIF2(IFF)
 163         ID1=IDIF1(IFF)
 164         ID2=IDIF2(IFF)
 165         SUM=MF1+MF2
 166         ALAM=SSDLAM(MH**2,MF1**2,MF2**2)
 167         IF(ALAM.LE.0.OR.SUM.GE.MH) GO TO 200
 168         IF(IFF.LE.3) THEN
 169           MF1RUN=MF1
 170           MF2RUN=MF2
 171           COLOR=1
 172         ELSE
 173           MF1RUN=SSMQCD(MF1,MH)
 174           MF2RUN=SSMQCD(MF2,MH)
 175           COLOR=3
 176         ENDIF
 177         TEMP1=MF1RUN**2*1./TAN(BETA)**2+MF2RUN**2*TAN(BETA)**2
 178         TEMP1=TEMP1*(MH**2-MF1**2-MF2**2)-4.*MF1**2*MF2**2
 179         IF (TEMP1.LT.0.0) GO TO 200
 180         DWID=G2*COLOR*SQRT(ALAM)*TEMP1/MH**3/(32.0*PI*AMW**2)
 181         WID=DWID
 182         CALL SSSAVE(ISHC,WID,ID1,-ID2,0,0,0)
 183 200   CONTINUE
 184 C
 185       RETURN
 186       END