ISAJET/isasusy/ssdhll.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE SSDHLL(DELHLL)
   3 C-----------------------------------------------------------------------
   4 C
   5 C          Calculates radiative correction to the
   6 C          H_h-H_l-H_l vertex.
   7 C          calculated by M. Bisset
   8 C
   9 C          This subroutine calculates the
  10 C          radiative correction to the
  11 C          H_h-H_l-H_l vertex which can be
  12 C          important in determining the
  13 C           H_h --> H_l H_l partial decay width.
  14 C
  15 C          Both top and bottom couplings are now
  16 C          included.  Non-degenerate mixed squark
  17 C          masses and A-terms are also included.
  18 C          The D-terms from the squark mass matrix
  19 C          (terms prop. to g**2 * Yukawa coupling)
  20 C          are included as an option:
  21 C                 INRAD = 1 ==> D-TERMS ON
  22 C                 INRAD = 2 ==> D-TERMS OFF    .
  23 C
  24 C         10/18/93 D-terms are now turned on.
  25 C                     INRAD = 1
  26 C
  27 C         There is an arbitrary mass scale that must
  28 C         chosen to avoid dimensionful logarithms.
  29 C         The choice does not matter if D-terms are
  30 C         not included, but it does matter if D-terms
  31 C         are included.
  32 C
  33 C         10/18/93 arbitrary mass scale set to H_h mass
  34 C                       QQQ = AMHH
  35 C
  36 C         It is assumed that the A-terms are real.
  37 C
  38 C-----------------------------------------------------------------------
  39 #if defined(CERNLIB_IMPNONE)
  40       IMPLICIT NONE
  41 #endif
  42 #include "isajet/sslun.inc"
  43 #include "isajet/sssm.inc"
  44 #include "isajet/sspar.inc"
  45 C
  46       REAL PI,PI2,SR2,G2,GP2,GGP,GG1,GG2
  47       REAL TANB,COTB,COSB,SINB,BE
  48       REAL SINB2,COSB2,COS2B,SIN2B
  49       REAL V2,VP2,V,VP,VVP,VPVM,VVPP,MT,MB
  50       REAL MT2,MB2,FT2,FB2,FT,FB,FT4,FB4
  51       REAL MW2,ZAP,QQQ2,EP,EP2,RR,MHP2
  52       REAL ATI,ABI,ATR,ABR,AT2,AB2
  53       REAL MSTL2,MSTR2,MSBL2,MSBR2
  54       REAL TLRM,BLRM
  55       REAL MST1SQ,MST2SQ,MSB1SQ,MSB2SQ
  56       REAL RTT,RBB
  57 C
  58       REAL A0,A1,A2,A1P,A2P,A3,A4
  59       REAL B0,B1,B2,B1P,B2P,B3,B4
  60       REAL MT1R,MT2R,MB1R,MB2R
  61       REAL MT1P,MT2P,MB1P,MB2P
  62       REAL MT1RR,MT2RR,MB1RR,MB2RR
  63       REAL MT1PP,MT2PP,MB1PP,MB2PP
  64       REAL MT1RP,MT2RP,MB1RP,MB2RP
  65       REAL MT1RRR,MT2RRR,MB1RRR,MB2RRR
  66       REAL MT1PRR,MT2PRR,MB1PRR,MB2PRR
  67       REAL MT1RPP,MT2RPP,MB1RPP,MB2RPP
  68       REAL MT1PPP,MT2PPP,MB1PPP,MB2PPP
  69 C
  70       REAL SQVT1,SQVT2,SQVB1,SQVB2
  71       REAL SQVRRR,SQVPPP,SQVPRR,SQVRPP
  72       REAL FVRRR,FVPPP
  73       REAL VRRR,VPPP,VPRR,VRPP
  74 C
  75       REAL ALPHAT,GGP1SQ,ALPHAB,GGP2SQ,TEMPSQ,BSQ
  76       REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS,HIGFRZ
  77       DOUBLE PRECISION SSMQCD
  78 C
  79       REAL CA2,SA2,DVHLL
  80       DOUBLE PRECISION DELHLL
  81 C
  82       INTEGER INRAD,ISPECT,ISPECB
  83 C
  84       MW2=AMW**2
  85       HIGFRZ=SQRT(AMTLSS*AMTRSS)
  86       QQQ2=HIGFRZ**2
  87       INRAD=1
  88       ZAP=1.0
  89 C
  90       PI=4.*ATAN(1.)
  91       PI2=PI**2
  92       SR2=SQRT(2.)
  93       G2=4.*PI*ALFAEM/SN2THW
  94       GP2=G2*SN2THW/(1.-SN2THW)
  95       ASMB=SUALFS(AMBT**2,.36,AMTP,3)
  96       MBMB=AMBT*(1.-4*ASMB/3./PI)
  97       MBQ=SSMQCD(DBLE(MBMB),DBLE(HIGFRZ))
  98       HIGFRZ=SQRT(AMTLSS*AMTRSS)
  99       ASMT=SUALFS(AMTP**2,.36,AMTP,3)
 100       MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
 101      $(ASMT/PI)**2)
 102       MTQ=SSMQCD(DBLE(MTMT),DBLE(HIGFRZ))
 103       MT=MTQ
 104       MB=MBQ
 105       MT2=MT**2
 106       MB2=MB**2
 107       EP=TWOM1
 108       EP2=EP**2
 109       MHP2=AMHA**2
 110       RR=RV2V1
 111       TANB=1.0/RR
 112       COTB=RR
 113       BE=ATAN(1./RV2V1)
 114       SINB=SIN(BE)
 115       COSB=COS(BE)
 116       SINB2=SINB**2
 117       COSB2=COSB**2
 118       SIN2B=SIN(2.*BE)
 119       COS2B=COS(2.*BE)
 120       V2=2.0*MW2*SINB2/G2
 121       VP2=2.0*MW2*COSB2/G2
 122       V=SQRT(V2)
 123       VP=SQRT(VP2)
 124       VVP=SQRT(V2*VP2)
 125       VPVM=VP2-V2
 126       GGP=G2+GP2
 127       GG1=G2-5.0*GP2/3.0
 128       GG2=G2-GP2/3.0
 129       VVPP=2.0*AMZ**2/GGP
 130       FT2=MT2/V2
 131       FB2=MB2/VP2
 132       FT=SQRT(FT2)
 133       FB=SQRT(FB2)
 134       FT4 = FT2**2
 135       FB4 = FB2**2
 136 C
 137 C      (AAT and AAB are also assumed to be real)
 138 C
 139       ATR=AAT
 140       ABR=AAB
 141       ATI=0.0
 142       ABI=0.0
 143       AT2=ATR**2+ATI**2
 144       AB2=ABR**2+ABI**2
 145 C
 146       MSTL2=AMTLSS**2
 147       MSTR2=AMTRSS**2
 148       MSBL2=AMBLSS**2
 149       MSBR2=AMBRSS**2
 150       TLRM=MSTL2-MSTR2
 151       BLRM=MSBL2-MSBR2
 152 C
 153 C      UNFORTUNATELY, I HAVE USED MY OLD CONVENTION
 154 C      FOR THE STOP AND SBOTTOM EIGENVALUES HERE
 155 C      (T1 <==> T2 OF NOTATION IN X. TATA'S AND OTHER
 156 C      PEOPLE'S NOTATION).  SO THE NEXT FOUR LINES ARE
 157 C      A FIX-UP UNTIL I GO THROUGH AND CHANGE THE
 158 C      NOTATION THROUGHOUT THIS SUBROUTINE.
 159 C
 160       MST2SQ=AMT1SS**2
 161       MST1SQ=AMT2SS**2
 162       MSB2SQ=AMB1SS**2
 163       MSB1SQ=AMB2SS**2
 164 C
 165 C
 166 C      Calculation of radiative correction to
 167 C      the H_H-H_l-H_l vertex
 168 C
 169 C
 170 C                  STOP TERMS
 171 C
 172       ISPECT=0
 173       RTT=(TLRM+VPVM*ZAP*GG1/4.0)**2
 174      $      +4.0*MT2*(EP*COTB+ATR)**2+4.0*MT2*ATI**2
 175 C
 176       IF(RTT.GT.0.0) THEN
 177         A0=SQRT(RTT)
 178         A1=-V*ZAP*GG1*(TLRM+ZAP*VPVM*GG1/4.0)/SR2
 179         A1=A1+4.0*SR2*FT*MT*(AT2+EP*ATR*COTB)
 180         A2=-ZAP*GG1*(TLRM+ZAP*VPVM*GG1/4.0)/2.0
 181         A2=A2 +V2*ZAP*GG1**2/4.0 +4.0*FT2*AT2
 182         A1P=VP*ZAP*GG1*(TLRM+ZAP*VPVM*GG1/4.0)/SR2
 183         A1P=A1P+4.0*SR2*FT*MT*EP*(ATR+EP*COTB)
 184         A2P=ZAP*GG1*(TLRM+ZAP*VPVM*GG1/4.0)/2.0
 185         A2P=A2P +VP2*ZAP*GG1**2/4.0 +4.0*FT2*EP2
 186         A3=SR2*ZAP*GG1**2/8.0
 187         A4=-VVP*ZAP*GG1**2/4.0 +4.0*FT2*EP*ATR
 188 C
 189         MT1R=SR2*FT*MT-SR2*V*ZAP*GGP/8.0 +A1/(4.0*A0)
 190         MT2R=SR2*FT*MT-SR2*V*ZAP*GGP/8.0 -A1/(4.0*A0)
 191         MT1P=SR2*VP*ZAP*GGP/8.0 +A1P/(4.0*A0)
 192         MT2P=SR2*VP*ZAP*GGP/8.0 -A1P/(4.0*A0)
 193         MT1RR=FT2 -ZAP*GGP/8.0 -A1**2/(8.0*A0**3) +A2/(4.0*A0)
 194         MT2RR=FT2 -ZAP*GGP/8.0 +A1**2/(8.0*A0**3) -A2/(4.0*A0)
 195         MT1PP=ZAP*GGP/8.0 -A1P**2/(8.0*A0**3) +A2P/(4.0*A0)
 196         MT2PP=ZAP*GGP/8.0 +A1P**2/(8.0*A0**3) -A2P/(4.0*A0)
 197         MT1RRR=3.0*A1**3/(16.0*A0**3)
 198         MT1RRR=MT1RRR/(A0**2) -3.0*A1*A2/(8.0*A0**3)
 199      $                            +3.0*V*A3/(4.0*A0)
 200         MT2RRR=-MT1RRR
 201         MT1PPP=3.0*A1P**3/(16.0*A0**3)
 202         MT1PPP=MT1PPP/(A0**2) -3.0*A1P*A2P/(8.0*A0**3)
 203      $                             +3.0*VP*A3/(4.0*A0)
 204         MT2PPP=-MT1PPP
 205         MT1RP=-A1*A1P/(8.0*A0**3) +A4/(4.0*A0)
 206         MT2RP=-MT1RP
 207         MT1PRR=3.0*A1P*A1**2/(16.0*A0**3)
 208         MT1PRR=MT1PRR/(A0**2)
 209      $          -(A2*A1P+2.0*A1*A4)/(8.0*A0**3) -VP*A3/(4.0*A0)
 210         MT2PRR=-MT1PRR
 211         MT1RPP=3.0*A1*A1P**2/(16.0*A0**3)
 212         MT1RPP=MT1RPP/(A0**2)
 213      $          -(A1*A2P+2.0*A1P*A4)/(8.0*A0**3) -V*A3/(4.0*A0)
 214         MT2RPP=-MT1RPP
 215       ELSEIF(RTT.EQ.0.0) THEN
 216         IF(INRAD.EQ.2.OR.TANB.EQ.1.0) THEN
 217           IF(EP.EQ.0.0.AND.TLRM.EQ.0.0) THEN
 218             IF(ATR.EQ.0.0.AND.ATI.EQ.0.0) THEN
 219               ISPECT=1
 220               MT1R=SR2*V*FT2
 221               MT2R=SR2*V*FT2
 222               MT1P=0.0
 223               MT2P=0.0
 224               MT1RR=FT2
 225               MT2RR=FT2
 226               MT1PP=0.0
 227               MT2PP=0.0
 228               MT1RRR=0.0
 229               MT2RRR=0.0
 230               MT1PPP=0.0
 231               MT2PPP=0.0
 232               MT1RP=0.0
 233               MT2RP=0.0
 234               MT1PRR=0.0
 235               MT2PRR=0.0
 236               MT1RPP=0.0
 237               MT2RPP=0.0
 238             ENDIF
 239           ENDIF
 240         ENDIF
 241       ENDIF
 242       IF(RTT.NE.0.0 .OR. ISPECT.EQ.1) THEN
 243         SQVT1=2.0*(3.0*MT1R*MT1RR+MST1SQ*MT1RRR)
 244         SQVT1=SQVT1*LOG(MST1SQ/QQQ2)
 245         SQVT1=SQVT1 +2.0*MT1R**3/MST1SQ +9.0*MT1R*MT1RR
 246         SQVT1=SQVT1+MST1SQ*MT1RRR
 247         SQVT2=2.0*(3.0*MT2R*MT2RR+MST2SQ*MT2RRR)
 248         SQVT2=SQVT2*LOG(MST2SQ/QQQ2)
 249         SQVT2=SQVT2 +2.0*MT2R**3/MST2SQ +9.0*MT2R*MT2RR
 250         SQVT2=SQVT2+MST2SQ*MT2RRR
 251         SQVRRR=SQVT1+SQVT2
 252 C
 253         SQVT1=2.0*(3.0*MT1P*MT1PP+MST1SQ*MT1PPP)
 254         SQVT1=SQVT1*LOG(MST1SQ/QQQ2)
 255         SQVT1=SQVT1 +2.0*MT1P**3/MST1SQ + 9.0*MT1P*MT1PP
 256         SQVT1=SQVT1+MST1SQ*MT1PPP
 257         SQVT2=2.0*(3.0*MT2P*MT2PP+MST2SQ*MT2PPP)
 258         SQVT2=SQVT2*LOG(MST2SQ/QQQ2)
 259         SQVT2=SQVT2 +2.0*MT2P**3/MST2SQ +9.0*MT2P*MT2PP
 260         SQVT2=SQVT2 +MST2SQ*MT2PPP
 261         SQVPPP = SQVT1 + SQVT2
 262 C
 263         SQVT1=2.0*MT1R*MT1RP+MT1P*MT1RR+MST1SQ*MT1PRR
 264         SQVT1=2.0*SQVT1*LOG(MST1SQ/QQQ2)
 265         SQVT1=SQVT1 +2.0*MT1P*MT1R**2/MST1SQ
 266         SQVT1=SQVT1+3.0*MT1P*MT1RR+6.0*MT1R*MT1RP
 267         SQVT1=SQVT1+MST1SQ*MT1PRR
 268         SQVT2=2.0*MT2R*MT2RP+MT2P*MT2RR+MST2SQ*MT2PRR
 269         SQVT2=2.0*SQVT2*LOG(MST2SQ/QQQ2)
 270         SQVT2=SQVT2 +2.0*MT2P*MT2R**2/MST2SQ
 271         SQVT2=SQVT2+3.0*MT2P*MT2RR+6.0*MT2R*MT2RP
 272         SQVT2=SQVT2+MST2SQ*MT2PRR
 273         SQVPRR=SQVT1+SQVT2
 274 C
 275         SQVT1=2.0*MT1P*MT1RP+MT1R*MT1PP+MST1SQ*MT1RPP
 276         SQVT1=2.0*SQVT1*LOG(MST1SQ/QQQ2)
 277         SQVT1=SQVT1 +2.0*MT1R*MT1P**2/MST1SQ
 278         SQVT1=SQVT1+3.0*MT1R*MT1PP+6.0*MT1P*MT1RP
 279         SQVT1=SQVT1+MST1SQ*MT1RPP
 280         SQVT2=2.0*MT2P*MT2RP+MT2R*MT2PP+MST2SQ*MT2RPP
 281         SQVT2=2.0*SQVT2*LOG(MST2SQ/QQQ2)
 282         SQVT2=SQVT2 +2.0*MT2R*MT2P**2/MST2SQ
 283         SQVT2=SQVT2+3.0*MT2R*MT2PP+6.0*MT2P*MT2RP
 284         SQVT2=SQVT2+MST2SQ*MT2RPP
 285         SQVRPP=SQVT1+SQVT2
 286 C
 287         FVRRR=-2.0*SR2*FT4*V*(6.0*LOG(MT2/QQQ2) + 13.0)
 288       ENDIF
 289 C
 290       IF(RTT.EQ.0.0 .AND. ISPECT.EQ.0) THEN
 291         ALPHAT=(MSTL2 + MSTR2)/2.0 + MT2
 292         ALPHAT=ALPHAT +VP2*(1.0-TANB**2)*ZAP*GGP/8.0
 293         GGP1SQ= ZAP*GGP**2 +ZAP*GG1**2
 294 C
 295         SQVRRR=12.0*FT4*LOG(ALPHAT/MT2)
 296         TEMPSQ=-FT2*ZAP*GGP +GGP1SQ/16.0
 297         SQVRRR=SQVRRR +3.0*TEMPSQ*LOG(ALPHAT/QQQ2)
 298         SQVRRR=SQVRRR -8.0*FT4 -9.0*FT2*ZAP*GGP/2.0
 299         SQVRRR=SQVRRR +9.0*GGP1SQ/32.0
 300         TEMPSQ=8.0*V2*(FT2-ZAP*GGP/8.0)**2
 301         TEMPSQ=TEMPSQ +3.0*V2*ZAP*GG1**2/8.0
 302         TEMPSQ=TEMPSQ +6.0*FT2*EP2*COTB**2
 303         SQVRRR=SQVRRR +TEMPSQ*(FT2-ZAP*GGP/8.0)/ALPHAT
 304         SQVRRR=SQVRRR*SR2*V
 305 C
 306         SQVPPP=3.0*GGP1SQ*(2.0*LOG(ALPHAT/QQQ2)+3.0)/32.0
 307         TEMPSQ=ZAP*GGP*(ZAP*GGP**2+3.0*GG1**2)*VP2/ALPHAT/64.0
 308         SQVPPP=SQVPPP+TEMPSQ
 309         TEMPSQ=3.0*FT2*EP2*ZAP*GGP/ALPHAT/4.0
 310         SQVPPP=(SQVPPP+TEMPSQ)*SR2*VP
 311 C
 312         TEMPSQ=FT2*ZAP*GGP -GGP1SQ/8.0
 313         SQVPRR=TEMPSQ*(2.0*LOG(ALPHAT/QQQ2)+3.0)
 314         TEMPSQ=4.0*ZAP*GGP*(FT2-ZAP*GGP/8.0)-FT2*ZAP*GG1**2
 315         TEMPSQ=TEMPSQ +3.0*ZAP*GGP*GG1**2/16.0
 316         TEMPSQ=V2*TEMPSQ+EP2*FT2*ZAP*GGP*(2.0+COTB**2)
 317         TEMPSQ=TEMPSQ-16.0*EP2*FT4
 318         SQVPRR=(SQVPRR+TEMPSQ/ALPHAT)*SR2*VP/4.0
 319 C
 320         TEMPSQ=FT2*ZAP*GGP -GGP1SQ/8.0
 321         SQVRPP=TEMPSQ*(2.0*LOG(ALPHAT/QQQ2)+3.0)
 322         TEMPSQ=GGP1SQ*(FT2-ZAP*GGP/8.0)-ZAP*GGP*GG1**2/4.0
 323         TEMPSQ=VP2*TEMPSQ/2.0 +8.0*EP2*FT4
 324         TEMPSQ=TEMPSQ+EP2*FT2*ZAP*GGP*(1.0+2.0*COTB**2)
 325         SQVRPP=(SQVRPP+TEMPSQ/ALPHAT)*SR2*V/4.0
 326 C
 327         FVRRR = 0.0
 328 C
 329 C        Fermion part (FRRR) is already combined
 330 C        with the squark part.
 331 C
 332       ENDIF
 333 C
 334 C
 335 C                  SBOTTOM TERMS
 336 C
 337       ISPECB=0
 338       RBB=(BLRM-VPVM*ZAP*GG2/4.0)**2
 339      $      +4.0*MB2*(EP*TANB+ABR)**2+4.0*MB2*ABI**2
 340 C
 341       IF(RBB.GT.0.0) THEN
 342         B0=SQRT(RBB)
 343         B1=V*ZAP*GG2*(BLRM-ZAP*VPVM*GG2/4.0)/SR2
 344         B1=B1+4.0*SR2*FB*MB*EP*(ABR+EP*TANB)
 345         B2=ZAP*GG2*(BLRM-ZAP*VPVM*GG2/4.0)/2.0
 346         B2=B2 +V2*ZAP*GG2**2/4.0 +4.0*FB2*EP2
 347         B1P=-VP*ZAP*GG2*(BLRM-ZAP*VPVM*GG2/4.0)/SR2
 348         B1P=B1P+4.0*SR2*FB*MB*(AB2+EP*ABR*TANB)
 349         B2P=-ZAP*GG2*(BLRM-ZAP*VPVM*GG2/4.0)/2.0
 350         B2P=B2P +VP2*ZAP*GG2**2/4.0 +4.0*FB2*AB2
 351         B3=SR2*ZAP*GG2**2/8.0
 352         B4=-VVP*ZAP*GG2**2/4.0 +4.0*FB2*EP*ABR
 353 C
 354         MB1R=SR2*V*ZAP*GGP/8.0 +B1/(4.0*B0)
 355         MB2R=SR2*V*ZAP*GGP/8.0 -B1/(4.0*B0)
 356         MB1P=SR2*FB*MB -SR2*VP*ZAP*GGP/8.0 +B1P/(4.0*B0)
 357         MB2P=SR2*FB*MB -SR2*VP*ZAP*GGP/8.0 -B1P/(4.0*B0)
 358         MB1RR=ZAP*GGP/8.0 -B1**2/(8.0*B0**3) +B2/(4.0*B0)
 359         MB2RR=ZAP*GGP/8.0 +B1**2/(8.0*B0**3) -B2/(4.0*B0)
 360         MB1PP=FB2 -ZAP*GGP/8.0
 361         MB2PP=MB1PP
 362         MB1PP=MB1PP -B1P**2/(8.0*B0**3) +B2P/(4.0*B0)
 363         MB2PP=MB2PP +B1P**2/(8.0*B0**3) -B2P/(4.0*B0)
 364         MB1RRR=3.0*B1**3/(16.0*B0**3)
 365         MB1RRR=MB1RRR/(B0**2) -3.0*B1*B2/(8.0*B0**3)
 366      $                            +3.0*V*B3/(4.0*B0)
 367         MB2RRR=-MB1RRR
 368         MB1PPP=3.0*B1P**3/(16.0*B0**3)
 369         MB1PPP=MB1PPP/(B0**2) -3.0*B1P*B2P/(8.0*B0**3)
 370         MB1PPP=MB1PPP +3.0*VP*B3/(4.0*B0)
 371         MB2PPP=-MB1PPP
 372         MB1RP=-B1*B1P/(8.0*B0**3) +B4/(4.0*B0)
 373         MB2RP=-MB1RP
 374         MB1PRR=3.0*B1P*B1**2/(16.0*B0**3)
 375         MB1PRR=MB1PRR/(B0**2) -(B2*B1P+2.0*B1*B4)/(8.0*B0**3)
 376         MB1PRR=MB1PRR -VP*B3/(4.0*B0)
 377         MB2PRR=-MB1PRR
 378         MB1RPP=3.0*B1*B1P**2/(16.0*B0**3)
 379         MB1RPP=MB1RPP/(B0**2) -(B1*B2P+2.0*B1P*B4)/(8.0*B0**3)
 380         MB1RPP=MB1RPP -V*B3/(4.0*B0)
 381         MB2RPP=-MB1RPP
 382       ELSEIF(RBB.EQ.0.0) THEN
 383         IF(INRAD.EQ.2.OR.TANB.EQ.1.0) THEN
 384           IF(EP.EQ.0.0.AND.BLRM.EQ.0.0) THEN
 385             IF(ABR.EQ.0.0.AND.ABI.EQ.0.0) THEN
 386               ISPECB=1
 387               MB1R=0.0
 388               MB2R=0.0
 389               MB1P=SR2*VP*FB2
 390               MB2P=SR2*VP*FB2
 391               MB1RR=0.0
 392               MB2RR=0.0
 393               MB1PP=FB2
 394               MB2PP=FB2
 395               MB1RRR=0.0
 396               MB2RRR=0.0
 397               MB1PPP=0.0
 398               MB2PPP=0.0
 399               MB1RP=0.0
 400               MB2RP=0.0
 401               MB1PRR=0.0
 402               MB1PRR=0.0
 403               MB2PRR=0.0
 404               MB1RPP=0.0
 405               MB2RPP=0.0
 406             ENDIF
 407           ENDIF
 408         ENDIF
 409       ENDIF
 410 C
 411       IF(RBB.NE.0.0 .OR. ISPECB.EQ.1) THEN
 412         SQVB1=2.0*(3.0*MB1R*MB1RR+MSB1SQ*MB1RRR)
 413         SQVB1=SQVB1*LOG(MSB1SQ/QQQ2)
 414         SQVB1=SQVB1 +2.0*MB1R**3/MSB1SQ +9.0*MB1R*MB1RR
 415         SQVB1=SQVB1+MSB1SQ*MB1RRR
 416         SQVB2=2.0*(3.0*MB2R*MB2RR+MSB2SQ*MB2RRR)
 417         SQVB2=SQVB2*LOG(MSB2SQ/QQQ2)
 418         SQVB2=SQVB2 +2.0*MB2R**3/MSB2SQ +9.0*MB2R*MB2RR
 419         SQVB2=SQVB2+MSB2SQ*MB2RRR
 420         SQVRRR = SQVRRR + SQVB1 + SQVB2
 421 C
 422         SQVB1=2.0*(3.0*MB1P*MB1PP+MSB1SQ*MB1PPP)
 423         SQVB1=SQVB1*LOG(MSB1SQ/QQQ2)
 424         SQVB1=SQVB1 +2.0*MB1P**3/MSB1SQ +9.0*MB1P*MB1PP
 425         SQVB1=SQVB1+MSB1SQ*MB1PPP
 426         SQVB2=2.0*(3.0*MB2P*MB2PP+MSB2SQ*MB2PPP)
 427         SQVB2=SQVB2*LOG(MSB2SQ/QQQ2)
 428         SQVB2=SQVB2 +2.0*MB2P**3/MSB2SQ +9.0*MB2P*MB2PP
 429         SQVB2=SQVB2+MSB2SQ*MB2PPP
 430         SQVPPP= SQVPPP+SQVB1+SQVB2
 431 C
 432         SQVB1=2.0*MB1R*MB1RP+MB1P*MB1RR+MSB1SQ*MB1PRR
 433         SQVB1=2.0*SQVB1*LOG(MSB1SQ/QQQ2)
 434         SQVB1=SQVB1 +2.0*MB1P*MB1R**2/MSB1SQ
 435         SQVB1=SQVB1 +3.0*MB1P*MB1RR +6.0*MB1R*MB1RP
 436         SQVB1=SQVB1+MSB1SQ*MB1PRR
 437         SQVB2=2.0*MB2R*MB2RP+MB2P*MB2RR+MSB2SQ*MB2PRR
 438         SQVB2=2.0*SQVB2*LOG(MSB2SQ/QQQ2)
 439         SQVB2=SQVB2 +2.0*MB2P*MB2R**2/MSB2SQ
 440         SQVB2=SQVB2 +3.0*MB2P*MB2RR +6.0*MB2R*MB2RP
 441         SQVB2=SQVB2+MSB2SQ*MB2PRR
 442         SQVPRR=SQVPRR+SQVB1+SQVB2
 443 C
 444         SQVB1=2.0*MB1P*MB1RP+MB1R*MB1PP+MSB1SQ*MB1RPP
 445         SQVB1=2.0*SQVB1*LOG(MSB1SQ/QQQ2)
 446         SQVB1=SQVB1 +2.0*MB1R*MB1P**2/MSB1SQ
 447         SQVB1=SQVB1+3.0*MB1R*MB1PP+6.0*MB1P*MB1RP
 448         SQVB1=SQVB1+MSB1SQ*MB1RPP
 449         SQVB2=2.0*MB2P*MB2RP+MB2R*MB2PP+MSB2SQ*MB2RPP
 450         SQVB2=2.0*SQVB2*LOG(MSB2SQ/QQQ2)
 451         SQVB2=SQVB2 +2.0*MB2R*MB2P**2/MSB2SQ
 452         SQVB2=SQVB2 +3.0*MB2R*MB2PP +6.0*MB2P*MB2RP
 453         SQVB2=SQVB2+MSB2SQ*MB2RPP
 454         SQVRPP=SQVRPP+SQVB1+SQVB2
 455 C
 456         IF(MB2.EQ.0.0) THEN
 457           FVPPP=0.0
 458         ELSE IF(MB2.NE.0.0) THEN
 459           FVPPP=-2.0*SR2*FB4*VP*(6.0*LOG(MB2/QQQ2)+13.0)
 460         ENDIF
 461 C
 462       ENDIF
 463 C
 464       IF(RBB.EQ.0.0 .AND. ISPECB.EQ.0) THEN
 465         ALPHAB=(MSBL2+MSBR2)/2.0 +MB2
 466         ALPHAB=ALPHAB -VP2*(1.0-TANB**2)*ZAP*GGP/8.0
 467         GGP2SQ=ZAP*GGP**2 +ZAP*GG2**2
 468 C
 469         BSQ=3.0*GGP2SQ*(2.0*LOG(ALPHAB/QQQ2)+3.0)/8.0
 470         TEMPSQ=V2*(ZAP*GGP**2 +3.0*ZAP*GG2**2)/16.0
 471      $                                   +3.0*FB2*EP2
 472         BSQ=(BSQ +ZAP*GGP*TEMPSQ/ALPHAB)*SR2*V/4.0
 473         SQVRRR=SQVRRR+BSQ
 474 C
 475         BSQ=12.0*FB4*LOG(ALPHAB/MB2) -8.0*FB4
 476         TEMPSQ=-FB2*ZAP*GGP +GGP2SQ/16.0
 477         BSQ=BSQ+3.0*TEMPSQ*(LOG(ALPHAB/QQQ2)+1.5)
 478         TEMPSQ=8.0*VP2*(FB2-ZAP*GGP/8.0)**2
 479      $     +3.0*VP2*ZAP*GG2**2/8.0 +6.0*FB2*EP2*TANB**2
 480         BSQ=BSQ +(FB2-ZAP*GGP/8.0)*TEMPSQ/ALPHAB
 481         BSQ=BSQ*SR2*VP
 482         SQVPPP=SQVPPP+BSQ
 483 C
 484         TEMPSQ=0.5*(FB2*ZAP*GGP -GGP2SQ/8.0)
 485         BSQ=TEMPSQ*(LOG(ALPHAB/QQQ2)+1.5)
 486         TEMPSQ=(FB2 -ZAP*GGP/8.0)*GGP2SQ -ZAP*GGP*GG2**2/4.0
 487         TEMPSQ=V2*TEMPSQ/4.0 +4.0*FB4*EP2 -FB2*EP2*ZAP*GGP/2.0
 488         TEMPSQ=(TEMPSQ-FB2*EP2*ZAP*GGP*TANB**2)/ALPHAB/2.0
 489         BSQ=(BSQ+TEMPSQ)*SR2*VP
 490         SQVPRR=SQVPRR+BSQ
 491 C
 492         TEMPSQ=0.5*(FB2*ZAP*GGP -GGP2SQ/8.0)
 493         BSQ=TEMPSQ*(LOG(ALPHAB/QQQ2)+1.5)
 494         TEMPSQ=4.0*ZAP*GGP*(FB2 -ZAP*GGP/8.0)**2
 495      $          -FB2*ZAP*GG2**2 +3.0*ZAP*GGP*GG2**2/16.0
 496         TEMPSQ=VP2*TEMPSQ-16.0*FB4*EP2
 497         TEMPSQ=TEMPSQ+FB2*EP2*ZAP*GGP*(TANB**2 +0.5)
 498         BSQ=(BSQ +TEMPSQ/ALPHAB/4.0)*SR2*V
 499         SQVRPP=SQVRPP+BSQ
 500 C
 501         FVPPP=0.0
 502 C
 503 C        Fermion part (FPPP) is already combined
 504 C        with the squark part.
 505 C
 506       ENDIF
 507 C
 508 C
 509       VRRR=3.0*(SQVRRR+FVRRR)/(32.0*PI2)
 510       VRRR=VRRR/6.0
 511 C
 512       VPPP=3.0*(SQVPPP+FVPPP)/(32.0*PI2)
 513       VPPP=VPPP/6.0
 514 C
 515       VPRR=3.0*(SQVPRR)/(32.0*PI2)
 516       VPRR=VPRR/2.0
 517 C
 518       VRPP=3.0*(SQVRPP)/(32.0*PI2)
 519       VRPP=VRPP/2.0
 520 C
 521 C
 522 C      Note in the following that the angle ALFAH
 523 C      calculated in the subroutine SSMHN must
 524 C      be input.
 525 C
 526       CA2=(COS(ALFAH))**2
 527       SA2=(SIN(ALFAH))**2
 528       DVHLL=-VRRR*CA2*SIN(ALFAH)
 529       DVHLL=DVHLL +VPRR*(CA2-2.0*SA2)*COS(ALFAH)
 530       DVHLL=DVHLL +VRPP*(2.0*CA2-SA2)*SIN(ALFAH)
 531       DVHLL=DVHLL +VPPP*SA2*COS(ALFAH)
 532 C
 533       DVHLL=3.0*DVHLL
 534       DVHLL=-DVHLL
 535 C
 536 C      Finally, multiply bt the coefficient of the
 537 C      tree-level Lagrangian level term (COEFF.)
 538 C      so that the answer may be written as:
 539 C        DW = (COEFF.)**2
 540 C            * (TREE-LEVEL ANGULAR DEPENDENCE + DVHLL)
 541 C
 542 C              *(LAMBDA KINEMATIC FCN)**0.5/(8*PI*MHH**3)
 543 C
 544 C
 545       DVHLL=4.0*SQRT((1.-SN2THW)/G2)*DVHLL/AMZ
 546 C
 547 C
 548 1000  DELHLL=DVHLL
 549       RETURN
 550       END