#include "isajet/pilot.h" C----------------------------------------------------------------- SUBROUTINE SURG26(T,G,F) C----------------------------------------------------------------- C C Right hand side of full renormalization group equations C dG_i/dT = F_i(G) C using the thresholds MSS for each mass calculated with the C couplings G0 frozen by SUGFRZ. C Added right neutrino RGE's on 9/24/99 #if defined(CERNLIB_IMPNONE) IMPLICIT NONE #endif #include "isajet/sslun.inc" #include "isajet/sugmg.inc" #include "isajet/sugpas.inc" REAL T,G(29),F(29) REAL FAC,COSB,TH2LP,B3,B2,XTAU,XT,XB,B1,THTOP,B12,B11,TANB,MT, $Q,B33,B32,PI,SINB,BETA,B21,B13,B22,B31,B23,XN REAL A11,A12,A13,A21,A22,A23,A31,A32,A33 REAL C11,C12,C13,C21,C22,C23,C31,C32,C33 REAL D11,D12,D13,D21,D22,D23,D31,D32,D33 REAL C1,C2,C3,CP1,CP2,CP3,CPP1,CPP2,CPP3,BSY1,BSY2,BSY3 REAL BY4,BY5,BY6 INTEGER THHL,THHH,THSHL,THSHH,THQ,THU,THD,THL,THE INTEGER TH1,TH2,TH3 INTEGER NSL,NSD,NSE,NSW,NSH,NSU,NSQ,NH,NSG,NU,NE,ND,NN C DATA ND/3/,NE/3/,NN/3/ DATA B11/7.96/,B12/5.4/,B13/17.6/,B21/1.8/,B22/25./,B23/24./ DATA B31/2.2/,B32/9./,B33/14./ DATA A11/5.2/,A12/2.8/,A13/3.6/,A21/6./,A22/6./,A23/2./ DATA A31/4./,A32/4./,A33/0./ DATA C11/1.7/,C12/.5/,C13/1.5/,C21/1.5/,C22/1.5/,C23/.5/ DATA C31/2./,C32/2./,C33/0./ DATA D11/3.98/,D12/2.7/,D13/8.8/,D21/.9/,D22/5.833/,D23/12./ DATA D31/1.1/,D32/4.5/,D33/-26./ DATA C1/.8666667/,C2/3./,C3/5.333333/ DATA CP1/.4666667/,CP2/3./,CP3/5.333333/ DATA CPP1/1.8/,CPP2/3./,CPP3/0./ DATA BSY1/6.6/,BSY2/1./,BSY3/-3./ C C-----THESE ARE VALID FROM MZ TO MGUT -------------------------------- PI=4.*ATAN(1.) TANB=XTANB MT=AMT Q=MGUT*EXP(T) BETA=ATAN(TANB) SINB=SIN(BETA) COSB=SQRT(1.-SINB**2) FAC=2./16./PI**2 C-----CALCULATE 1-LOOP THRESHOLD EFFECTS ----------------------------- IF (Q.GT.MSUSY) THEN TH2LP=1. ELSE TH2LP=0. END IF IF (Q.GT.MSS(2)) THEN NSQ=3 NSU=3 NSD=3 ELSE NSQ=0 NSU=0 NSD=0 END IF IF (Q.GT.MSS(17)) THEN NSL=3 NSE=3 ELSE NSL=0 NSE=0 END IF IF (Q.GT.ABS(MU)) THEN NSH=2 ELSE NSH=0 END IF IF (Q.GT.ABS(MSS(27))) THEN NSW=1 ELSE NSW=0 END IF IF (Q.GT.MSS(1)) THEN NSG=1 ELSE NSG=0 END IF IF (Q.GT.MSS(31)) THEN NH=2 ELSE NH=1 END IF IF (Q.GT.GSS(7)) THEN TH1=1 ELSE TH1=0 END IF IF (Q.GT.GSS(8)) THEN TH2=1 ELSE TH2=0 END IF IF (Q.GT.MT) THEN NU=3 THTOP=1. ELSE NU=2 THTOP=0. END IF TH2LP=1. THTOP=1. THHL=1 THHH=NH/2 THSHL=NSH/2 THSHH=NSH/2 THQ=NSQ/3 THU=NSQ/3 THD=NSQ/3 THL=NSL/3 THE=NSE/3 TH3=NSG B1=2.*(17*NU/12.+5*ND/12.+5*NE/4.+NN/4.)/5.+ $ NSQ/30.+4*NSU/15.+NSD/15.+NSL/10.+NSE/5.+ $ 1.*NSH/5.+1.*NH/10. B2=-22./3.+.5*(NU+ND)+1.*(NE+NN)/6.+ $ 1.*NSQ/2.+1.*NSL/6.+1.*NSH/3.+1.*NH/6.+4.*NSW/3. B3=2.*(NU+ND)/3.+1.*NSQ/3.+1.*NSU/6.+1.*NSD/6.+2.*NSG-11. IF (Q.GT.MSUSY) THEN F(1)=G(1)/16./PI**2*(B1*G(1)**2+TH2LP/16./PI**2*G(1)**2* $(B11*G(1)**2+B12*G(2)**2+B13*G(3)**2-A11*G(6)**2-A12*G(5)**2 $-A13*G(4)**2)) F(2)=G(2)/16./PI**2*(B2*G(2)**2+TH2LP/16./PI**2*G(2)**2* $(B21*G(1)**2+B22*G(2)**2+B23*G(3)**2-A21*G(6)**2-A22*G(5)**2 $-A23*G(4)**2)) F(3)=G(3)/16./PI**2*(B3*G(3)**2+TH2LP/16./PI**2*G(3)**2* $(B31*G(1)**2+B32*G(2)**2+B33*G(3)**2-A31*G(6)**2-A32*G(5)**2 $-A33*G(4)**2)) ELSE F(1)=G(1)/16./PI**2*(B1*G(1)**2+TH2LP/16./PI**2*G(1)**2* $(D11*G(1)**2+D12*G(2)**2+D13*G(3)**2-C11*G(6)**2-C12*G(5)**2 $-C13*G(4)**2)) F(2)=G(2)/16./PI**2*(B2*G(2)**2+TH2LP/16./PI**2*G(2)**2* $(D21*G(1)**2+D22*G(2)**2+D23*G(3)**2-C21*G(6)**2-C22*G(5)**2 $-C23*G(4)**2)) F(3)=G(3)/16./PI**2*(B3*G(3)**2+TH2LP/16./PI**2*G(3)**2* $(D31*G(1)**2+D32*G(2)**2+D33*G(3)**2-C31*G(6)**2-C32*G(5)**2 $-C33*G(4)**2)) ENDIF IF (Q.LT.MSUSY) THEN F(4)=G(4)/16./PI**2*(5*G(4)**2*COSB**2/2.+3*G(6)**2*SINB**2* $ THTOP+3*G(5)**2*COSB**2-9*G(1)**2/4.-9*G(2)**2/4. $ -SINB**2*(3*G(6)**2*THTOP-3*G(5)**2-G(4)**2)) F(5)=G(5)/16./PI**2*(9*G(5)**2*COSB**2/2.+3*G(6)**2*SINB**2* $ THTOP/2.+G(4)**2*COSB**2-G(1)**2/4.-9*G(2)**2/4. $ -8*G(3)**2-SINB**2*(3*G(6)**2*THTOP-3*G(5)**2-G(4)**2)) F(6)=G(6)/16./PI**2*(9*G(6)**2*SINB**2/2.*THTOP+ $ 3*G(5)**2*COSB**2/2.+G(4)**2*COSB**2-17.*G(1)**2/20. $ -9*G(2)**2/4.-8*G(3)**2+COSB**2* $ (3*G(6)**2*THTOP-3*G(5)**2-G(4)**2)) ELSE C BY4=4*G(4)**2+3*G(5)**2-9*G(1)**2/5.-3*G(2)**2+G(27)**2 C BY5=6*G(5)**2+G(6)**2*THTOP+G(4)**2-7*G(1)**2/15.-3*G(2)**2- C $16*G(3)**2/3.+ C BY6=6*G(6)**2*THTOP+G(5)**2-13*G(1)**2/15.-3*G(2)**2- C $16*G(3)**2/3.+G(27)**2 BY4=2.5*(COSB**2*THHL+SINB**2*THHH)*G(4)**2+ $1.5*(COSB**2*THSHL+SINB**2*THSHH)*G(4)**2+ $3*SINB**2*(THHL-THHH)*G(6)**2+3*(COSB**2*THHL+SINB**2*THHH)* $G(5)**2- $(3*G(1)**2/5.*(15./4.+3*THSHL/4.-TH1*(THL/4.+THE+THSHL/4.)) $+G(2)**2*(9./4.+9*THSHL/4.-3*(THL+THSHL)*TH2/4.)) BY4=BY4+G(27)**2 BY5=4.5*(COSB**2*THHL+SINB**2*THHH)*G(5)**2+ $1.5*(COSB**2*THSHL+SINB**2*THSHH)*G(5)**2+ $.5*(SINB**2*THHL+COSB**2*THHH- $4*SINB**2*(THHL-THHH))*G(6)**2+.5*(SINB**2*THSHL+COSB**2*THSHH)* $G(6)**2+3*SINB**2*(THHL-THHH)*G(6)**2+ $(COSB**2*THHL+SINB**2*THHH)*G(4)**2- $(3*G(1)**2/5.*(5./12.+3*THSHL/4.-TH1*(THQ/36.+THD/9.+THSHL/4.)) $+G(2)**2*(9./4.+9*THSHL/4.-3*(THQ+THSHL)*TH2/4.)+ $G(3)**2*(8.-4*(THQ+THD)*TH3/3.)) BY6=4.5*(SINB**2*THHL+COSB**2*THHH)*G(6)**2+ $1.5*(SINB**2*THSHL+COSB**2*THSHH)*G(6)**2+ $.5*(COSB**2*THHL+SINB**2*THHH- $4*COSB**2*(THHL-THHH))*G(5)**2+.5*(COSB**2*THSHL+SINB**2*THSHH)* $G(5)**2+COSB**2*(THHL-THHH)*(3*G(5)**2+G(4)**2)- $(3*G(1)**2/5.*(17./12.+3*THSHL/4.-TH1*(THQ/36.+4*THU/9.+THSHL/4.)) $+G(2)**2*(9./4.+9*THSHL/4.-3*(THQ+THSHL)*TH2/4.)+ $G(3)**2*(8.-4*(THQ+THU)*TH3/3.)) BY6=BY6+G(27)**2 F(4)=G(4)/16./PI**2*(BY4+ $ TH2LP/16./PI**2*((CPP1*BSY1+CPP1**2/2.)*G(1)**4+ $ (CPP2*BSY2+CPP2**2/2.)*G(2)**4+(CPP3*BSY3+CPP3**2/2.)* $ G(3)**4+9*G(1)**2*G(2)**2/5.+ $ G(5)**2*(-.4*G(1)**2+16*G(3)**2)+ $ G(4)**2*(1.2*G(1)**2+6*G(2)**2)- $ (3*G(6)**2*G(5)**2+9*G(5)**4+9*G(5)**2*G(4)**2+ $ 10*G(4)**4))) F(5)=G(5)/16./PI**2*(BY5+ $ TH2LP/16./PI**2*((CP1*BSY1+CP1**2/2.)*G(1)**4+ $ (CP2*BSY2+CP2**2/2.)*G(2)**4+(CP3*BSY3+CP3**2/2.)* $ G(3)**4+G(1)**2*G(2)**2+8*G(1)**2*G(3)**2/9.+ $ 8*G(2)**2*G(3)**2+.8*G(6)**2*G(1)**2+ $ G(5)**2*(.4*G(1)**2+6*G(2)**2+16*G(3)**2)+ $ 1.2*G(4)**2*G(1)**2- $ (22*G(5)**4+5*G(6)**2*G(5)**2+3*G(5)**2*G(4)**2+ $ 3*G(4)**4+5*G(6)**4))) F(6)=G(6)/16./PI**2*(BY6+ $ TH2LP/16./PI**2*((C1*BSY1+C1**2/2.)*G(1)**4+ $ (C2*BSY2+C2**2/2.)*G(2)**4+(C3*BSY3+C3**2/2.)* $ G(3)**4+G(1)**2*G(2)**2+136*G(1)**2*G(3)**2/45.+ $ 8*G(2)**2*G(3)**2+G(6)**2*(1.2*G(1)**2+6*G(2)**2+ $ 16*G(3)**2)+.4*G(5)**2*G(1)**2- $ (22*G(6)**4+5*G(6)**2*G(5)**2+5*G(5)**4+ $ G(5)**2*G(4)**2))) END IF F(7)=FAC*B1*G(1)**2*G(7) F(8)=FAC*B2*G(2)**2*G(8) F(9)=FAC*B3*G(3)**2*G(9) F(10)=FAC*(9*G(1)**2*G(7)/5.+3*G(2)**2*G(8)+3*G(5)**2*G(11)+ $ 4*G(4)**2*G(10)+G(27)**2*G(29)) F(11)=FAC*(7*G(1)**2*G(7)/15.+3*G(2)**2*G(8)+16*G(3)**2*G(9)/3.+ $ 6*G(5)**2*G(11)+G(6)**2*G(12)+G(4)**2*G(10)) F(12)=FAC*(13*G(1)**2*G(7)/15.+3*G(2)**2*G(8)+16*G(3)**2*G(9)/3.+ $ G(5)**2*G(11)+6*G(6)**2*G(12)+G(27)**2*G(29)) XTAU=G(21)+G(20)+G(13)+G(10)**2 XB=G(24)+G(22)+G(13)+G(11)**2 XT=G(24)+G(23)+G(14)+G(12)**2 XN=G(21)+G(28)+G(14)+G(29)**2 F(13)=FAC*(-.6*G(1)**2*G(7)**2-3*G(2)**2*G(8)**2+ $ 3*G(5)**2*XB+G(4)**2*XTAU) F(14)=FAC*(-.6*G(1)**2*G(7)**2-3*G(2)**2*G(8)**2+3*G(6)**2*XT $+G(27)**2*XN) F(15)=-2.4*FAC*G(1)**2*G(7)**2 F(16)=FAC*(-.6*G(1)**2*G(7)**2-3*G(2)**2*G(8)**2) F(17)=FAC*(-4*G(1)**2*G(7)**2/15.-16*G(3)**2*G(9)**2/3.) F(18)=FAC*(-16*G(1)**2*G(7)**2/15.-16*G(3)**2*G(9)**2/3.) F(19)=FAC*(-G(1)**2*G(7)**2/15.-3*G(2)**2*G(8)**2- $ 16*G(3)**2*G(9)**2/3.) F(20)=FAC*(-2.4*G(1)**2*G(7)**2+2*G(4)**2*XTAU) F(21)=FAC*(-.6*G(1)**2*G(7)**2-3*G(2)**2*G(8)**2+G(4)**2*XTAU $+G(27)**2*XN) F(22)=FAC*(-4*G(1)**2*G(7)**2/15.-16*G(3)**2*G(9)**2/3.+ $ 2*G(5)**2*XB) F(23)=FAC*(-16*G(1)**2*G(7)**2/15.-16*G(3)**2*G(9)**2/3.+ $ 2*G(6)**2*XT) F(24)=FAC*(-G(1)**2*G(7)**2/15.-3*G(2)**2*G(8)**2- $ 16*G(3)**2*G(9)**2/3.+G(6)**2*XT+G(5)**2*XB) F(25)=FAC*G(25)/2.*(-.6*G(1)**2-3*G(2)**2+3*G(6)**2+ $ 3*G(5)**2+G(4)**2) F(26)=FAC*(.6*G(1)**2*G(7)+3*G(2)**2*G(8)+3*G(6)**2*G(12)+ $ 3*G(5)**2*G(11)+G(4)**2*G(10)) IF (Q.GT.AMNRMJ) THEN F(27)=G(27)/16./PI**2*(3*G(6)**2+G(4)**2+4*G(27)**2- $ 3*G(2)**2-3*G(1)**2/5.) F(28)=FAC*2*G(27)**2*XN F(29)=FAC*(3*G(1)**2*G(7)/5.+3*G(2)**2*G(8)+3*G(6)**2*G(12)+ $4*G(27)**2*G(29)+G(4)**2*G(10)) ELSE F(27)=0. F(28)=0. F(29)=0. END IF RETURN END