1 #include "isajet/pilot.h"
2 SUBROUTINE SSMASS(XM1,XM2,IALLOW,ILOOP,MHLNEG,MHCNEG,IMODEL)
3 C-----------------------------------------------------------------------
5 C Diagonalize neutralino, chargino, and Higgs mass matrices
6 C and save results in /SSPAR/.
8 C If XM1, XM2 < 1E19, use them for the U(1) and SU(2) mass
9 C terms. Otherwise calculate them from AMGLSS and unification.
11 C Return IALLOW = 1 if Z1SS is not LSP
12 C IALLOW = 0 otherwise
14 C-----------------------------------------------------------------------
15 #if defined(CERNLIB_IMPNONE)
18 #include "isajet/sslun.inc"
19 #include "isajet/sssm.inc"
20 #include "isajet/sspar.inc"
21 #include "isajet/ssinf.inc"
24 INTEGER IALLOW,MHLNEG,MHCNEG,IMODEL
25 REAL AR(4,4),WORK(4),WR(4)
26 REAL ZETA,ZETAS,YM,XM,COS2A,SINA,AL,SIN2A,COSA,MU2,GP,G,
27 $TEMP,VS,VP,V,MTAMTA,MTAMB,MTAMZ,ASMB,MBMB,
28 $ASMT,MTMT,SUALFE,SUALFS
29 REAL MW1,MW2,THX,THY,MU1
30 REAL COSB,SINB,BE,COS2B,SIN2B,PI,SR2,HIGFRZ,MTQ,MBQ
31 REAL TERM1,TERM2,TERM3,TANTHT,AMGLMZ,SSPOLE,TANTHB,TANTHL
32 REAL CS2THW,DELCHI,AM2
33 DOUBLE PRECISION SSMQCD
34 COMPLEX*16 SSB0,SSB1,ZZZ
36 INTEGER I,J,K,IERR,ILOOP
42 G=SQRT(4.*PI*ALFAEM/SN2THW)
43 GP=G*SQRT(SN2THW/(1.-SN2THW))
51 HIGFRZ=MAX(AMZ,SQRT(AMTLSS*AMTRSS))
53 C Compute m(tau), m(b) at z scale using qcd, qed
55 MTAMTA=AMTAU*(1.-SUALFE(AMTAU**2)/PI)
56 MTAMB=MTAMTA*(SUALFE(AMBT**2)/SUALFE(AMTAU**2))**(-27./76.)
57 MTAMZ=MTAMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.)
58 ASMB=SUALFS(AMBT**2,.36,AMTP,3)
59 MBMB=AMBT*(1.-4*ASMB/3./PI)
60 MBQ=SSMQCD(DBLE(MBMB),DBLE(HIGFRZ))
61 ASMT=SUALFS(AMTP**2,.36,AMTP,3)
62 MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
64 MTQ=SSMQCD(DBLE(MTMT),DBLE(HIGFRZ))
66 C Light/heavy stop states and mixing angle
68 TERM1=(AMTLSS**2+AMTRSS**2)/2.+AMZ**2*COS2B/4.+MTQ**2
69 TERM2=((AMTLSS**2-AMTRSS**2)/2.+COS2B*(8.*AMW**2-5.*AMZ**2)
71 TERM3=SQRT(TERM2+MTQ**2*(TWOM1*COSB/SINB+AAT)**2)
72 IF (TERM1.GT.TERM3) THEN
73 AMT1SS=SQRT(TERM1-TERM3)
77 AMT2SS=SQRT(TERM1+TERM3)
78 IF (AAT.NE.TWOM1*COSB/SINB) THEN
79 TANTHT=(AMT1SS**2-MTQ**2+AMZ**2*COS2B*(-.5+2*SN2THW/3.)-
80 $ AMTLSS**2)/MTQ/(TWOM1*COSB/SINB+AAT)
86 C Light/heavy sbottom states and mixing angle
88 TERM1=(AMBLSS**2+AMBRSS**2)/2.-AMZ**2*COS2B/4.+MBQ**2
89 TERM2=((AMBLSS**2-AMBRSS**2)/2.-COS2B*(4.*AMW**2-AMZ**2)
91 TERM3=SQRT(TERM2+MBQ**2*(TWOM1*SINB/COSB+AAB)**2)
92 IF (TERM1.GT.TERM3) THEN
93 AMB1SS=SQRT(TERM1-TERM3)
97 AMB2SS=SQRT(TERM1+TERM3)
98 TANTHB=(AMB1SS**2-MBQ**2+AMZ**2*COS2B*(.5-SN2THW/3.)-
99 $AMBLSS**2)/MBQ/(TWOM1*SINB/COSB+AAB)
102 C Light/heavy stau states and mixing angle
104 TERM1=(AMLLSS**2+AMLRSS**2)/2.-AMZ**2*COS2B/4.+MTAMZ**2
105 TERM2=((AMLLSS**2-AMLRSS**2)/2.-COS2B*(4.*AMW**2-3*AMZ**2)
107 TERM3=SQRT(TERM2+MTAMZ**2*(TWOM1*SINB/COSB+AAL)**2)
108 C if stau mass^2<0, then set to tiny mass so point is excluded
109 IF (TERM1.GT.TERM3) THEN
110 AML1SS=SQRT(TERM1-TERM3)
114 AML2SS=SQRT(TERM1+TERM3)
115 TANTHL=(AML1SS**2-MTAMZ**2+AMZ**2*COS2B*(.5-SN2THW)-
116 $AMLLSS**2)/MTAMZ/(TWOM1*SINB/COSB+AAL)
119 C define msbar gluino mass at mz from physical gluino mass
120 AMGLMZ=SSPOLE(AMGLSS,AMZ**2,-ALFA3)
121 VS=2.*AMW**2/G**2/(1.+RV2V1**2)
125 C Use either explicit values or scaling to determine SU(2)
126 C and U(1) mass terms. NOTE SIGN CONVENTION!
128 IF(ABS(XM2).LT.1.E19.AND.ABS(XM1).LT.1.E19) THEN
132 MU2=-ALFA2*AMGLMZ/ALFA3
133 MU1=5*SN2THW/3./(1.-SN2THW)*MU2
136 C Neutralino mass matrix
155 CALL EISRS1(4,4,AR,WR,ZMIXSS,IERR,WORK)
157 WRITE(LOUT,*) 'EISRS1 ERROR IN SSMASS, IERR=',IERR
161 C Sort eigenvectors and eigenvalues according to masses
165 IF (ABS(WR(I)).GT.ABS(WR(J))) THEN
171 ZMIXSS(K,J)=ZMIXSS(K,I)
183 C Chargino mass matrix
190 ZETAS=(TWOM1**2-MU2**2)**2
191 $+4*AMW**2*(AMW**2*COS2A**2+TWOM1**2+MU2**2+2*TWOM1*MU2*SIN2A)
193 XM=-(TWOM1**2-MU2**2-2*AMW**2*COS2A-ZETA)
194 $/(2*SR2*AMW*(MU2*SINA+TWOM1*COSA))
195 YM=-(TWOM1**2-MU2**2+2*AMW**2*COS2A-ZETA)
196 $/(2*SR2*AMW*(MU2*COSA+TWOM1*SINA))
207 IF (GAMMAL.LT.0.) GAMMAL=GAMMAL+PI
208 IF (GAMMAR.LT.0.) GAMMAR=GAMMAR+PI
211 AMW2SS=THX*THY*(COS(GAMMAR)*(MU2*COS(GAMMAL)+G*VP*SIN(GAMMAL))
212 $-SIN(GAMMAR)*(-G*V*COS(GAMMAL)-TWOM1*SIN(GAMMAL)))
213 AMW1SS=SIN(GAMMAR)*(MU2*SIN(GAMMAL)-G*VP*COS(GAMMAL))
214 $+COS(GAMMAR)*(-G*V*SIN(GAMMAL)+TWOM1*COS(GAMMAL))
215 C IMPLEMENT INO MASS SPLITTING FOR AMSB MODELS
219 DELCHI=G**2*MW1/8./PI**2*(2*CS2THW*REAL8(SSB0(MW1**2,MW1,AMZ))+
220 $2*SN2THW*REAL8(SSB0(MW1**2,MW1,0.))-2*REAL8(SSB0(MW1**2,MW1,AMW))
221 $-CS2THW*REAL8(SSB1(MW1**2,MW1,AMZ))-SN2THW*
222 $REAL8(SSB1(MW1**2,MW1,0.))+REAL8(SSB1(MW1**2,MW1,AMW)))
223 AMW1SS=AMW1SS+SIGN(1.,AMW1SS)*DELCHI
227 C Check validity of parameters
229 IF (IMODEL.EQ.1.OR.IMODEL.EQ.7) THEN
230 IF(MW1.LE.ABS(AMZ1SS)) IALLOW=1
231 IF(AMT1SS.LE.ABS(AMZ1SS)) IALLOW=1
232 IF(AMB1SS.LE.ABS(AMZ1SS)) IALLOW=1
233 IF(AML1SS.LE.ABS(AMZ1SS)) IALLOW=1
235 C IF(IALLOW.NE.0) RETURN