1 #include "isajet/pilot.h"
3 C-----------------------------------------------------------------------
5 C Calculate all decays higgs -> f fbar, including QCD radiative
6 C corrections for quarks.
8 C Bisset's SETFAC, WDHFFN, QCDRAD
10 C-----------------------------------------------------------------------
11 #if defined(CERNLIB_IMPNONE)
14 #include "isajet/sslun.inc"
15 #include "isajet/sssm.inc"
16 #include "isajet/sspar.inc"
17 #include "isajet/sstype.inc"
19 DOUBLE PRECISION PI,SR2,G2,DWID,MHIH,BETA,BEFAC,ALFAC,MH,MF
20 $,MFRUN,FACTOR,ALAM,MF1,MF2,SUM,MF1RUN,MF2RUN,COLOR,TEMP1
22 DOUBLE PRECISION MFIFF(9),MFIF1(6),MFIF2(6)
23 DOUBLE PRECISION SSDLAM,SSMQCD,SSHFF1
25 INTEGER IH,IDIH,IFF,IDF,ID1,ID2
26 INTEGER IDIFF(9),IDIF1(6),IDIF2(6)
31 G2=4.0*PI*ALFAEM/SN2THW
33 C Loop over HL, HH, HA and fermions
77 FACTOR=1.-4.*MF**2/MH**2
78 IF(FACTOR.LE.0) GO TO 110
83 QCDFAC=SSHFF1(MH,MF,IH)
89 DWID=G2*MFRUN**2*MH*ALFAC**2/(32.*PI*AMW**2*BEFAC**2)
90 IF(IH.EQ.1.OR.IH.EQ.2) THEN
95 DWID=DWID*COLOR*QCDFAC
97 CALL SSSAVE(IDIH,WID,IDF,-IDF,0,0,0)
115 FACTOR=1.-4.*MF**2/MH**2
116 IF(FACTOR.LE.0) GO TO 120
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
128 CALL SSSAVE(IDIH,WID,IDF,-IDF,0,0,0)
132 C HC decays. F1 has Iz=+1/2, F2 has Iz=-1/2
166 ALAM=SSDLAM(MH**2,MF1**2,MF2**2)
167 IF(ALAM.LE.0.OR.SUM.GE.MH) GO TO 200
173 MF1RUN=SSMQCD(MF1,MH)
174 MF2RUN=SSMQCD(MF2,MH)
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)
182 CALL SSSAVE(ISHC,WID,ID1,-ID2,0,0,0)