]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | C----------------------------------------------------------------- | |
3 | SUBROUTINE SUGEFF(G0,SIG1,SIG2) | |
4 | C----------------------------------------------------------------- | |
5 | C | |
6 | C Compute Higgs mass shift due to 1-loop effective potential | |
7 | C | |
8 | #if defined(CERNLIB_IMPNONE) | |
9 | IMPLICIT NONE | |
10 | #endif | |
11 | #include "isajet/sslun.inc" | |
12 | #include "isajet/sugpas.inc" | |
13 | #include "isajet/sugmg.inc" | |
14 | REAL G0(29),SIG1,SIG2 | |
15 | REAL DT1,DELT1S,SIG1T1,DMSDV2,FB,FT,MST2,MSB2,MSB1,SIG2B1, | |
16 | $SIG1B1,SIG2B2,SIG1B2,DB1,SIG1T2,SIG2T1,DELB1S,SIG2T2,MST1,COS2W, | |
17 | $MT,PI,COTB,TANB,MB,SIG2B,SIG1B,G,GGP,SIG2T,E,FAC,SIG1T,QS, | |
18 | $BETA,SINB,COSB,FAC4,FL,ML,SIG1L,SIG2L,MSL1,MSL2, | |
19 | $DELL1S,DL1,SIG1L1,SIG2L1,SIG1L2,SIG2L2 | |
20 | C | |
21 | G=G2 | |
22 | TANB=XTANB | |
23 | COS2W=1.-XW | |
24 | COTB=1./TANB | |
25 | BETA=ATAN(TANB) | |
26 | SINB=SIN(BETA) | |
27 | COSB=COS(BETA) | |
28 | PI=4.*ATAN(1.) | |
29 | FAC=3./8./PI**2 | |
30 | FAC4=FAC/3. | |
31 | E=EXP(1.) | |
32 | QS=HIGFRZ**2 | |
33 | C-----CALCULATE TOP AND BOTTOM CONTRIBUTIONS; USE RUNNING MASSES-- | |
34 | FL=G0(4) | |
35 | FB=G0(5) | |
36 | FT=G0(6) | |
37 | ML=FL*VEV*COSB | |
38 | MT=FT*VEV*SINB | |
39 | MB=FB*VEV*COSB | |
40 | SIG1T=0. | |
41 | SIG2T=-FAC*MT**2*G0(6)**2*LOG(MT**2/E/QS) | |
42 | SIG1B=-FAC*MB**2*G0(5)**2*LOG(MB**2/E/QS) | |
43 | SIG2B=0. | |
44 | SIG1L=-FAC4*ML**2*G0(4)**2*LOG(ML**2/E/QS) | |
45 | SIG2L=0. | |
46 | GGP=(G**2+GP**2)/2. | |
47 | MST1=MSS(12) | |
48 | MST2=MSS(13) | |
49 | MSB1=MSS(10) | |
50 | MSB2=MSS(11) | |
51 | MSL1=MSS(21) | |
52 | MSL2=MSS(22) | |
53 | C-----CALCULATE STOP_1 CONTRIBUTION ------------------------------- | |
54 | DELT1S=(.5*(G0(24)-G0(23))+(8*COS2W-5.)*GGP* | |
55 | $ (VP*VP-V*V)/12.)**2+G0(6)**2*V*V*(G0(12)-MU*COTB)**2 | |
56 | DT1=.5*(G0(24)-G0(23))+(8*COS2W-5.)*GGP* | |
57 | $ (VP*VP-V*V)/12. | |
58 | DMSDV2=GGP/4.-(2*DT1*(8*COS2W-5.)*GGP/12.- | |
59 | $ FT**2*MU*(G0(12)*TANB-MU))/2./SQRT(DELT1S) | |
60 | SIG1T1=FAC/2.*MST1**2*LOG(MST1**2/E/QS)*DMSDV2 | |
61 | DMSDV2=-GGP/4.+FT**2-(-2*DT1*(8*COS2W-5.)*GGP/12.+ | |
62 | $ FT**2*G0(12)*(G0(12)-MU*COTB))/2./SQRT(DELT1S) | |
63 | SIG2T1=FAC/2.*MST1**2*LOG(MST1**2/E/QS)*DMSDV2 | |
64 | C-----CALCULATE STOP_2 CONTRIBUTION ------------------------------- | |
65 | DMSDV2=GGP/4.+(2*DT1*(8*COS2W-5.)*GGP/12.- | |
66 | $ FT**2*MU*(G0(12)*TANB-MU))/2./SQRT(DELT1S) | |
67 | SIG1T2=FAC/2.*MST2**2*LOG(MST2**2/E/QS)*DMSDV2 | |
68 | DMSDV2=-GGP/4.+FT**2+(-2*DT1*(8*COS2W-5.)*GGP/12.+ | |
69 | $ FT**2*G0(12)*(G0(12)-MU*COTB))/2./SQRT(DELT1S) | |
70 | SIG2T2=FAC/2.*MST2**2*LOG(MST2**2/E/QS)*DMSDV2 | |
71 | C-----CALCULATE SBOT_1 CONTRIBUTION ------------------------------- | |
72 | DELB1S=(.5*(G0(24)-G0(22))-(4*COS2W-1.)*GGP* | |
73 | $ (VP*VP-V*V)/12.)**2+G0(5)**2*VP*VP*(G0(11)-MU*TANB)**2 | |
74 | DB1=.5*(G0(24)-G0(22))-(4*COS2W-1.)*GGP* | |
75 | $ (VP*VP-V*V)/12. | |
76 | DMSDV2=-GGP/4.+FB**2-(-2*DB1*(4*COS2W-1.)*GGP/12.+ | |
77 | $ FB**2*G0(11)*(G0(11)-MU*TANB))/2./SQRT(DELB1S) | |
78 | SIG1B1=FAC/2.*MSB1**2*LOG(MSB1**2/E/QS)*DMSDV2 | |
79 | DMSDV2=GGP/4.-(2*DB1*(4*COS2W-1.)*GGP/12.- | |
80 | $ FB**2*MU*(G0(11)*COTB-MU))/2./SQRT(DELB1S) | |
81 | SIG2B1=FAC/2.*MSB1**2*LOG(MSB1**2/E/QS)*DMSDV2 | |
82 | C-----CALCULATE SBOT_2 CONTRIBUTION ------------------------------- | |
83 | DMSDV2=-GGP/4.+FB**2+(-2*DB1*(4*COS2W-1.)*GGP/12.+ | |
84 | $ FB**2*G0(11)*(G0(11)-MU*TANB))/2./SQRT(DELB1S) | |
85 | SIG1B2=FAC/2.*MSB2**2*LOG(MSB2**2/E/QS)*DMSDV2 | |
86 | DMSDV2=GGP/4.+(2*DB1*(4*COS2W-1.)*GGP/12.- | |
87 | $ FB**2*MU*(G0(11)*COTB-MU))/2./SQRT(DELB1S) | |
88 | SIG2B2=FAC/2.*MSB2**2*LOG(MSB2**2/E/QS)*DMSDV2 | |
89 | C-----CALCULATE STAU_1 CONTRIBUTION ------------------------------- | |
90 | DELL1S=(.5*(G0(21)-G0(20))-(4*COS2W-3.)*GGP* | |
91 | $ (VP*VP-V*V)/4.)**2+G0(4)**2*VP*VP*(G0(10)-MU*TANB)**2 | |
92 | DL1=.5*(G0(21)-G0(20))-(4*COS2W-3.)*GGP* | |
93 | $ (VP*VP-V*V)/4. | |
94 | DMSDV2=-GGP/4.+FL**2-(-2*DL1*(4*COS2W-3.)*GGP/4.+ | |
95 | $ FL**2*G0(10)*(G0(10)-MU*TANB))/2./SQRT(DELL1S) | |
96 | SIG1L1=FAC4/2.*MSL1**2*LOG(MSL1**2/E/QS)*DMSDV2 | |
97 | DMSDV2=GGP/4.-(2*DL1*(4*COS2W-3.)*GGP/4.- | |
98 | $ FL**2*MU*(G0(10)*COTB-MU))/2./SQRT(DELL1S) | |
99 | SIG2L1=FAC4/2.*MSL1**2*LOG(MSL1**2/E/QS)*DMSDV2 | |
100 | C-----CALCULATE STAU_2 CONTRIBUTION ------------------------------- | |
101 | DMSDV2=-GGP/4.+FL**2+(-2*DL1*(4*COS2W-3.)*GGP/4.+ | |
102 | $ FL**2*G0(10)*(G0(10)-MU*TANB))/2./SQRT(DELL1S) | |
103 | SIG1L2=FAC4/2.*MSL2**2*LOG(MSL2**2/E/QS)*DMSDV2 | |
104 | DMSDV2=GGP/4.+(2*DL1*(4*COS2W-3.)*GGP/4.- | |
105 | $ FL**2*MU*(G0(10)*COTB-MU))/2./SQRT(DELL1S) | |
106 | SIG2L2=FAC4/2.*MSL2**2*LOG(MSL2**2/E/QS)*DMSDV2 | |
107 | C-----ADD ALL TERMS ------------------------------------------------ | |
108 | SIG1=SIG1B+SIG1B1+SIG1B2+SIG1T+SIG1T1+SIG1T2+ | |
109 | $SIG1L+SIG1L1+SIG1L2 | |
110 | SIG2=SIG2B+SIG2B1+SIG2B2+SIG2T+SIG2T1+SIG2T2+ | |
111 | $SIG2L+SIG2L1+SIG2L2 | |
112 | RETURN | |
113 | END |