ISAJET/isasusy/surg26.F

   1 #include "isajet/pilot.h"
   2 C-----------------------------------------------------------------
   3       SUBROUTINE SURG26(T,G,F)
   4 C-----------------------------------------------------------------
   5 C
   6 C     Right hand side of full renormalization group equations
   7 C          dG_i/dT = F_i(G)
   8 C     using the thresholds MSS for each mass calculated with the
   9 C     couplings G0 frozen by SUGFRZ.
  10 C     Added right neutrino RGE's on 9/24/99
  11 C     Upgrade to 2-loop RGE's for MSSM on 2/11/00
  12 #if defined(CERNLIB_IMPNONE)
  13       IMPLICIT NONE
  14 #endif
  15 #include "isajet/sslun.inc"
  16 #include "isajet/sugmg.inc"
  17 #include "isajet/sugpas.inc"
  18       REAL T,G(29),F(29)
  19       REAL FAC,COSB,TH2LP,B3,B2,XTAU,XT,XB,B1,THTOP,B12,B11,TANB,MT,
  20      $Q,B33,B32,PI,SINB,BETA,B21,B13,B22,B31,B23,XN,SMV,SPMV
  21       REAL A11,A12,A13,A21,A22,A23,A31,A32,A33,A14,A24
  22       REAL C11,C12,C13,C21,C22,C23,C31,C32,C33
  23       REAL D11,D12,D13,D21,D22,D23,D31,D32,D33
  24       REAL C1,C2,C3,CP1,CP2,CP3,CPP1,CPP2,CPP3,BSY1,BSY2,BSY3
  25       REAL BY4,BY5,BY6,SIG1,SIG2,SIG3
  26       REAL G72LP,G82LP,G92LP,FAC2LP,G102lP,G112lP,G122lP
  27       REAL G132LP,G142LP,G152LP,G162LP,G172LP,G182LP,G192LP,G262LP
  28       REAL G202LP,G212LP,G222LP,G232LP,G242LP,G252LP,G292LP,G282LP
  29       INTEGER THHL,THHH,THSHL,THSHH,THQ,THU,THD,THL,THE
  30       INTEGER TH1,TH2,TH3
  31       INTEGER NSL,NSD,NSE,NSW,NSH,NSU,NSQ,NH,NSG,NU,NE,ND,NN
  32 C
  33       DATA ND/3/,NE/3/,NN/3/
  34       DATA B11/7.96/,B12/5.4/,B13/17.6/,B21/1.8/,B22/25./,B23/24./
  35       DATA B31/2.2/,B32/9./,B33/14./
  36       DATA A11/5.2/,A12/2.8/,A13/3.6/,A21/6./,A22/6./,A23/2./
  37       DATA A31/4./,A32/4./,A33/0./,A14/1.2/,A24/2./
  38       DATA C11/1.7/,C12/.5/,C13/1.5/,C21/1.5/,C22/1.5/,C23/.5/
  39       DATA C31/2./,C32/2./,C33/0./
  40       DATA D11/3.98/,D12/2.7/,D13/8.8/,D21/.9/,D22/5.833/,D23/12./
  41       DATA D31/1.1/,D32/4.5/,D33/-26./
  42       DATA C1/.8666667/,C2/3./,C3/5.333333/
  43       DATA CP1/.4666667/,CP2/3./,CP3/5.333333/
  44       DATA CPP1/1.8/,CPP2/3./,CPP3/0./
  45       DATA BSY1/6.6/,BSY2/1./,BSY3/-3./
  46 C
  47 C-----THESE ARE VALID FROM MZ TO MGUT --------------------------------
  48       PI=4.*ATAN(1.)
  49       TANB=XTANB
  50       MT=AMT
  51       Q=MGUT*EXP(T)
  52       BETA=ATAN(TANB)
  53       SINB=SIN(BETA)
  54       COSB=SQRT(1.-SINB**2)
  55       FAC=2./16./PI**2
  56 C-----CALCULATE 1-LOOP THRESHOLD EFFECTS -----------------------------
  57       IF (Q.GT.MSUSY) THEN
  58         TH2LP=1.
  59       ELSE
  60         TH2LP=0.
  61       END IF
  62       IF (Q.GT.MSS(2)) THEN
  63         NSQ=3
  64         NSU=3
  65         NSD=3
  66       ELSE
  67         NSQ=0
  68         NSU=0
  69         NSD=0
  70       END IF
  71       IF (Q.GT.MSS(17)) THEN
  72         NSL=3
  73         NSE=3
  74       ELSE
  75         NSL=0
  76         NSE=0
  77       END IF
  78       IF (Q.GT.ABS(MU)) THEN
  79         NSH=2
  80       ELSE
  81         NSH=0
  82       END IF
  83       IF (Q.GT.ABS(MSS(27))) THEN
  84         NSW=1
  85       ELSE
  86         NSW=0
  87       END IF
  88       IF (Q.GT.MSS(1)) THEN
  89         NSG=1
  90       ELSE
  91         NSG=0
  92       END IF
  93       IF (Q.GT.MSS(31)) THEN
  94         NH=2
  95       ELSE
  96         NH=1
  97       END IF
  98       IF (Q.GT.GSS(7)) THEN
  99         TH1=1
 100       ELSE
 101         TH1=0
 102       END IF
 103       IF (Q.GT.GSS(8)) THEN
 104         TH2=1
 105       ELSE
 106         TH2=0
 107       END IF
 108       IF (Q.GT.MT) THEN
 109         NU=3
 110         THTOP=1.
 111       ELSE
 112         NU=2
 113         THTOP=0.
 114       END IF
 115       TH2LP=1.
 116       THTOP=1.
 117       THHL=1
 118       THHH=NH/2
 119       THSHL=NSH/2
 120       THSHH=NSH/2
 121       THQ=NSQ/3
 122       THU=NSQ/3
 123       THD=NSQ/3
 124       THL=NSL/3
 125       THE=NSE/3
 126       TH3=NSG
 127       B1=2.*(17*NU/12.+5*ND/12.+5*NE/4.+NN/4.)/5.+
 128      $ NSQ/30.+4*NSU/15.+NSD/15.+NSL/10.+NSE/5.+
 129      $ 1.*NSH/5.+1.*NH/10.
 130       B2=-22./3.+.5*(NU+ND)+1.*(NE+NN)/6.+
 131      $ 1.*NSQ/2.+1.*NSL/6.+1.*NSH/3.+1.*NH/6.+4.*NSW/3.
 132       B3=2.*(NU+ND)/3.+1.*NSQ/3.+1.*NSU/6.+1.*NSD/6.+2.*NSG-11.
 133       IF (Q.GT.MSUSY) THEN
 134       F(1)=G(1)/16./PI**2*(B1*G(1)**2+TH2LP/16./PI**2*G(1)**2*
 135      $(B11*G(1)**2+B12*G(2)**2+B13*G(3)**2-A11*G(6)**2-A12*G(5)**2
 136      $-A13*G(4)**2-A14*G(27)**2))
 137       F(2)=G(2)/16./PI**2*(B2*G(2)**2+TH2LP/16./PI**2*G(2)**2*
 138      $(B21*G(1)**2+B22*G(2)**2+B23*G(3)**2-A21*G(6)**2-A22*G(5)**2
 139      $-A23*G(4)**2-A24*G(27)**2))
 140       F(3)=G(3)/16./PI**2*(B3*G(3)**2+TH2LP/16./PI**2*G(3)**2*
 141      $(B31*G(1)**2+B32*G(2)**2+B33*G(3)**2-A31*G(6)**2-A32*G(5)**2
 142      $-A33*G(4)**2))
 143       ELSE
 144       F(1)=G(1)/16./PI**2*(B1*G(1)**2+TH2LP/16./PI**2*G(1)**2*
 145      $(D11*G(1)**2+D12*G(2)**2+D13*G(3)**2-C11*G(6)**2-C12*G(5)**2
 146      $-C13*G(4)**2))
 147       F(2)=G(2)/16./PI**2*(B2*G(2)**2+TH2LP/16./PI**2*G(2)**2*
 148      $(D21*G(1)**2+D22*G(2)**2+D23*G(3)**2-C21*G(6)**2-C22*G(5)**2
 149      $-C23*G(4)**2))
 150       F(3)=G(3)/16./PI**2*(B3*G(3)**2+TH2LP/16./PI**2*G(3)**2*
 151      $(D31*G(1)**2+D32*G(2)**2+D33*G(3)**2-C31*G(6)**2-C32*G(5)**2
 152      $-C33*G(4)**2))
 153       ENDIF
 154       IF (Q.LT.MSUSY) THEN
 155         F(4)=G(4)/16./PI**2*(5*G(4)**2*COSB**2/2.+3*G(6)**2*SINB**2*
 156      $   THTOP+3*G(5)**2*COSB**2-9*G(1)**2/4.-9*G(2)**2/4.
 157      $   -SINB**2*(3*G(6)**2*THTOP-3*G(5)**2-G(4)**2))
 158         F(5)=G(5)/16./PI**2*(9*G(5)**2*COSB**2/2.+3*G(6)**2*SINB**2*
 159      $   THTOP/2.+G(4)**2*COSB**2-G(1)**2/4.-9*G(2)**2/4.
 160      $   -8*G(3)**2-SINB**2*(3*G(6)**2*THTOP-3*G(5)**2-G(4)**2))
 161         F(6)=G(6)/16./PI**2*(9*G(6)**2*SINB**2/2.*THTOP+
 162      $   3*G(5)**2*COSB**2/2.+G(4)**2*COSB**2-17.*G(1)**2/20.
 163      $   -9*G(2)**2/4.-8*G(3)**2+COSB**2*
 164      $   (3*G(6)**2*THTOP-3*G(5)**2-G(4)**2))
 165       ELSE
 166 C        BY4=4*G(4)**2+3*G(5)**2-9*G(1)**2/5.-3*G(2)**2+G(27)**2
 167 C        BY5=6*G(5)**2+G(6)**2*THTOP+G(4)**2-7*G(1)**2/15.-3*G(2)**2-
 168 C     $16*G(3)**2/3.+
 169 C        BY6=6*G(6)**2*THTOP+G(5)**2-13*G(1)**2/15.-3*G(2)**2-
 170 C     $16*G(3)**2/3.+G(27)**2
 171         BY4=2.5*(COSB**2*THHL+SINB**2*THHH)*G(4)**2+
 172      $1.5*(COSB**2*THSHL+SINB**2*THSHH)*G(4)**2+
 173      $3*SINB**2*(THHL-THHH)*G(6)**2+3*(COSB**2*THHL+SINB**2*THHH)*
 174      $G(5)**2-
 175      $(3*G(1)**2/5.*(15./4.+3*THSHL/4.-TH1*(THL/4.+THE+THSHL/4.))
 176      $+G(2)**2*(9./4.+9*THSHL/4.-3*(THL+THSHL)*TH2/4.))
 177         BY4=BY4+G(27)**2
 178         BY5=4.5*(COSB**2*THHL+SINB**2*THHH)*G(5)**2+
 179      $1.5*(COSB**2*THSHL+SINB**2*THSHH)*G(5)**2+
 180      $.5*(SINB**2*THHL+COSB**2*THHH-
 181      $4*SINB**2*(THHL-THHH))*G(6)**2+.5*(SINB**2*THSHL+COSB**2*THSHH)*
 182      $G(6)**2+3*SINB**2*(THHL-THHH)*G(6)**2+
 183      $(COSB**2*THHL+SINB**2*THHH)*G(4)**2-
 184      $(3*G(1)**2/5.*(5./12.+3*THSHL/4.-TH1*(THQ/36.+THD/9.+THSHL/4.))
 185      $+G(2)**2*(9./4.+9*THSHL/4.-3*(THQ+THSHL)*TH2/4.)+
 186      $G(3)**2*(8.-4*(THQ+THD)*TH3/3.))
 187         BY6=4.5*(SINB**2*THHL+COSB**2*THHH)*G(6)**2+
 188      $1.5*(SINB**2*THSHL+COSB**2*THSHH)*G(6)**2+
 189      $.5*(COSB**2*THHL+SINB**2*THHH-
 190      $4*COSB**2*(THHL-THHH))*G(5)**2+.5*(COSB**2*THSHL+SINB**2*THSHH)*
 191      $G(5)**2+COSB**2*(THHL-THHH)*(3*G(5)**2+G(4)**2)-
 192      $(3*G(1)**2/5.*(17./12.+3*THSHL/4.-TH1*(THQ/36.+4*THU/9.+THSHL/4.))
 193      $+G(2)**2*(9./4.+9*THSHL/4.-3*(THQ+THSHL)*TH2/4.)+
 194      $G(3)**2*(8.-4*(THQ+THU)*TH3/3.))
 195         BY6=BY6+G(27)**2
 196         F(4)=G(4)/16./PI**2*(BY4+
 197      $       TH2LP/16./PI**2*((CPP1*BSY1+CPP1**2/2.)*G(1)**4+
 198      $       (CPP2*BSY2+CPP2**2/2.)*G(2)**4+(CPP3*BSY3+CPP3**2/2.)*
 199      $       G(3)**4+9*G(1)**2*G(2)**2/5.+
 200      $       G(5)**2*(-.4*G(1)**2+16*G(3)**2)+
 201      $       G(4)**2*(1.2*G(1)**2+6*G(2)**2)-
 202      $       (3*G(6)**2*G(5)**2+9*G(5)**4+9*G(5)**2*G(4)**2+
 203      $       10*G(4)**4)-3*G(27)**4-3*G(27)**2*G(6)**2
 204      $       -3*G(27)**2*G(4)**2))
 205         F(5)=G(5)/16./PI**2*(BY5+
 206      $       TH2LP/16./PI**2*((CP1*BSY1+CP1**2/2.)*G(1)**4+
 207      $       (CP2*BSY2+CP2**2/2.)*G(2)**4+(CP3*BSY3+CP3**2/2.)*
 208      $       G(3)**4+G(1)**2*G(2)**2+8*G(1)**2*G(3)**2/9.+
 209      $       8*G(2)**2*G(3)**2+.8*G(6)**2*G(1)**2+
 210      $       G(5)**2*(.4*G(1)**2+6*G(2)**2+16*G(3)**2)+
 211      $       1.2*G(4)**2*G(1)**2-
 212      $       (22*G(5)**4+5*G(6)**2*G(5)**2+3*G(5)**2*G(4)**2+
 213      $       3*G(4)**4+5*G(6)**4)-G(27)**2*G(6)**2-G(27)**2*G(4)**2))
 214         F(6)=G(6)/16./PI**2*(BY6+
 215      $       TH2LP/16./PI**2*((C1*BSY1+C1**2/2.)*G(1)**4+
 216      $       (C2*BSY2+C2**2/2.)*G(2)**4+(C3*BSY3+C3**2/2.)*
 217      $       G(3)**4+G(1)**2*G(2)**2+136*G(1)**2*G(3)**2/45.+
 218      $       8*G(2)**2*G(3)**2+G(6)**2*(1.2*G(1)**2+6*G(2)**2+
 219      $       16*G(3)**2)+.4*G(5)**2*G(1)**2-
 220      $       (22*G(6)**4+5*G(6)**2*G(5)**2+5*G(5)**4+
 221      $       G(5)**2*G(4)**2)-3*G(27)**4-3*G(27)**2*G(6)**2-
 222      $       G(27)**2*G(4)**2))
 223       END IF
 224       FAC2LP=(FAC/2.)**2
 225       G72LP=2*G(1)**2*(B11*G(1)**2*(G(7)+G(7))+
 226      $B12*G(2)**2*(G(7)+G(8))+B13*G(3)**2*(G(7)+G(9))+
 227      $A11*G(6)**2*(G(12)-G(7))+A12*G(5)**2*(G(11)-G(7))+
 228      $A13*G(4)**2*(G(10)-G(7))+A14*G(27)**2*(G(29)-G(7)))
 229       G82LP=2*G(2)**2*(B21*G(1)**2*(G(7)+G(8))+
 230      $B22*G(2)**2*(G(8)+G(8))+B23*G(3)**2*(G(8)+G(9))+
 231      $A21*G(6)**2*(G(12)-G(8))+A22*G(5)**2*(G(11)-G(8))+
 232      $A23*G(4)**2*(G(10)-G(8))+A24*G(27)**2*(G(29)-G(8)))
 233       G92LP=2*G(3)**2*(B31*G(1)**2*(G(9)+G(7))+
 234      $B32*G(2)**2*(G(9)+G(8))+B33*G(3)**2*(G(9)+G(9))+
 235      $A31*G(6)**2*(G(12)-G(9))+A32*G(5)**2*(G(11)-G(9))+
 236      $A33*G(4)**2*(G(10)-G(9)))
 237       F(7)=FAC*B1*G(1)**2*G(7)+TH2LP*FAC2LP*G72LP
 238       F(8)=FAC*B2*G(2)**2*G(8)+TH2LP*FAC2LP*G82LP
 239       F(9)=FAC*B3*G(3)**2*G(9)+TH2LP*FAC2LP*G92LP
 240       XTAU=G(21)+G(20)+G(13)+G(10)**2
 241       XB=G(24)+G(22)+G(13)+G(11)**2
 242       XT=G(24)+G(23)+G(14)+G(12)**2
 243       XN=G(21)+G(28)+G(14)+G(29)**2
 244       SMV=G(14)-G(13)+(2*G(19)+G(24))-(2*G(16)+G(21))
 245      $    -2*(2*G(18)+G(23))+(2*G(17)+G(22))+(2*G(15)+G(20))
 246       SPMV=-2*G(4)**2*G(20)+(6*G(1)**2*(2*G(15)+G(20)))/5.+
 247      $G(4)**2*(G(13)+G(21))-(3*(G(1)**2+5*G(2)**2)*
 248      $(G(13)-G(14)+2*G(16)+G(21)))/10.-2*G(5)**2*G(22)+
 249      $(2*(G(1)**2+20*G(3)**2)*(2*G(17)+G(22)))/15.+
 250      $4*G(6)**2*G(23)-(16*(G(1)**2+5*G(3)**2)*(2*G(18)+G(23)))/15.+
 251      $G(5)**2*(3*G(13)-G(24))-G(6)**2*(3*G(14)+G(24))+
 252      $((G(1)**2+45*G(2)**2+80*G(3)**2)*(2*G(19)+G(24)))/30.
 253      $+G(27)**2*(G(21)-G(14))
 254       SIG1=.2*G(1)**2*(3*(G(14)+G(13))+(2*G(19)+G(24))+
 255      $3*(2*G(16)+G(21))+8*(2*G(18)+G(23))+2*(2*G(17)+G(22))+
 256      $6*(2*G(15)+G(20)))
 257       SIG2=G(2)**2*(G(14)+G(13)+3*(2*G(19)+G(24))+(2*G(16)+G(21)))
 258       SIG3=G(3)**2*(2*(2*G(19)+G(24))+(2*G(18)+G(23))+
 259      $(2*G(17)+G(22)))
 260 C
 261       G102LP=(-54*G(1)**4-(18*G(1)**2*G(2)**2)/5.-
 262      $(12*G(1)**2*G(4)**2)/5.+(4*G(1)**2*G(5)**2)/5.)*G(7)+
 263      $((-18*G(1)**2*G(2)**2)/5.-30*G(2)**4-
 264      $12*G(2)**2*G(4)**2)*G(8)-32*G(3)**2*G(5)**2*G(9)+
 265      $(12*G(1)**2*G(4)**2*G(10))/5.+12*G(2)**2*G(4)**2*G(10)-
 266      $40*G(4)**4*G(10)-36*G(5)**4*G(11)+G(5)**2*
 267      $((-4*G(1)**2)/5.+32*G(3)**2-18*G(4)**2-6*G(6)**2)*G(11)+
 268      $G(5)**2*(-18*G(4)**2*G(10)-6*G(6)**2*G(12))
 269      $-6*G(27)**2*(G(12)*G(6)**2+G(10)*G(4)**2+G(29)*
 270      $(G(6)**2+2*G(27)**2+G(4)**2))
 271       G112LP=((-574*G(1)**4)/45.-2*G(1)**2*G(2)**2-
 272      $(16*G(1)**2*G(3)**2)/9.-(12*G(1)**2*G(4)**2)/5.-
 273      $(4*G(1)**2*G(5)**2)/5.-(8*G(1)**2*G(6)**2)/5.)*G(7)+
 274      $(-2*G(1)**2*G(2)**2-30*G(2)**4-16*G(2)**2*G(3)**2-
 275      $12*G(2)**2*G(5)**2)*G(8)+((-16*G(1)**2*G(3)**2)/9.-
 276      $16*G(2)**2*G(3)**2+(64*G(3)**4)/9.-
 277      $32*G(3)**2*G(5)**2)*G(9)+(12*G(1)**2*G(4)**2*G(10))/5.-
 278      $12*G(4)**4*G(10)-88*G(5)**4*G(11)+G(5)**2*((4*G(1)**2)/5.+
 279      $12*G(2)**2+32*G(3)**2-6*G(4)**2-10*G(6)**2)*G(11)+
 280      $(8*G(1)**2*G(6)**2*G(12))/5.-20*G(6)**4*G(12)+
 281      $G(5)**2*(-6*G(4)**2*G(10)-10*G(6)**2*G(12))
 282      $-2*G(27)**2*(G(12)*G(6)**2+G(10)*G(4)**2+G(29)*
 283      $(G(6)**2+G(4)**2))
 284       G122LP=((-5486*G(1)**4)/225.-2*G(1)**2*G(2)**2-
 285      $(272*G(1)**2*G(3)**2)/45.-(4*G(1)**2*G(5)**2)/5.-
 286      $(12*G(1)**2*G(6)**2)/5.)*G(7)+(-2*G(1)**2*G(2)**2-
 287      $30*G(2)**4-16*G(2)**2*G(3)**2-12*G(2)**2*G(6)**2)*G(8)+
 288      $((-272*G(1)**2*G(3)**2)/45.-16*G(2)**2*G(3)**2+
 289      $(64*G(3)**4)/9.-32*G(3)**2*G(6)**2)*G(9)-
 290      $2*G(4)**2*G(5)**2*G(10)+(4*G(1)**2*G(5)**2*G(11))/5.-
 291      $2*G(4)**2*G(5)**2*G(11)-20*G(5)**4*G(11)-
 292      $10*G(5)**2*G(6)**2*G(11)+((12*G(1)**2)/5.+12*G(2)**2+
 293      $32*G(3)**2-10*G(5)**2)*G(6)**2*G(12)-88*G(6)**4*G(12)
 294      $-2*G(27)**2*(3*G(12)*G(6)**2+G(10)*G(4)**2+G(29)*
 295      $(3*G(6)**2+6*G(27)**2+G(4)**2))
 296       G132LP=XB*((-4*G(1)**2)/5.+32*G(3)**2)*G(5)**2-
 297      $(12*G(1)**2*G(4)**2*(-XTAU-2*G(7)**2+2*G(7)*G(10)))/5.-
 298      $12*G(4)**4*(XTAU+G(10)**2)-36*G(5)**4*(XB+G(11)**2)+
 299      $G(5)**2*((8*G(1)**2*G(7)*(-G(7)+G(11)))/5.-
 300      $64*G(3)**2*G(9)*(-G(9)+G(11)))-6*G(5)**2*G(6)**2*
 301      $(XB+XT+2*G(11)*G(12))-1.2*G(1)**2*SPMV+33*G(2)**4*G(8)**2+
 302      $3.6*G(2)**2*G(1)**2*(G(8)**2+G(7)**2+G(7)*G(8))+
 303      $24.84*G(1)**4*G(7)**2+3*G(2)**2*SIG2+.6*G(1)**2*SIG1
 304      $-2*G(4)**2*G(27)**2*(2*G(29)*G(10)+XN+XTAU)
 305       G142LP=XT*((8*G(1)**2)/5.+32*G(3)**2)*G(6)**2-
 306      $6*G(5)**2*G(6)**2*(XB+XT+2*G(11)*G(12))-
 307      $36*G(6)**4*(XT+G(12)**2)+G(6)**2*((-16*G(1)**2*G(7)*
 308      $(-G(7)+G(12)))/5.-64*G(3)**2*G(9)*(-G(9)+G(12)))+
 309      $1.2*G(1)**2*SPMV+33*G(2)**4*G(8)**2+
 310      $3.6*G(2)**2*G(1)**2*(G(8)**2+G(7)**2+G(7)*G(8))+
 311      $24.84*G(1)**4*G(7)**2+3*G(2)**2*SIG2+.6*G(1)**2*SIG1
 312      $-12*G(27)**4*(G(29)**2+XN)-2*G(4)**2*G(27)**2*
 313      $(2*G(29)*G(10)+XN+XTAU)
 314       G152LP=2.4*G(1)**2*SPMV+112.32*G(1)**4*G(7)**2+2.4*G(1)**2*SIG1
 315       G162LP=-1.2*G(1)**2*SPMV+33*G(2)**4*G(8)**2+3.6*G(2)**2*G(1)**2*
 316      $(G(8)**2+G(7)**2+G(8)*G(7))+24.84*G(1)**4*G(7)**2+
 317      $3*G(2)**2*SIG2+.6*G(1)**2*SIG1
 318       G172LP=.8*G(1)**2*SPMV-42.66667*G(3)**4*G(9)**2+
 319      $128*G(3)**2*G(1)**2*(G(9)**2+G(7)**2+G(7)*G(9))/45.+
 320      $808*G(1)**4*G(7)**2/75.+16*G(3)**2*SIG3/3.+4*G(1)**2*SIG1/15.
 321       G182LP=-1.6*G(1)**2*SPMV-42.66667*G(3)**4*G(9)**2+
 322      $512*G(3)**2*G(1)**2*(G(9)**2+G(7)**2+G(7)*G(9))/45.+
 323      $3424*G(1)**4*G(7)**2/75.+16*G(3)**2*SIG3/3.+16*G(1)**2*SIG1/15.
 324       G192LP=.4*G(1)**2*SPMV-128.*G(3)**4*G(9)**2/3.+
 325      $32*G(3)**2*G(2)**2*(G(9)**2+G(8)**2+G(9)*G(8))+
 326      $32*G(3)**2*G(1)**2*(G(9)**2+G(7)**2+G(9)*G(7))/45.+
 327      $33*G(2)**4*G(8)**2+
 328      $.4*G(2)**2*G(1)**2*(G(8)**2+G(7)**2+G(8)*G(7))+
 329      $199*G(1)**4*G(7)**2/75.+16*G(3)**2*SIG3/3.+
 330      $3*G(2)**2*SIG2+G(1)**2*SIG1/15.
 331       G202LP=XTAU*((-12*G(1)**2)/5.+12*G(2)**2)*G(4)**2-
 332      $16*G(4)**4*(XTAU+G(10)**2)+G(4)**2*((24*G(1)**2*G(7)*
 333      $(-G(7)+G(10)))/5.-24*G(2)**2*G(8)*(-G(8)+G(10)))-
 334      $12*G(4)**2*G(5)**2*(XB+XTAU+2*G(10)*G(11))+
 335      $2.4*G(1)**2*SPMV+112.32*G(1)**4*G(7)**2+2.4*G(1)**2*SIG1
 336      $-4*G(27)**2*G(4)**2*(2*G(29)*G(10)+XN+XTAU)
 337       G212LP=(-12*G(1)**2*G(4)**2*(-XTAU-2*G(7)**2+2*G(7)*G(10)))/5.-
 338      $12*G(4)**4*(XTAU+G(10)**2)-6*G(4)**2*G(5)**2*
 339      $(XB+XTAU+2*G(10)*G(11))
 340      $-1.2*G(1)**2*SPMV+33*G(2)**4*G(8)**2+3.6*G(2)**2*G(1)**2*
 341      $(G(8)**2+G(7)**2+G(8)*G(7))+24.84*G(1)**4*G(7)**2+
 342      $3*G(2)**2*SIG2+.6*G(1)**2*SIG1
 343      $-6*G(27)**2*(2*G(12)*G(29)*G(6)**2+G(6)**2*XT+
 344      $G(6)**2*XN+2*G(27)**2*(XN+G(29)**2))
 345       G222LP=XB*((4*G(1)**2)/5.+12*G(2)**2)*G(5)**2-
 346      $4*G(4)**2*G(5)**2*(XB+XTAU+2*G(10)*G(11))-32*G(5)**4*
 347      $(XB+G(11)**2)+G(5)**2*((-8*G(1)**2*G(7)*(-G(7)+G(11)))/5.-
 348      $24*G(2)**2*G(8)*(-G(8)+G(11)))-4*G(5)**2*G(6)**2*(XB+XT+
 349      $2*G(11)*G(12))+.8*G(1)**2*SPMV-42.66667*G(3)**4*G(9)**2+
 350      $128*G(3)**2*G(1)**2*(G(9)**2+G(7)**2+G(7)*G(9))/45.+
 351      $808*G(1)**4*G(7)**2/75.+16*G(3)**2*SIG3/3.+4*G(1)**2*SIG1/15.
 352       G232LP=XT*((-4*G(1)**2)/5.+12*G(2)**2)*G(6)**2-
 353      $4*G(5)**2*G(6)**2*(XB+XT+2*G(11)*G(12))-32*G(6)**4*
 354      $(XT+G(12)**2)+G(6)**2*((8*G(1)**2*G(7)*(-G(7)+G(12)))/5.-
 355      $24*G(2)**2*G(8)*(-G(8)+G(12)))-1.6*G(1)**2*SPMV-42.66667*
 356      $G(3)**4*G(9)**2+
 357      $512*G(3)**2*G(1)**2*(G(9)**2+G(7)**2+G(7)*G(9))/45.+
 358      $3424*G(1)**4*G(7)**2/75.+16*G(3)**2*SIG3/3.+16*G(1)**2*SIG1/15.
 359      $-4*G(6)**2*G(27)**2*(2*G(12)*G(29)+XT+XN)
 360       G242LP=(-4*G(1)**2*G(5)**2*(-XB-2*G(7)**2+2*G(7)*G(11)))/5.-
 361      $2*G(4)**2*G(5)**2*(XB+XTAU+2*G(10)*G(11))-20*G(5)**4*
 362      $(XB+G(11)**2)-(8*G(1)**2*G(6)**2*(-XT-2*G(7)**2+
 363      $2*G(7)*G(12)))/5.-20*G(6)**4*(XT+G(12)**2)
 364      $+.4*G(1)**2*SPMV-128.*G(3)**4*G(9)**2/3.+
 365      $32*G(3)**2*G(2)**2*(G(9)**2+G(8)**2+G(9)*G(8))+
 366      $32*G(3)**2*G(1)**2*(G(9)**2+G(7)**2+G(9)*G(7))/45.+
 367      $33*G(2)**4*G(8)**2+
 368      $.4*G(2)**2*G(1)**2*(G(8)**2+G(7)**2+G(8)*G(7))+
 369      $199*G(1)**4*G(7)**2/75.+16*G(3)**2*SIG3/3.+
 370      $3*G(2)**2*SIG2+G(1)**2*SIG1/15.
 371      $-2*G(6)**2*G(27)**2*(2*G(12)*G(29)+XT+XN)
 372 C     ADD IN MU 2-LOOP TERM SOMETIME...
 373       G252LP=0.
 374       G282LP=-16*G(27)**4*(G(29)**2+XN)-12*G(6)**2*G(27)**2*
 375      $(2*G(12)*G(29)+XT+XN)-4*G(4)**2*G(27)**2*(2*G(29)*G(10)+
 376      $XN+XTAU)+G(27)**2*G(29)*(-24*G(7)*G(1)**2/5.-24*G(2)**2*G(8))
 377      $+G(27)**2*(2.4*G(1)**2*(2*G(7)**2+XN)+12*G(2)**2*
 378      $(2*G(8)**2+XN))
 379       G292LP=-32*G(3)**2*G(6)**2*G(9)+(12*G(1)**2*G(4)**2*G(10))/5.
 380      $-12*G(4)**4*G(10)-6*G(4)**2*G(5)**2*G(10)-
 381      $6*G(4)**2*G(5)**2*G(11)-6*G(5)**2*G(6)**2*G(11)+
 382      $(8*G(1)**2*G(6)**2*G(12))/5.+32*G(3)**2*G(6)**2*G(12)-
 383      $6*G(5)**2*G(6)**2*G(12)-36*G(6)**4*G(12)+
 384      $(-6*G(4)**2*G(10)-18*G(6)**2*G(12))*G(27)**2+
 385      $((12*G(1)**2)/5.+12*G(2)**2-6*G(4)**2-
 386      $18*G(6)**2)*G(29)*G(27)**2-40*G(29)*G(27)**4+
 387      $G(7)*((-414*G(1)**4)/25.-(18*G(1)**2*G(2)**2)/5.-
 388      $(12*G(1)**2*G(4)**2)/5.-(8*G(1)**2*G(6)**2)/5.-
 389      $(12*G(1)**2*G(27)**2)/5.)+G(8)*((-18*G(1)**2*G(2)**2)/5.-
 390      $30*G(2)**4-12*G(2)**2*G(27)**2)
 391 C
 392       F(10)=FAC*(9*G(1)**2*G(7)/5.+3*G(2)**2*G(8)+3*G(5)**2*G(11)+
 393      $4*G(4)**2*G(10)+G(27)**2*G(29))+FAC2LP*TH2LP*G102LP
 394       F(11)=FAC*(7*G(1)**2*G(7)/15.+3*G(2)**2*G(8)+16*G(3)**2*G(9)/3.+
 395      $6*G(5)**2*G(11)+G(6)**2*G(12)+G(4)**2*G(10))
 396      $+FAC2LP*TH2LP*G112LP
 397       F(12)=FAC*(13*G(1)**2*G(7)/15.+3*G(2)**2*G(8)+16*G(3)**2*G(9)/3.+
 398      $G(5)**2*G(11)+6*G(6)**2*G(12)+G(27)**2*G(29))
 399      $+FAC2LP*TH2LP*G122LP
 400       F(13)=FAC*(-.6*G(1)**2*G(7)**2-3*G(2)**2*G(8)**2+
 401      $3*G(5)**2*XB+G(4)**2*XTAU-.3*G(1)**2*SMV)+FAC2LP*TH2LP*G132LP
 402       F(14)=FAC*(-.6*G(1)**2*G(7)**2-3*G(2)**2*G(8)**2+3*G(6)**2*XT
 403      $+G(27)**2*XN+.3*G(1)**2*SMV)+FAC2LP*TH2LP*G142LP
 404       F(15)=FAC*(-2.4*G(1)**2*G(7)**2+.6*G(1)**2*SMV)
 405      $+FAC2LP*TH2LP*G152LP
 406       F(16)=FAC*(-.6*G(1)**2*G(7)**2-3*G(2)**2*G(8)**2-.3*G(1)**2*SMV)
 407      $+FAC2LP*TH2LP*G162LP
 408       F(17)=FAC*(-4*G(1)**2*G(7)**2/15.-16*G(3)**2*G(9)**2/3.+
 409      $.2*G(1)**2*SMV)+FAC2LP*TH2LP*G172LP
 410       F(18)=FAC*(-16*G(1)**2*G(7)**2/15.-16*G(3)**2*G(9)**2/3.-
 411      $.4*G(1)**2*SMV)+FAC2LP*TH2LP*G182LP
 412       F(19)=FAC*(-G(1)**2*G(7)**2/15.-3*G(2)**2*G(8)**2-
 413      $16*G(3)**2*G(9)**2/3.+.1*G(1)**2*SMV)+FAC2LP*TH2LP*G192LP
 414       F(20)=FAC*(-2.4*G(1)**2*G(7)**2+2*G(4)**2*XTAU+.6*G(1)**2*SMV)+
 415      $FAC2LP*TH2LP*G202LP
 416       F(21)=FAC*(-.6*G(1)**2*G(7)**2-3*G(2)**2*G(8)**2+G(4)**2*XTAU
 417      $+G(27)**2*XN-.3*G(1)**2*SMV)+FAC2LP*TH2LP*G212LP
 418       F(22)=FAC*(-4*G(1)**2*G(7)**2/15.-16*G(3)**2*G(9)**2/3.+
 419      $2*G(5)**2*XB+.2*G(1)**2*SMV)+FAC2LP*TH2LP*G222LP
 420       F(23)=FAC*(-16*G(1)**2*G(7)**2/15.-16*G(3)**2*G(9)**2/3.+
 421      $2*G(6)**2*XT-.4*G(1)**2*SMV)+FAC2LP*TH2LP*G232LP
 422       F(24)=FAC*(-G(1)**2*G(7)**2/15.-3*G(2)**2*G(8)**2-
 423      $16*G(3)**2*G(9)**2/3.+G(6)**2*XT+G(5)**2*XB+.1*G(1)**2*SMV)
 424      $+FAC2LP*TH2LP*G242LP
 425       F(25)=FAC*G(25)/2.*(-.6*G(1)**2-3*G(2)**2+3*G(6)**2+
 426      $       3*G(5)**2+G(4)**2+G(27)**2)+FAC2LP*TH2LP*G252LP
 427       F(26)=FAC*(-.6*G(1)**2*G(7)-3*G(2)**2*G(8)+3*G(6)**2*G(12)+
 428      $       3*G(5)**2*G(11)+G(4)**2*G(10))
 429       IF (Q.GT.AMNRMJ) THEN
 430       F(27)=G(27)/16./PI**2*(3*G(6)**2+G(4)**2+4*G(27)**2-
 431      $       3*G(2)**2-3*G(1)**2/5.+TH2LP/16./PI**2*(-10*G(27)**4
 432      $       +G(27)**2*(1.2*G(1)**2+3*(2*G(2)**2-3*G(6)**2-
 433      $       G(4)**2))+207*G(1)**4/50.+.2*G(1)**2*(9*G(2)**2+
 434      $       4*G(6)**2+6*G(4)**2)+.5*(15*G(2)**4-18*G(6)**4-
 435      $       6*G(4)**4+32*G(3)**2*G(6)**2-6*G(5)**2*G(6)**2-
 436      $       6*G(5)**2*G(4)**2)))
 437       F(28)=FAC*2*G(27)**2*XN+FAC2LP*TH2LP*G282LP
 438       F(29)=FAC*(3*G(1)**2*G(7)/5.+3*G(2)**2*G(8)+3*G(6)**2*G(12)+
 439      $4*G(27)**2*G(29)+G(4)**2*G(10))+FAC2LP*TH2LP*G292LP
 440       ELSE
 441       F(27)=0.
 442       F(28)=0.
 443       F(29)=0.
 444       END IF
 445       RETURN
 446       END
 447
 448
 449