1 #include "isajet/pilot.h"
3 C-----------------------------------------------------------------------
5 C This program gives stop squark branching fractions to gauginos
6 C according to Baer and Tata.
7 C If no other modes are allowed, stop -> c z_i through loops is
10 C-----------------------------------------------------------------------
11 #if defined(CERNLIB_IMPNONE)
14 #include "isajet/sslun.inc"
15 #include "isajet/ssmode.inc"
16 #include "isajet/sssm.inc"
17 #include "isajet/sspar.inc"
18 #include "isajet/sstype.inc"
20 COMPLEX ZI,ZONE,ZA,ZB,ZPP,ZPM,ZAUIZ,ZBUIZ
21 DOUBLE PRECISION SSALFS,SSMQCD
23 REAL WID,AWD(2),BW(2),FB,FT,XM,YM,THX,THY,AU1,MZ1,WIDC1
24 REAL PI,SR2,G,GP,TANB,COTB,MPL,MMI,AH
25 REAL AUIZ,MZIZ,SINT,COST,AS,BS,SNZI,THIZ
26 INTEGER IZ,ISTOP,IDSTOP
28 REAL MW1,MW2,SNW1,SNW2,CS2THW,BETA,TN2THW,SINB,COSB
29 REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS
31 DATA ZONE/(1.,0.)/,ZI/(0.,1.)/
33 C Partly duplicated from SSMASS
39 G=SQRT(4*PI*ALFAEM/SN2THW)
40 GP=G*SQRT(SN2THW/(1.-SN2THW))
44 C Reconstruct masses from SSMASS
45 ASMB=SUALFS(AMBT**2,.36,AMTP,3)
46 MBMB=AMBT*(1.-4*ASMB/3./PI)
47 MBQ=SSMQCD(DBLE(MBMB),DBLE(AMT1SS))
48 ASMT=SUALFS(AMTP**2,.36,AMTP,3)
49 MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))*
51 MTQ=SSMQCD(DBLE(MTMT),DBLE(AMT1SS))
52 FB=G*MBQ/SR2/AMW/COS(BETA)
53 FT=G*MTQ/SR2/AMW/SIN(BETA)
63 AWD(1)=-G*SNW1*SIN(GAMMAR)
64 AWD(2)=-G*SNW2*THY*COS(GAMMAR)
65 BW(1)=-FT*SNW1*COS(GAMMAR)
66 BW(2)=FT*SNW2*THY*SIN(GAMMAR)
67 BWP(1)=-FB*COS(GAMMAL)
68 BWP(2)=FB*THX*SIN(GAMMAL)
76 C Compute stop_i branching fractions to charm + zi if no other
77 C modes are allowed. WIDC1 is an unknown amplitude arbitrarily
84 AU1=-G/SR2*ZMIXSS(3,1)-GP/3./SR2*ZMIXSS(4,1)
96 IF(AMSTOP.LT.(MW1+AMBT).AND.AMSTOP.GT.(AMCH+MZ1)) THEN
99 AUIZ=-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ)
100 IF (AMT1SS.GT.(AMCH+MZIZ)) THEN
101 WID=AUIZ**2*(AMSTOP**2-MZIZ**2)/AU1**2
102 $ /(AMSTOP**2-MZ1**2)*WIDC1
103 CALL SSSAVE(IDSTOP,WID,ISZIZ(IZ),IDCH,0,0,0)
106 ELSEIF(AMSTOP.LT.(MW1+AMBT).AND.AMSTOP.LE.(AMCH+MZ1)) THEN
107 WRITE(LOUT,1000) ISTOP
108 1000 FORMAT(' ERROR IN SSSTBF: NO ALLOWED MODE FOR STOP',I2)
112 C stop_i -> gluino + top
114 IF (AMT1SS.GT.(AMGLSS+AMTP)) THEN
115 WID=2*SSALFS(DBLE(AMT1SS**2))*AMT1SS*((1.-AMGLSS**2/AMT1SS**2-
116 $ AMTP**2/AMT1SS**2)-2*SIN(2*THETAT)*AMTP*AMGLSS/AMT1SS**2)
117 $ *SQRT(SSXLAM(1.,AMGLSS**2/AMT1SS**2,AMTP**2/AMT1SS**2))/3.
118 CALL SSSAVE(ISTP1,WID,ISGL,IDTP,0,0,0)
121 IF (AMT2SS.GT.(AMGLSS+AMTP)) THEN
122 WID=2*SSALFS(DBLE(AMT2SS**2))*AMT2SS*((1.-AMGLSS**2/AMT2SS**2-
123 $ AMTP**2/AMT2SS**2)+2*SIN(2*THETAT)*AMTP*AMGLSS/AMT2SS**2)
124 $ *SQRT(SSXLAM(1.,AMGLSS**2/AMT2SS**2,AMTP**2/AMT2SS**2))/3.
125 CALL SSSAVE(ISTP2,WID,ISGL,IDTP,0,0,0)
128 C stop_1 -> top + zino_i
132 SNZI=SIGN(1.,AMZISS(IZ))
138 ZAUIZ=ZI**(THIZ-1.)*SNZI
139 $ *(-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
140 ZBUIZ=ZI**(THIZ-1.)*4*GP*ZMIXSS(4,IZ)/3./SR2
143 ZA=((ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST
144 $ -(ZI*ZBUIZ-ZPM*FT*ZMIXSS(1,IZ))*SINT)/2.
145 ZB=((-ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST
146 $ -(ZI*ZBUIZ+ZPM*FT*ZMIXSS(1,IZ))*SINT)/2.
149 IF (AMT1SS.GT.(AMTP+MZIZ)) THEN
150 WID=(AS*(AMT1SS**2-(AMTP+MZIZ)**2)+BS
151 $ *(AMT1SS**2-(AMTP-MZIZ)**2))/8./PI/AMT1SS
152 $ *SQRT(SSXLAM(1.,AMTP**2/AMT1SS**2,MZIZ**2/AMT1SS**2))
153 CALL SSSAVE(ISTP1,WID,ISZIZ(IZ),IDTP,0,0,0)
159 IF (AMT1SS.GT.(AMBT+MW1)) THEN
160 A=AWD(1)*COST-BW(1)*SINT
162 WID=AMT1SS*((AS+BWP(1)**2*COST**2)*(1.-MW1**2/AMT1SS**2-
163 $ AMBT**2/AMT1SS**2)-4*MW1*AMBT*BWP(1)*COST*A/AMT1SS**2)
164 $ *SQRT(SSXLAM(1.,MW1**2/AMT1SS**2,AMBT**2/AMT1SS**2))/16./PI
165 CALL SSSAVE(ISTP1,WID,ISW1,IDBT,0,0,0)
167 IF (AMT1SS.GT.(AMBT+MW2)) THEN
168 A=AWD(2)*COST-BW(2)*SINT
170 WID=AMT1SS*((AS+BWP(2)**2*COST**2)*(1.-MW2**2/AMT1SS**2-
171 $ AMBT**2/AMT1SS**2)-4*MW2*AMBT*BWP(2)*COST*A/AMT1SS**2)
172 $ *SQRT(SSXLAM(1.,MW2**2/AMT1SS**2,AMBT**2/AMT1SS**2))/16./PI
173 CALL SSSAVE(ISTP1,WID,ISW2,IDBT,0,0,0)
176 IF (AMT2SS.GT.(AMBT+MW1)) THEN
177 A=AWD(1)*SINT+BW(1)*COST
179 WID=AMT2SS*((AS+BWP(1)**2*SINT**2)*(1.-MW1**2/AMT2SS**2-
180 $ AMBT**2/AMT2SS**2)-4*MW1*AMBT*BWP(1)*SINT*A/AMT2SS**2)
181 $ *SQRT(SSXLAM(1.,MW1**2/AMT2SS**2,AMBT**2/AMT2SS**2))/16./PI
182 CALL SSSAVE(ISTP2,WID,ISW1,IDBT,0,0,0)
184 IF (AMT2SS.GT.(AMBT+MW2)) THEN
185 A=AWD(2)*SINT+BW(2)*COST
187 WID=AMT2SS*((AS+BWP(2)**2*SINT**2)*(1.-MW2**2/AMT2SS**2-
188 $ AMBT**2/AMT2SS**2)-4*MW2*AMBT*BWP(2)*SINT*A/AMT2SS**2)
189 $ *SQRT(SSXLAM(1.,MW2**2/AMT2SS**2,AMBT**2/AMT2SS**2))/16./PI
190 CALL SSSAVE(ISTP2,WID,ISW2,IDBT,0,0,0)
193 C stop_2 -> stop_1 + X modes
195 IF (AMT2SS.GT.(AMT1SS+AMZ)) THEN
196 WID=G**2*COST**2*SINT**2
197 $ *(SQRT(SSXLAM(AMT2SS**2,AMZ**2,AMT1SS**2)))**3
198 $ /64./PI/CS2THW/AMT2SS**3/AMZ**2
199 CALL SSSAVE(ISTP2,WID,IDZ,ISTP1,0,0,0)
202 IF (AMT2SS.GT.(AMT1SS+AMHL)) THEN
203 AH=G*AMW*SIN(BETA-ALFAH)*(1.-5.*TN2THW/3.)*SINT*COST/2.
204 $ +G*AMTP*COS(2.*THETAT)*(TWOM1*SIN(ALFAH)+AAT*COS(ALFAH))/2.
206 WID=AH**2/16./PI/AMT2SS**3
207 $ *SQRT(SSXLAM(AMT2SS**2,AMHL**2,AMT1SS**2))
208 CALL SSSAVE(ISTP2,WID,ISHL,ISTP1,0,0,0)
211 IF (AMT2SS.GT.(AMT1SS+AMHH)) THEN
212 AH=-G*AMW*COS(BETA-ALFAH)*(1.-5.*TN2THW/3.)*SINT*COST/2.
213 $ +G*AMTP*COS(2.*THETAT)*(TWOM1*COS(ALFAH)-AAT*SIN(ALFAH))/2.
215 WID=AH**2/16./PI/AMT2SS**3
216 $ *SQRT(SSXLAM(AMT2SS**2,AMHH**2,AMT1SS**2))
217 CALL SSSAVE(ISTP2,WID,ISHH,ISTP1,0,0,0)
220 IF (AMT2SS.GT.(AMT1SS+AMHA)) THEN
221 AH=G*AMTP*(TWOM1-AAT/TANB)/2./AMW
222 WID=AH**2/16./PI/AMT2SS**3
223 $ *SQRT(SSXLAM(AMT2SS**2,AMHA**2,AMT1SS**2))
224 CALL SSSAVE(ISTP2,WID,ISHA,ISTP1,0,0,0)
227 C t_i --> b_i + W decays
229 IF (AMT1SS.GT.(AMB1SS+AMW)) THEN
230 WID=G**2*COST**2*COSB**2*(SSXLAM(AMT1SS**2,AMB1SS**2,
231 $AMW**2))**1.5/32./PI/AMT1SS**3/AMW**2
232 CALL SSSAVE(ISTP1,WID,IDW,ISBT1,0,0,0)
235 IF (AMT1SS.GT.(AMB2SS+AMW)) THEN
236 WID=G**2*COST**2*SINB**2*(SSXLAM(AMT1SS**2,AMB2SS**2,
237 $AMW**2))**1.5/32./PI/AMT1SS**3/AMW**2
238 CALL SSSAVE(ISTP1,WID,IDW,ISBT2,0,0,0)
241 IF (AMT2SS.GT.(AMB1SS+AMW)) THEN
242 WID=G**2*SINT**2*COSB**2*(SSXLAM(AMT2SS**2,AMB1SS**2,
243 $AMW**2))**1.5/32./PI/AMT2SS**3/AMW**2
244 CALL SSSAVE(ISTP2,WID,IDW,ISBT1,0,0,0)
247 IF (AMT2SS.GT.(AMB2SS+AMW)) THEN
248 WID=G**2*SINT**2*SINB**2*(SSXLAM(AMT2SS**2,AMB2SS**2,
249 $AMW**2))**1.5/32./PI/AMT2SS**3/AMW**2
250 CALL SSSAVE(ISTP2,WID,IDW,ISBT2,0,0,0)
253 C t_i --> b_i + H+ decays
255 IF (AMT1SS.GT.(AMB1SS+AMHC)) THEN
256 A=G/SR2/AMW*(AMTP*AMBT*(COTB+TANB)*SINT*SINB+
257 $(AMBT**2*TANB+AMTP**2*COTB-AMW**2*SIN(2*BETA))*
258 $COST*COSB-AMTP*(TWOM1-AAT*COTB)*SINT*COSB-AMBT*
259 $(TWOM1-AAB*TANB)*SINB*COST)
261 WID=AS*SQRT(SSXLAM(AMT1SS**2,AMB1SS**2,AMHC**2))/
263 CALL SSSAVE(ISTP1,WID,ISHC,ISBT1,0,0,0)
266 IF (AMT1SS.GT.(AMB2SS+AMHC)) THEN
267 A=G/SR2/AMW*(-AMTP*AMBT*(COTB+TANB)*SINT*COSB+
268 $(AMBT**2*TANB+AMTP**2*COTB-AMW**2*SIN(2*BETA))*
269 $COST*SINB-AMTP*(TWOM1-AAT*COTB)*SINT*SINB+AMBT*
270 $(TWOM1-AAB*TANB)*COST*COSB)
272 WID=AS*SQRT(SSXLAM(AMT1SS**2,AMB2SS**2,AMHC**2))/
274 CALL SSSAVE(ISTP1,WID,ISHC,ISBT2,0,0,0)
277 IF (AMT2SS.GT.(AMB1SS+AMHC)) THEN
278 A=G/SR2/AMW*(-AMTP*AMBT*(COTB+TANB)*COST*SINT+
279 $(AMBT**2*TANB+AMTP**2*COTB-AMW**2*SIN(2*BETA))*
280 $SINT*COSB+AMTP*(TWOM1-AAT*COTB)*COST*COSB-AMBT*
281 $(TWOM1-AAB*TANB)*SINT*SINB)
283 WID=AS*SQRT(SSXLAM(AMT2SS**2,AMB1SS**2,AMHC**2))/
285 CALL SSSAVE(ISTP2,WID,ISHC,ISBT1,0,0,0)
288 IF (AMT2SS.GT.(AMB2SS+AMHC)) THEN
289 A=G/SR2/AMW*(AMTP*AMBT*(COTB+TANB)*COST*COSB+
290 $(AMBT**2*TANB+AMTP**2*COTB-AMW**2*SIN(2*BETA))*
291 $SINT*SINB+AMTP*(TWOM1-AAT*COTB)*SINB*COST+AMBT*
292 $(TWOM1-AAB*TANB)*COSB*SINT)
294 WID=AS*SQRT(SSXLAM(AMT2SS**2,AMB2SS**2,AMHC**2))/
296 CALL SSSAVE(ISTP2,WID,ISHC,ISBT2,0,0,0)
300 C stop_2 -> top + zino_i
304 SNZI=SIGN(1.,AMZISS(IZ))
310 ZAUIZ=ZI**(THIZ-1.)*SNZI
311 $ *(-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ))
312 ZBUIZ=ZI**(THIZ-1.)*4*GP*ZMIXSS(4,IZ)/3./SR2
315 ZA=((ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*SINT
316 $ +(ZI*ZBUIZ-ZPM*FT*ZMIXSS(1,IZ))*COST)/2.
317 ZB=((-ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*SINT
318 $ +(ZI*ZBUIZ+ZPM*FT*ZMIXSS(1,IZ))*COST)/2.
321 IF (AMT2SS.GT.(AMTP+MZIZ)) THEN
322 WID=(AS*(AMT2SS**2-(AMTP+MZIZ)**2)+BS
323 $ *(AMT2SS**2-(AMTP-MZIZ)**2))/8./PI/AMT2SS
324 $ *SQRT(SSXLAM(1.,AMTP**2/AMT2SS**2,MZIZ**2/AMT2SS**2))
325 CALL SSSAVE(ISTP2,WID,ISZIZ(IZ),IDTP,0,0,0)
329 C Normalize branching ratios