[u/mrichter/AliRoot.git] / ISAJET / isasusy / surg06.F

#include "isajet/pilot.h"
C-----------------------------------------------------------------
      SUBROUTINE SURG06(T,GY,F)
C-----------------------------------------------------------------
C
C     Right hand side of truncated renormalization group equations
C          dG_i/dT = F_i(G)
C     using a single common scale MSUSY for SUSY mass thresholds.
C     Added GY(7) for neutrino Yukawa coupling on 9/24/99
C     THIS INCLUDES TWO-LOOP YUKAWAS FOR MSSM ONLY 11/18/99
C
#if defined(CERNLIB_IMPNONE)
      IMPLICIT NONE
#endif
#include "isajet/sslun.inc"
#include "isajet/sugpas.inc"
      REAL T,GY(7),F(7)
      REAL SINB,BETA,COSB,B1,THTOP,B3,B2,TH2LP,B12,B11,B21,B13,
     $MT,TANB,PI,Q,B33,B23,B22,B32,B31
      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
      INTEGER NSL,NSD,NSH,NSE,NSU,NSQ,NU,NSG,NSW,NH,NN,NE,ND
c      DATA MZ/91.17/
      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./
      TANB=XTANB
      MT=AMT
C-----THESE ARE VALID FROM MZ TO MGUT 
      Q=MGUT*EXP(T)
      PI=4.*ATAN(1.)
      BETA=ATAN(TANB)
      SINB=SIN(BETA)
      COSB=SQRT(1.-SINB**2)
      IF (Q.GT.MSUSY) THEN
        NSQ=3
        NSU=3
        NSD=3
        NSL=3
        NSE=3
        NSH=2
        NSW=1
        NSG=1
        NH=2
        TH2LP=1.
      ELSE
        NSQ=0
        NSU=0
        NSD=0
        NSL=0
        NSE=0
        NSH=0
        NSW=0
        NSG=0
        NH=1
        TH2LP=0.
      END IF
      IF (Q.GT.MT) THEN
        NU=3
        THTOP=1.
      ELSE
        NU=2
        THTOP=0.
      END IF
      TH2LP=1.
      THTOP=1.
      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)=GY(1)/16./PI**2*(B1*GY(1)**2+TH2LP/16./PI**2*GY(1)**2*
     $(B11*GY(1)**2+B12*GY(2)**2+B13*GY(3)**2-A11*GY(6)**2-A12*GY(5)**2
     $-A13*GY(4)**2))
      F(2)=GY(2)/16./PI**2*(B2*GY(2)**2+TH2LP/16./PI**2*GY(2)**2*
     $(B21*GY(1)**2+B22*GY(2)**2+B23*GY(3)**2-A21*GY(6)**2-A22*GY(5)**2
     $-A23*GY(4)**2))
      F(3)=GY(3)/16./PI**2*(B3*GY(3)**2+TH2LP/16./PI**2*GY(3)**2*
     $(B31*GY(1)**2+B32*GY(2)**2+B33*GY(3)**2-A31*GY(6)**2-A32*GY(5)**2
     $-A33*GY(4)**2))
      ELSE
      F(1)=GY(1)/16./PI**2*(B1*GY(1)**2+TH2LP/16./PI**2*GY(1)**2*
     $(D11*GY(1)**2+D12*GY(2)**2+D13*GY(3)**2-C11*GY(6)**2-C12*GY(5)**2
     $-C13*GY(4)**2))
      F(2)=GY(2)/16./PI**2*(B2*GY(2)**2+TH2LP/16./PI**2*GY(2)**2*
     $(D21*GY(1)**2+D22*GY(2)**2+D23*GY(3)**2-C21*GY(6)**2-C22*GY(5)**2
     $-C23*GY(4)**2))
      F(3)=GY(3)/16./PI**2*(B3*GY(3)**2+TH2LP/16./PI**2*GY(3)**2*
     $(D31*GY(1)**2+D32*GY(2)**2+D33*GY(3)**2-C31*GY(6)**2-C32*GY(5)**2
     $-C33*GY(4)**2))
      ENDIF
      IF (Q.LT.MSUSY) THEN
        F(4)=GY(4)/16./PI**2*(5*GY(4)**2*COSB**2/2.+3*GY(6)**2*SINB**2*
     $   THTOP+3*GY(5)**2*COSB**2-9*GY(1)**2/4.-9*GY(2)**2/4.
     $   -SINB**2*(3*GY(6)**2*THTOP-3*GY(5)**2-GY(4)**2))
        F(5)=GY(5)/16./PI**2*(9*GY(5)**2*COSB**2/2.+3*GY(6)**2*SINB**2*
     $   THTOP/2.+GY(4)**2*COSB**2-GY(1)**2/4.-9*GY(2)**2/4.
     $   -8*GY(3)**2-SINB**2*(3*GY(6)**2*THTOP-3*GY(5)**2-GY(4)**2))
        F(6)=GY(6)/16./PI**2*(9*GY(6)**2*SINB**2/2.*THTOP+
     $   3*GY(5)**2*COSB**2/2.+GY(4)**2*COSB**2-17.*GY(1)**2/20.
     $   -9*GY(2)**2/4.-8*GY(3)**2+COSB**2*
     $   (3*GY(6)**2*THTOP-3*GY(5)**2-GY(4)**2))
      ELSE
        F(4)=GY(4)/16./PI**2*(4*GY(4)**2+3*GY(5)**2+GY(7)**2-
     $        9*GY(1)**2/5.-3*GY(2)**2+
     $       TH2LP/16./PI**2*((CPP1*BSY1+CPP1**2/2.)*GY(1)**4+
     $       (CPP2*BSY2+CPP2**2/2.)*GY(2)**4+(CPP3*BSY3+CPP3**2/2.)*
     $       GY(3)**4+9*GY(1)**2*GY(2)**2/5.+
     $       GY(5)**2*(-.4*GY(1)**2+16*GY(3)**2)+
     $       GY(4)**2*(1.2*GY(1)**2+6*GY(2)**2)-
     $       (3*GY(6)**2*GY(5)**2+9*GY(5)**4+9*GY(5)**2*GY(4)**2+
     $       10*GY(4)**4)))
        F(5)=GY(5)/16./PI**2*(6*GY(5)**2+GY(6)**2*THTOP+
     $        GY(4)**2-7*GY(1)**2/15.-3*GY(2)**2-16*GY(3)**2/3.+
     $       TH2LP/16./PI**2*((CP1*BSY1+CP1**2/2.)*GY(1)**4+
     $       (CP2*BSY2+CP2**2/2.)*GY(2)**4+(CP3*BSY3+CP3**2/2.)*
     $       GY(3)**4+GY(1)**2*GY(2)**2+8*GY(1)**2*GY(3)**2/9.+
     $       8*GY(2)**2*GY(3)**2+.8*GY(6)**2*GY(1)**2+
     $       GY(5)**2*(.4*GY(1)**2+6*GY(2)**2+16*GY(3)**2)+
     $       1.2*GY(4)**2*GY(1)**2-
     $       (22*GY(5)**4+5*GY(6)**2*GY(5)**2+3*GY(5)**2*GY(4)**2+
     $       3*GY(4)**4+5*GY(6)**4)))
        F(6)=GY(6)/16./PI**2*(6*GY(6)**2*THTOP+GY(5)**2+GY(7)**2-
     $        13*GY(1)**2/15.-3*GY(2)**2-16*GY(3)**2/3.+
     $       TH2LP/16./PI**2*((C1*BSY1+C1**2/2.)*GY(1)**4+
     $       (C2*BSY2+C2**2/2.)*GY(2)**4+(C3*BSY3+C3**2/2.)*
     $       GY(3)**4+GY(1)**2*GY(2)**2+136*GY(1)**2*GY(3)**2/45.+
     $       8*GY(2)**2*GY(3)**2+GY(6)**2*(1.2*GY(1)**2+6*GY(2)**2+
     $       16*GY(3)**2)+.4*GY(5)**2*GY(1)**2-
     $       (22*GY(6)**4+5*GY(6)**2*GY(5)**2+5*GY(5)**4+
     $       GY(5)**2*GY(4)**2)))
      END IF
