]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | SUBROUTINE SSZIBF | |
3 | C----------------------------------------------------------------------- | |
4 | C This subroutine calculates the neutralino (zi) subset of | |
5 | C SSWZBF, which was too long. | |
6 | C Valid for all scalar masses (functions in double precision) | |
7 | C Includes Higgs sector radiative corrections (Aug. 31) | |
8 | C | |
9 | C Auxiliary functions are called SSWxyi, SSZxyi, where normally | |
10 | C x indicates the SUSY particle, y the SM particle(s), and i is | |
11 | C a counter. | |
12 | C | |
13 | C Part of Baer's GAUGBF | |
14 | C | |
15 | C----------------------------------------------------------------------- | |
16 | #if defined(CERNLIB_IMPNONE) | |
17 | IMPLICIT NONE | |
18 | #endif | |
19 | #include "isajet/sslun.inc" | |
20 | #include "isajet/ssmode.inc" | |
21 | #include "isajet/sssm.inc" | |
22 | #include "isajet/sspar.inc" | |
23 | #include "isajet/sstype.inc" | |
24 | #include "isajet/sstmp.inc" | |
25 | #include "isajet/sspols.inc" | |
26 | C | |
27 | EXTERNAL SSZWF1,SSZZF1,SSZZF2,SSZZF3,SSWZF1,SSWZF2,SSWZF3 | |
28 | $,SSWZF4,SSWZF5,SSWZF6,SSWZF7,SSWWF1,SSZZF4,SSZZF5,SSGX1 | |
29 | $,SSGX2,SSGX8 | |
30 | C | |
31 | REAL AUI(4),BUI(4),ADI(4),BDI(4),ALI(4),BLI(4),ANI(4),BNI(4) | |
32 | $,WIJ(4,4),AUWI(4),ADWI(4),ANWI(4),ALWI(4),XIM(4),YIM(4) | |
33 | $,XIP(4),YIP(4),SNIJ(4,4),XLIJ(4,4),HIJ(4,4) | |
34 | $,V1I(4),V2I(4),V3I(4),V4I(4),XHIJ(4,4),XPIJ(4,4),AMWISS(2) | |
35 | C | |
36 | INTEGER ISZ(4),THIZ,THJZ | |
37 | C | |
38 | REAL MWIW,SL,PP,SP,PL,MZIZ,PH,SH,COSB,SINB,COSA,SINA | |
39 | $,MZJZ,FAC,COSBE,SINBE | |
40 | REAL STHW,CTHW | |
41 | REAL EF,A,B,TANB,FB,FT,SR2,G,PI,GP,FL | |
42 | $,MW2,SNW1,MW1,YM,BE,SNW2,XM,THX,THY | |
43 | $,BTN,APD,APL,APU,BTD,BTL,APN,BTU | |
44 | REAL TANW,COTW,XWINO,YWINO,SN,SNIW,SNJZ | |
45 | REAL SSXINT,SSXLAM | |
46 | REAL WID,WIDEE,E, | |
47 | $TERMLL,TERMRR,TERMZZ,TERMLZ,TERMRZ,TERMLH | |
48 | REAL FACT,ALIZ1,ALJZ1,ALIZ2,ALJZ2,BEIZ1,BEJZ1,BEIZ2,BEJZ2, | |
49 | $SGNIJ,XUPPER,GLLF1,GRRF1,GLRF1,GF1,GLLF2,GRRF2,GLRF2,GF2, | |
50 | $GLGL,GRGR,GLGR,GRGL,GF12,GF,TERMHL,TERMHH,TERMHA, | |
51 | $TERM1Z,TERM2Z,TERM1L,TERM2L,TERM1H,TERM2H,TERM1A,TERM2A | |
52 | REAL COSL,SINL,BPWI(2),BPLWI(2) | |
53 | REAL BWI(2),AS,BS,COST,SINT,SNIZ | |
54 | REAL GLLF1L,GRRF1R,GLGLL,GRGRR,GLLF2L,GRRF2R,TMZZRL,TMZZLR, | |
55 | $TM1ZRL,TM1ZLR,TM2ZRL,TM2ZLR | |
56 | REAL ASMB,MBMB,MBQ,ASMT,MTMT,MTQ,SUALFS | |
57 | REAL SUALFE,MTAMTA,MTAMB,MTAMZ,AMPL | |
58 | REAL FUDGE | |
59 | DOUBLE PRECISION SSMQCD | |
60 | COMPLEX ZI,ZONE,ZA,ZB,ZPP,ZPM,ZAUIZ,ZBUIZ,ZADIZ,ZBDIZ, | |
61 | $ZALIZ,ZBLIZ,Z1(2),Z2(2) | |
62 | INTEGER IW,JZ,IZ,ISZIZ,ISWIW | |
63 | DATA FUDGE/1.0/ | |
64 | DATA ZONE/(1.,0.)/,ZI/(0.,1.)/ | |
65 | C | |
66 | C Constants from neutralino mass matrix | |
67 | C | |
68 | AMPL=2.4E18 | |
69 | PI=4.*ATAN(1.) | |
70 | SR2=SQRT(2.) | |
71 | G=SQRT(4*PI*ALFAEM/SN2THW) | |
72 | GP=G*SQRT(SN2THW/(1.-SN2THW)) | |
73 | E=SQRT(4*PI/128.) | |
74 | C | |
75 | TANW=SQRT(SN2THW/(1.-SN2THW)) | |
76 | COTW=1./TANW | |
77 | STHW=SQRT(SN2THW) | |
78 | CTHW=SQRT(1.-SN2THW) | |
79 | APL=.25*(3*TANW-COTW) | |
80 | BTL=.25*(COTW+TANW) | |
81 | APN=.25*(TANW+COTW) | |
82 | BTN=-.25*(COTW+TANW) | |
83 | APU=-5*TANW/12.+COTW/4. | |
84 | BTU=-.25*(COTW+TANW) | |
85 | APD=-COTW/4.+TANW/12. | |
86 | BTD=.25*(COTW+TANW) | |
87 | C | |
88 | TANB=1./RV2V1 | |
89 | BE=ATAN(1./RV2V1) | |
90 | SINBE=SIN(BE) | |
91 | COSBE=COS(BE) | |
92 | XM=1./TAN(GAMMAL) | |
93 | YM=1./TAN(GAMMAR) | |
94 | THX=SIGN(1.,XM) | |
95 | THY=SIGN(1.,YM) | |
96 | ASMB=SUALFS(AMBT**2,.36,AMTP,3) | |
97 | MBMB=AMBT*(1.-4*ASMB/3./PI) | |
98 | MBQ=SSMQCD(DBLE(MBMB),DBLE(AMTP)) | |
99 | ASMT=SUALFS(AMTP**2,.36,AMTP,3) | |
100 | MTMT=AMTP/(1.+4*ASMT/3./PI+(16.11-1.04*(5.-6.63/AMTP))* | |
101 | $(ASMT/PI)**2) | |
102 | MTQ=SSMQCD(DBLE(MTMT),DBLE(AMTP)) | |
103 | FB=G*MBQ/SR2/AMW/COS(BE) | |
104 | FT=G*MTQ/SR2/AMW/SIN(BE) | |
105 | MTAMTA=AMTAU*(1.-SUALFE(AMTAU**2)/PI) | |
106 | MTAMB=MTAMTA*(SUALFE(AMBT**2)/SUALFE(AMTAU**2))**(-27./76.) | |
107 | MTAMZ=MTAMB*(SUALFE(AMZ**2)/SUALFE(AMBT**2))**(-27./80.) | |
108 | FL=G*MTAMZ/SR2/AMW/COS(BE) | |
109 | SNW1=SIGN(1.,AMW1SS) | |
110 | SNW2=SIGN(1.,AMW2SS) | |
111 | AMWISS(1)=AMW1SS | |
112 | AMWISS(2)=AMW2SS | |
113 | BWI(1)=-FT*SNW1*COS(GAMMAR) | |
114 | BWI(2)=FT*SNW2*THY*SIN(GAMMAR) | |
115 | BPWI(1)=-FB*COS(GAMMAL) | |
116 | BPWI(2)=FB*THX*SIN(GAMMAL) | |
117 | BPLWI(1)=-FL*COS(GAMMAL) | |
118 | BPLWI(2)=FL*THX*SIN(GAMMAL) | |
119 | MW1=ABS(AMW1SS) | |
120 | MW2=ABS(AMW2SS) | |
121 | XWINO=.5*(THX*SIN(GAMMAL)*COS(GAMMAL) | |
122 | $-THY*SIN(GAMMAR)*COS(GAMMAR)) | |
123 | YWINO=.5*(THX*SIN(GAMMAL)*COS(GAMMAL) | |
124 | $+THY*SIN(GAMMAR)*COS(GAMMAR)) | |
125 | COST=COS(THETAT) | |
126 | SINT=SIN(THETAT) | |
127 | COSB=COS(THETAB) | |
128 | SINB=SIN(THETAB) | |
129 | COSL=COS(THETAL) | |
130 | SINL=SIN(THETAL) | |
131 | C | |
132 | C Constants from Higgs mass matrix | |
133 | C | |
134 | SINA=SIN(ALFAH) | |
135 | COSA=COS(ALFAH) | |
136 | C | |
137 | C Gaugino couplings | |
138 | C | |
139 | DO 100 IZ=1,4 | |
140 | AUI(IZ)=G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ) | |
141 | BUI(IZ)=4.*GP/3./SR2*ZMIXSS(4,IZ) | |
142 | ADI(IZ)=-G/SR2*ZMIXSS(3,IZ)+GP/3./SR2*ZMIXSS(4,IZ) | |
143 | BDI(IZ)=-2.*GP/3./SR2*ZMIXSS(4,IZ) | |
144 | ALI(IZ)=G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ) | |
145 | BLI(IZ)=-SR2*GP*ZMIXSS(4,IZ) | |
146 | ANI(IZ)=G/SR2*ZMIXSS(3,IZ)-GP/SR2*ZMIXSS(4,IZ) | |
147 | BNI(IZ)=0.0 | |
148 | 100 CONTINUE | |
149 | C | |
150 | DO 110 IZ=1,4 | |
151 | DO 110 JZ=1,4 | |
152 | IF(IZ.LT.JZ) THEN | |
153 | WIJ(IZ,JZ)=SQRT(G**2+GP**2) | |
154 | $ *(ZMIXSS(1,IZ)*ZMIXSS(1,JZ)-ZMIXSS(2,IZ)*ZMIXSS(2,JZ))/4. | |
155 | ELSEIF(IZ.GT.JZ) THEN | |
156 | WIJ(IZ,JZ)=-SQRT(G**2+GP**2) | |
157 | $ *(ZMIXSS(1,IZ)*ZMIXSS(1,JZ)-ZMIXSS(2,IZ)*ZMIXSS(2,JZ))/4. | |
158 | ENDIF | |
159 | 110 CONTINUE | |
160 | C | |
161 | AUWI(2)=G*THX*COS(GAMMAL) | |
162 | ADWI(2)=SNW2*G*THY*COS(GAMMAR) | |
163 | ALWI(2)=ADWI(2) | |
164 | ANWI(2)=AUWI(2) | |
165 | AUWI(1)=G*SIN(GAMMAL) | |
166 | ADWI(1)=SNW1*G*SIN(GAMMAR) | |
167 | ALWI(1)=ADWI(1) | |
168 | ANWI(1)=AUWI(1) | |
169 | C | |
170 | DO 120 IZ=1,4 | |
171 | XIM(IZ)=.5*(SNW1*SIGN(1.,AMZISS(IZ))*(COS(GAMMAR) | |
172 | $ *ZMIXSS(1,IZ)/SR2+SIN(GAMMAR)*ZMIXSS(3,IZ))-COS(GAMMAL) | |
173 | $ *ZMIXSS(2,IZ)/SR2+SIN(GAMMAL)*ZMIXSS(3,IZ)) | |
174 | YIM(IZ)=.5*(-SNW1*SIGN(1.,AMZISS(IZ))*(COS(GAMMAR) | |
175 | $ *ZMIXSS(1,IZ)/SR2+SIN(GAMMAR)*ZMIXSS(3,IZ))-COS(GAMMAL) | |
176 | $ *ZMIXSS(2,IZ)/SR2+SIN(GAMMAL)*ZMIXSS(3,IZ)) | |
177 | XIP(IZ)=.5*(SNW2*SIGN(1.,AMZISS(IZ))*THY*(-SIN(GAMMAR) | |
178 | $ *ZMIXSS(1,IZ)/SR2+COS(GAMMAR)*ZMIXSS(3,IZ))+THX*(SIN(GAMMAL) | |
179 | $ *ZMIXSS(2,IZ)/SR2+COS(GAMMAL)*ZMIXSS(3,IZ))) | |
180 | YIP(IZ)=.5*(-SNW2*SIGN(1.,AMZISS(IZ))*THY*(-SIN(GAMMAR) | |
181 | $ *ZMIXSS(1,IZ)/SR2+COS(GAMMAR)*ZMIXSS(3,IZ))+THX*(SIN(GAMMAL) | |
182 | $ *ZMIXSS(2,IZ)/SR2+COS(GAMMAL)*ZMIXSS(3,IZ))) | |
183 | 120 CONTINUE | |
184 | C | |
185 | DO 130 IZ=1,4 | |
186 | DO 130 JZ=1,4 | |
187 | IF(IZ.NE.JZ) THEN | |
188 | SNIJ(IZ,JZ)=-1.*SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ)) | |
189 | XLIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ)) | |
190 | $ *(ZMIXSS(2,IZ)*SINA-ZMIXSS(1,IZ)*COSA) | |
191 | $ *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2. | |
192 | XHIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ)) | |
193 | $ *(ZMIXSS(2,IZ)*COSA+ZMIXSS(1,IZ)*SINA) | |
194 | $ *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2. | |
195 | XPIJ(IZ,JZ)=SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ)) | |
196 | $ *(ZMIXSS(2,IZ)*SINBE-ZMIXSS(1,IZ)*COSBE) | |
197 | $ *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2. | |
198 | HIJ(IZ,JZ)=-SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ)) | |
199 | $ *(ZMIXSS(2,IZ)*COSA+ZMIXSS(1,IZ)*SINA) | |
200 | $ *(G*ZMIXSS(3,JZ)-GP*ZMIXSS(4,JZ))/2. | |
201 | ENDIF | |
202 | 130 CONTINUE | |
203 | C | |
204 | SP=-.5*(-THY*SNW2*COSBE*SIN(GAMMAL)*SIN(GAMMAR)+ | |
205 | $THY*SNW2*SINBE*COS(GAMMAL)*COS(GAMMAR)- | |
206 | $THX*SNW1*COSBE*COS(GAMMAL)*COS(GAMMAR)+ | |
207 | $THX*SNW1*SINBE*SIN(GAMMAL)*SIN(GAMMAR)) | |
208 | PP=-.5*(-THY*SNW2*COSBE*SIN(GAMMAL)*SIN(GAMMAR)+ | |
209 | $THY*SNW2*SINBE*COS(GAMMAL)*COS(GAMMAR)+ | |
210 | $THX*SNW1*COSBE*COS(GAMMAL)*COS(GAMMAR)- | |
211 | $THX*SNW1*SINBE*SIN(GAMMAL)*SIN(GAMMAR)) | |
212 | C | |
213 | SL=.5*(THY*SNW2*SINA*COS(GAMMAL)*COS(GAMMAR)- | |
214 | $THY*SNW2*COSA*SIN(GAMMAL)*SIN(GAMMAR)+ | |
215 | $THX*SNW1*COSA*COS(GAMMAL)*COS(GAMMAR)- | |
216 | $THX*SNW1*SINA*SIN(GAMMAL)*SIN(GAMMAR)) | |
217 | PL=.5*(THY*SNW2*SINA*COS(GAMMAL)*COS(GAMMAR)- | |
218 | $THY*SNW2*COSA*SIN(GAMMAL)*SIN(GAMMAR)- | |
219 | $THX*SNW1*COSA*COS(GAMMAL)*COS(GAMMAR)+ | |
220 | $THX*SNW1*SINA*SIN(GAMMAL)*SIN(GAMMAR)) | |
221 | C | |
222 | SH=.5*(THY*SNW2*COSA*COS(GAMMAL)*COS(GAMMAR)+ | |
223 | $THY*SNW2*SINA*SIN(GAMMAL)*SIN(GAMMAR)- | |
224 | $THX*SNW1*SINA*COS(GAMMAL)*COS(GAMMAR)- | |
225 | $THX*SNW1*COSA*SIN(GAMMAL)*SIN(GAMMAR)) | |
226 | PH=.5*(THY*SNW2*COSA*COS(GAMMAL)*COS(GAMMAR)+ | |
227 | $THY*SNW2*SINA*SIN(GAMMAL)*SIN(GAMMAR)+ | |
228 | $THX*SNW1*SINA*COS(GAMMAL)*COS(GAMMAR)+ | |
229 | $THX*SNW1*COSA*SIN(GAMMAL)*SIN(GAMMAR)) | |
230 | C | |
231 | DO 140 IZ=1,4 | |
232 | V1I(IZ)=-SIN(GAMMAR)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ)) | |
233 | $ -G*COS(GAMMAR)*ZMIXSS(1,IZ) | |
234 | V2I(IZ)=COS(GAMMAR)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ)) | |
235 | $ -G*SIN(GAMMAR)*ZMIXSS(1,IZ) | |
236 | V3I(IZ)=-SIN(GAMMAL)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ)) | |
237 | $ +G*COS(GAMMAL)*ZMIXSS(2,IZ) | |
238 | V4I(IZ)=COS(GAMMAL)/SR2*(G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ)) | |
239 | $ +G*SIN(GAMMAL)*ZMIXSS(2,IZ) | |
240 | 140 CONTINUE | |
241 | C----------------------------------------------------------------------- | |
242 | C Generate Neutralino zi Branching Fractions | |
243 | C----------------------------------------------------------------------- | |
244 | ISZ(1)=ISZ1 | |
245 | ISZ(2)=ISZ2 | |
246 | ISZ(3)=ISZ3 | |
247 | ISZ(4)=ISZ4 | |
248 | C | |
249 | DO 200 IZ=2,4 | |
250 | MZIZ=ABS(AMZISS(IZ)) | |
251 | SNIZ=SIGN(1.,AMZISS(IZ)) | |
252 | IF (SNIZ.EQ.1.) THEN | |
253 | THIZ=0 | |
254 | ELSE | |
255 | THIZ=1 | |
256 | END IF | |
257 | ISZIZ=ISZ(IZ) | |
258 | C z2 --> z1 + photon | |
259 | C !!! NEEDS UPDATING | |
260 | IF (IZ.EQ.2.AND.(ABS(ZMIXSS(3,1)).LE.1.E-4).AND. | |
261 | $ (ABS(ZMIXSS(4,1)).LE.1.E-4)) THEN | |
262 | CALL SSSAVE(ISZ2,1.,ISZ1,IDGM,0,0,0) | |
263 | GOTO 200 | |
264 | END IF | |
265 | C | |
266 | C zi --> wi + x partial widths | |
267 | C | |
268 | DO 205 IW=1,2 | |
269 | C Loop over w1, w2 | |
270 | IF(IW.EQ.1) THEN | |
271 | MWIW=MW1 | |
272 | SNIW=SNW1 | |
273 | ISWIW=ISW1 | |
274 | ELSE | |
275 | MWIW=MW2 | |
276 | SNIW=SNW2 | |
277 | ISWIW=ISW2 | |
278 | ENDIF | |
279 | C | |
280 | C zi -> wj + w | |
281 | C | |
282 | IF(MZIZ.GT.MWIW+AMW) THEN | |
283 | EF=MZIZ**2+MWIW**2-AMW**2+((MZIZ**2-MWIW**2)**2 | |
284 | $ -AMW**4)/AMW/AMW | |
285 | WID=G*G*SQRT(SSXLAM(MZIZ**2,AMW**2,MWIW**2))/32./PI/ | |
286 | $ MZIZ**3*(2.*EF*(XIM(IZ)**2+YIM(IZ)**2)-12.* | |
287 | $ MZIZ*MWIW*(XIM(IZ)**2-YIM(IZ)**2)) | |
288 | CALL SSSAVE(ISZIZ,WID,+ISWIW,-IDW,0,0,0) | |
289 | CALL SSSAVE(ISZIZ,WID,-ISWIW,+IDW,0,0,0) | |
290 | C | |
291 | C zi -> wj + f + fbar (w forbidden) ONLY W CONTRIBUTION INCLUDED! | |
292 | C | |
293 | ELSEIF(MZIZ.GT.FUDGE*MWIW) THEN | |
294 | TMP(1)=MWIW | |
295 | TMP(2)=MZIZ | |
296 | IF(IW.EQ.1) THEN | |
297 | TMP(3)=XIM(IZ) | |
298 | TMP(4)=YIM(IZ) | |
299 | ELSE | |
300 | TMP(3)=XIP(IZ) | |
301 | TMP(4)=YIP(IZ) | |
302 | ENDIF | |
303 | WID=G**4/96./PI**3/MZIZ | |
304 | $ *SSXINT(MWIW,SSZWF1,(MWIW**2+MZIZ**2)/2./MZIZ) | |
305 | CALL SSSAVE(ISZIZ,3.*WID,-ISWIW,IDUP,-IDDN,0,0) | |
306 | Z1(1)=(-ZI)**THIZ*G*TMP(3) | |
307 | Z1(2)=(-ZI)**THIZ*G*TMP(4) | |
308 | Z2(1)=G/2./SR2 | |
309 | Z2(2)=-Z2(1) | |
310 | CALL SSME3(1,AMW,Z1,Z2) | |
311 | CALL SSSAVE(ISZIZ,WID,-ISWIW,IDNE,-IDE,0,0) | |
312 | CALL SSME3(1,AMW,Z1,Z2) | |
313 | CALL SSSAVE(ISZIZ,WID,-ISWIW,IDNM,-IDMU,0,0) | |
314 | CALL SSME3(1,AMW,Z1,Z2) | |
315 | CALL SSSAVE(ISZIZ,3.*WID,ISWIW,IDDN,-IDUP,0,0) | |
316 | CALL SSME3(1,AMW,Z1,Z2) | |
317 | CALL SSSAVE(ISZIZ,WID,ISWIW,IDE,-IDNE,0,0) | |
318 | CALL SSME3(1,AMW,Z1,Z2) | |
319 | CALL SSSAVE(ISZIZ,WID,ISWIW,IDMU,-IDNM,0,0) | |
320 | CALL SSME3(1,AMW,Z1,Z2) | |
321 | IF (MZIZ.GT.(MWIW+AMCH+AMST)) THEN | |
322 | CALL SSSAVE(ISZIZ,3.*WID,-ISWIW,IDCH,-IDST,0,0) | |
323 | CALL SSME3(1,AMW,Z1,Z2) | |
324 | CALL SSSAVE(ISZIZ,WID,-ISWIW,IDNT,-IDTAU,0,0) | |
325 | CALL SSME3(1,AMW,Z1,Z2) | |
326 | CALL SSSAVE(ISZIZ,3.*WID,ISWIW,IDST,-IDCH,0,0) | |
327 | CALL SSME3(1,AMW,Z1,Z2) | |
328 | CALL SSSAVE(ISZIZ,WID,ISWIW,IDTAU,-IDNT,0,0) | |
329 | CALL SSME3(1,AMW,Z1,Z2) | |
330 | ENDIF | |
331 | ENDIF | |
332 | C | |
333 | C zi --> wj + hc | |
334 | C | |
335 | IF (MZIZ.GT.(MWIW+AMHC)) THEN | |
336 | IF (IW.EQ.1) THEN | |
337 | A=(SNW1*COSBE*V2I(IZ)-SIGN(1.,AMZISS(IZ)) | |
338 | $ *SINBE*V4I(IZ))/2. | |
339 | B=(SNW1*COSBE*V2I(IZ)+SIGN(1.,AMZISS(IZ)) | |
340 | $ *SINBE*V4I(IZ))/2. | |
341 | ELSE | |
342 | A=(THY*SNW2*COSBE*V1I(IZ)-SIGN(1.,AMZISS(IZ)) | |
343 | $ *THX*SINBE*V3I(IZ))/2. | |
344 | B=(THY*SNW2*COSBE*V1I(IZ)+SIGN(1.,AMZISS(IZ)) | |
345 | $ *THX*SINBE*V3I(IZ))/2. | |
346 | END IF | |
347 | WID=SQRT(MWIW**4+MZIZ**4+AMHC**4-2.*(MWIW*MZIZ)**2 | |
348 | $ -2*(MWIW*AMHC)**2-2*(MZIZ*AMHC)**2)/8./PI/MZIZ**3 | |
349 | $ *((A*A+B*B)*(MWIW*MWIW+MZIZ*MZIZ-AMHC*AMHC)/2. | |
350 | $ +(A*A-B*B)*MWIW*MZIZ) | |
351 | CALL SSSAVE(ISZIZ,WID,+ISWIW,-ISHC,0,0,0) | |
352 | CALL SSSAVE(ISZIZ,WID,-ISWIW,+ISHC,0,0,0) | |
353 | ENDIF | |
354 | 205 CONTINUE | |
355 | C | |
356 | C zi --> zj + z | |
357 | C | |
358 | C Note that if m(zi) > m(zj) + m(z), then the z terms are | |
359 | C omitted from the zi -> zj + f + fbar calculation, so there | |
360 | C is no double counting with zi -> zj + z. | |
361 | C | |
362 | DO 210 JZ=1,IZ-1 | |
363 | MZJZ=ABS(AMZISS(JZ)) | |
364 | IF(MZIZ.GT.(MZJZ+AMZ)) THEN | |
365 | SN=-1.*SIGN(1.,AMZISS(IZ))*SIGN(1.,AMZISS(JZ)) | |
366 | WID=WIJ(JZ,IZ)**2/(2.*PI)/(MZIZ)**3*SQRT(MZIZ**4+MZJZ**4 | |
367 | $ +AMZ**4-2.*(MZIZ*MZJZ)**2-2.*(MZIZ*AMZ)**2-2.*(MZJZ*AMZ)**2) | |
368 | $ *((MZIZ**2+MZJZ**2-AMZ**2)/2.+((MZIZ**2-MZJZ**2)**2-AMZ**4)/ | |
369 | $ 2./AMZ**2-3.*SN*MZIZ*MZJZ) | |
370 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDZ,0,0,0) | |
371 | END IF | |
372 | 210 CONTINUE | |
373 | C | |
374 | C zi --> zj + f + fbar | |
375 | C | |
376 | DO 220 JZ=1,IZ-1 | |
377 | MZJZ=ABS(AMZISS(JZ)) | |
378 | SNJZ=SIGN(1.,AMZISS(JZ)) | |
379 | IF (SNJZ.EQ.1.) THEN | |
380 | THJZ=0 | |
381 | ELSE | |
382 | THJZ=1 | |
383 | END IF | |
384 | IF(MZIZ.LT.FUDGE*MZJZ) GO TO 220 | |
385 | FAC=1./2./MZIZ/(2.*PI)**5*PI**2*MZIZ**2 | |
386 | C Leptons -- Z decay allowed, so omit Z | |
387 | TMP(1)=MZIZ | |
388 | TMP(2)=MZJZ | |
389 | TMP(4)=-SNIJ(JZ,IZ) | |
390 | TMP(6)=0. | |
391 | C zi -> zj + u + ubar | |
392 | IF (MZIZ.GT.(MZJZ+2*AMUP)) THEN | |
393 | IF (MZIZ.LT.AMULSS) THEN | |
394 | TMP(3)=AMULSS | |
395 | TERMLL=3*2*AUI(JZ)**2*AUI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
396 | ELSE | |
397 | TERMLL=0. | |
398 | END IF | |
399 | IF (MZIZ.LT.AMURSS) THEN | |
400 | TMP(3)=AMURSS | |
401 | TERMRR=3*2*BUI(JZ)**2*BUI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
402 | ELSE | |
403 | TERMRR=0. | |
404 | END IF | |
405 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
406 | TMP(4)=+SNIJ(JZ,IZ) | |
407 | TERMZZ=3*(APU**2+BTU**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
408 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ) | |
409 | ELSE | |
410 | TERMZZ=0. | |
411 | END IF | |
412 | IF (MZIZ.LT.AMULSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
413 | TMP(4)=+SNIJ(JZ,IZ) | |
414 | TMP(3)=AMULSS | |
415 | TERMLZ=3*8*E*(APU-BTU)*WIJ(JZ,IZ)*AUI(JZ)*AUI(IZ)/MZIZ | |
416 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
417 | ELSE | |
418 | TERMLZ=0. | |
419 | END IF | |
420 | IF (MZIZ.LT.AMURSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
421 | TMP(4)=+SNIJ(JZ,IZ) | |
422 | TMP(3)=AMURSS | |
423 | TERMRZ=-3*8*E*(APU+BTU)*WIJ(JZ,IZ)*BUI(JZ)*BUI(IZ)/MZIZ | |
424 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
425 | ELSE | |
426 | TERMRZ=0. | |
427 | END IF | |
428 | WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ | |
429 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDUP,-IDUP,0,0) | |
430 | C Enter information for decay matrix element | |
431 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
432 | Z1(1)=0. | |
433 | Z1(2)=2*E*WIJ(JZ,IZ) | |
434 | ELSE | |
435 | Z1(1)=2*E*WIJ(JZ,IZ) | |
436 | Z1(2)=0. | |
437 | END IF | |
438 | Z2(1)=APU | |
439 | Z2(2)=BTU | |
440 | CALL SSME3(1,AMZ,Z1,Z2) | |
441 | Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*AUI(IZ)/2. | |
442 | Z1(2)=-Z1(1) | |
443 | Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.) | |
444 | Z2(2)=Z2(1) | |
445 | CALL SSME3(2,AMULSS,Z1,Z2) | |
446 | Z1(1)=ZI**(THIZ-1)*BUI(IZ)/2. | |
447 | Z1(2)=Z1(1) | |
448 | Z2(1)=CONJG(ZI**(THJZ-1)*BUI(JZ)/2.) | |
449 | Z2(2)=-Z2(1) | |
450 | CALL SSME3(2,AMURSS,Z1,Z2) | |
451 | Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2. | |
452 | Z1(2)=-Z1(1) | |
453 | Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*AUI(IZ)/2.) | |
454 | Z2(2)=Z2(1) | |
455 | CALL SSME3(3,AMULSS,Z1,Z2) | |
456 | Z1(1)=ZI**(THJZ-1)*BUI(JZ)/2. | |
457 | Z1(2)=Z1(1) | |
458 | Z2(1)=-CONJG(ZI**(THIZ-1)*BUI(IZ)/2.) | |
459 | Z2(2)=-Z2(1) | |
460 | CALL SSME3(3,AMURSS,Z1,Z2) | |
461 | END IF | |
462 | C zi -> zj + d + dbar | |
463 | TMP(4)=-SNIJ(JZ,IZ) | |
464 | IF (MZIZ.GT.(MZJZ+2*AMDN)) THEN | |
465 | IF (MZIZ.LT.AMDLSS) THEN | |
466 | TMP(3)=AMDLSS | |
467 | TERMLL=3*2*ADI(JZ)**2*ADI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
468 | ELSE | |
469 | TERMLL=0. | |
470 | END IF | |
471 | IF (MZIZ.LT.AMDRSS) THEN | |
472 | TMP(3)=AMDRSS | |
473 | TERMRR=3*2*BDI(JZ)**2*BDI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
474 | ELSE | |
475 | TERMRR=0. | |
476 | END IF | |
477 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
478 | TMP(4)=+SNIJ(JZ,IZ) | |
479 | TERMZZ=3*(APD**2+BTD**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
480 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ) | |
481 | ELSE | |
482 | TERMZZ=0. | |
483 | END IF | |
484 | IF (MZIZ.LT.AMDLSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
485 | TMP(4)=+SNIJ(JZ,IZ) | |
486 | TMP(3)=AMDLSS | |
487 | TERMLZ=3*8*E*(APD-BTD)*WIJ(JZ,IZ)*ADI(JZ)*ADI(IZ)/MZIZ | |
488 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
489 | ELSE | |
490 | TERMLZ=0. | |
491 | END IF | |
492 | IF (MZIZ.LT.AMDRSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
493 | TMP(4)=+SNIJ(JZ,IZ) | |
494 | TMP(3)=AMDRSS | |
495 | TERMRZ=-3*8*E*(APD+BTD)*WIJ(JZ,IZ)*BDI(JZ)*BDI(IZ)/MZIZ | |
496 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
497 | ELSE | |
498 | TERMRZ=0. | |
499 | END IF | |
500 | WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ | |
501 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDDN,-IDDN,0,0) | |
502 | C Enter information for decay matrix element | |
503 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
504 | Z1(1)=0. | |
505 | Z1(2)=2*E*WIJ(JZ,IZ) | |
506 | ELSE | |
507 | Z1(1)=2*E*WIJ(JZ,IZ) | |
508 | Z1(2)=0. | |
509 | END IF | |
510 | Z2(1)=APD | |
511 | Z2(2)=BTD | |
512 | CALL SSME3(1,AMZ,Z1,Z2) | |
513 | Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ADI(IZ)/2. | |
514 | Z1(2)=-Z1(1) | |
515 | Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ)/2.) | |
516 | Z2(2)=Z2(1) | |
517 | CALL SSME3(2,AMDLSS,Z1,Z2) | |
518 | Z1(1)=ZI**(THIZ-1)*BDI(IZ)/2. | |
519 | Z1(2)=Z1(1) | |
520 | Z2(1)=CONJG(ZI**(THJZ-1)*BDI(JZ)/2.) | |
521 | Z2(2)=-Z2(1) | |
522 | CALL SSME3(2,AMDRSS,Z1,Z2) | |
523 | Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ)/2. | |
524 | Z1(2)=-Z1(1) | |
525 | Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ADI(IZ)/2.) | |
526 | Z2(2)=Z2(1) | |
527 | CALL SSME3(3,AMDLSS,Z1,Z2) | |
528 | Z1(1)=ZI**(THJZ-1)*BDI(JZ)/2. | |
529 | Z1(2)=Z1(1) | |
530 | Z2(1)=-CONJG(ZI**(THIZ-1)*BDI(IZ)/2.) | |
531 | Z2(2)=-Z2(1) | |
532 | CALL SSME3(3,AMDRSS,Z1,Z2) | |
533 | END IF | |
534 | C zi -> zj + s + sbar | |
535 | TMP(4)=-SNIJ(JZ,IZ) | |
536 | IF (MZIZ.GT.(MZJZ+2*AMST)) THEN | |
537 | IF (MZIZ.LT.AMSLSS) THEN | |
538 | TMP(3)=AMSLSS | |
539 | TERMLL=3*2*ADI(JZ)**2*ADI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
540 | ELSE | |
541 | TERMLL=0. | |
542 | END IF | |
543 | IF (MZIZ.LT.AMSRSS) THEN | |
544 | TMP(3)=AMSRSS | |
545 | TERMRR=3*2*BDI(JZ)**2*BDI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
546 | ELSE | |
547 | TERMRR=0. | |
548 | END IF | |
549 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
550 | TMP(4)=+SNIJ(JZ,IZ) | |
551 | TERMZZ=3*(APD**2+BTD**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
552 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ) | |
553 | ELSE | |
554 | TERMZZ=0. | |
555 | END IF | |
556 | IF (MZIZ.LT.AMSLSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
557 | TMP(4)=+SNIJ(JZ,IZ) | |
558 | TMP(3)=AMSLSS | |
559 | TERMLZ=3*8*E*(APD-BTD)*WIJ(JZ,IZ)*ADI(JZ)*ADI(IZ)/MZIZ | |
560 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
561 | ELSE | |
562 | TERMLZ=0. | |
563 | END IF | |
564 | IF (MZIZ.LT.AMSRSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
565 | TMP(4)=+SNIJ(JZ,IZ) | |
566 | TMP(3)=AMSRSS | |
567 | TERMRZ=-3*8*E*(APD+BTD)*WIJ(JZ,IZ)*BDI(JZ)*BDI(IZ)/MZIZ | |
568 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
569 | ELSE | |
570 | TERMRZ=0. | |
571 | END IF | |
572 | WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ | |
573 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDST,-IDST,0,0) | |
574 | C Enter information for decay matrix element | |
575 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
576 | Z1(1)=0. | |
577 | Z1(2)=2*E*WIJ(JZ,IZ) | |
578 | ELSE | |
579 | Z1(1)=2*E*WIJ(JZ,IZ) | |
580 | Z1(2)=0. | |
581 | END IF | |
582 | Z2(1)=APD | |
583 | Z2(2)=BTD | |
584 | CALL SSME3(1,AMZ,Z1,Z2) | |
585 | Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ADI(IZ)/2. | |
586 | Z1(2)=-Z1(1) | |
587 | Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ)/2.) | |
588 | Z2(2)=Z2(1) | |
589 | CALL SSME3(2,AMSLSS,Z1,Z2) | |
590 | Z1(1)=ZI**(THIZ-1)*BDI(IZ)/2. | |
591 | Z1(2)=Z1(1) | |
592 | Z2(1)=CONJG(ZI**(THJZ-1)*BDI(JZ)/2.) | |
593 | Z2(2)=-Z2(1) | |
594 | CALL SSME3(2,AMSRSS,Z1,Z2) | |
595 | Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ADI(JZ)/2. | |
596 | Z1(2)=-Z1(1) | |
597 | Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ADI(IZ)/2.) | |
598 | Z2(2)=Z2(1) | |
599 | CALL SSME3(3,AMSLSS,Z1,Z2) | |
600 | Z1(1)=ZI**(THJZ-1)*BDI(JZ)/2. | |
601 | Z1(2)=Z1(1) | |
602 | Z2(1)=-CONJG(ZI**(THIZ-1)*BDI(IZ)/2.) | |
603 | Z2(2)=-Z2(1) | |
604 | CALL SSME3(3,AMSRSS,Z1,Z2) | |
605 | END IF | |
606 | C zi -> zj + c + cbar | |
607 | TMP(4)=-SNIJ(JZ,IZ) | |
608 | IF (MZIZ.GT.(MZJZ+2*AMCH)) THEN | |
609 | IF (MZIZ.LT.AMCLSS) THEN | |
610 | TMP(3)=AMCLSS | |
611 | TERMLL=3*2*AUI(JZ)**2*AUI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
612 | ELSE | |
613 | TERMLL=0. | |
614 | END IF | |
615 | IF (MZIZ.LT.AMCRSS) THEN | |
616 | TMP(3)=AMCRSS | |
617 | TERMRR=3*2*BUI(JZ)**2*BUI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
618 | ELSE | |
619 | TERMRR=0. | |
620 | END IF | |
621 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
622 | TMP(4)=+SNIJ(JZ,IZ) | |
623 | TERMZZ=3*(APU**2+BTU**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
624 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ) | |
625 | ELSE | |
626 | TERMZZ=0. | |
627 | END IF | |
628 | IF (MZIZ.LT.AMCLSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
629 | TMP(4)=+SNIJ(JZ,IZ) | |
630 | TMP(3)=AMCLSS | |
631 | TERMLZ=3*8*E*(APU-BTU)*WIJ(JZ,IZ)*AUI(JZ)*AUI(IZ)/MZIZ | |
632 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
633 | ELSE | |
634 | TERMLZ=0. | |
635 | END IF | |
636 | IF (MZIZ.LT.AMCRSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
637 | TMP(4)=+SNIJ(JZ,IZ) | |
638 | TMP(3)=AMCRSS | |
639 | TERMRZ=-3*8*E*(APU+BTU)*WIJ(JZ,IZ)*BUI(JZ)*BUI(IZ)/MZIZ | |
640 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
641 | ELSE | |
642 | TERMRZ=0. | |
643 | END IF | |
644 | WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ | |
645 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDCH,-IDCH,0,0) | |
646 | C Enter information for decay matrix element | |
647 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
648 | Z1(1)=0. | |
649 | Z1(2)=2*E*WIJ(JZ,IZ) | |
650 | ELSE | |
651 | Z1(1)=2*E*WIJ(JZ,IZ) | |
652 | Z1(2)=0. | |
653 | END IF | |
654 | Z2(1)=APU | |
655 | Z2(2)=BTU | |
656 | CALL SSME3(1,AMZ,Z1,Z2) | |
657 | Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*AUI(IZ)/2. | |
658 | Z1(2)=-Z1(1) | |
659 | Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2.) | |
660 | Z2(2)=Z2(1) | |
661 | CALL SSME3(2,AMCLSS,Z1,Z2) | |
662 | Z1(1)=ZI**(THIZ-1)*BUI(IZ)/2. | |
663 | Z1(2)=Z1(1) | |
664 | Z2(1)=CONJG(ZI**(THJZ-1)*BUI(JZ)/2.) | |
665 | Z2(2)=-Z2(1) | |
666 | CALL SSME3(2,AMCRSS,Z1,Z2) | |
667 | Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*AUI(JZ)/2. | |
668 | Z1(2)=-Z1(1) | |
669 | Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*AUI(IZ)/2.) | |
670 | Z2(2)=Z2(1) | |
671 | CALL SSME3(3,AMCLSS,Z1,Z2) | |
672 | Z1(1)=ZI**(THJZ-1)*BUI(JZ)/2. | |
673 | Z1(2)=Z1(1) | |
674 | Z2(1)=-CONJG(ZI**(THIZ-1)*BUI(IZ)/2.) | |
675 | Z2(2)=-Z2(1) | |
676 | CALL SSME3(3,AMCRSS,Z1,Z2) | |
677 | END IF | |
678 | C zi -> zj + b+ bbar ; mixing/yukawa effects now included! | |
679 | C thanks to M. Drees | |
680 | FACT=1./2./(2*PI)**5/2./MZIZ | |
681 | ALIZ1=ADI(IZ)*COSB-FB*ZMIXSS(2,IZ)*SINB | |
682 | ALJZ1=ADI(JZ)*COSB-FB*ZMIXSS(2,JZ)*SINB | |
683 | ALIZ2=ADI(IZ)*SINB+FB*ZMIXSS(2,IZ)*COSB | |
684 | ALJZ2=ADI(JZ)*SINB+FB*ZMIXSS(2,JZ)*COSB | |
685 | BEIZ1=BDI(IZ)*SINB+FB*ZMIXSS(2,IZ)*COSB | |
686 | BEJZ1=BDI(JZ)*SINB+FB*ZMIXSS(2,JZ)*COSB | |
687 | BEIZ2=-BDI(IZ)*COSB+FB*ZMIXSS(2,IZ)*SINB | |
688 | BEJZ2=-BDI(JZ)*COSB+FB*ZMIXSS(2,JZ)*SINB | |
689 | SGNIJ=-SNIJ(JZ,IZ) | |
690 | XUPPER=(MZIZ**2+AMBT**2-(AMBT+MZJZ)**2)/2./MZIZ | |
691 | IF (MZIZ.GT.(MZJZ+2*AMBT)) THEN | |
692 | TMP(1)=MZIZ | |
693 | TMP(2)=AMBT | |
694 | TMP(3)=MZJZ | |
695 | TMP(4)=AMB1SS | |
696 | TMP(5)=AMB1SS | |
697 | TMP(6)=AMBT | |
698 | IF (MZIZ.LT.AMB1SS) THEN | |
699 | GLLF1=4*ALIZ1**2*((ALJZ1**2+BEJZ1**2)* | |
700 | , SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*ALJZ1**2* | |
701 | , SSXINT(AMBT,SSGX2,XUPPER)) | |
702 | GRRF1=4*BEIZ1**2*((ALJZ1**2+BEJZ1**2)* | |
703 | , SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*BEJZ1**2* | |
704 | , SSXINT(AMBT,SSGX2,XUPPER)) | |
705 | GLRF1=-8*ALIZ1*BEIZ1*ALJZ1*BEJZ1* | |
706 | , SSXINT(AMBT,SSGX8,XUPPER) | |
707 | GF1=GLLF1+GRRF1+GLRF1 | |
708 | ELSE | |
709 | GF1=0. | |
710 | END IF | |
711 | IF (MZIZ.LT.AMB1SS) THEN | |
712 | TMP(4)=AMB1SS | |
713 | TMP(5)=AMB2SS | |
714 | GLGL=8*ALIZ1*ALIZ2*(ALJZ1*ALJZ2+BEJZ1*BEJZ2)* | |
715 | , SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*8*ALIZ1*ALIZ2* | |
716 | , ALJZ1*ALJZ2*SSXINT(AMBT,SSGX2,XUPPER) | |
717 | GRGR=8*BEIZ1*BEIZ2*(ALJZ1*ALJZ2+BEJZ1*BEJZ2)* | |
718 | , SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*8*BEIZ1*BEIZ2* | |
719 | , BEJZ1*BEJZ2*SSXINT(AMBT,SSGX2,XUPPER) | |
720 | GLGR=-8*ALIZ1*BEIZ2*ALJZ2*BEJZ1* | |
721 | , SSXINT(AMBT,SSGX8,XUPPER) | |
722 | GRGL=-8*ALIZ2*BEIZ1*ALJZ1*BEJZ2* | |
723 | , SSXINT(AMBT,SSGX8,XUPPER) | |
724 | GF12=GLGL+GRGR+GLGR+GRGL | |
725 | ELSE | |
726 | GF12=0. | |
727 | END IF | |
728 | IF (MZIZ.LT.AMB2SS) THEN | |
729 | TMP(4)=AMB2SS | |
730 | TMP(5)=AMB2SS | |
731 | GLLF2=4*ALIZ2**2*((ALJZ2**2+BEJZ2**2)* | |
732 | , SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*ALJZ2**2* | |
733 | , SSXINT(AMBT,SSGX2,XUPPER)) | |
734 | GRRF2=4*BEIZ2**2*((ALJZ2**2+BEJZ2**2)* | |
735 | , SSXINT(AMBT,SSGX1,XUPPER)+SGNIJ*BEJZ2**2* | |
736 | , SSXINT(AMBT,SSGX2,XUPPER)) | |
737 | GLRF2=-8*ALIZ2*BEIZ2*ALJZ2*BEJZ2* | |
738 | , SSXINT(AMBT,SSGX8,XUPPER) | |
739 | GF2=GLLF2+GRRF2+GLRF2 | |
740 | ELSE | |
741 | GF2=0. | |
742 | END IF | |
743 | GF=FACT*(GF1+GF2+GF12) | |
744 | TMP(2)=MZJZ | |
745 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
746 | TMP(4)=+SNIJ(JZ,IZ) | |
747 | TERMZZ=(APD**2+BTD**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
748 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ) | |
749 | ELSE | |
750 | TERMZZ=0. | |
751 | END IF | |
752 | IF (MZIZ.LE.(MZJZ+AMHL)) THEN | |
753 | TMP(3)=AMHL | |
754 | TMP(5)=AMHL | |
755 | TMP(4)=-SNIJ(JZ,IZ) | |
756 | TERMHL=G**2/64./PI**3/MZIZ*(MBQ*SINA*(XLIJ(JZ,IZ)+ | |
757 | $ XLIJ(IZ,JZ))/AMW/COSBE)**2* | |
758 | $ SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ) | |
759 | ELSE | |
760 | TERMHL=0. | |
761 | END IF | |
762 | IF (MZIZ.LE.(MZJZ+AMHH)) THEN | |
763 | TMP(3)=AMHH | |
764 | TMP(5)=AMHH | |
765 | TMP(4)=-SNIJ(JZ,IZ) | |
766 | TERMHH=G**2/64./PI**3/MZIZ*(MBQ*COSA*(XHIJ(JZ,IZ)+ | |
767 | $ XHIJ(IZ,JZ))/AMW/COSBE)**2* | |
768 | $ SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ) | |
769 | ELSE | |
770 | TERMHH=0. | |
771 | END IF | |
772 | IF (MZIZ.LE.(MZJZ+AMHH).AND.MZIZ.LE.(MZJZ+AMHL)) THEN | |
773 | TMP(3)=AMHL | |
774 | TMP(5)=AMHH | |
775 | TMP(4)=-SNIJ(JZ,IZ) | |
776 | TERMLH=2*G**2/64./PI**3/MZIZ*(MBQ/AMW/COSBE)**2* | |
777 | $ (COSA*(XHIJ(JZ,IZ)+XHIJ(IZ,JZ))*SINA* | |
778 | $ (XLIJ(JZ,IZ)+XLIJ(IZ,JZ)))* | |
779 | $ SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ) | |
780 | ELSE | |
781 | TERMLH=0. | |
782 | END IF | |
783 | IF (MZIZ.LE.(MZJZ+AMHA)) THEN | |
784 | TMP(3)=AMHA | |
785 | TMP(5)=AMHA | |
786 | TMP(4)=SNIJ(JZ,IZ) | |
787 | TERMHA=G**2*TANB**2/64./PI**3/MZIZ*(MBQ*(XPIJ(JZ,IZ)+ | |
788 | $ XPIJ(IZ,JZ))/AMW)**2* | |
789 | $ SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMBT**2)/2./MZIZ) | |
790 | ELSE | |
791 | TERMHA=0. | |
792 | END IF | |
793 | IF (MZIZ.LT.AMB1SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
794 | TMP(4)=-SGNIJ | |
795 | TMP(3)=AMB1SS | |
796 | TERM1Z=8*E*WIJ(JZ,IZ)*(ALJZ1*ALIZ1*(APD-BTD)- | |
797 | $ BEJZ1*BEIZ1*(APD+BTD))/MZIZ | |
798 | $ /(2*PI)**5*SSXINT(4*AMBT**2,SSZZF3,(MZIZ-MZJZ)**2) | |
799 | ELSE | |
800 | TERM1Z=0. | |
801 | END IF | |
802 | IF (MZIZ.LT.AMB2SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
803 | TMP(4)=-SGNIJ | |
804 | TMP(3)=AMB2SS | |
805 | TERM2Z=8*E*WIJ(JZ,IZ)*(ALJZ2*ALIZ2*(APD-BTD)- | |
806 | $ BEJZ2*BEIZ2*(APD+BTD))/MZIZ | |
807 | $ /(2*PI)**5*SSXINT(4*AMBT**2,SSZZF3,(MZIZ-MZJZ)**2) | |
808 | ELSE | |
809 | TERM2Z=0. | |
810 | END IF | |
811 | IF (MZIZ.LT.AMB1SS.AND.MZIZ.LT.(MZJZ+AMHL)) THEN | |
812 | TMP(3)=AMB1SS | |
813 | TMP(4)=SGNIJ | |
814 | TMP(5)=AMHL | |
815 | TERM1L=2*PI*PI*G*MBQ*SINA*SGNIJ*(XLIJ(IZ,JZ)+ | |
816 | $ XLIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1) | |
817 | $ /MZIZ/AMW/COSBE*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2) | |
818 | TERM1L=FACT*TERM1L | |
819 | ELSE | |
820 | TERM1L=0. | |
821 | END IF | |
822 | IF (MZIZ.LT.AMB2SS.AND.MZIZ.LT.(MZJZ+AMHL)) THEN | |
823 | TMP(3)=AMB2SS | |
824 | TMP(4)=SGNIJ | |
825 | TMP(5)=AMHL | |
826 | TERM2L=2*PI*PI*G*MBQ*SINA*SGNIJ*(XLIJ(IZ,JZ)+ | |
827 | $ XLIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2) | |
828 | $ /MZIZ/AMW/COSBE*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2) | |
829 | TERM2L=FACT*TERM2L | |
830 | ELSE | |
831 | TERM2L=0. | |
832 | END IF | |
833 | IF (MZIZ.LT.AMB1SS.AND.MZIZ.LT.(MZJZ+AMHH)) THEN | |
834 | TMP(3)=AMB1SS | |
835 | TMP(4)=SGNIJ | |
836 | TMP(5)=AMHH | |
837 | TERM1H=2*PI*PI*G*MBQ*COSA*SGNIJ*(XHIJ(IZ,JZ)+ | |
838 | $ XHIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1) | |
839 | $ /MZIZ/AMW/COSBE*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2) | |
840 | TERM1H=FACT*TERM1H | |
841 | ELSE | |
842 | TERM1H=0. | |
843 | END IF | |
844 | IF (MZIZ.LT.AMB2SS.AND.MZIZ.LT.(MZJZ+AMHH)) THEN | |
845 | TMP(3)=AMB2SS | |
846 | TMP(4)=SGNIJ | |
847 | TMP(5)=AMHH | |
848 | TERM2H=2*PI*PI*G*MBQ*COSA*SGNIJ*(XHIJ(IZ,JZ)+ | |
849 | $ XHIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2) | |
850 | $ /MZIZ/AMW/COSBE*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2) | |
851 | TERM2H=FACT*TERM2H | |
852 | ELSE | |
853 | TERM2H=0. | |
854 | END IF | |
855 | IF (MZIZ.LT.AMB1SS.AND.MZIZ.LT.(MZJZ+AMHA)) THEN | |
856 | TMP(3)=AMB1SS | |
857 | TMP(4)=-SGNIJ | |
858 | TMP(5)=AMHA | |
859 | TERM1A=-2*PI*PI*G*MBQ*TANB*SGNIJ*(XPIJ(IZ,JZ)+ | |
860 | $ XPIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1) | |
861 | $ /MZIZ/AMW*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2) | |
862 | TERM1A=FACT*TERM1A | |
863 | ELSE | |
864 | TERM1A=0. | |
865 | END IF | |
866 | IF (MZIZ.LT.AMB2SS.AND.MZIZ.LT.(MZJZ+AMHA)) THEN | |
867 | TMP(3)=AMB2SS | |
868 | TMP(4)=-SGNIJ | |
869 | TMP(5)=AMHA | |
870 | TERM2A=-2*PI*PI*G*MBQ*TANB*SGNIJ*(XPIJ(IZ,JZ)+ | |
871 | $ XPIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2) | |
872 | $ /MZIZ/AMW*SSXINT(4*AMBT**2,SSZZF5,(MZIZ-MZJZ)**2) | |
873 | TERM2A=FACT*TERM2A | |
874 | ELSE | |
875 | TERM2A=0. | |
876 | END IF | |
877 | WID=3*(GF+TERMZZ+TERMHL+TERMHH+TERMLH+TERMHA+TERM1Z+ | |
878 | $ TERM2Z+TERM1L+TERM2L+TERM1H+TERM2H+TERM1A+TERM2A) | |
879 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDBT,-IDBT,0,0) | |
880 | C Enter information for decay matrix element | |
881 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
882 | Z1(1)=0. | |
883 | Z1(2)=2*E*WIJ(JZ,IZ) | |
884 | ELSE | |
885 | Z1(1)=2*E*WIJ(JZ,IZ) | |
886 | Z1(2)=0. | |
887 | END IF | |
888 | Z2(1)=APD | |
889 | Z2(2)=BTD | |
890 | CALL SSME3(1,AMZ,Z1,Z2) | |
891 | Z1(1)=((ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)- | |
892 | $ FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-(ZI**THIZ*BDI(IZ)- | |
893 | $ FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2. | |
894 | Z1(2)=((-ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)- | |
895 | $ FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-(ZI**THIZ*BDI(IZ)+ | |
896 | $ FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2. | |
897 | Z2(1)=CONJG(((ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)- | |
898 | $ FB*ZMIXSS(2,JZ)*ZI**THJZ)*COSB-(ZI**THJZ*BDI(JZ)- | |
899 | $ FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINB)/2.) | |
900 | Z2(2)=-CONJG(((-ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)- | |
901 | $ FB*ZMIXSS(2,JZ)*ZI**THJZ)*COSB-(ZI**THJZ*BDI(JZ)+ | |
902 | $ FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINB)/2.) | |
903 | CALL SSME3(2,AMB1SS,Z1,Z2) | |
904 | Z1(1)=((ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)- | |
905 | $ FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+(ZI**THIZ*BDI(IZ)- | |
906 | $ FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2. | |
907 | Z1(2)=((-ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)- | |
908 | $ FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+(ZI**THIZ*BDI(IZ)+ | |
909 | $ FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2. | |
910 | Z2(1)=CONJG(((ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)- | |
911 | $ FB*ZMIXSS(2,JZ)*ZI**THJZ)*SINB+(ZI**THJZ*BDI(JZ)- | |
912 | $ FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSB)/2.) | |
913 | Z2(2)=-CONJG(((-ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)- | |
914 | $ FB*ZMIXSS(2,JZ)*ZI**THJZ)*SINB+(ZI**THJZ*BDI(JZ)+ | |
915 | $ FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSB)/2.) | |
916 | CALL SSME3(2,AMB2SS,Z1,Z2) | |
917 | Z1(1)=((ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)- | |
918 | $ FB*ZMIXSS(2,JZ)*ZI**THJZ)*COSB-(ZI**THJZ*BDI(JZ)- | |
919 | $ FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINB)/2. | |
920 | Z1(2)=((-ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)- | |
921 | $ FB*ZMIXSS(2,JZ)*ZI**THJZ)*COSB-(ZI**THJZ*BDI(JZ)+ | |
922 | $ FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINB)/2. | |
923 | Z2(1)=-CONJG(((ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)- | |
924 | $ FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-(ZI**THIZ*BDI(IZ)- | |
925 | $ FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.) | |
926 | Z2(2)=CONJG(((-ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)- | |
927 | $ FB*ZMIXSS(2,IZ)*ZI**THIZ)*COSB-(ZI**THIZ*BDI(IZ)+ | |
928 | $ FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINB)/2.) | |
929 | CALL SSME3(3,AMB1SS,Z1,Z2) | |
930 | Z1(1)=((ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)- | |
931 | $ FB*ZMIXSS(2,JZ)*ZI**THJZ)*SINB+(ZI**THJZ*BDI(JZ)- | |
932 | $ FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSB)/2. | |
933 | Z1(2)=((-ZI**(THJZ)*(-1.)**(THJZ+1)*ADI(JZ)- | |
934 | $ FB*ZMIXSS(2,JZ)*ZI**THJZ)*SINB+(ZI**THJZ*BDI(JZ)+ | |
935 | $ FB*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSB)/2. | |
936 | Z2(1)=-CONJG(((ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)- | |
937 | $ FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+(ZI**THIZ*BDI(IZ)- | |
938 | $ FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.) | |
939 | Z2(2)=CONJG(((-ZI**(THIZ)*(-1.)**(THIZ+1)*ADI(IZ)- | |
940 | $ FB*ZMIXSS(2,IZ)*ZI**THIZ)*SINB+(ZI**THIZ*BDI(IZ)+ | |
941 | $ FB*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSB)/2.) | |
942 | CALL SSME3(3,AMB2SS,Z1,Z2) | |
943 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
944 | Z1(1)=(-ZI)**(THIZ+THJZ)*(XLIJ(IZ,JZ)+XLIJ(JZ,IZ)) | |
945 | Z1(2)=0. | |
946 | ELSE | |
947 | Z1(1)=0. | |
948 | Z1(2)=(-ZI)**(THIZ+THJZ)*(XLIJ(IZ,JZ)+XLIJ(JZ,IZ)) | |
949 | END IF | |
950 | Z2(1)=-G*MBQ*SINA/2./AMW/COSBE | |
951 | Z2(2)=0. | |
952 | CALL SSME3(4,AMHL,Z1,Z2) | |
953 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
954 | Z1(1)=(-ZI)**(THIZ+THJZ)*(XHIJ(IZ,JZ)+XHIJ(JZ,IZ)) | |
955 | Z1(2)=0. | |
956 | ELSE | |
957 | Z1(1)=0. | |
958 | Z1(2)=(-ZI)**(THIZ+THJZ)*(XHIJ(IZ,JZ)+XHIJ(JZ,IZ)) | |
959 | END IF | |
960 | Z2(1)=-G*MBQ*COSA/2./AMW/COSBE | |
961 | Z2(2)=0. | |
962 | CALL SSME3(4,AMHH,Z1,Z2) | |
963 | IF (-SNIJ(JZ,IZ).GT.0.) THEN | |
964 | Z1(1)=0. | |
965 | Z1(2)=(-ZI)**(THIZ+THJZ+1)*(XPIJ(IZ,JZ)+XPIJ(JZ,IZ)) | |
966 | ELSE | |
967 | Z1(1)=(-ZI)**(THIZ+THJZ+1)*(XPIJ(IZ,JZ)+XPIJ(JZ,IZ)) | |
968 | Z1(2)=0. | |
969 | END IF | |
970 | Z2(1)=0. | |
971 | Z2(2)=ZI*G*MBQ*TANB/2./AMW | |
972 | CALL SSME3(4,AMHA,Z1,Z2) | |
973 | END IF | |
974 | C zi -> zj + e + ebar | |
975 | TMP(2)=MZJZ | |
976 | TMP(4)=-SNIJ(JZ,IZ) | |
977 | TMP(6)=0. | |
978 | IF (MZIZ.GT.(MZJZ+2*AME)) THEN | |
979 | IF (MZIZ.LT.AMELSS) THEN | |
980 | TMP(3)=AMELSS | |
981 | TERMLL=2*ALI(JZ)**2*ALI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
982 | ELSE | |
983 | TERMLL=0. | |
984 | END IF | |
985 | IF (MZIZ.LT.AMERSS) THEN | |
986 | TMP(3)=AMERSS | |
987 | TERMRR=2*BLI(JZ)**2*BLI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
988 | ELSE | |
989 | TERMRR=0. | |
990 | END IF | |
991 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
992 | TMP(4)=+SNIJ(JZ,IZ) | |
993 | TERMZZ=(APL**2+BTL**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
994 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ) | |
995 | ELSE | |
996 | TERMZZ=0. | |
997 | END IF | |
998 | IF (MZIZ.LT.AMELSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
999 | TMP(4)=+SNIJ(JZ,IZ) | |
1000 | TMP(3)=AMELSS | |
1001 | TERMLZ=8*E*(APL-BTL)*WIJ(JZ,IZ)*ALI(JZ)*ALI(IZ)/MZIZ | |
1002 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
1003 | ELSE | |
1004 | TERMLZ=0. | |
1005 | END IF | |
1006 | IF (MZIZ.LT.AMERSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
1007 | TMP(4)=+SNIJ(JZ,IZ) | |
1008 | TMP(3)=AMERSS | |
1009 | TERMRZ=-8*E*(APL+BTL)*WIJ(JZ,IZ)*BLI(JZ)*BLI(IZ)/MZIZ | |
1010 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
1011 | ELSE | |
1012 | TERMRZ=0. | |
1013 | END IF | |
1014 | WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ | |
1015 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDE,-IDE,0,0) | |
1016 | C Enter information for decay matrix element | |
1017 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
1018 | Z1(1)=0. | |
1019 | Z1(2)=2*E*WIJ(JZ,IZ) | |
1020 | ELSE | |
1021 | Z1(1)=2*E*WIJ(JZ,IZ) | |
1022 | Z1(2)=0. | |
1023 | END IF | |
1024 | Z2(1)=APL | |
1025 | Z2(2)=BTL | |
1026 | CALL SSME3(1,AMZ,Z1,Z2) | |
1027 | Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ALI(IZ)/2. | |
1028 | Z1(2)=-Z1(1) | |
1029 | Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ)/2.) | |
1030 | Z2(2)=Z2(1) | |
1031 | CALL SSME3(2,AMELSS,Z1,Z2) | |
1032 | Z1(1)=ZI**(THIZ-1)*BLI(IZ)/2. | |
1033 | Z1(2)=Z1(1) | |
1034 | Z2(1)=CONJG(ZI**(THJZ-1)*BLI(JZ)/2.) | |
1035 | Z2(2)=-Z2(1) | |
1036 | CALL SSME3(2,AMERSS,Z1,Z2) | |
1037 | Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ)/2. | |
1038 | Z1(2)=-Z1(1) | |
1039 | Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ALI(IZ)/2.) | |
1040 | Z2(2)=Z2(1) | |
1041 | CALL SSME3(3,AMELSS,Z1,Z2) | |
1042 | Z1(1)=ZI**(THJZ-1)*BLI(JZ)/2. | |
1043 | Z1(2)=Z1(1) | |
1044 | Z2(1)=-CONJG(ZI**(THIZ-1)*BLI(IZ)/2.) | |
1045 | Z2(2)=-Z2(1) | |
1046 | CALL SSME3(3,AMERSS,Z1,Z2) | |
1047 | END IF | |
1048 | C zi -> zj + mu + mubar | |
1049 | TMP(4)=-SNIJ(JZ,IZ) | |
1050 | IF (MZIZ.GT.(MZJZ+2*AMMU)) THEN | |
1051 | IF (MZIZ.LT.AMMLSS) THEN | |
1052 | TMP(3)=AMMLSS | |
1053 | TERMLL=2*ALI(JZ)**2*ALI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
1054 | ELSE | |
1055 | TERMLL=0. | |
1056 | END IF | |
1057 | IF (MZIZ.LT.AMMRSS) THEN | |
1058 | TMP(3)=AMMRSS | |
1059 | TERMRR=2*BLI(JZ)**2*BLI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
1060 | ELSE | |
1061 | TERMRR=0. | |
1062 | END IF | |
1063 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
1064 | TMP(4)=+SNIJ(JZ,IZ) | |
1065 | TERMZZ=(APL**2+BTL**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
1066 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ) | |
1067 | ELSE | |
1068 | TERMZZ=0. | |
1069 | END IF | |
1070 | IF (MZIZ.LT.AMMLSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
1071 | TMP(4)=+SNIJ(JZ,IZ) | |
1072 | TMP(3)=AMMLSS | |
1073 | TERMLZ=8*E*(APL-BTL)*WIJ(JZ,IZ)*ALI(JZ)*ALI(IZ)/MZIZ | |
1074 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
1075 | ELSE | |
1076 | TERMLZ=0. | |
1077 | END IF | |
1078 | IF (MZIZ.LT.AMMRSS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
1079 | TMP(4)=+SNIJ(JZ,IZ) | |
1080 | TMP(3)=AMMRSS | |
1081 | TERMRZ=-8*E*(APL+BTL)*WIJ(JZ,IZ)*BLI(JZ)*BLI(IZ)/MZIZ | |
1082 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
1083 | ELSE | |
1084 | TERMRZ=0. | |
1085 | END IF | |
1086 | WID=TERMLL+TERMRR+TERMZZ+TERMLZ+TERMRZ | |
1087 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDMU,-IDMU,0,0) | |
1088 | C Enter information for decay matrix element | |
1089 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
1090 | Z1(1)=0. | |
1091 | Z1(2)=2*E*WIJ(JZ,IZ) | |
1092 | ELSE | |
1093 | Z1(1)=2*E*WIJ(JZ,IZ) | |
1094 | Z1(2)=0. | |
1095 | END IF | |
1096 | Z2(1)=APL | |
1097 | Z2(2)=BTL | |
1098 | CALL SSME3(1,AMZ,Z1,Z2) | |
1099 | Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ALI(IZ)/2. | |
1100 | Z1(2)=-Z1(1) | |
1101 | Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ)/2.) | |
1102 | Z2(2)=Z2(1) | |
1103 | CALL SSME3(2,AMMLSS,Z1,Z2) | |
1104 | Z1(1)=ZI**(THIZ-1)*BLI(IZ)/2. | |
1105 | Z1(2)=Z1(1) | |
1106 | Z2(1)=CONJG(ZI**(THJZ-1)*BLI(JZ)/2.) | |
1107 | Z2(2)=-Z2(1) | |
1108 | CALL SSME3(2,AMMRSS,Z1,Z2) | |
1109 | Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ALI(JZ)/2. | |
1110 | Z1(2)=-Z1(1) | |
1111 | Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ALI(IZ)/2.) | |
1112 | Z2(2)=Z2(1) | |
1113 | CALL SSME3(3,AMMLSS,Z1,Z2) | |
1114 | Z1(1)=ZI**(THJZ-1)*BLI(JZ)/2. | |
1115 | Z1(2)=Z1(1) | |
1116 | Z2(1)=-CONJG(ZI**(THIZ-1)*BLI(IZ)/2.) | |
1117 | Z2(2)=-Z2(1) | |
1118 | CALL SSME3(3,AMMRSS,Z1,Z2) | |
1119 | END IF | |
1120 | C zi -> zj + tau + taubar. | |
1121 | C Mixing/yukawa effects now included thanks to M. Drees | |
1122 | ALIZ1=-ALI(IZ)*COSL-FL*ZMIXSS(2,IZ)*SINL | |
1123 | ALJZ1=-ALI(JZ)*COSL-FL*ZMIXSS(2,JZ)*SINL | |
1124 | ALIZ2=-ALI(IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL | |
1125 | ALJZ2=-ALI(JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL | |
1126 | BEIZ1=BLI(IZ)*SINL+FL*ZMIXSS(2,IZ)*COSL | |
1127 | BEJZ1=BLI(JZ)*SINL+FL*ZMIXSS(2,JZ)*COSL | |
1128 | BEIZ2=-BLI(IZ)*COSL+FL*ZMIXSS(2,IZ)*SINL | |
1129 | BEJZ2=-BLI(JZ)*COSL+FL*ZMIXSS(2,JZ)*SINL | |
1130 | SGNIJ=-SNIJ(JZ,IZ) | |
1131 | XUPPER=(MZIZ**2+AMTAU**2-(AMTAU+MZJZ)**2)/2./MZIZ | |
1132 | C Polarization for stau_i -> z2ss+tau, z3ss+tau, z4ss+tau. | |
1133 | C See below for z1ss+tau. | |
1134 | IF(JZ.EQ.1) THEN | |
1135 | PTAU1(IZ)=(BEIZ1**2-ALIZ1**2)/(BEIZ1**2+ALIZ1**2) | |
1136 | PTAU2(IZ)=(BEIZ2**2-ALIZ2**2)/(BEIZ2**2+ALIZ2**2) | |
1137 | ENDIF | |
1138 | IF (MZIZ.GT.(MZJZ+2*AMTAU)) THEN | |
1139 | TMP(1)=MZIZ | |
1140 | TMP(2)=AMTAU | |
1141 | TMP(3)=MZJZ | |
1142 | TMP(4)=AML1SS | |
1143 | TMP(5)=AML1SS | |
1144 | TMP(6)=AMTAU | |
1145 | IF (MZIZ.LT.AML1SS) THEN | |
1146 | GLLF1=4*ALIZ1**2*((ALJZ1**2+BEJZ1**2)* | |
1147 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*ALJZ1**2* | |
1148 | , SSXINT(AMTAU,SSGX2,XUPPER)) | |
1149 | GLLF1L=FACT*4*ALIZ1**2*(ALJZ1**2* | |
1150 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*ALJZ1**2* | |
1151 | , SSXINT(AMTAU,SSGX2,XUPPER)) | |
1152 | GRRF1=4*BEIZ1**2*((ALJZ1**2+BEJZ1**2)* | |
1153 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*BEJZ1**2* | |
1154 | , SSXINT(AMTAU,SSGX2,XUPPER)) | |
1155 | GRRF1R=FACT*4*BEIZ1**2*(BEJZ1**2* | |
1156 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*BEJZ1**2* | |
1157 | , SSXINT(AMTAU,SSGX2,XUPPER)) | |
1158 | GLRF1=-8*ALIZ1*BEIZ1*ALJZ1*BEJZ1* | |
1159 | , SSXINT(AMTAU,SSGX8,XUPPER) | |
1160 | GF1=GLLF1+GRRF1+GLRF1 | |
1161 | ELSE | |
1162 | GLLF1L=0. | |
1163 | GRRF1R=0. | |
1164 | GF1=0. | |
1165 | END IF | |
1166 | IF (MZIZ.LT.AML1SS) THEN | |
1167 | TMP(4)=AML1SS | |
1168 | TMP(5)=AML2SS | |
1169 | GLGL=8*ALIZ1*ALIZ2*(ALJZ1*ALJZ2+BEJZ1*BEJZ2)* | |
1170 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*8*ALIZ1*ALIZ2* | |
1171 | , ALJZ1*ALJZ2*SSXINT(AMTAU,SSGX2,XUPPER) | |
1172 | GLGLL=FACT*(8*ALIZ1*ALIZ2*ALJZ1*ALJZ2* | |
1173 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*8*ALIZ1*ALIZ2* | |
1174 | , ALJZ1*ALJZ2*SSXINT(AMTAU,SSGX2,XUPPER)) | |
1175 | GRGR=8*BEIZ1*BEIZ2*(ALJZ1*ALJZ2+BEJZ1*BEJZ2)* | |
1176 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*8*BEIZ1*BEIZ2* | |
1177 | , BEJZ1*BEJZ2*SSXINT(AMTAU,SSGX2,XUPPER) | |
1178 | GRGRR=FACT*(8*BEIZ1*BEIZ2*BEJZ1*BEJZ2* | |
1179 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*8*BEIZ1*BEIZ2* | |
1180 | , BEJZ1*BEJZ2*SSXINT(AMTAU,SSGX2,XUPPER)) | |
1181 | GLGR=-8*ALIZ1*BEIZ2*ALJZ2*BEJZ1* | |
1182 | , SSXINT(AMTAU,SSGX8,XUPPER) | |
1183 | GRGL=-8*ALIZ2*BEIZ1*ALJZ1*BEJZ2* | |
1184 | , SSXINT(AMTAU,SSGX8,XUPPER) | |
1185 | GF12=GLGL+GRGR+GLGR+GRGL | |
1186 | ELSE | |
1187 | GLGLL=0. | |
1188 | GRGRR=0. | |
1189 | GF12=0. | |
1190 | END IF | |
1191 | IF (MZIZ.LT.AML2SS) THEN | |
1192 | TMP(4)=AML2SS | |
1193 | TMP(5)=AML2SS | |
1194 | GLLF2=4*ALIZ2**2*((ALJZ2**2+BEJZ2**2)* | |
1195 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*ALJZ2**2* | |
1196 | , SSXINT(AMTAU,SSGX2,XUPPER)) | |
1197 | GLLF2L=FACT*4*ALIZ2**2*(ALJZ2**2* | |
1198 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*ALJZ2**2* | |
1199 | , SSXINT(AMTAU,SSGX2,XUPPER)) | |
1200 | GRRF2=4*BEIZ2**2*((ALJZ2**2+BEJZ2**2)* | |
1201 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*BEJZ2**2* | |
1202 | , SSXINT(AMTAU,SSGX2,XUPPER)) | |
1203 | GRRF2R=FACT*4*BEIZ2**2*(BEJZ2**2* | |
1204 | , SSXINT(AMTAU,SSGX1,XUPPER)+SGNIJ*BEJZ2**2* | |
1205 | , SSXINT(AMTAU,SSGX2,XUPPER)) | |
1206 | GLRF2=-8*ALIZ2*BEIZ2*ALJZ2*BEJZ2* | |
1207 | , SSXINT(AMTAU,SSGX8,XUPPER) | |
1208 | GF2=GLLF2+GRRF2+GLRF2 | |
1209 | ELSE | |
1210 | GLLF2L=0. | |
1211 | GRRF2R=0. | |
1212 | GF2=0. | |
1213 | END IF | |
1214 | GF=FACT*(GF1+GF2+GF12) | |
1215 | TMP(2)=MZJZ | |
1216 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
1217 | TMP(4)=+SNIJ(JZ,IZ) | |
1218 | TERMZZ=(APL**2+BTL**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
1219 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ) | |
1220 | TMZZRL=TERMZZ*(APL+BTL)**2/2./(APL**2+BTL**2) | |
1221 | TMZZLR=TERMZZ*(APL-BTL)**2/2./(APL**2+BTL**2) | |
1222 | ELSE | |
1223 | TERMZZ=0. | |
1224 | TMZZRL=0. | |
1225 | TMZZLR=0. | |
1226 | END IF | |
1227 | IF (MZIZ.LE.(MZJZ+AMHL)) THEN | |
1228 | TMP(3)=AMHL | |
1229 | TMP(5)=AMHL | |
1230 | TMP(4)=-SNIJ(JZ,IZ) | |
1231 | TERMHL=G**2/64./PI**3/MZIZ*(MTAMZ*SINA*(XLIJ(JZ,IZ)+ | |
1232 | $ XLIJ(IZ,JZ))/AMW/COSBE)**2* | |
1233 | $ SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ) | |
1234 | ELSE | |
1235 | TERMHL=0. | |
1236 | END IF | |
1237 | IF (MZIZ.LE.(MZJZ+AMHH)) THEN | |
1238 | TMP(3)=AMHH | |
1239 | TMP(5)=AMHH | |
1240 | TMP(4)=-SNIJ(JZ,IZ) | |
1241 | TERMHH=G**2/64./PI**3/MZIZ*(MTAMZ*COSA*(XHIJ(JZ,IZ)+ | |
1242 | $ XHIJ(IZ,JZ))/AMW/COSBE)**2* | |
1243 | $ SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ) | |
1244 | ELSE | |
1245 | TERMHH=0. | |
1246 | END IF | |
1247 | IF (MZIZ.LE.(MZJZ+AMHH).AND.MZIZ.LE.(MZJZ+AMHL)) THEN | |
1248 | TMP(3)=AMHL | |
1249 | TMP(5)=AMHH | |
1250 | TMP(4)=-SNIJ(JZ,IZ) | |
1251 | TERMLH=2*G**2/64./PI**3/MZIZ*(MTAMZ/AMW/COSBE)**2* | |
1252 | $ (COSA*(XHIJ(JZ,IZ)+XHIJ(IZ,JZ))*SINA* | |
1253 | $ (XLIJ(JZ,IZ)+XLIJ(IZ,JZ)))* | |
1254 | $ SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ) | |
1255 | ELSE | |
1256 | TERMLH=0. | |
1257 | END IF | |
1258 | IF (MZIZ.LE.(MZJZ+AMHA)) THEN | |
1259 | TMP(3)=AMHA | |
1260 | TMP(5)=AMHA | |
1261 | TMP(4)=SNIJ(JZ,IZ) | |
1262 | TERMHA=G**2*TANB**2/64./PI**3/MZIZ*(MTAMZ*(XPIJ(JZ,IZ)+ | |
1263 | $ XPIJ(IZ,JZ))/AMW)**2* | |
1264 | $ SSXINT(MZJZ,SSZZF4,(MZIZ**2+MZJZ**2-4*AMTAU**2)/2./MZIZ) | |
1265 | ELSE | |
1266 | TERMHA=0. | |
1267 | END IF | |
1268 | IF (MZIZ.LT.AML1SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
1269 | TMP(4)=-SGNIJ | |
1270 | TMP(3)=AML1SS | |
1271 | TERM1Z=8*E*WIJ(JZ,IZ)*(ALJZ1*ALIZ1*(APL-BTL)- | |
1272 | $ BEJZ1*BEIZ1*(APL+BTL))/MZIZ | |
1273 | $ /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2) | |
1274 | TM1ZRL=-8*E*WIJ(JZ,IZ)*BEJZ1*BEIZ1*(APL+BTL)/MZIZ | |
1275 | $ /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2) | |
1276 | TM1ZLR=8*E*WIJ(JZ,IZ)*ALJZ1*ALIZ1*(APL-BTL)/MZIZ | |
1277 | $ /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2) | |
1278 | ELSE | |
1279 | TERM1Z=0. | |
1280 | TM1ZRL=0. | |
1281 | TM1ZLR=0. | |
1282 | END IF | |
1283 | IF (MZIZ.LT.AML2SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
1284 | TMP(4)=-SGNIJ | |
1285 | TMP(3)=AML2SS | |
1286 | TERM2Z=8*E*WIJ(JZ,IZ)*(ALJZ2*ALIZ2*(APL-BTL)- | |
1287 | $ BEJZ2*BEIZ2*(APL+BTL))/MZIZ | |
1288 | $ /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2) | |
1289 | TM2ZRL=-8*E*WIJ(JZ,IZ)*BEJZ2*BEIZ2*(APL+BTL)/MZIZ | |
1290 | $ /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2) | |
1291 | TM2ZLR=8*E*WIJ(JZ,IZ)*ALJZ2*ALIZ2*(APL-BTL)/MZIZ | |
1292 | $ /(2*PI)**5*SSXINT(4*AMTAU**2,SSZZF3,(MZIZ-MZJZ)**2) | |
1293 | ELSE | |
1294 | TERM2Z=0. | |
1295 | TM2ZRL=0. | |
1296 | TM2ZLR=0. | |
1297 | END IF | |
1298 | IF (MZIZ.LT.AML1SS.AND.MZIZ.LT.(MZJZ+AMHL)) THEN | |
1299 | TMP(3)=AML1SS | |
1300 | TMP(4)=SGNIJ | |
1301 | TMP(5)=AMHL | |
1302 | TERM1L=2*PI*PI*G*MTAMZ*SINA*SGNIJ*(XLIJ(IZ,JZ)+ | |
1303 | $ XLIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1) | |
1304 | $ /MZIZ/AMW/COSBE*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2) | |
1305 | TERM1L=FACT*TERM1L | |
1306 | ELSE | |
1307 | TERM1L=0. | |
1308 | END IF | |
1309 | IF (MZIZ.LT.AML2SS.AND.MZIZ.LT.(MZJZ+AMHL)) THEN | |
1310 | TMP(3)=AML2SS | |
1311 | TMP(4)=SGNIJ | |
1312 | TMP(5)=AMHL | |
1313 | TERM2L=2*PI*PI*G*MTAMZ*SINA*SGNIJ*(XLIJ(IZ,JZ)+ | |
1314 | $ XLIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2) | |
1315 | $ /MZIZ/AMW/COSBE*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2) | |
1316 | TERM2L=FACT*TERM2L | |
1317 | ELSE | |
1318 | TERM2L=0. | |
1319 | END IF | |
1320 | IF (MZIZ.LT.AML1SS.AND.MZIZ.LT.(MZJZ+AMHH)) THEN | |
1321 | TMP(3)=AML1SS | |
1322 | TMP(4)=SGNIJ | |
1323 | TMP(5)=AMHH | |
1324 | TERM1H=2*PI*PI*G*MTAMZ*COSA*SGNIJ*(XHIJ(IZ,JZ)+ | |
1325 | $ XHIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1) | |
1326 | $ /MZIZ/AMW/COSBE*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2) | |
1327 | TERM1H=FACT*TERM1H | |
1328 | ELSE | |
1329 | TERM1H=0. | |
1330 | END IF | |
1331 | IF (MZIZ.LT.AML2SS.AND.MZIZ.LT.(MZJZ+AMHH)) THEN | |
1332 | TMP(3)=AML2SS | |
1333 | TMP(4)=SGNIJ | |
1334 | TMP(5)=AMHH | |
1335 | TERM2H=2*PI*PI*G*MTAMZ*COSA*SGNIJ*(XHIJ(IZ,JZ)+ | |
1336 | $ XHIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2) | |
1337 | $ /MZIZ/AMW/COSBE*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2) | |
1338 | TERM2H=FACT*TERM2H | |
1339 | ELSE | |
1340 | TERM2H=0. | |
1341 | END IF | |
1342 | IF (MZIZ.LT.AML1SS.AND.MZIZ.LT.(MZJZ+AMHA)) THEN | |
1343 | TMP(3)=AML1SS | |
1344 | TMP(4)=-SGNIJ | |
1345 | TMP(5)=AMHA | |
1346 | TERM1A=-2*PI*PI*G*MTAMZ*TANB*SGNIJ*(XPIJ(IZ,JZ)+ | |
1347 | $ XPIJ(JZ,IZ))*(ALIZ1*BEJZ1+ALJZ1*BEIZ1) | |
1348 | $ /MZIZ/AMW*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2) | |
1349 | TERM1A=FACT*TERM1A | |
1350 | ELSE | |
1351 | TERM1A=0. | |
1352 | END IF | |
1353 | IF (MZIZ.LT.AML2SS.AND.MZIZ.LT.(MZJZ+AMHA)) THEN | |
1354 | TMP(3)=AML2SS | |
1355 | TMP(4)=-SGNIJ | |
1356 | TMP(5)=AMHA | |
1357 | TERM2A=-2*PI*PI*G*MTAMZ*TANB*SGNIJ*(XPIJ(IZ,JZ)+ | |
1358 | $ XPIJ(JZ,IZ))*(ALIZ2*BEJZ2+ALJZ2*BEIZ2) | |
1359 | $ /MZIZ/AMW*SSXINT(4*AMTAU**2,SSZZF5,(MZIZ-MZJZ)**2) | |
1360 | TERM2A=FACT*TERM2A | |
1361 | ELSE | |
1362 | TERM2A=0. | |
1363 | END IF | |
1364 | WID=GF+TERMZZ+TERMHL+TERMHH+TERMLH+TERMHA+TERM1Z+ | |
1365 | $ TERM2Z+TERM1L+TERM2L+TERM1H+TERM2H+TERM1A+TERM2A | |
1366 | C tau polarization for 3-body z2 -> z1 tau tau | |
1367 | IF (IZ.EQ.2.AND.JZ.EQ.1.AND.WID.GT.0.) THEN | |
1368 | PTAUZZ=(GRRF1R+GRGRR+GRRF2R+TMZZRL+TM1ZRL+TM2ZRL- | |
1369 | $ (GLLF1L+GLGLL+GLLF2L+TMZZLR+TM1ZLR+TM2ZLR)) | |
1370 | $ /WID | |
1371 | END IF | |
1372 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDTAU,-IDTAU,0,0) | |
1373 | C Enter information for decay matrix element | |
1374 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
1375 | Z1(1)=0. | |
1376 | Z1(2)=2*E*WIJ(JZ,IZ) | |
1377 | ELSE | |
1378 | Z1(1)=2*E*WIJ(JZ,IZ) | |
1379 | Z1(2)=0. | |
1380 | END IF | |
1381 | Z2(1)=APL | |
1382 | Z2(2)=BTL | |
1383 | CALL SSME3(1,AMZ,Z1,Z2) | |
1384 | Z1(1)=((ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)- | |
1385 | $ FL*ZMIXSS(2,IZ)*ZI**THIZ)*COSL-(ZI**THIZ*BLI(IZ)- | |
1386 | $ FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINL)/2. | |
1387 | Z1(2)=((-ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)- | |
1388 | $ FL*ZMIXSS(2,IZ)*ZI**THIZ)*COSL-(ZI**THIZ*BLI(IZ)+ | |
1389 | $ FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINL)/2. | |
1390 | Z2(1)=CONJG(((ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)- | |
1391 | $ FL*ZMIXSS(2,JZ)*ZI**THJZ)*COSL-(ZI**THJZ*BLI(JZ)- | |
1392 | $ FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINL)/2.) | |
1393 | Z2(2)=-CONJG(((-ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)- | |
1394 | $ FL*ZMIXSS(2,JZ)*ZI**THJZ)*COSL-(ZI**THJZ*BLI(JZ)+ | |
1395 | $ FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINL)/2.) | |
1396 | CALL SSME3(2,AML1SS,Z1,Z2) | |
1397 | Z1(1)=((ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)- | |
1398 | $ FL*ZMIXSS(2,IZ)*ZI**THIZ)*SINL+(ZI**THIZ*BLI(IZ)- | |
1399 | $ FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSL)/2. | |
1400 | Z1(2)=((-ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)- | |
1401 | $ FL*ZMIXSS(2,IZ)*ZI**THIZ)*SINL+(ZI**THIZ*BLI(IZ)+ | |
1402 | $ FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSL)/2. | |
1403 | Z2(1)=CONJG(((ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)- | |
1404 | $ FL*ZMIXSS(2,JZ)*ZI**THJZ)*SINL+(ZI**THJZ*BLI(JZ)- | |
1405 | $ FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSL)/2.) | |
1406 | Z2(2)=-CONJG(((-ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)- | |
1407 | $ FL*ZMIXSS(2,JZ)*ZI**THJZ)*SINL+(ZI**THJZ*BLI(JZ)+ | |
1408 | $ FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSL)/2.) | |
1409 | CALL SSME3(2,AML2SS,Z1,Z2) | |
1410 | Z1(1)=((ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)- | |
1411 | $ FL*ZMIXSS(2,JZ)*ZI**THJZ)*COSL-(ZI**THJZ*BLI(JZ)- | |
1412 | $ FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINL)/2. | |
1413 | Z1(2)=((-ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)- | |
1414 | $ FL*ZMIXSS(2,JZ)*ZI**THJZ)*COSL-(ZI**THJZ*BLI(JZ)+ | |
1415 | $ FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*SINL)/2. | |
1416 | Z2(1)=-CONJG(((ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)- | |
1417 | $ FL*ZMIXSS(2,IZ)*ZI**THIZ)*COSL-(ZI**THIZ*BLI(IZ)- | |
1418 | $ FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINL)/2.) | |
1419 | Z2(2)=CONJG(((-ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)- | |
1420 | $ FL*ZMIXSS(2,IZ)*ZI**THIZ)*COSL-(ZI**THIZ*BLI(IZ)+ | |
1421 | $ FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*SINL)/2.) | |
1422 | CALL SSME3(3,AML1SS,Z1,Z2) | |
1423 | Z1(1)=((ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)- | |
1424 | $ FL*ZMIXSS(2,JZ)*ZI**THJZ)*SINL+(ZI**THJZ*BLI(JZ)- | |
1425 | $ FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSL)/2. | |
1426 | Z1(2)=((-ZI**(THJZ)*(-1.)**(THJZ+1)*ALI(JZ)- | |
1427 | $ FL*ZMIXSS(2,JZ)*ZI**THJZ)*SINL+(ZI**THJZ*BLI(JZ)+ | |
1428 | $ FL*ZMIXSS(2,JZ)*(-ZI)**THJZ)*COSL)/2. | |
1429 | Z2(1)=-CONJG(((ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)- | |
1430 | $ FL*ZMIXSS(2,IZ)*ZI**THIZ)*SINL+(ZI**THIZ*BLI(IZ)- | |
1431 | $ FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSL)/2.) | |
1432 | Z2(2)=CONJG(((-ZI**(THIZ)*(-1.)**(THIZ+1)*ALI(IZ)- | |
1433 | $ FL*ZMIXSS(2,IZ)*ZI**THIZ)*SINL+(ZI**THIZ*BLI(IZ)+ | |
1434 | $ FL*ZMIXSS(2,IZ)*(-ZI)**THIZ)*COSL)/2.) | |
1435 | CALL SSME3(3,AML2SS,Z1,Z2) | |
1436 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
1437 | Z1(1)=(-ZI)**(THIZ+THJZ)*(XLIJ(IZ,JZ)+XLIJ(JZ,IZ)) | |
1438 | Z1(2)=0. | |
1439 | ELSE | |
1440 | Z1(1)=0. | |
1441 | Z1(2)=(-ZI)**(THIZ+THJZ)*(XLIJ(IZ,JZ)+XLIJ(JZ,IZ)) | |
1442 | END IF | |
1443 | Z2(1)=-G*MTAMZ*SINA/2./AMW/COSBE | |
1444 | Z2(2)=0. | |
1445 | CALL SSME3(4,AMHL,Z1,Z2) | |
1446 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
1447 | Z1(1)=(-ZI)**(THIZ+THJZ)*(XHIJ(IZ,JZ)+XHIJ(JZ,IZ)) | |
1448 | Z1(2)=0. | |
1449 | ELSE | |
1450 | Z1(1)=0. | |
1451 | Z1(2)=(-ZI)**(THIZ+THJZ)*(XHIJ(IZ,JZ)+XHIJ(JZ,IZ)) | |
1452 | END IF | |
1453 | Z2(1)=-G*MTAMZ*COSA/2./AMW/COSBE | |
1454 | Z2(2)=0. | |
1455 | CALL SSME3(4,AMHH,Z1,Z2) | |
1456 | IF (-SNIJ(JZ,IZ).GT.0.) THEN | |
1457 | Z1(1)=0. | |
1458 | Z1(2)=(-ZI)**(THIZ+THJZ+1)*(XPIJ(IZ,JZ)+XPIJ(JZ,IZ)) | |
1459 | ELSE | |
1460 | Z1(1)=(-ZI)**(THIZ+THJZ+1)*(XPIJ(IZ,JZ)+XPIJ(JZ,IZ)) | |
1461 | Z1(2)=0. | |
1462 | END IF | |
1463 | Z2(1)=0. | |
1464 | Z2(2)=ZI*G*MTAMZ*TANB/2./AMW | |
1465 | CALL SSME3(4,AMHA,Z1,Z2) | |
1466 | END IF | |
1467 | C zi -> zj + nu_e + nu_e bar | |
1468 | TMP(6)=0. | |
1469 | IF (MZIZ.GT.MZJZ) THEN | |
1470 | IF (MZIZ.LT.AMN1SS) THEN | |
1471 | TMP(4)=-SNIJ(JZ,IZ) | |
1472 | TMP(3)=AMN1SS | |
1473 | TERMLL=2*ANI(JZ)**2*ANI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
1474 | ELSE | |
1475 | TERMLL=0. | |
1476 | END IF | |
1477 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
1478 | TMP(4)=+SNIJ(JZ,IZ) | |
1479 | TERMZZ=(APN**2+BTN**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
1480 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ) | |
1481 | ELSE | |
1482 | TERMZZ=0. | |
1483 | END IF | |
1484 | IF (MZIZ.LT.AMN1SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
1485 | TMP(4)=+SNIJ(JZ,IZ) | |
1486 | TMP(3)=AMN1SS | |
1487 | TERMLZ=8*E*(APN-BTN)*WIJ(JZ,IZ)*ANI(JZ)*ANI(IZ)/MZIZ | |
1488 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
1489 | ELSE | |
1490 | TERMLZ=0. | |
1491 | END IF | |
1492 | WID=TERMLL+TERMZZ+TERMLZ | |
1493 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDNE,-IDNE,0,0) | |
1494 | C Enter information for decay matrix element | |
1495 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
1496 | Z1(1)=0. | |
1497 | Z1(2)=2*E*WIJ(JZ,IZ) | |
1498 | ELSE | |
1499 | Z1(1)=2*E*WIJ(JZ,IZ) | |
1500 | Z1(2)=0. | |
1501 | END IF | |
1502 | Z2(1)=APN | |
1503 | Z2(2)=BTN | |
1504 | CALL SSME3(1,AMZ,Z1,Z2) | |
1505 | Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2. | |
1506 | Z1(2)=-Z1(1) | |
1507 | Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.) | |
1508 | Z2(2)=Z2(1) | |
1509 | CALL SSME3(2,AMN1SS,Z1,Z2) | |
1510 | Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2. | |
1511 | Z1(2)=-Z1(1) | |
1512 | Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2.) | |
1513 | Z2(2)=Z2(1) | |
1514 | CALL SSME3(3,AMN1SS,Z1,Z2) | |
1515 | END IF | |
1516 | C zi -> zj + nu_mu + nu_mu bar | |
1517 | IF (MZIZ.GT.MZJZ) THEN | |
1518 | IF (MZIZ.LT.AMN2SS) THEN | |
1519 | TMP(4)=-SNIJ(JZ,IZ) | |
1520 | TMP(3)=AMN2SS | |
1521 | TERMLL=2*ANI(JZ)**2*ANI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
1522 | ELSE | |
1523 | TERMLL=0. | |
1524 | END IF | |
1525 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
1526 | TMP(4)=+SNIJ(JZ,IZ) | |
1527 | TERMZZ=(APN**2+BTN**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
1528 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ) | |
1529 | ELSE | |
1530 | TERMZZ=0. | |
1531 | END IF | |
1532 | IF (MZIZ.LT.AMN2SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
1533 | TMP(4)=+SNIJ(JZ,IZ) | |
1534 | TMP(3)=AMN2SS | |
1535 | TERMLZ=8*E*(APN-BTN)*WIJ(JZ,IZ)*ANI(JZ)*ANI(IZ)/MZIZ | |
1536 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
1537 | ELSE | |
1538 | TERMLZ=0. | |
1539 | END IF | |
1540 | WID=TERMLL+TERMZZ+TERMLZ | |
1541 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDNM,-IDNM,0,0) | |
1542 | C Enter information for decay matrix element | |
1543 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
1544 | Z1(1)=0. | |
1545 | Z1(2)=2*E*WIJ(JZ,IZ) | |
1546 | ELSE | |
1547 | Z1(1)=2*E*WIJ(JZ,IZ) | |
1548 | Z1(2)=0. | |
1549 | END IF | |
1550 | Z2(1)=APN | |
1551 | Z2(2)=BTN | |
1552 | CALL SSME3(1,AMZ,Z1,Z2) | |
1553 | Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2. | |
1554 | Z1(2)=-Z1(1) | |
1555 | Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.) | |
1556 | Z2(2)=Z2(1) | |
1557 | CALL SSME3(2,AMN2SS,Z1,Z2) | |
1558 | Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2. | |
1559 | Z1(2)=-Z1(1) | |
1560 | Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2.) | |
1561 | Z2(2)=Z2(1) | |
1562 | CALL SSME3(3,AMN2SS,Z1,Z2) | |
1563 | END IF | |
1564 | C zi -> zj + nu_tau + nu_tau bar | |
1565 | IF (MZIZ.GT.MZJZ) THEN | |
1566 | IF (MZIZ.LT.AMN3SS) THEN | |
1567 | TMP(4)=-SNIJ(JZ,IZ) | |
1568 | TMP(3)=AMN3SS | |
1569 | TERMLL=2*ANI(JZ)**2*ANI(IZ)**2*FAC*SSXINT(0.,SSZZF1,1.) | |
1570 | ELSE | |
1571 | TERMLL=0. | |
1572 | END IF | |
1573 | IF (MZIZ.LT.(MZJZ+AMZ)) THEN | |
1574 | TMP(4)=+SNIJ(JZ,IZ) | |
1575 | TERMZZ=(APN**2+BTN**2)*E**2*WIJ(JZ,IZ)**2/MZIZ | |
1576 | $ *SSXINT(MZJZ,SSZZF2,(MZIZ**2+MZJZ**2)/2./MZIZ) | |
1577 | ELSE | |
1578 | TERMZZ=0. | |
1579 | END IF | |
1580 | IF (MZIZ.LT.AMN3SS.AND.MZIZ.LT.(MZJZ+AMZ)) THEN | |
1581 | TMP(4)=+SNIJ(JZ,IZ) | |
1582 | TMP(3)=AMN3SS | |
1583 | TERMLZ=8*E*(APN-BTN)*WIJ(JZ,IZ)*ANI(JZ)*ANI(IZ)/MZIZ | |
1584 | $ /(2*PI)**5*SSXINT(0.,SSZZF3,(MZIZ-MZJZ)**2) | |
1585 | ELSE | |
1586 | TERMLZ=0. | |
1587 | END IF | |
1588 | WID=TERMLL+TERMZZ+TERMLZ | |
1589 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),IDNT,-IDNT,0,0) | |
1590 | C Enter information for decay matrix element | |
1591 | IF (-1.*SNIJ(JZ,IZ).GT.0.) THEN | |
1592 | Z1(1)=0. | |
1593 | Z1(2)=2*E*WIJ(JZ,IZ) | |
1594 | ELSE | |
1595 | Z1(1)=2*E*WIJ(JZ,IZ) | |
1596 | Z1(2)=0. | |
1597 | END IF | |
1598 | Z2(1)=APN | |
1599 | Z2(2)=BTN | |
1600 | CALL SSME3(1,AMZ,Z1,Z2) | |
1601 | Z1(1)=ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2. | |
1602 | Z1(2)=-Z1(1) | |
1603 | Z2(1)=CONJG(ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2.) | |
1604 | Z2(2)=Z2(1) | |
1605 | CALL SSME3(2,AMN3SS,Z1,Z2) | |
1606 | Z1(1)=ZI**(THJZ-1)*(-1.)**(THJZ+1)*ANI(JZ)/2. | |
1607 | Z1(2)=-Z1(1) | |
1608 | Z2(1)=-CONJG(ZI**(THIZ-1)*(-1.)**(THIZ+1)*ANI(IZ)/2.) | |
1609 | Z2(2)=Z2(1) | |
1610 | CALL SSME3(3,AMN3SS,Z1,Z2) | |
1611 | END IF | |
1612 | 220 CONTINUE | |
1613 | C | |
1614 | C zi --> zj + higgs | |
1615 | C | |
1616 | DO 230 JZ=1,IZ-1 | |
1617 | C zi --> zj + hl | |
1618 | MZJZ=ABS(AMZISS(JZ)) | |
1619 | IF (MZIZ.GT.(MZJZ+AMHL)) THEN | |
1620 | SN=SIGN(1.,AMZISS(JZ))*SIGN(1.,AMZISS(IZ)) | |
1621 | WID=(XLIJ(JZ,IZ)+XLIJ(IZ,JZ))**2/8./PI/(MZIZ)**3 | |
1622 | $ *SQRT(MZIZ**4+MZJZ**4+AMHL**4-2.*(MZIZ*MZJZ)**2 | |
1623 | $ -2.*(MZIZ*AMHL)**2-2.*(MZJZ*AMHL)**2)*((MZIZ**2+MZJZ**2 | |
1624 | $ -AMHL**2)/2.+SN*MZIZ*MZJZ) | |
1625 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),ISHL,0,0,0) | |
1626 | END IF | |
1627 | C zi --> zj + hh | |
1628 | IF (MZIZ.GT.(MZJZ+AMHH)) THEN | |
1629 | SN=SIGN(1.,AMZISS(JZ))*SIGN(1.,AMZISS(IZ)) | |
1630 | WID=(HIJ(JZ,IZ)+HIJ(IZ,JZ))**2/8./PI/(MZIZ)**3 | |
1631 | $ *SQRT(MZIZ**4+MZJZ**4+AMHH**4-2.*(MZIZ*MZJZ)**2 | |
1632 | $ -2.*(MZIZ*AMHH)**2-2.*(MZJZ*AMHH)**2) | |
1633 | $ *((MZIZ**2+MZJZ**2-AMHH**2)/2.+SN*MZIZ*MZJZ) | |
1634 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),ISHH,0,0,0) | |
1635 | END IF | |
1636 | C zi --> zj + ha | |
1637 | IF (MZIZ.GT.(MZJZ+AMHA)) THEN | |
1638 | SN=-SIGN(1.,AMZISS(JZ))*SIGN(1.,AMZISS(IZ)) | |
1639 | WID=(XPIJ(IZ,JZ)+XPIJ(JZ,IZ))**2/8./PI/(MZIZ)**3 | |
1640 | $ *SQRT(MZIZ**4+MZJZ**4+AMHA**4-2.*(MZIZ*MZJZ)**2 | |
1641 | $ -2.*(MZIZ*AMHA)**2-2.*(MZJZ*AMHA)**2)*((MZIZ**2+MZJZ**2 | |
1642 | $ -AMHA**2)/2.+SN*MZIZ*MZJZ) | |
1643 | CALL SSSAVE(ISZIZ,WID,ISZ(JZ),ISHA,0,0,0) | |
1644 | END IF | |
1645 | 230 CONTINUE | |
1646 | 200 CONTINUE | |
1647 | C | |
1648 | C zi --> squark + qbar; enlarge to include Z1 decays | |
1649 | C in case of models with light gravitino | |
1650 | C | |
1651 | DO 245 IZ=1,4 | |
1652 | MZIZ=ABS(AMZISS(IZ)) | |
1653 | SNIZ=SIGN(1.,AMZISS(IZ)) | |
1654 | IF (SNIZ.EQ.1.) THEN | |
1655 | THIZ=0 | |
1656 | ELSE | |
1657 | THIZ=1 | |
1658 | END IF | |
1659 | ISZIZ=ISZ(IZ) | |
1660 | C left squarks | |
1661 | IF (MZIZ.GT.(AMULSS+AMUP)) THEN | |
1662 | WID=3*AUI(IZ)**2*(MZIZ**2+AMUP**2-AMULSS**2)/MZIZ**3/ | |
1663 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMUP**2,AMULSS**2)) | |
1664 | CALL SSSAVE(ISZIZ,WID,ISUPL,-IDUP,0,0,0) | |
1665 | CALL SSSAVE(ISZIZ,WID,-ISUPL,IDUP,0,0,0) | |
1666 | ENDIF | |
1667 | IF (MZIZ.GT.(AMDLSS+AMDN)) THEN | |
1668 | WID=3*ADI(IZ)**2*(MZIZ**2+AMDN**2-AMDLSS**2)/MZIZ**3/ | |
1669 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMDN**2,AMDLSS**2)) | |
1670 | CALL SSSAVE(ISZIZ,WID,ISDNL,-IDDN,0,0,0) | |
1671 | CALL SSSAVE(ISZIZ,WID,-ISDNL,IDDN,0,0,0) | |
1672 | END IF | |
1673 | IF (MZIZ.GT.(AMSLSS+AMST)) THEN | |
1674 | WID=3*ADI(IZ)**2*(MZIZ**2+AMST**2-AMSLSS**2)/MZIZ**3/ | |
1675 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMST**2,AMSLSS**2)) | |
1676 | CALL SSSAVE(ISZIZ,WID,ISSTL,-IDST,0,0,0) | |
1677 | CALL SSSAVE(ISZIZ,WID,-ISSTL,IDST,0,0,0) | |
1678 | END IF | |
1679 | IF (MZIZ.GT.(AMCLSS+AMCH)) THEN | |
1680 | WID=3*AUI(IZ)**2*(MZIZ**2+AMCH**2-AMCLSS**2)/MZIZ**3/ | |
1681 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMCH**2,AMCLSS**2)) | |
1682 | CALL SSSAVE(ISZIZ,WID,ISCHL,-IDCH,0,0,0) | |
1683 | CALL SSSAVE(ISZIZ,WID,-ISCHL,IDCH,0,0,0) | |
1684 | ENDIF | |
1685 | C right squarks | |
1686 | IF (MZIZ.GT.(AMURSS+AMUP)) THEN | |
1687 | WID=3*BUI(IZ)**2*(MZIZ**2+AMUP**2-AMURSS**2)/MZIZ**3/ | |
1688 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMUP**2,AMURSS**2)) | |
1689 | CALL SSSAVE(ISZIZ,WID,ISUPR,-IDUP,0,0,0) | |
1690 | CALL SSSAVE(ISZIZ,WID,-ISUPR,IDUP,0,0,0) | |
1691 | END IF | |
1692 | IF (MZIZ.GT.(AMDRSS+AMDN)) THEN | |
1693 | WID=3*BDI(IZ)**2*(MZIZ**2+AMDN**2-AMDRSS**2)/MZIZ**3/ | |
1694 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMDN**2,AMDRSS**2)) | |
1695 | CALL SSSAVE(ISZIZ,WID,ISDNR,-IDDN,0,0,0) | |
1696 | CALL SSSAVE(ISZIZ,WID,-ISDNR,IDDN,0,0,0) | |
1697 | END IF | |
1698 | IF (MZIZ.GT.(AMSRSS+AMST)) THEN | |
1699 | WID=3*BDI(IZ)**2*(MZIZ**2+AMST**2-AMSRSS**2)/MZIZ**3/ | |
1700 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMST**2,AMSRSS**2)) | |
1701 | CALL SSSAVE(ISZIZ,WID,ISSTR,-IDST,0,0,0) | |
1702 | CALL SSSAVE(ISZIZ,WID,-ISSTR,IDST,0,0,0) | |
1703 | END IF | |
1704 | IF(MZIZ.GT.(AMCRSS+AMCH)) THEN | |
1705 | WID=3*BUI(IZ)**2*(MZIZ**2+AMCH**2-AMCRSS**2)/MZIZ**3/ | |
1706 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMCH**2,AMCRSS**2)) | |
1707 | CALL SSSAVE(ISZIZ,WID,ISCHR,-IDCH,0,0,0) | |
1708 | CALL SSSAVE(ISZIZ,WID,-ISCHR,IDCH,0,0,0) | |
1709 | ENDIF | |
1710 | C z_i --> sbottom + bottom | |
1711 | C | |
1712 | ZAUIZ=ZI**(THIZ-1)*SNIZ | |
1713 | $ *(-G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ)) | |
1714 | ZBUIZ=ZI**(THIZ-1)*4*GP*ZMIXSS(4,IZ)/3./SR2 | |
1715 | ZADIZ=ZI**(THIZ-1)*SNIZ | |
1716 | $ *(G/SR2*ZMIXSS(3,IZ)-GP/3./SR2*ZMIXSS(4,IZ)) | |
1717 | ZBDIZ=-2*ZI**(THIZ-1)*GP*ZMIXSS(4,IZ)/3./SR2 | |
1718 | ZALIZ=ZI**(THIZ-1)*SNIZ | |
1719 | $ *(G/SR2*ZMIXSS(3,IZ)+GP/SR2*ZMIXSS(4,IZ)) | |
1720 | ZBLIZ=-1*ZI**(THIZ-1)*SR2*GP*ZMIXSS(4,IZ) | |
1721 | ZPP=ZI**THIZ | |
1722 | ZPM=(-ZI)**THIZ | |
1723 | IF(MZIZ.GT.(AMB1SS+AMBT)) THEN | |
1724 | ZA=((ZI*ZADIZ-ZPP*FB*ZMIXSS(2,IZ))*COSB | |
1725 | $ -(ZI*ZBDIZ-ZPM*FB*ZMIXSS(2,IZ))*SINB)/2. | |
1726 | ZB=((-ZI*ZADIZ-ZPP*FB*ZMIXSS(2,IZ))*COSB | |
1727 | $ -(ZI*ZBDIZ+ZPM*FB*ZMIXSS(2,IZ))*SINB)/2. | |
1728 | AS=ZA*CONJG(ZA) | |
1729 | BS=ZB*CONJG(ZB) | |
1730 | WID=3*(AS*((AMBT+MZIZ)**2-AMB1SS**2)+BS*((MZIZ-AMBT)**2- | |
1731 | $ AMB1SS**2))/16./PI/MZIZ**3* | |
1732 | $ SQRT(SSXLAM(MZIZ**2,AMB1SS**2,AMBT**2)) | |
1733 | CALL SSSAVE(ISZIZ,WID,ISBT1,-IDBT,0,0,0) | |
1734 | CALL SSSAVE(ISZIZ,WID,-ISBT1,IDBT,0,0,0) | |
1735 | ENDIF | |
1736 | IF(MZIZ.GT.(AMB2SS+AMBT)) THEN | |
1737 | ZA=((ZI*ZADIZ-ZPP*FB*ZMIXSS(2,IZ))*SINB | |
1738 | $ +(ZI*ZBDIZ-ZPM*FB*ZMIXSS(2,IZ))*COSB)/2. | |
1739 | ZB=((-ZI*ZADIZ-ZPP*FB*ZMIXSS(2,IZ))*SINB | |
1740 | $ +(ZI*ZBDIZ+ZPM*FB*ZMIXSS(2,IZ))*COSB)/2. | |
1741 | AS=ZA*CONJG(ZA) | |
1742 | BS=ZB*CONJG(ZB) | |
1743 | WID=3*(AS*((AMBT+MZIZ)**2-AMB2SS**2)+BS*((MZIZ-AMBT)**2- | |
1744 | $ AMB2SS**2))/16./PI/MZIZ**3* | |
1745 | $ SQRT(SSXLAM(MZIZ**2,AMB2SS**2,AMBT**2)) | |
1746 | CALL SSSAVE(ISZIZ,WID,ISBT2,-IDBT,0,0,0) | |
1747 | CALL SSSAVE(ISZIZ,WID,-ISBT2,IDBT,0,0,0) | |
1748 | ENDIF | |
1749 | C z_i --> stop + top | |
1750 | C | |
1751 | IF(MZIZ.GT.AMT1SS+AMTP) THEN | |
1752 | ZA=((ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST | |
1753 | $ -(ZI*ZBUIZ-ZPM*FT*ZMIXSS(1,IZ))*SINT)/2. | |
1754 | ZB=((-ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*COST | |
1755 | $ -(ZI*ZBUIZ+ZPM*FT*ZMIXSS(1,IZ))*SINT)/2. | |
1756 | AS=ZA*CONJG(ZA) | |
1757 | BS=ZB*CONJG(ZB) | |
1758 | WID=3*(AS*((AMTP+MZIZ)**2-AMT1SS**2)+BS*((MZIZ-AMTP)**2- | |
1759 | $ AMT1SS**2))/16./PI/MZIZ**3* | |
1760 | $ SQRT(SSXLAM(MZIZ**2,AMT1SS**2,AMTP**2)) | |
1761 | CALL SSSAVE(ISZIZ,WID,ISTP1,-IDTP,0,0,0) | |
1762 | CALL SSSAVE(ISZIZ,WID,-ISTP1,IDTP,0,0,0) | |
1763 | ENDIF | |
1764 | IF(MZIZ.GT.AMT2SS+AMTP) THEN | |
1765 | ZA=((ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*SINT | |
1766 | $ +(ZI*ZBUIZ-ZPM*FT*ZMIXSS(1,IZ))*COST)/2. | |
1767 | ZB=((-ZI*ZAUIZ-ZPP*FT*ZMIXSS(1,IZ))*SINT | |
1768 | $ +(ZI*ZBUIZ+ZPM*FT*ZMIXSS(1,IZ))*COST)/2. | |
1769 | AS=ZA*CONJG(ZA) | |
1770 | BS=ZB*CONJG(ZB) | |
1771 | WID=3*(AS*((AMTP+MZIZ)**2-AMT2SS**2)+BS*((MZIZ-AMTP)**2- | |
1772 | $ AMT2SS**2))/16./PI/MZIZ**3* | |
1773 | $ SQRT(SSXLAM(MZIZ**2,AMT2SS**2,AMTP**2)) | |
1774 | CALL SSSAVE(ISZIZ,WID,ISTP2,-IDTP,0,0,0) | |
1775 | CALL SSSAVE(ISZIZ,WID,-ISTP2,IDTP,0,0,0) | |
1776 | ENDIF | |
1777 | C | |
1778 | C zi --> slepton + lepton | |
1779 | C | |
1780 | IF(MZIZ.GT.(AMELSS+AME)) THEN | |
1781 | WID=ALI(IZ)**2*(MZIZ**2+AME**2-AMELSS**2)/MZIZ**3/ | |
1782 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AME**2,AMELSS**2)) | |
1783 | CALL SSSAVE(ISZIZ,WID,ISEL,-IDE,0,0,0) | |
1784 | CALL SSSAVE(ISZIZ,WID,-ISEL,IDE,0,0,0) | |
1785 | END IF | |
1786 | IF(MZIZ.GT.(AMMLSS+AMMU)) THEN | |
1787 | WID=ALI(IZ)**2*(MZIZ**2+AMMU**2-AMMLSS**2)/MZIZ**3/ | |
1788 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMMU**2,AMMLSS**2)) | |
1789 | CALL SSSAVE(ISZIZ,WID,ISMUL,-IDMU,0,0,0) | |
1790 | CALL SSSAVE(ISZIZ,WID,-ISMUL,IDMU,0,0,0) | |
1791 | END IF | |
1792 | IF(MZIZ.GT.(AMERSS+AME)) THEN | |
1793 | WID=BLI(IZ)**2*(MZIZ**2+AME**2-AMERSS**2)/MZIZ**3/ | |
1794 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AME**2,AMERSS**2)) | |
1795 | CALL SSSAVE(ISZIZ,WID,ISER,-IDE,0,0,0) | |
1796 | CALL SSSAVE(ISZIZ,WID,-ISER,IDE,0,0,0) | |
1797 | END IF | |
1798 | IF(MZIZ.GT.(AMMRSS+AMMU)) THEN | |
1799 | WID=BLI(IZ)**2*(MZIZ**2+AMMU**2-AMMRSS**2)/MZIZ**3/ | |
1800 | $ 32./PI*SQRT(SSXLAM(MZIZ**2,AMMU**2,AMMRSS**2)) | |
1801 | CALL SSSAVE(ISZIZ,WID,ISMUR,-IDMU,0,0,0) | |
1802 | CALL SSSAVE(ISZIZ,WID,-ISMUR,IDMU,0,0,0) | |
1803 | END IF | |
1804 | IF(MZIZ.GT.(AML1SS+AMTAU)) THEN | |
1805 | ZA=((ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*COSL | |
1806 | $ -(ZI*ZBLIZ-ZPM*FL*ZMIXSS(2,IZ))*SINL)/2. | |
1807 | ZB=((-ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*COSL | |
1808 | $ -(ZI*ZBLIZ+ZPM*FL*ZMIXSS(2,IZ))*SINL)/2. | |
1809 | AS=ZA*CONJG(ZA) | |
1810 | BS=ZB*CONJG(ZB) | |
1811 | WID=(AS*((AMTAU+MZIZ)**2-AML1SS**2)+BS*((MZIZ-AMTAU)**2- | |
1812 | $ AML1SS**2))/16./PI/MZIZ**3* | |
1813 | $ SQRT(SSXLAM(MZIZ**2,AML1SS**2,AMTAU**2)) | |
1814 | CALL SSSAVE(ISZIZ,WID,ISTAU1,-IDTAU,0,0,0) | |
1815 | CALL SSSAVE(ISZIZ,WID,-ISTAU1,IDTAU,0,0,0) | |
1816 | END IF | |
1817 | IF(MZIZ.GT.(AML2SS+AMTAU)) THEN | |
1818 | ZA=((ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*SINL | |
1819 | $ +(ZI*ZBLIZ-ZPM*FL*ZMIXSS(2,IZ))*COSL)/2. | |
1820 | ZB=((-ZI*ZALIZ-ZPP*FL*ZMIXSS(2,IZ))*SINL | |
1821 | $ +(ZI*ZBLIZ+ZPM*FL*ZMIXSS(2,IZ))*COSL)/2. | |
1822 | AS=ZA*CONJG(ZA) | |
1823 | BS=ZB*CONJG(ZB) | |
1824 | WID=(AS*((AMTAU+MZIZ)**2-AML2SS**2)+BS*((MZIZ-AMTAU)**2- | |
1825 | $ AML2SS**2))/16./PI/MZIZ**3* | |
1826 | $ SQRT(SSXLAM(MZIZ**2,AML2SS**2,AMTAU**2)) | |
1827 | CALL SSSAVE(ISZIZ,WID,ISTAU2,-IDTAU,0,0,0) | |
1828 | CALL SSSAVE(ISZIZ,WID,-ISTAU2,IDTAU,0,0,0) | |
1829 | ENDIF | |
1830 | IF(MZIZ.GT.AMN1SS) THEN | |
1831 | WID=(MZIZ**2-AMN1SS**2)**2/32./PI/MZIZ**3 | |
1832 | CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,ISNEL,-IDNE,0,0,0) | |
1833 | CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,-ISNEL,IDNE,0,0,0) | |
1834 | END IF | |
1835 | IF(MZIZ.GT.AMN2SS) THEN | |
1836 | WID=(MZIZ**2-AMN2SS**2)**2/32./PI/MZIZ**3 | |
1837 | CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,ISNML,-IDNM,0,0,0) | |
1838 | CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,-ISNML,IDNM,0,0,0) | |
1839 | END IF | |
1840 | IF(MZIZ.GT.AMN3SS) THEN | |
1841 | WID=(MZIZ**2-AMN3SS**2)**2/32./PI/MZIZ**3 | |
1842 | CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,ISNTL,-IDNT,0,0,0) | |
1843 | CALL SSSAVE(ISZIZ,ANI(IZ)**2*WID,-ISNTL,IDNT,0,0,0) | |
1844 | END IF | |
1845 | 245 CONTINUE | |
1846 | C | |
1847 | C Compute decays to gravitino for GMSB models | |
1848 | C | |
1849 | DO 250 IZ=1,4 | |
1850 | MZIZ=ABS(AMZISS(IZ)) | |
1851 | ISZIZ=ISZ(IZ) | |
1852 | IF (MZIZ.GT.AMGVSS) THEN | |
1853 | WID=(ZMIXSS(4,IZ)*CTHW+ZMIXSS(3,IZ)*STHW)**2*MZIZ**5/ | |
1854 | , 48./PI/(AMGVSS*AMPL)**2 | |
1855 | CALL SSSAVE(ISZIZ,WID,ISGRAV,IDGM,0,0,0) | |
1856 | C Dalitz decay | |
1857 | IF (MZIZ.GT.(AMGVSS+2*AME)) THEN | |
1858 | WIDEE=WID*2*ALFAEM/(3*PI)*LOG(MZIZ/AME) | |
1859 | CALL SSSAVE(ISZIZ,WIDEE,ISGRAV,IDE,-IDE,0,0) | |
1860 | END IF | |
1861 | END IF | |
1862 | IF (MZIZ.GT.(AMZ+AMGVSS)) THEN | |
1863 | WID=(2*(ZMIXSS(4,IZ)*STHW-ZMIXSS(3,IZ)*CTHW)**2+ | |
1864 | , (ZMIXSS(1,IZ)*SINBE-ZMIXSS(2,IZ)*COSBE)**2)* | |
1865 | , (MZIZ**2-AMZ**2)**4/96./PI/MZIZ**3/(AMGVSS*AMPL)**2 | |
1866 | CALL SSSAVE(ISZIZ,WID,ISGRAV,IDZ,0,0,0) | |
1867 | END IF | |
1868 | IF (MZIZ.GT.(AMHL+AMGVSS)) THEN | |
1869 | WID=(ZMIXSS(1,IZ)*COSA+ZMIXSS(2,IZ)*SINA)**2/6./ | |
1870 | , (AMGVSS*AMPL)**2*(MZIZ**2-AMHL**2)**4/16./PI/MZIZ**3 | |
1871 | CALL SSSAVE(ISZIZ,WID,ISGRAV,ISHL,0,0,0) | |
1872 | END IF | |
1873 | IF (MZIZ.GT.(AMHH+AMGVSS)) THEN | |
1874 | WID=(-ZMIXSS(1,IZ)*SINA+ZMIXSS(2,IZ)*COSA)**2/6./ | |
1875 | , (AMGVSS*AMPL)**2*(MZIZ**2-AMHH**2)**4/16./PI/MZIZ**3 | |
1876 | CALL SSSAVE(ISZIZ,WID,ISGRAV,ISHH,0,0,0) | |
1877 | END IF | |
1878 | IF (MZIZ.GT.(AMHA+AMGVSS)) THEN | |
1879 | WID=(ZMIXSS(1,IZ)*COSBE+ZMIXSS(2,IZ)*SINBE)**2/6./ | |
1880 | , (AMGVSS*AMPL)**2*(MZIZ**2-AMHA**2)**4/16./PI/MZIZ**3 | |
1881 | CALL SSSAVE(ISZIZ,WID,ISGRAV,ISHA,0,0,0) | |
1882 | END IF | |
1883 | 250 CONTINUE | |
1884 | C | |
1885 | C Normalize zi branching ratios | |
1886 | C | |
1887 | CALL SSNORM(ISZ1) | |
1888 | CALL SSNORM(ISZ2) | |
1889 | CALL SSNORM(ISZ3) | |
1890 | CALL SSNORM(ISZ4) | |
1891 | C | |
1892 | RETURN | |
1893 | END |