ISAJET/isasusy/sugeff.F

   1 #include "isajet/pilot.h"
   2 C-----------------------------------------------------------------
   3       SUBROUTINE SUGEFF(G0,SIG1,SIG2)
   4 C-----------------------------------------------------------------
   5 C
   6 C     Compute Higgs mass shift due to 1-loop effective potential
   7 C
   8 #if defined(CERNLIB_IMPNONE)
   9       IMPLICIT NONE
  10 #endif
  11 #include "isajet/sslun.inc"
  12 #include "isajet/sugpas.inc"
  13 #include "isajet/sugmg.inc"
  14       REAL G0(29),SIG1,SIG2
  15       REAL DT1,DELT1S,SIG1T1,DMSDV2,FB,FT,MST2,MSB2,MSB1,SIG2B1,
  16      $SIG1B1,SIG2B2,SIG1B2,DB1,SIG1T2,SIG2T1,DELB1S,SIG2T2,MST1,COS2W,
  17      $MT,PI,COTB,TANB,MB,SIG2B,SIG1B,G,GGP,SIG2T,E,FAC,SIG1T,QS,
  18      $BETA,SINB,COSB,FAC4,FL,ML,SIG1L,SIG2L,MSL1,MSL2,
  19      $DELL1S,DL1,SIG1L1,SIG2L1,SIG1L2,SIG2L2
  20 C
  21       G=G2
  22       TANB=XTANB
  23       COS2W=1.-XW
  24       COTB=1./TANB
  25       BETA=ATAN(TANB)
  26       SINB=SIN(BETA)
  27       COSB=COS(BETA)
  28       PI=4.*ATAN(1.)
  29       FAC=3./8./PI**2
  30       FAC4=FAC/3.
  31       E=EXP(1.)
  32       QS=HIGFRZ**2
  33 C-----CALCULATE TOP AND BOTTOM CONTRIBUTIONS; USE RUNNING MASSES--
  34       FL=G0(4)
  35       FB=G0(5)
  36       FT=G0(6)
  37       ML=FL*VEV*COSB
  38       MT=FT*VEV*SINB
  39       MB=FB*VEV*COSB
  40       SIG1T=0.
  41       SIG2T=-FAC*MT**2*G0(6)**2*LOG(MT**2/E/QS)
  42       SIG1B=-FAC*MB**2*G0(5)**2*LOG(MB**2/E/QS)
  43       SIG2B=0.
  44       SIG1L=-FAC4*ML**2*G0(4)**2*LOG(ML**2/E/QS)
  45       SIG2L=0.
  46       GGP=(G**2+GP**2)/2.
  47       MST1=MSS(12)
  48       MST2=MSS(13)
  49       MSB1=MSS(10)
  50       MSB2=MSS(11)
  51       MSL1=MSS(21)
  52       MSL2=MSS(22)
  53 C-----CALCULATE STOP_1 CONTRIBUTION -------------------------------
  54       DELT1S=(.5*(G0(24)-G0(23))+(8*COS2W-5.)*GGP*
  55      $        (VP*VP-V*V)/12.)**2+G0(6)**2*V*V*(G0(12)-MU*COTB)**2
  56       DT1=.5*(G0(24)-G0(23))+(8*COS2W-5.)*GGP*
  57      $       (VP*VP-V*V)/12.
  58       DMSDV2=GGP/4.-(2*DT1*(8*COS2W-5.)*GGP/12.-
  59      $        FT**2*MU*(G0(12)*TANB-MU))/2./SQRT(DELT1S)
  60       SIG1T1=FAC/2.*MST1**2*LOG(MST1**2/E/QS)*DMSDV2
  61       DMSDV2=-GGP/4.+FT**2-(-2*DT1*(8*COS2W-5.)*GGP/12.+
  62      $        FT**2*G0(12)*(G0(12)-MU*COTB))/2./SQRT(DELT1S)
  63       SIG2T1=FAC/2.*MST1**2*LOG(MST1**2/E/QS)*DMSDV2
  64 C-----CALCULATE STOP_2 CONTRIBUTION -------------------------------
  65       DMSDV2=GGP/4.+(2*DT1*(8*COS2W-5.)*GGP/12.-
  66      $        FT**2*MU*(G0(12)*TANB-MU))/2./SQRT(DELT1S)
  67       SIG1T2=FAC/2.*MST2**2*LOG(MST2**2/E/QS)*DMSDV2
  68       DMSDV2=-GGP/4.+FT**2+(-2*DT1*(8*COS2W-5.)*GGP/12.+
  69      $        FT**2*G0(12)*(G0(12)-MU*COTB))/2./SQRT(DELT1S)
  70       SIG2T2=FAC/2.*MST2**2*LOG(MST2**2/E/QS)*DMSDV2
  71 C-----CALCULATE SBOT_1 CONTRIBUTION -------------------------------
  72       DELB1S=(.5*(G0(24)-G0(22))-(4*COS2W-1.)*GGP*
  73      $  (VP*VP-V*V)/12.)**2+G0(5)**2*VP*VP*(G0(11)-MU*TANB)**2
  74       DB1=.5*(G0(24)-G0(22))-(4*COS2W-1.)*GGP*
  75      $       (VP*VP-V*V)/12.
  76       DMSDV2=-GGP/4.+FB**2-(-2*DB1*(4*COS2W-1.)*GGP/12.+
  77      $        FB**2*G0(11)*(G0(11)-MU*TANB))/2./SQRT(DELB1S)
  78       SIG1B1=FAC/2.*MSB1**2*LOG(MSB1**2/E/QS)*DMSDV2
  79       DMSDV2=GGP/4.-(2*DB1*(4*COS2W-1.)*GGP/12.-
  80      $        FB**2*MU*(G0(11)*COTB-MU))/2./SQRT(DELB1S)
  81       SIG2B1=FAC/2.*MSB1**2*LOG(MSB1**2/E/QS)*DMSDV2
  82 C-----CALCULATE SBOT_2 CONTRIBUTION -------------------------------
  83       DMSDV2=-GGP/4.+FB**2+(-2*DB1*(4*COS2W-1.)*GGP/12.+
  84      $        FB**2*G0(11)*(G0(11)-MU*TANB))/2./SQRT(DELB1S)
  85       SIG1B2=FAC/2.*MSB2**2*LOG(MSB2**2/E/QS)*DMSDV2
  86       DMSDV2=GGP/4.+(2*DB1*(4*COS2W-1.)*GGP/12.-
  87      $        FB**2*MU*(G0(11)*COTB-MU))/2./SQRT(DELB1S)
  88       SIG2B2=FAC/2.*MSB2**2*LOG(MSB2**2/E/QS)*DMSDV2
  89 C-----CALCULATE STAU_1 CONTRIBUTION -------------------------------
  90       DELL1S=(.5*(G0(21)-G0(20))-(4*COS2W-3.)*GGP*
  91      $  (VP*VP-V*V)/4.)**2+G0(4)**2*VP*VP*(G0(10)-MU*TANB)**2
  92       DL1=.5*(G0(21)-G0(20))-(4*COS2W-3.)*GGP*
  93      $       (VP*VP-V*V)/4.
  94       DMSDV2=-GGP/4.+FL**2-(-2*DL1*(4*COS2W-3.)*GGP/4.+
  95      $        FL**2*G0(10)*(G0(10)-MU*TANB))/2./SQRT(DELL1S)
  96       SIG1L1=FAC4/2.*MSL1**2*LOG(MSL1**2/E/QS)*DMSDV2
  97       DMSDV2=GGP/4.-(2*DL1*(4*COS2W-3.)*GGP/4.-
  98      $        FL**2*MU*(G0(10)*COTB-MU))/2./SQRT(DELL1S)
  99       SIG2L1=FAC4/2.*MSL1**2*LOG(MSL1**2/E/QS)*DMSDV2
 100 C-----CALCULATE STAU_2 CONTRIBUTION -------------------------------
 101       DMSDV2=-GGP/4.+FL**2+(-2*DL1*(4*COS2W-3.)*GGP/4.+
 102      $        FL**2*G0(10)*(G0(10)-MU*TANB))/2./SQRT(DELL1S)
 103       SIG1L2=FAC4/2.*MSL2**2*LOG(MSL2**2/E/QS)*DMSDV2
 104       DMSDV2=GGP/4.+(2*DL1*(4*COS2W-3.)*GGP/4.-
 105      $        FL**2*MU*(G0(10)*COTB-MU))/2./SQRT(DELL1S)
 106       SIG2L2=FAC4/2.*MSL2**2*LOG(MSL2**2/E/QS)*DMSDV2
 107 C-----ADD ALL TERMS ------------------------------------------------
 108       SIG1=SIG1B+SIG1B1+SIG1B2+SIG1T+SIG1T1+SIG1T2+
 109      $SIG1L+SIG1L1+SIG1L2
 110       SIG2=SIG2B+SIG2B1+SIG2B2+SIG2T+SIG2T1+SIG2T2+
 111      $SIG2L+SIG2L1+SIG2L2
 112       RETURN
 113       END