1 #include "isajet/pilot.h"
2 C-----------------------------------------------------------------
3 SUBROUTINE SURG26(T,G,F)
4 C-----------------------------------------------------------------
6 C Right hand side of full renormalization group equations
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)
15 #include "isajet/sslun.inc"
16 #include "isajet/sugmg.inc"
17 #include "isajet/sugpas.inc"
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
31 INTEGER NSL,NSD,NSE,NSW,NSH,NSU,NSQ,NH,NSG,NU,NE,ND,NN
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./
47 C-----THESE ARE VALID FROM MZ TO MGUT --------------------------------
56 C-----CALCULATE 1-LOOP THRESHOLD EFFECTS -----------------------------
71 IF (Q.GT.MSS(17)) THEN
78 IF (Q.GT.ABS(MU)) THEN
83 IF (Q.GT.ABS(MSS(27))) THEN
93 IF (Q.GT.MSS(31)) THEN
103 IF (Q.GT.GSS(8)) THEN
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.
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
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
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
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
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))
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-
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)*
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.))
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.))
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-
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))+
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))+
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)*
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.+
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*
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.+
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...
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*
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)
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)+
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