1 #include "isajet/pilot.h"
2 SUBROUTINE SSTEST(IALLOW)
4 C Test MSSM parameters against existing bounds on SUSY from
6 C IALLOW = 1 Z1 is not LSP
7 C IALLOW = 2 Gamma(Z -> Z1SS Z1SS) < GAMINV
8 C IALLOW = 4 Z -> charginos allowed
9 C IALLOW = 8 BF(Z -> Z1SS Z2SS)>10^5
10 C IALLOW = 16 Z -> squarks, sleptons
11 C IALLOW = 32 BR(Z -> Z* HL0) < B(Z -> Z* H(M=MHSM))
12 C IALLOW = 64 BR(Z -> HL0 HA0) > 0
13 C IALLOW = 128 M(H+) > M(Z)/2
14 C where GAMINV is the present bound on the invisible width,
15 C and MHSM is the lower bound on the standard Higgs mass.
17 C Bounds on the other modes are only approximate, but the
18 C error in the allowed region of masses must be tiny.
19 C Updated by H. Baer on 5/25/95
21 #if defined(CERNLIB_IMPNONE)
24 #include "isajet/sslun.inc"
25 #include "isajet/ssmode.inc"
26 #include "isajet/sssm.inc"
27 #include "isajet/sspar.inc"
28 #include "isajet/sstmp.inc"
29 #include "isajet/sstype.inc"
30 #include "isajet/sugmg.inc"
31 #include "isajet/xmssm.inc"
35 REAL MHSM,GAMINV,PI,SR2,G,GP,MZ,MZ1,MZ2,MZ3,MZ4,MW1,MW2,
36 $TANB,BETA,COS2B,SIN2B,VS,V,VP,FT,MHL,ALPHA,SUSYCC,
37 $GAMSS,W11,GZ1Z1,GAMSM,SSXINT,SSXLAM,COS2W,
38 $W12,GZ1Z2,DGAMZ,BFZ,BFZ1Z2
41 DATA MHSM/64./,GAMINV/.0043/,DGAMZ/.0115/,BFZ/1.E-5/
48 G=SQRT(4*PI*ALFAEM/SN2THW)
49 GP=G*SQRT(SN2THW/(1.-SN2THW))
59 C Check that Z1SS is LSP
61 IF(MZ1.GT.MW1.OR.MZ1.GT.AMGLSS.OR.MZ1.GT.AMULSS
62 $.OR.MZ1.GT.AMERSS.OR.MZ1.GT.AMELSS.OR.MZ1.GT.AMN1SS
63 $.OR.MZ1.GT.AMB1SS.OR.MZ1.GT.AMT1SS.OR.MZ1.GT.AML1SS) THEN
69 IF (AMZ.GT.2*MZ1) THEN
71 $ *(ZMIXSS(1,1)*ZMIXSS(1,1)-ZMIXSS(2,1)*ZMIXSS(2,1))/4.
72 GZ1Z1=SQRT(SSXLAM(MZ**2,MZ1**2,MZ1**2))/12./PI/MZ**3*W11**2
73 $ *(2*MZ**2-MZ1**2-MZ1**2-(MZ1**2-MZ1**2)**2/MZ**2
74 $ -6*MZ1*MZ1*SIGN(1.,AMZ1SS*AMZ1SS))
75 IF(GZ1Z1.GT.GAMINV) THEN
80 C Check for other allowed visible modes modes
86 C Check funny Z branching fractions
89 IF (AMZ.GT.MZ1+MZ2) THEN
91 $ *(ZMIXSS(1,1)*ZMIXSS(1,2)-ZMIXSS(2,1)*ZMIXSS(2,2))/4.
92 GZ1Z2=SQRT(SSXLAM(MZ**2,MZ1**2,MZ2**2))/6./PI/MZ**3*W12**2
93 $ *(2*MZ**2-MZ1**2-MZ2**2-(MZ1**2-MZ2**2)**2/MZ**2
94 $ -6*MZ1*MZ2*SIGN(1.,AMZ1SS*AMZ2SS))
97 IF(BFZ1Z2.GT.BFZ) THEN
101 IF(AMZ.GT.2*AMULSS.OR.AMZ.GT.2*AMELSS.OR.AMZ.GT.2*AMERSS
102 $.OR.AMZ.GT.2*AMN1SS.OR.AMZ.GT.2*AMB1SS.OR.AMZ.GT.2*AMT1SS)THEN
109 GAMSM=SSXINT(2*MHSM/MZ,SSZHX,(1.+MHSM**2/MZ**2))
116 VS=2*AMW**2/G**2/(1.+RV2V1**2)
119 FT=G*AMTP/SR2/AMW/V*SQRT(V**2+VP**2)
122 SUSYCC=SIN(ALPHA+BETA)
124 GAMSS=SSXINT(2*MHL/AMZ,SSZHX,(1.+MHL**2/AMZ**2))*SUSYCC**2
125 IF(GAMSS.GE.GAMSM) IALLOW=IALLOW+32
128 IF (AMZ.GT.(AMHL+AMHA)) THEN
132 IF(AMZ.GT.2*AMHC) THEN