]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ISAJET/isasusy/ssstbf.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Separated TOF libraries (base,rec,sim)
[u/mrichter/AliRoot.git] / ISAJET / isasusy / ssstbf.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE SSSTBF
3C-----------------------------------------------------------------------
4C
5C This program gives stop squark branching fractions to gauginos
6C according to Baer and Tata.
7C If no other modes are allowed, stop -> c z_i through loops is
8C used as the default.
9C
10C-----------------------------------------------------------------------
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"
19C
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.)/
32C
33C Partly duplicated from SSMASS
34C
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)
44C 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)
62C
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)
75C
76C Compute stop_i branching fractions to charm + zi if no other
77C modes are allowed. WIDC1 is an unknown amplitude arbitrarily
78C set equal to 1.0.
79C
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
105110 CONTINUE
106 ELSEIF(AMSTOP.LT.(MW1+AMBT).AND.AMSTOP.LE.(AMCH+MZ1)) THEN
107 WRITE(LOUT,1000) ISTOP
1081000 FORMAT(' ERROR IN SSSTBF: NO ALLOWED MODE FOR STOP',I2)
109 END IF
110100 CONTINUE
111C
112C stop_i -> gluino + top
113C
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
120C
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
127C
128C stop_1 -> top + zino_i
129C
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
155200 CONTINUE
156C
157C Wino decays
158C
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
175C
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
192C
193C stop_2 -> stop_1 + X modes
194C
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
201C
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
210C
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
219C
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
226C
227C t_i --> b_i + W decays
228C
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
234C
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
240C
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
246C
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
252C
253C t_i --> b_i + H+ decays
254C
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
265C
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
276C
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
287C
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
298C
299C
300C stop_2 -> top + zino_i
301C
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
327500 CONTINUE
328C
329C Normalize branching ratios
330C
331 CALL SSNORM(ISTP1)
332 CALL SSNORM(ISTP2)
333C
334 RETURN
335 END