C     THE NEUTRINO YUKAWA IS TO 1-LOOP ONLY FOR NOW...
      IF (Q.GT.AMNRMJ) THEN
        F(7)=GY(7)/16./PI**2*(3*GY(6)**2+GY(4)**2+4*GY(7)**2-
     $        3*GY(2)**2-3*GY(1)**2/5.)
      ELSE
        F(7)=0.
      END IF
      RETURN
      END
Commit	Line	Data
0795afa3	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.(17NU/12.+5ND/12.+5NE/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(B1GY(1)2+TH2LP/16./PI2GY(1)*2
87	$(B11GY(1)2+B12GY(2)*2+B13GY(3)*2-A11GY(6)*2-A12GY(5)**2
88	$-A13GY(4)*2))
89	F(2)=GY(2)/16./PI*2(B2GY(2)2+TH2LP/16./PI2GY(2)*2
90	$(B21GY(1)2+B22GY(2)*2+B23GY(3)*2-A21GY(6)*2-A22GY(5)**2
91	$-A23GY(4)*2))
92	F(3)=GY(3)/16./PI*2(B3GY(3)2+TH2LP/16./PI2GY(3)*2
93	$(B31GY(1)2+B32GY(2)*2+B33GY(3)*2-A31GY(6)*2-A32GY(5)**2
94	$-A33GY(4)*2))
95	ELSE
96	F(1)=GY(1)/16./PI*2(B1GY(1)2+TH2LP/16./PI2GY(1)*2
97	$(D11GY(1)2+D12GY(2)*2+D13GY(3)*2-C11GY(6)*2-C12GY(5)**2
98	$-C13GY(4)*2))
99	F(2)=GY(2)/16./PI*2(B2GY(2)2+TH2LP/16./PI2GY(2)*2
100	$(D21GY(1)2+D22GY(2)*2+D23GY(3)*2-C21GY(6)*2-C22GY(5)**2
101	$-C23GY(4)*2))
102	F(3)=GY(3)/16./PI*2(B3GY(3)2+TH2LP/16./PI2GY(3)*2
103	$(D31GY(1)2+D32GY(2)*2+D33GY(3)*2-C31GY(6)*2-C32GY(5)**2
104	$-C33GY(4)*2))
105	ENDIF
106	IF (Q.LT.MSUSY) THEN
107	F(4)=GY(4)/16./PI*2(5GY(4)2COSB*2/2.+3GY(6)*2SINB*2
108	$ THTOP+3GY(5)2COSB*2-9GY(1)*2/4.-9GY(2)**2/4.
109	$ -SINB*2(3GY(6)2THTOP-3GY(5)2-GY(4)*2))
110	F(5)=GY(5)/16./PI*2(9GY(5)2COSB*2/2.+3GY(6)*2SINB*2
111	$ THTOP/2.+GY(4)*2COSB2-GY(1)2/4.-9GY(2)*2/4.
112	$ -8GY(3)2-SINB2(3GY(6)2THTOP-3GY(5)2-GY(4)*2))
113	F(6)=GY(6)/16./PI*2(9GY(6)2SINB*2/2.THTOP+
114	$ 3GY(5)2COSB2/2.+GY(4)2COSB2-17.GY(1)**2/20.
115	$ -9GY(2)2/4.-8GY(3)2+COSB2*
116	$ (3GY(6)2THTOP-3GY(5)2-GY(4)*2))
117	ELSE
118	F(4)=GY(4)/16./PI*2(4GY(4)2+3GY(5)2+GY(7)2-
119	$ 9GY(1)2/5.-3GY(2)**2+
120	$ TH2LP/16./PI*2((CPP1BSY1+CPP12/2.)GY(1)**4+
121	$ (CPP2BSY2+CPP22/2.)GY(2)*4+(CPP3BSY3+CPP3*2/2.)
122	$ GY(3)*4+9GY(1)*2GY(2)**2/5.+
123	$ GY(5)*2(-.4GY(1)2+16GY(3)**2)+
124	$ GY(4)*2(1.2GY(1)2+6GY(2)**2)-
125	$ (3GY(6)2GY(5)*2+9GY(5)*4+9GY(5)*2GY(4)**2+
126	$ 10GY(4)*4)))
127	F(5)=GY(5)/16./PI*2(6GY(5)2+GY(6)2THTOP+
128	$ GY(4)*2-7GY(1)*2/15.-3GY(2)*2-16GY(3)**2/3.+
129	$ TH2LP/16./PI*2((CP1BSY1+CP12/2.)GY(1)**4+
130	$ (CP2BSY2+CP22/2.)GY(2)*4+(CP3BSY3+CP3*2/2.)
131	$ GY(3)4+GY(1)2GY(2)2+8GY(1)*2GY(3)**2/9.+
132	$ 8GY(2)2GY(3)*2+.8GY(6)*2GY(1)**2+
133	$ GY(5)*2(.4GY(1)2+6GY(2)*2+16GY(3)**2)+
134	$ 1.2GY(4)2GY(1)**2-
135	$ (22GY(5)4+5GY(6)*2GY(5)*2+3GY(5)*2GY(4)**2+
136	$ 3GY(4)4+5GY(6)**4)))
137	F(6)=GY(6)/16./PI*2(6GY(6)2THTOP+GY(5)2+GY(7)2-
138	$ 13GY(1)2/15.-3GY(2)*2-16GY(3)**2/3.+
139	$ TH2LP/16./PI*2((C1BSY1+C12/2.)GY(1)**4+
140	$ (C2BSY2+C22/2.)GY(2)*4+(C3BSY3+C3*2/2.)
141	$ GY(3)4+GY(1)2GY(2)2+136GY(1)*2GY(3)**2/45.+
142	$ 8GY(2)2GY(3)2+GY(6)2(1.2GY(1)*2+6GY(2)**2+
143	$ 16GY(3)2)+.4GY(5)*2GY(1)**2-
144	$ (22GY(6)4+5GY(6)*2GY(5)*2+5GY(5)**4+
145	$ GY(5)*2GY(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(3GY(6)2+GY(4)2+4GY(7)**2-
150	$ 3GY(2)2-3GY(1)**2/5.)
151	ELSE
152	F(7)=0.
153	END IF
154	RETURN
155	END