]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | SUBROUTINE SSSTBF | |
3 | C----------------------------------------------------------------------- | |
4 | 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 | |
8 | C used as the default. | |
9 | C | |
10 | C----------------------------------------------------------------------- | |
11 | #if defined(CERNLIB_IMPNONE) | |
12 | IMPLICIT NONE | |
13 | #endif | |
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" | |
19 | C | |
20 | COMPLEX ZI,ZONE,ZA,ZB,ZPP,ZPM,ZAUIZ,ZBUIZ | |
21 | DOUBLE PRECISION SSALFS,SSMQCD | |
22 | REAL SSXLAM | |
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 | |
27 | REAL AMSTOP,BWP(2),A | |
28 | REAL MW1,MW2,SNW1,SNW2,CS2THW,BETA,TN2THW,SINB,COSB | |
29 | REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS | |
30 | INTEGER ISZIZ(4) | |
31 | DATA ZONE/(1.,0.)/,ZI/(0.,1.)/ | |
32 | C | |
33 | C Partly duplicated from SSMASS | |
34 | C | |
35 | CS2THW=1.-SN2THW | |
36 | TN2THW=SN2THW/CS2THW | |
37 | PI=4.*ATAN(1.) | |
38 | SR2=SQRT(2.) | |
39 | G=SQRT(4*PI*ALFAEM/SN2THW) | |
40 | GP=G*SQRT(SN2THW/(1.-SN2THW)) | |
41 | TANB=1./RV2V1 | |
42 | COTB=RV2V1 | |
43 | BETA=ATAN(TANB) | |
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))* | |
50 | $(ASMT/PI)**2) | |
51 | MTQ=SSMQCD(DBLE(MTMT),DBLE(AMT1SS)) | |
52 | FB=G*MBQ/SR2/AMW/COS(BETA) | |
53 | FT=G*MTQ/SR2/AMW/SIN(BETA) | |
54 | MW1=ABS(AMW1SS) | |
55 | MW2=ABS(AMW2SS) | |
56 | SNW1=SIGN(1.,AMW1SS) | |
57 | SNW2=SIGN(1.,AMW2SS) | |
58 | XM=1./TAN(GAMMAL) | |
59 | YM=1./TAN(GAMMAR) | |
60 | THX=SIGN(1.,XM) | |
61 | THY=SIGN(1.,YM) | |
62 | C | |
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) | |
69 | MMI=AMW1SS | |
70 | MPL=AMW2SS | |
71 | COST=COS(THETAT) | |
72 | SINT=SIN(THETAT) | |
73 | COSB=COS(THETAB) | |
74 | SINB=SIN(THETAB) | |
75 | C | |
76 | C Compute stop_i branching fractions to charm + zi if no other | |
77 | C modes are allowed. WIDC1 is an unknown amplitude arbitrarily | |
78 | C set equal to 1.0. | |
79 | C | |
80 | ISZIZ(1)=ISZ1 | |
81 | ISZIZ(2)=ISZ2 | |
82 | ISZIZ(3)=ISZ3 | |
83 | ISZIZ(4)=ISZ4 | |
84 | AU1=-G/SR2*ZMIXSS(3,1)-GP/3./SR2*ZMIXSS(4,1) | |
85 | MZ1=ABS(AMZISS(1)) | |
86 | DO 100 ISTOP=1,2 | |
87 | IF(ISTOP.EQ.1) THEN | |
88 | AMSTOP=AMT1SS | |
89 | IDSTOP=ISTP1 | |
90 | WIDC1=1.0E-6 | |
91 | ELSE | |
92 | AMSTOP=AMT2SS | |
93 | IDSTOP=ISTP2 | |
94 | WIDC1=1.0E-6 | |
95 | ENDIF | |
96 | IF(AMSTOP.LT.(MW1+AMBT).AND.AMSTOP.GT.(AMCH+MZ1)) THEN | |
97 | DO 110 IZ=1,4 | |
98 | MZIZ=ABS(AMZISS(IZ)) | |
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) | |
104 | END IF | |
105 | 110 CONTINUE | |
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) | |
109 | END IF | |
110 | 100 CONTINUE | |
111 | C | |
112 | C stop_i -> gluino + top | |
113 | C | |
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) | |
119 | END IF | |
120 | C | |
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) | |
126 | END IF | |
127 | C | |
128 | C stop_1 -> top + zino_i | |
129 | C | |
130 | DO 200 IZ=1,4 | |
131 | MZIZ=ABS(AMZISS(IZ)) | |
132 | SNZI=SIGN(1.,AMZISS(IZ)) | |
133 | IF (SNZI.EQ.1.) THEN | |
134 | THIZ=0. | |
135 | ELSE | |
136 | THIZ=1. | |
137 | END IF | |
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 | |
141 | ZPP=ZI**THIZ | |
142 | ZPM=(-ZI)**THIZ | |
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. | |
147 | AS=ZA*CONJG(ZA) | |
148 | BS=ZB*CONJG(ZB) | |
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) | |
154 | END IF | |
155 | 200 CONTINUE | |
156 | C | |
157 | C Wino decays | |
158 | C | |
159 | IF (AMT1SS.GT.(AMBT+MW1)) THEN | |
160 | A=AWD(1)*COST-BW(1)*SINT | |
161 | AS=A*A | |
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) | |
166 | END IF | |
167 | IF (AMT1SS.GT.(AMBT+MW2)) THEN | |
168 | A=AWD(2)*COST-BW(2)*SINT | |
169 | AS=A*A | |
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) | |
174 | END IF | |
175 | C | |
176 | IF (AMT2SS.GT.(AMBT+MW1)) THEN | |
177 | A=AWD(1)*SINT+BW(1)*COST | |
178 | AS=A*A | |
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) | |
183 | END IF | |
184 | IF (AMT2SS.GT.(AMBT+MW2)) THEN | |
185 | A=AWD(2)*SINT+BW(2)*COST | |
186 | AS=A*A | |
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) | |
191 | END IF | |
192 | C | |
193 | C stop_2 -> stop_1 + X modes | |
194 | C | |
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) | |
200 | END IF | |
201 | C | |
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. | |
205 | $ /AMW/SIN(BETA) | |
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) | |
209 | END IF | |
210 | C | |
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. | |
214 | $ /AMW/SIN(BETA) | |
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) | |
218 | END IF | |
219 | C | |
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) | |
225 | END IF | |
226 | C | |
227 | C t_i --> b_i + W decays | |
228 | C | |
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) | |
233 | END IF | |
234 | C | |
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) | |
239 | END IF | |
240 | C | |
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) | |
245 | END IF | |
246 | C | |
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) | |
251 | END IF | |
252 | C | |
253 | C t_i --> b_i + H+ decays | |
254 | C | |
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) | |
260 | AS=A*A | |
261 | WID=AS*SQRT(SSXLAM(AMT1SS**2,AMB1SS**2,AMHC**2))/ | |
262 | $ 16./PI/AMT1SS**3 | |
263 | CALL SSSAVE(ISTP1,WID,ISHC,ISBT1,0,0,0) | |
264 | END IF | |
265 | C | |
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) | |
271 | AS=A*A | |
272 | WID=AS*SQRT(SSXLAM(AMT1SS**2,AMB2SS**2,AMHC**2))/ | |
273 | $ 16./PI/AMT1SS**3 | |
274 | CALL SSSAVE(ISTP1,WID,ISHC,ISBT2,0,0,0) | |
275 | END IF | |
276 | C | |
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) | |
282 | AS=A*A | |
283 | WID=AS*SQRT(SSXLAM(AMT2SS**2,AMB1SS**2,AMHC**2))/ | |
284 | $ 16./PI/AMT2SS**3 | |
285 | CALL SSSAVE(ISTP2,WID,ISHC,ISBT1,0,0,0) | |
286 | END IF | |
287 | C | |
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) | |
293 | AS=A*A | |
294 | WID=AS*SQRT(SSXLAM(AMT2SS**2,AMB2SS**2,AMHC**2))/ | |
295 | $ 16./PI/AMT2SS**3 | |
296 | CALL SSSAVE(ISTP2,WID,ISHC,ISBT2,0,0,0) | |
297 | END IF | |
298 | C | |
299 | C | |
300 | C stop_2 -> top + zino_i | |
301 | C | |
302 | DO 500 IZ=1,4 | |
303 | MZIZ=ABS(AMZISS(IZ)) | |
304 | SNZI=SIGN(1.,AMZISS(IZ)) | |
305 | IF (SNZI.EQ.1.) THEN | |
306 | THIZ=0. | |
307 | ELSE | |
308 | THIZ=1. | |
309 | END IF | |
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 | |
313 | ZPP=ZI**THIZ | |
314 | ZPM=(-ZI)**THIZ | |
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. | |
319 | AS=ZA*CONJG(ZA) | |
320 | BS=ZB*CONJG(ZB) | |
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) | |
326 | END IF | |
327 | 500 CONTINUE | |
328 | C | |
329 | C Normalize branching ratios | |
330 | C | |
331 | CALL SSNORM(ISTP1) | |
332 | CALL SSNORM(ISTP2) | |
333 | C | |
334 | RETURN | |
335 | END |