ISAJET/isasusy/surg06.F

   1 #include "isajet/pilot.h"
   2 C-----------------------------------------------------------------
   3       SUBROUTINE SURG06(T,GY,F)
   4 C-----------------------------------------------------------------
   5 C
   6 C     Right hand side of truncated renormalization group equations
   7 C          dG_i/dT = F_i(G)
   8 C     using a single common scale MSUSY for SUSY mass thresholds.
   9 C     Added GY(7) for neutrino Yukawa coupling on 9/24/99
  10 C     THIS INCLUDES TWO-LOOP YUKAWAS FOR MSSM ONLY 11/18/99
  11 C
  12 #if defined(CERNLIB_IMPNONE)
  13       IMPLICIT NONE
  14 #endif
  15 #include "isajet/sslun.inc"
  16 #include "isajet/sugpas.inc"
  17       REAL T,GY(7),F(7)
  18       REAL SINB,BETA,COSB,B1,THTOP,B3,B2,TH2LP,B12,B11,B21,B13,
  19      $MT,TANB,PI,Q,B33,B23,B22,B32,B31
  20       REAL A11,A12,A13,A21,A22,A23,A31,A32,A33
  21       REAL C11,C12,C13,C21,C22,C23,C31,C32,C33
  22       REAL D11,D12,D13,D21,D22,D23,D31,D32,D33
  23       REAL C1,C2,C3,CP1,CP2,CP3,CPP1,CPP2,CPP3,BSY1,BSY2,BSY3
  24       INTEGER NSL,NSD,NSH,NSE,NSU,NSQ,NU,NSG,NSW,NH,NN,NE,ND
  25 c      DATA MZ/91.17/
  26       DATA ND/3/,NE/3/,NN/3/
  27       DATA B11/7.96/,B12/5.4/,B13/17.6/,B21/1.8/,B22/25./,B23/24./
  28       DATA B31/2.2/,B32/9./,B33/14./
  29       DATA A11/5.2/,A12/2.8/,A13/3.6/,A21/6./,A22/6./,A23/2./
  30       DATA A31/4./,A32/4./,A33/0./
  31       DATA C11/1.7/,C12/.5/,C13/1.5/,C21/1.5/,C22/1.5/,C23/.5/
  32       DATA C31/2./,C32/2./,C33/0./
  33       DATA D11/3.98/,D12/2.7/,D13/8.8/,D21/.9/,D22/5.833/,D23/12./
  34       DATA D31/1.1/,D32/4.5/,D33/-26./
  35       DATA C1/.8666667/,C2/3./,C3/5.333333/
  36       DATA CP1/.4666667/,CP2/3./,CP3/5.333333/
  37       DATA CPP1/1.8/,CPP2/3./,CPP3/0./
  38       DATA BSY1/6.6/,BSY2/1./,BSY3/-3./
  39       TANB=XTANB
  40       MT=AMT
  41 C-----THESE ARE VALID FROM MZ TO MGUT
  42       Q=MGUT*EXP(T)
  43       PI=4.*ATAN(1.)
  44       BETA=ATAN(TANB)
  45       SINB=SIN(BETA)
  46       COSB=SQRT(1.-SINB**2)
  47       IF (Q.GT.MSUSY) THEN
  48         NSQ=3
  49         NSU=3
  50         NSD=3
  51         NSL=3
  52         NSE=3
  53         NSH=2
  54         NSW=1
  55         NSG=1
  56         NH=2
  57         TH2LP=1.
  58       ELSE
  59         NSQ=0
  60         NSU=0
  61         NSD=0
  62         NSL=0
  63         NSE=0
  64         NSH=0
  65         NSW=0
  66         NSG=0
  67         NH=1
  68         TH2LP=0.
  69       END IF
  70       IF (Q.GT.MT) THEN
  71         NU=3
  72         THTOP=1.
  73       ELSE
  74         NU=2
  75         THTOP=0.
  76       END IF
  77       TH2LP=1.
  78       THTOP=1.
  79       B1=2.*(17*NU/12.+5*ND/12.+5*NE/4.+NN/4.)/5.+
  80      $ NSQ/30.+4*NSU/15.+NSD/15.+NSL/10.+NSE/5.+
  81      $ 1.*NSH/5.+1.*NH/10.
  82       B2=-22./3.+.5*(NU+ND)+1.*(NE+NN)/6.+
  83      $ 1.*NSQ/2.+1.*NSL/6.+1.*NSH/3.+1.*NH/6.+4.*NSW/3.
  84       B3=2.*(NU+ND)/3.+1.*NSQ/3.+1.*NSU/6.+1.*NSD/6.+2.*NSG-11.
  85       IF (Q.GT.MSUSY) THEN
  86       F(1)=GY(1)/16./PI**2*(B1*GY(1)**2+TH2LP/16./PI**2*GY(1)**2*
  87      $(B11*GY(1)**2+B12*GY(2)**2+B13*GY(3)**2-A11*GY(6)**2-A12*GY(5)**2
  88      $-A13*GY(4)**2))
  89       F(2)=GY(2)/16./PI**2*(B2*GY(2)**2+TH2LP/16./PI**2*GY(2)**2*
  90      $(B21*GY(1)**2+B22*GY(2)**2+B23*GY(3)**2-A21*GY(6)**2-A22*GY(5)**2
  91      $-A23*GY(4)**2))
  92       F(3)=GY(3)/16./PI**2*(B3*GY(3)**2+TH2LP/16./PI**2*GY(3)**2*
  93      $(B31*GY(1)**2+B32*GY(2)**2+B33*GY(3)**2-A31*GY(6)**2-A32*GY(5)**2
  94      $-A33*GY(4)**2))
  95       ELSE
  96       F(1)=GY(1)/16./PI**2*(B1*GY(1)**2+TH2LP/16./PI**2*GY(1)**2*
  97      $(D11*GY(1)**2+D12*GY(2)**2+D13*GY(3)**2-C11*GY(6)**2-C12*GY(5)**2
  98      $-C13*GY(4)**2))
  99       F(2)=GY(2)/16./PI**2*(B2*GY(2)**2+TH2LP/16./PI**2*GY(2)**2*
 100      $(D21*GY(1)**2+D22*GY(2)**2+D23*GY(3)**2-C21*GY(6)**2-C22*GY(5)**2
 101      $-C23*GY(4)**2))
 102       F(3)=GY(3)/16./PI**2*(B3*GY(3)**2+TH2LP/16./PI**2*GY(3)**2*
 103      $(D31*GY(1)**2+D32*GY(2)**2+D33*GY(3)**2-C31*GY(6)**2-C32*GY(5)**2
 104      $-C33*GY(4)**2))
 105       ENDIF
 106       IF (Q.LT.MSUSY) THEN
 107         F(4)=GY(4)/16./PI**2*(5*GY(4)**2*COSB**2/2.+3*GY(6)**2*SINB**2*
 108      $   THTOP+3*GY(5)**2*COSB**2-9*GY(1)**2/4.-9*GY(2)**2/4.
 109      $   -SINB**2*(3*GY(6)**2*THTOP-3*GY(5)**2-GY(4)**2))
 110         F(5)=GY(5)/16./PI**2*(9*GY(5)**2*COSB**2/2.+3*GY(6)**2*SINB**2*
 111      $   THTOP/2.+GY(4)**2*COSB**2-GY(1)**2/4.-9*GY(2)**2/4.
 112      $   -8*GY(3)**2-SINB**2*(3*GY(6)**2*THTOP-3*GY(5)**2-GY(4)**2))
 113         F(6)=GY(6)/16./PI**2*(9*GY(6)**2*SINB**2/2.*THTOP+
 114      $   3*GY(5)**2*COSB**2/2.+GY(4)**2*COSB**2-17.*GY(1)**2/20.
 115      $   -9*GY(2)**2/4.-8*GY(3)**2+COSB**2*
 116      $   (3*GY(6)**2*THTOP-3*GY(5)**2-GY(4)**2))
 117       ELSE
 118         F(4)=GY(4)/16./PI**2*(4*GY(4)**2+3*GY(5)**2+GY(7)**2-
 119      $        9*GY(1)**2/5.-3*GY(2)**2+
 120      $       TH2LP/16./PI**2*((CPP1*BSY1+CPP1**2/2.)*GY(1)**4+
 121      $       (CPP2*BSY2+CPP2**2/2.)*GY(2)**4+(CPP3*BSY3+CPP3**2/2.)*
 122      $       GY(3)**4+9*GY(1)**2*GY(2)**2/5.+
 123      $       GY(5)**2*(-.4*GY(1)**2+16*GY(3)**2)+
 124      $       GY(4)**2*(1.2*GY(1)**2+6*GY(2)**2)-
 125      $       (3*GY(6)**2*GY(5)**2+9*GY(5)**4+9*GY(5)**2*GY(4)**2+
 126      $       10*GY(4)**4)))
 127         F(5)=GY(5)/16./PI**2*(6*GY(5)**2+GY(6)**2*THTOP+
 128      $        GY(4)**2-7*GY(1)**2/15.-3*GY(2)**2-16*GY(3)**2/3.+
 129      $       TH2LP/16./PI**2*((CP1*BSY1+CP1**2/2.)*GY(1)**4+
 130      $       (CP2*BSY2+CP2**2/2.)*GY(2)**4+(CP3*BSY3+CP3**2/2.)*
 131      $       GY(3)**4+GY(1)**2*GY(2)**2+8*GY(1)**2*GY(3)**2/9.+
 132      $       8*GY(2)**2*GY(3)**2+.8*GY(6)**2*GY(1)**2+
 133      $       GY(5)**2*(.4*GY(1)**2+6*GY(2)**2+16*GY(3)**2)+
 134      $       1.2*GY(4)**2*GY(1)**2-
 135      $       (22*GY(5)**4+5*GY(6)**2*GY(5)**2+3*GY(5)**2*GY(4)**2+
 136      $       3*GY(4)**4+5*GY(6)**4)))
 137         F(6)=GY(6)/16./PI**2*(6*GY(6)**2*THTOP+GY(5)**2+GY(7)**2-
 138      $        13*GY(1)**2/15.-3*GY(2)**2-16*GY(3)**2/3.+
 139      $       TH2LP/16./PI**2*((C1*BSY1+C1**2/2.)*GY(1)**4+
 140      $       (C2*BSY2+C2**2/2.)*GY(2)**4+(C3*BSY3+C3**2/2.)*
 141      $       GY(3)**4+GY(1)**2*GY(2)**2+136*GY(1)**2*GY(3)**2/45.+
 142      $       8*GY(2)**2*GY(3)**2+GY(6)**2*(1.2*GY(1)**2+6*GY(2)**2+
 143      $       16*GY(3)**2)+.4*GY(5)**2*GY(1)**2-
 144      $       (22*GY(6)**4+5*GY(6)**2*GY(5)**2+5*GY(5)**4+
 145      $       GY(5)**2*GY(4)**2)))
 146       END IF
 147 C     THE NEUTRINO YUKAWA IS TO 1-LOOP ONLY FOR NOW...
 148       IF (Q.GT.AMNRMJ) THEN
 149         F(7)=GY(7)/16./PI**2*(3*GY(6)**2+GY(4)**2+4*GY(7)**2-
 150      $        3*GY(2)**2-3*GY(1)**2/5.)
 151       ELSE
 152         F(7)=0.
 153       END IF
 154       RETURN
 155       END