]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | SUBROUTINE SIGSSE | |
3 | C | |
4 | C Compute d(sigma)/d(cos theta) for | |
5 | C e+ e- ----> SUSY particles | |
6 | C See Baer et. al., IJMP A4, 4111 (1989) for sigma's | |
7 | C Polarized cross sections added 9/18/95 hb | |
8 | C Mixed sbottoms and staus included 10/23/96 hb | |
9 | C | |
10 | C SIGMA = cross section summed over quark types allowed by | |
11 | C JETTYPE and WTYPE cards. | |
12 | C SIGS(I) = partial cross section for I1 + I2 --> I3 + I4. | |
13 | C INOUT(I) = IOPAK**3*I4 + IOPAK**2*I3 + IOPAK*I2 + I1 | |
14 | C using JETTYPE code. | |
15 | C | |
16 | C Extra factor of 1/2 needed because all jets are treated | |
17 | C as identical. | |
18 | C | |
19 | #if defined(CERNLIB_IMPNONE) | |
20 | IMPLICIT NONE | |
21 | #endif | |
22 | #include "isajet/itapes.inc" | |
23 | #include "isajet/jetsig.inc" | |
24 | #include "isajet/eepar.inc" | |
25 | #include "isajet/primar.inc" | |
26 | #include "isajet/jetpar.inc" | |
27 | #include "isajet/q1q2.inc" | |
28 | #include "isajet/wcon.inc" | |
29 | #include "isajet/const.inc" | |
30 | #include "isajet/sspar.inc" | |
31 | #include "isajet/sssm.inc" | |
32 | #include "isajet/sstype.inc" | |
33 | #include "isajet/brembm.inc" | |
34 | C | |
35 | REAL ALQ(2),BEQ(2),E,CS2THW,TNTHW,CTTHW,AE,BE,AM1,AM2, | |
36 | $EQ,ALR,Z,PHIZ,PROPZ,SIG,PCM,AMASS,ALL(2),BEL(2), | |
37 | $G,MSNE,TM2,TM3,TM4,TM5,TM6,AZJ,AZI,MEL,MER, | |
38 | $AEZS,BEZS,SR2,GP,AN,BN,AEZJS,BEZJS,SSXLAM, | |
39 | $TGG,TNN,TGN,TZN,AMWI,XS,YS,XC,YC,SINGL,SINGR, | |
40 | $COSGL,COSGR,XM,YM,THX,THY,XI,DEL,AMWISS(2),KK, | |
41 | $AMZIZ1,AMZIZ2,SIGLL,SIGRR,SIGLZ,SIGRZ,SSGT,SSGST, | |
42 | $FAC1,EZ0,BETA,EEL,EER, | |
43 | $FLEP,FLEM,FREP,FREM,SIGLR,SIGRL,PHIZLR,PHIZRL, | |
44 | $TM1LR,TM1RL,TZZRL,TZZLR,TGZLR,TGZRL,SIGZZL,SIGZZR, | |
45 | $FACLR,FACRL,RSH,JAC,ESTRUC,SH,SSFEL | |
46 | COMPLEX AEZ(4),BEZ(4),ZI,ZONE,WIJ | |
47 | INTEGER IS2UD(25),IUD(13),JS2JT(25),IQ1,IQ2,IFL1,IFL2, | |
48 | $IFLQ,IFM,I,IDQSS(25),MATCHL(18),IL2JS(18),IS2LN(18), | |
49 | $I1,I2,IL1,IL2,IDL1,IDL2,IZ,IZ1,IP,ITHZ(4),IDLSS(18), | |
50 | $IW2JS(4),IW1,JW1,JTW1,JTW2,IZ2JS(4), | |
51 | $IZ2,JTYPZ1,JTYPZ2 | |
52 | INTEGER MSUPL,MSDNL,MSSTL,MSCHL,MSBT1,MSTP1, | |
53 | $MSUPR,MSDNR,MSSTR,MSCHR,MSBT2,MSTP2,MSW1,MSW2, | |
54 | $MSNEL,MSEL,MSNML,MSMUL,MSNTL,MSTAU1,MSER,MSMUR,MSTAU2 | |
55 | PARAMETER (MSUPL=-ISUPL) | |
56 | PARAMETER (MSDNL=-ISDNL) | |
57 | PARAMETER (MSSTL=-ISSTL) | |
58 | PARAMETER (MSCHL=-ISCHL) | |
59 | PARAMETER (MSBT1=-ISBT1) | |
60 | PARAMETER (MSTP1=-ISTP1) | |
61 | PARAMETER (MSUPR=-ISUPR) | |
62 | PARAMETER (MSDNR=-ISDNR) | |
63 | PARAMETER (MSSTR=-ISSTR) | |
64 | PARAMETER (MSCHR=-ISCHR) | |
65 | PARAMETER (MSBT2=-ISBT2) | |
66 | PARAMETER (MSTP2=-ISTP2) | |
67 | PARAMETER (MSW1=-ISW1) | |
68 | PARAMETER (MSW2=-ISW2) | |
69 | PARAMETER (MSNEL=-ISNEL) | |
70 | PARAMETER (MSEL=-ISEL) | |
71 | PARAMETER (MSNML=-ISNML) | |
72 | PARAMETER (MSMUL=-ISMUL) | |
73 | PARAMETER (MSNTL=-ISNTL) | |
74 | PARAMETER (MSTAU1=-ISTAU1) | |
75 | PARAMETER (MSER=-ISER) | |
76 | PARAMETER (MSMUR=-ISMUR) | |
77 | PARAMETER (MSTAU2=-ISTAU2) | |
78 | DATA IDQSS/0, | |
79 | $ISUPL,MSUPL,ISDNL,MSDNL,ISSTL,MSSTL,ISCHL,MSCHL,ISBT1,MSBT1, | |
80 | $ISTP1,MSTP1, | |
81 | $ISUPR,MSUPR,ISDNR,MSDNR,ISSTR,MSSTR,ISCHR,MSCHR,ISBT2,MSBT2, | |
82 | $ISTP2,MSTP2/ | |
83 | DATA IDLSS/ISNEL,MSNEL,ISEL,MSEL,ISNML,MSNML,ISMUL,MSMUL, | |
84 | $ISNTL,MSNTL,ISTAU1,MSTAU1,ISER,MSER,ISMUR,MSMUR, | |
85 | $ISTAU2,MSTAU2/ | |
86 | DATA IS2UD/0,1,1,2,2,2,2,1,1,2,2,1,1,1,1,2,2,2,2,1,1,2,2,1,1/ | |
87 | DATA IUD/0,1,-1,2,-2,2,-2,1,-1,2,-2,1,-1/ | |
88 | DATA JS2JT/1, | |
89 | $2,3,4,5,6,7,8,9,10,11,12,13,2,3,4,5,6,7,8,9,10,11,12,13/ | |
90 | DATA MATCHL/2,1,4,3,6,5,8,7,10,9,12,11,14,13,16,15,18,17/ | |
91 | DATA IL2JS/34,35,36,37,38,39,40,41,42,43,44,45,46,47, | |
92 | $48,49,50,51/ | |
93 | DATA IS2LN/1,1,2,2,1,1,2,2,1,1,2,2,2,2,2,2,2,2/ | |
94 | DATA IW2JS/26,27,28,29/ | |
95 | DATA IZ2JS/30,31,32,33/ | |
96 | DATA ZONE,ZI/(1.,0.),(0.,1.)/ | |
97 | C | |
98 | C FUNCTIONS | |
99 | IF (IBREM) THEN | |
100 | SH=SHAT | |
101 | JAC=2*(1.-SHAT/SCM)*2*SQRT(SHAT)*(RSHMAX-RSHMIN)/SCM/(X1+X2) | |
102 | ELSE | |
103 | SH=SCM | |
104 | END IF | |
105 | PROPZ=(SH-AMZ**2)**2+AMZ**2*GAMZ**2 | |
106 | C | |
107 | C CONSTANTS | |
108 | RSH=SQRT(SH) | |
109 | EB=RSH/2. | |
110 | QSQBM=QSQ | |
111 | E=SQRT(4*PI*ALFAEM) | |
112 | G=SQRT(4*PI*ALFAEM/SN2THW) | |
113 | GP=G*SQRT(SN2THW/(1.-SN2THW)) | |
114 | BETA=ATAN(1./RV2V1) | |
115 | SR2=SQRT(2.) | |
116 | CS2THW=1.-SN2THW | |
117 | TNTHW=SQRT(SN2THW/CS2THW) | |
118 | CTTHW=1./TNTHW | |
119 | ALQ(1)=CTTHW/4.-5*TNTHW/12. | |
120 | BEQ(1)=-(CTTHW+TNTHW)/4. | |
121 | ALQ(2)=TNTHW/12.-CTTHW/4. | |
122 | BEQ(2)=-BEQ(1) | |
123 | ALL(1)=(CTTHW+TNTHW)/4. | |
124 | BEL(1)=-(CTTHW+TNTHW)/4. | |
125 | ALL(2)=(3*TNTHW-CTTHW)/4. | |
126 | BEL(2)=-BEL(1) | |
127 | AE=ALL(2) | |
128 | BE=BEL(2) | |
129 | AN=ALL(1) | |
130 | BN=BEL(1) | |
131 | FLEP=(1.+PLEP)/2. | |
132 | FLEM=(1.+PLEM)/2. | |
133 | FREP=(1.-PLEP)/2. | |
134 | FREM=(1.-PLEM)/2. | |
135 | MEL=AMASS(ISEL) | |
136 | MER=AMASS(ISER) | |
137 | MSNE=AMASS(ISNEL) | |
138 | XM=1./TAN(GAMMAL) | |
139 | YM=1./TAN(GAMMAR) | |
140 | THX=SIGN(1.,XM) | |
141 | THY=SIGN(1.,YM) | |
142 | AMWISS(1)=ABS(AMW1SS) | |
143 | AMWISS(2)=ABS(AMW2SS) | |
144 | DO 5 IZ=1,4 | |
145 | ITHZ(IZ)=0 | |
146 | IF (AMZISS(IZ).LT.0.) ITHZ(IZ)=1 | |
147 | AEZ(IZ)=-1*ZI**(ITHZ(IZ)-1)*(-1)**(ITHZ(IZ)+1)* | |
148 | $ (G*ZMIXSS(3,IZ)+GP*ZMIXSS(4,IZ))/SR2 | |
149 | BEZ(IZ)=-1*ZI**(ITHZ(IZ)-1)*SR2*GP*ZMIXSS(4,IZ) | |
150 | 5 CONTINUE | |
151 | C | |
152 | C ENTRY | |
153 | SIG=0. | |
154 | SIGMA=0. | |
155 | NSIGS=0 | |
156 | DO 10 I=1,MXSIGS | |
157 | SIGS(I)=0. | |
158 | 10 CONTINUE | |
159 | C | |
160 | C First do squark pairs: IQ1 labels JETTYPE1. | |
161 | C | |
162 | DO 100 IQ1=2,25 | |
163 | IQ2=MATCH(IQ1,4) | |
164 | IF(.NOT.(GOQ(IQ1,1).AND.GOQ(IQ2,2))) GO TO 100 | |
165 | IFL1=IDQSS(IQ1) | |
166 | IFL2=IDQSS(IQ2) | |
167 | AM1=AMASS(IFL1) | |
168 | AM2=AMASS(IFL2) | |
169 | IF((AM1+AM2).GE.RSH) GO TO 100 | |
170 | IFLQ=IS2UD(IQ1) | |
171 | IF (IFLQ.EQ.1) THEN | |
172 | EQ=2./3. | |
173 | ELSE | |
174 | EQ=-1./3. | |
175 | END IF | |
176 | C Left squarks | |
177 | IF(IQ1.LE.9) THEN | |
178 | ALR=2*(ALQ(IFLQ)-BEQ(IFLQ)) | |
179 | C Right squarks | |
180 | ELSEIF(IQ1.GE.14.AND.IQ1.LE.21) THEN | |
181 | ALR=2*(ALQ(IFLQ)+BEQ(IFLQ)) | |
182 | C Mixed stops and sbottoms | |
183 | ELSEIF(IQ1.EQ.10.OR.IQ1.EQ.11) THEN | |
184 | ALR=2*(ALQ(IFLQ)-BEQ(IFLQ)*COS(2*THETAB)) | |
185 | ELSEIF(IQ1.EQ.12.OR.IQ1.EQ.13) THEN | |
186 | ALR=2*(ALQ(IFLQ)-BEQ(IFLQ)*COS(2*THETAT)) | |
187 | ELSEIF(IQ1.EQ.22.OR.IQ1.EQ.23) THEN | |
188 | ALR=2*(ALQ(IFLQ)+BEQ(IFLQ)*COS(2*THETAB)) | |
189 | ELSEIF(IQ1.EQ.24.OR.IQ1.EQ.25) THEN | |
190 | ALR=2*(ALQ(IFLQ)+BEQ(IFLQ)*COS(2*THETAT)) | |
191 | END IF | |
192 | PCM=.5*SQRT(SH-4.*AM1**2) | |
193 | IFM=ISIGN(1,IUD(JS2JT(IQ1))) | |
194 | IF (IFM.GT.0) THEN | |
195 | Z=CTH(1) | |
196 | ELSE | |
197 | Z=-CTH(1) | |
198 | END IF | |
199 | C Calculate d(sigma)/d(cos theta) in mb | |
200 | PHIZLR=2*E**4*(1.-Z**2)*(8*EQ**2/SH+(2*ALR**2*(AE-BE)**2* | |
201 | $ SH-8*(AE-BE)*EQ*ALR*(SH-AMZ**2))/PROPZ) | |
202 | PHIZRL=2*E**4*(1.-Z**2)*(8*EQ**2/SH+(2*ALR**2*(AE+BE)**2* | |
203 | $ SH-8*(AE+BE)*EQ*ALR*(SH-AMZ**2))/PROPZ) | |
204 | SIGLR=3*PCM**3/512./PI/EB**3*PHIZLR | |
205 | SIGRL=3*PCM**3/512./PI/EB**3*PHIZRL | |
206 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
207 | IF (IBREM.AND..NOT.IBEAM) THEN | |
208 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
209 | ELSE IF (IBEAM) THEN | |
210 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
211 | END IF | |
212 | CALL SIGFIL(SIG,0,0,IQ1,IQ2) | |
213 | 100 CONTINUE | |
214 | C Mixed sbottom_1 and sbottom_2 production | |
215 | IF ((AMB1SS+AMB2SS).LT.RSH) THEN | |
216 | Z=CTH(1) | |
217 | PCM=SQRT(SSXLAM(SH,AMB1SS**2,AMB2SS**2))/2./RSH | |
218 | SIGLR=2*3*8*PI*ALFAEM**2*BEQ(2)**2*COS(THETAB)**2* | |
219 | $ SIN(THETAB)**2*(AE-BE)**2*PCM**3*(1.-Z**2)/RSH/PROPZ | |
220 | SIGRL=2*3*8*PI*ALFAEM**2*BEQ(2)**2*COS(THETAB)**2* | |
221 | $ SIN(THETAB)**2*(AE+BE)**2*PCM**3*(1.-Z**2)/RSH/PROPZ | |
222 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
223 | IF (IBREM.AND..NOT.IBEAM) THEN | |
224 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
225 | ELSE IF (IBEAM) THEN | |
226 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
227 | END IF | |
228 | IF(GOQ(10,1).AND.GOQ(23,2)) THEN | |
229 | CALL SIGFIL(SIG,0,0,10,23) | |
230 | END IF | |
231 | IF(GOQ(23,1).AND.GOQ(10,2)) THEN | |
232 | CALL SIGFIL(SIG,0,0,23,10) | |
233 | END IF | |
234 | IF(GOQ(11,1).AND.GOQ(22,2)) THEN | |
235 | CALL SIGFIL(SIG,0,0,11,22) | |
236 | END IF | |
237 | IF(GOQ(22,1).AND.GOQ(11,2)) THEN | |
238 | CALL SIGFIL(SIG,0,0,22,11) | |
239 | END IF | |
240 | ENDIF | |
241 | C Mixed stop_1 and stop_2 production | |
242 | IF ((AMT1SS+AMT2SS).LT.RSH) THEN | |
243 | Z=CTH(1) | |
244 | PCM=SQRT(SSXLAM(SH,AMT1SS**2,AMT2SS**2))/2./RSH | |
245 | SIGLR=2*3*8*PI*ALFAEM**2*BEQ(1)**2*COS(THETAT)**2* | |
246 | $ SIN(THETAT)**2*(AE-BE)**2*PCM**3*(1.-Z**2)/RSH/PROPZ | |
247 | SIGRL=2*3*8*PI*ALFAEM**2*BEQ(1)**2*COS(THETAT)**2* | |
248 | $ SIN(THETAT)**2*(AE+BE)**2*PCM**3*(1.-Z**2)/RSH/PROPZ | |
249 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
250 | IF (IBREM.AND..NOT.IBEAM) THEN | |
251 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
252 | ELSE IF (IBEAM) THEN | |
253 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
254 | END IF | |
255 | IF(GOQ(12,1).AND.GOQ(25,2)) THEN | |
256 | CALL SIGFIL(SIG,0,0,12,25) | |
257 | END IF | |
258 | IF(GOQ(25,1).AND.GOQ(12,2)) THEN | |
259 | CALL SIGFIL(SIG,0,0,25,12) | |
260 | END IF | |
261 | IF(GOQ(13,1).AND.GOQ(24,2)) THEN | |
262 | CALL SIGFIL(SIG,0,0,13,24) | |
263 | END IF | |
264 | IF(GOQ(24,1).AND.GOQ(13,2)) THEN | |
265 | CALL SIGFIL(SIG,0,0,24,13) | |
266 | END IF | |
267 | ENDIF | |
268 | C | |
269 | C 2nd and 3rd generation sleptons: IL1 labels JETTYPE1. | |
270 | C | |
271 | DO 200 I=5,16 | |
272 | I1=I | |
273 | IF (I1.GE.13) I1=I1+2 | |
274 | I2=MATCHL(I1) | |
275 | IL1=IL2JS(I1) | |
276 | IL2=IL2JS(I2) | |
277 | IF(.NOT.(GOQ(IL1,1).AND.GOQ(IL2,2))) GO TO 200 | |
278 | IDL1=IDLSS(I1) | |
279 | IDL2=IDLSS(I2) | |
280 | AM1=AMASS(IDL1) | |
281 | AM2=AMASS(IDL2) | |
282 | IF((AM1+AM2).GE.RSH) GO TO 200 | |
283 | IFL1=IS2LN(I1) | |
284 | IFL2=IS2LN(I2) | |
285 | IF (IFL1.EQ.1) THEN | |
286 | EQ=0. | |
287 | ELSE | |
288 | EQ=-1. | |
289 | END IF | |
290 | IF (I1.EQ.15.OR.I1.EQ.16) THEN | |
291 | ALR=2*(ALL(IFL1)+BEL(IFL1)) | |
292 | ELSE IF (I1.GE.5.AND.I1.LE.10) THEN | |
293 | ALR=2*(ALL(IFL1)-BEL(IFL1)) | |
294 | ELSE IF (I1.EQ.11.OR.I1.EQ.12) THEN | |
295 | ALR=2*(ALL(IFL1)-BEL(IFL1)*COS(2*THETAL)) | |
296 | ELSE IF (I1.EQ.17.OR.I1.EQ.18) THEN | |
297 | ALR=2*(ALL(IFL1)+BEL(IFL1)*COS(2*THETAL)) | |
298 | END IF | |
299 | PCM=.5*SQRT(SH-4.*AM1**2) | |
300 | IFM=ISIGN(1,IDL1) | |
301 | IF (IFM.GT.0) THEN | |
302 | Z=CTH(1) | |
303 | ELSE | |
304 | Z=-CTH(1) | |
305 | END IF | |
306 | C Calculate d(sigma)/d(cos theta) in mb | |
307 | PHIZLR=2*E**4*(1.-Z**2)*(8*EQ**2/SH+(2*ALR**2*(AE-BE)**2* | |
308 | $ SH-8*(AE-BE)*EQ*ALR*(SH-AMZ**2))/PROPZ) | |
309 | PHIZRL=2*E**4*(1.-Z**2)*(8*EQ**2/SH+(2*ALR**2*(AE+BE)**2* | |
310 | $ SH-8*(AE+BE)*EQ*ALR*(SH-AMZ**2))/PROPZ) | |
311 | SIGLR=PCM**3/512./PI/EB**3*PHIZLR | |
312 | SIGRL=PCM**3/512./PI/EB**3*PHIZRL | |
313 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
314 | IF (IBREM.AND..NOT.IBEAM) THEN | |
315 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
316 | ELSE IF (IBEAM) THEN | |
317 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
318 | END IF | |
319 | CALL SIGFIL(SIG,0,0,IL1,IL2) | |
320 | 200 CONTINUE | |
321 | C Mixed stau_1 and stau_2 production | |
322 | IF ((AML1SS+AML2SS).LT.RSH) THEN | |
323 | Z=CTH(1) | |
324 | PCM=SQRT(SSXLAM(SH,AML1SS**2,AML2SS**2))/2./RSH | |
325 | SIGLR=2*8*PI*ALFAEM**2*BEL(2)**2*COS(THETAL)**2* | |
326 | $ SIN(THETAL)**2*(AE-BE)**2*PCM**3*(1.-Z**2)/RSH/PROPZ | |
327 | SIGRL=2*8*PI*ALFAEM**2*BEL(2)**2*COS(THETAL)**2* | |
328 | $ SIN(THETAL)**2*(AE+BE)**2*PCM**3*(1.-Z**2)/RSH/PROPZ | |
329 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
330 | IF (IBREM.AND..NOT.IBEAM) THEN | |
331 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
332 | ELSE IF (IBEAM) THEN | |
333 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
334 | END IF | |
335 | IF(GOQ(44,1).AND.GOQ(51,2)) THEN | |
336 | CALL SIGFIL(SIG,0,0,44,51) | |
337 | END IF | |
338 | IF(GOQ(51,1).AND.GOQ(44,2)) THEN | |
339 | CALL SIGFIL(SIG,0,0,51,44) | |
340 | END IF | |
341 | IF(GOQ(45,1).AND.GOQ(50,2)) THEN | |
342 | CALL SIGFIL(SIG,0,0,45,50) | |
343 | END IF | |
344 | IF(GOQ(50,1).AND.GOQ(45,2)) THEN | |
345 | CALL SIGFIL(SIG,0,0,50,45) | |
346 | END IF | |
347 | ENDIF | |
348 | C | |
349 | C Next do 1st generation sleptons | |
350 | C | |
351 | C Sneutrino_e pairs | |
352 | DO 210 I1=1,2 | |
353 | I2=MATCHL(I1) | |
354 | IL1=IL2JS(I1) | |
355 | IL2=IL2JS(I2) | |
356 | IF(.NOT.(GOQ(IL1,1).AND.GOQ(IL2,2))) GO TO 210 | |
357 | MSNE=AMASS(ISNEL) | |
358 | IF((2*MSNE).GE.RSH) GO TO 210 | |
359 | IF (I1.EQ.1) THEN | |
360 | Z=CTH(1) | |
361 | ELSE | |
362 | Z=-CTH(1) | |
363 | END IF | |
364 | PCM=.5*SQRT(SH-4*MSNE**2) | |
365 | TM1LR=32*E**4*(AN-BN)**2*(AE-BE)**2/PROPZ | |
366 | TM1RL=32*E**4*(AN-BN)**2*(AE+BE)**2/PROPZ | |
367 | TM2=8*G**4*SIN(GAMMAR)**4/(2*EB*(EB-PCM*Z)+AMW1SS**2-MSNE**2)**2 | |
368 | TM3=8*G**4*COS(GAMMAR)**4/(2*EB*(EB-PCM*Z)+AMW2SS**2-MSNE**2)**2 | |
369 | TM4=-32*E**2*(AN-BN)*G**2*SIN(GAMMAR)**2*(SH-AMZ**2)*(AE-BE)/ | |
370 | $ PROPZ/(2*EB*(EB-PCM*Z)+AMW1SS**2-MSNE**2) | |
371 | TM5=-32*E**2*(AN-BN)*G**2*COS(GAMMAR)**2*(SH-AMZ**2)*(AE-BE)/ | |
372 | $ PROPZ/(2*EB*(EB-PCM*Z)+AMW2SS**2-MSNE**2) | |
373 | TM6=16*G**4*SIN(GAMMAR)**2*COS(GAMMAR)**2/ | |
374 | $ (2*EB*(EB-PCM*Z)+AMW1SS**2-MSNE**2)/ | |
375 | $ (2*EB*(EB-PCM*Z)+AMW2SS**2-MSNE**2) | |
376 | SIGLR=2*PCM**3*EB*(1.-Z**2)/128./PI/SH* | |
377 | $ (TM1LR+TM2+TM3+TM4+TM5+TM6) | |
378 | SIGRL=2*PCM**3*EB*(1.-Z**2)/128./PI/SH*TM1RL | |
379 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
380 | IF (IBREM.AND..NOT.IBEAM) THEN | |
381 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
382 | ELSE IF (IBEAM) THEN | |
383 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
384 | END IF | |
385 | CALL SIGFIL(SIG,0,0,IL1,IL2) | |
386 | 210 CONTINUE | |
387 | C E_L~ pairs | |
388 | DO 220 I1=3,4 | |
389 | I2=MATCHL(I1) | |
390 | IL1=IL2JS(I1) | |
391 | IL2=IL2JS(I2) | |
392 | IF(.NOT.(GOQ(IL1,1).AND.GOQ(IL2,2))) GO TO 220 | |
393 | IF(2*MEL.GE.RSH) GO TO 220 | |
394 | PCM=.5*SQRT(SH-4.*MEL**2) | |
395 | EQ=-1. | |
396 | ALR=2*(AE-BE) | |
397 | IF (I1.EQ.3) THEN | |
398 | Z=CTH(1) | |
399 | ELSE | |
400 | Z=-CTH(1) | |
401 | END IF | |
402 | PHIZLR=E**4*(1.-Z**2)*(8*EQ**2/SH+(2*ALR**2*(AE-BE)**2* | |
403 | $ SH-8*(AE-BE)*EQ*ALR*(SH-AMZ**2))/PROPZ) | |
404 | PHIZRL=E**4*(1.-Z**2)*(8*EQ**2/SH+(2*ALR**2*(AE+BE)**2* | |
405 | $ SH-8*(AE+BE)*EQ*ALR*(SH-AMZ**2))/PROPZ) | |
406 | DO 221 IZ1=1,4 | |
407 | AEZS=AEZ(IZ1)*CONJG(AEZ(IZ1)) | |
408 | PHIZLR=PHIZLR+2*AEZS**2*SH*(1.-Z**2)/(2*EB*(EB-PCM*Z)- | |
409 | $ MEL**2+AMZISS(IZ1)**2)**2-4*E**2*(1.-Z**2)*AEZS/ | |
410 | $ (2*EB*(EB-PCM*Z)-MEL**2+AMZISS(IZ1)**2)*(2.+(AE-BE)*ALR* | |
411 | $ SH*(SH-AMZ**2)/PROPZ) | |
412 | IF (IZ1.LE.3) THEN | |
413 | DO 222 IP=IZ1+1,4 | |
414 | AEZJS=AEZ(IP)*CONJG(AEZ(IP)) | |
415 | PHIZLR=PHIZLR+4*AEZS*AEZJS*SH*(1.-Z**2)/ | |
416 | $ (2*EB*(EB-PCM*Z)-MEL**2+AMZISS(IZ1)**2)/ | |
417 | $ (2*EB*(EB-PCM*Z)-MEL**2+AMZISS(IP)**2) | |
418 | 222 CONTINUE | |
419 | END IF | |
420 | 221 CONTINUE | |
421 | SIGLR=2*PCM**3/512./PI/EB**3*PHIZLR | |
422 | SIGRL=2*PCM**3/512./PI/EB**3*PHIZRL | |
423 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
424 | IF (IBREM.AND..NOT.IBEAM) THEN | |
425 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
426 | ELSE IF (IBEAM) THEN | |
427 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
428 | END IF | |
429 | CALL SIGFIL(SIG,0,0,IL1,IL2) | |
430 | 220 CONTINUE | |
431 | C E_R~ pairs | |
432 | DO 230 I1=13,14 | |
433 | I2=MATCHL(I1) | |
434 | IL1=IL2JS(I1) | |
435 | IL2=IL2JS(I2) | |
436 | IF(.NOT.(GOQ(IL1,1).AND.GOQ(IL2,2))) GO TO 230 | |
437 | IF(2*MER.GE.RSH) GO TO 230 | |
438 | PCM=.5*SQRT(SH-4.*MER**2) | |
439 | EQ=-1. | |
440 | ALR=2*(AE+BE) | |
441 | IF (I1.EQ.13) THEN | |
442 | Z=CTH(1) | |
443 | ELSE | |
444 | Z=-CTH(1) | |
445 | END IF | |
446 | PHIZLR=E**4*(1.-Z**2)*(8*EQ**2/SH+(2*ALR**2*(AE-BE)**2* | |
447 | $ SH-8*(AE-BE)*EQ*ALR*(SH-AMZ**2))/PROPZ) | |
448 | PHIZRL=E**4*(1.-Z**2)*(8*EQ**2/SH+(2*ALR**2*(AE+BE)**2* | |
449 | $ SH-8*(AE+BE)*EQ*ALR*(SH-AMZ**2))/PROPZ) | |
450 | DO 231 IZ1=1,4 | |
451 | BEZS=BEZ(IZ1)*CONJG(BEZ(IZ1)) | |
452 | PHIZRL=PHIZRL+2*BEZS**2*SH*(1.-Z**2)/(2*EB*(EB-PCM*Z)- | |
453 | $ MER**2+AMZISS(IZ1)**2)**2-4*E**2*(1.-Z**2)*BEZS/ | |
454 | $ (2*EB*(EB-PCM*Z)-MER**2+AMZISS(IZ1)**2)*(2.+(AE+BE)*ALR* | |
455 | $ SH*(SH-AMZ**2)/PROPZ) | |
456 | IF (IZ1.LE.3) THEN | |
457 | DO 232 IP=IZ1+1,4 | |
458 | BEZJS=BEZ(IP)*CONJG(BEZ(IP)) | |
459 | PHIZRL=PHIZRL+4*BEZS*BEZJS*SH*(1.-Z**2)/ | |
460 | $ (2*EB*(EB-PCM*Z)-MER**2+AMZISS(IZ1)**2)/ | |
461 | $ (2*EB*(EB-PCM*Z)-MER**2+AMZISS(IP)**2) | |
462 | 232 CONTINUE | |
463 | END IF | |
464 | 231 CONTINUE | |
465 | SIGLR=2*PCM**3/512./PI/EB**3*PHIZLR | |
466 | SIGRL=2*PCM**3/512./PI/EB**3*PHIZRL | |
467 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
468 | IF (IBREM.AND..NOT.IBEAM) THEN | |
469 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
470 | ELSE IF (IBEAM) THEN | |
471 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
472 | END IF | |
473 | CALL SIGFIL(SIG,0,0,IL1,IL2) | |
474 | 230 CONTINUE | |
475 | C E_L~+E_R~bar and E_R~+E_L~bar pairs; now has MEL =/ MER ! | |
476 | IF((MEL+MER).GE.RSH) GO TO 270 | |
477 | IF(GOQ(36,1).AND.GOQ(47,2)) THEN | |
478 | PCM=SQRT(SSXLAM(SH,MEL**2,MER**2))/4./EB | |
479 | EEL=SQRT(PCM**2+MEL**2) | |
480 | Z=CTH(1) | |
481 | PHIZ=0. | |
482 | DO 241 IZ1=1,4 | |
483 | BEZS=BEZ(IZ1)*CONJG(BEZ(IZ1)) | |
484 | AEZS=AEZ(IZ1)*CONJG(AEZ(IZ1)) | |
485 | AZI=(AMZISS(IZ1)**2-MEL**2)/2./EB | |
486 | PHIZ=PHIZ+AEZS*BEZS*AMZISS(IZ1)**2/(EEL-PCM*Z+AZI)**2 | |
487 | IF (IZ1.LE.3) THEN | |
488 | DO 242 IP=IZ1+1,4 | |
489 | AZJ=(AMZISS(IP)**2-MEL**2)/2./EB | |
490 | PHIZ=PHIZ+2*ABS(AMZISS(IZ1)*AMZISS(IP))* | |
491 | $ REAL(AEZ(IZ1)*CONJG(AEZ(IP))*CONJG(BEZ(IZ1))*BEZ(IP))/ | |
492 | $ (EEL-PCM*Z+AZI)/(EEL-PCM*Z+AZJ) | |
493 | 242 CONTINUE | |
494 | END IF | |
495 | 241 CONTINUE | |
496 | SIG=4*PCM/128./PI/SH/EB*PHIZ | |
497 | SIG=FLEM*FLEP*SIG*UNITS/2. | |
498 | IF (IBREM.AND..NOT.IBEAM) THEN | |
499 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
500 | ELSE IF (IBEAM) THEN | |
501 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
502 | END IF | |
503 | CALL SIGFIL(SIG,0,0,36,47) | |
504 | ENDIF | |
505 | IF(GOQ(46,1).AND.GOQ(37,2)) THEN | |
506 | PCM=SQRT(SSXLAM(SH,MEL**2,MER**2))/4./EB | |
507 | EER=SQRT(PCM**2+MER**2) | |
508 | Z=CTH(1) | |
509 | PHIZ=0. | |
510 | DO 243 IZ1=1,4 | |
511 | BEZS=BEZ(IZ1)*CONJG(BEZ(IZ1)) | |
512 | AEZS=AEZ(IZ1)*CONJG(AEZ(IZ1)) | |
513 | AZI=(AMZISS(IZ1)**2-MER**2)/2./EB | |
514 | PHIZ=PHIZ+AEZS*BEZS*AMZISS(IZ1)**2/(EER-PCM*Z+AZI)**2 | |
515 | IF (IZ1.LE.3) THEN | |
516 | DO 244 IP=IZ1+1,4 | |
517 | AZJ=(AMZISS(IP)**2-MER**2)/2./EB | |
518 | PHIZ=PHIZ+2*ABS(AMZISS(IZ1)*AMZISS(IP))* | |
519 | $ REAL(AEZ(IZ1)*CONJG(AEZ(IP))*CONJG(BEZ(IZ1))*BEZ(IP))/ | |
520 | $ (EER-PCM*Z+AZI)/(EER-PCM*Z+AZJ) | |
521 | 244 CONTINUE | |
522 | END IF | |
523 | 243 CONTINUE | |
524 | SIG=4*PCM/128./PI/SH/EB*PHIZ | |
525 | SIG=FREM*FREP*SIG*UNITS/2. | |
526 | IF (IBREM.AND..NOT.IBEAM) THEN | |
527 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
528 | ELSE IF (IBEAM) THEN | |
529 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
530 | END IF | |
531 | CALL SIGFIL(SIG,0,0,46,37) | |
532 | ENDIF | |
533 | C E_R~bar+E_L~ and E_L~bar+E_R~ pairs; now assumes MEL =/ MER ! | |
534 | IF(GOQ(47,1).AND.GOQ(36,2)) THEN | |
535 | PCM=SQRT(SSXLAM(SH,MEL**2,MER**2))/4./EB | |
536 | EEL=SQRT(PCM**2+MEL**2) | |
537 | Z=-CTH(1) | |
538 | PHIZ=0. | |
539 | DO 251 IZ1=1,4 | |
540 | BEZS=BEZ(IZ1)*CONJG(BEZ(IZ1)) | |
541 | AEZS=AEZ(IZ1)*CONJG(AEZ(IZ1)) | |
542 | AZI=(AMZISS(IZ1)**2-MEL**2)/2./EB | |
543 | PHIZ=PHIZ+AEZS*BEZS*AMZISS(IZ1)**2/(EEL-PCM*Z+AZI)**2 | |
544 | IF (IZ1.LE.3) THEN | |
545 | DO 252 IP=IZ1+1,4 | |
546 | AZJ=(AMZISS(IP)**2-MEL**2)/2./EB | |
547 | PHIZ=PHIZ+2*ABS(AMZISS(IZ1)*AMZISS(IP))* | |
548 | $ REAL(AEZ(IZ1)*CONJG(AEZ(IP))*CONJG(BEZ(IZ1))*BEZ(IP))/ | |
549 | $ (EEL-PCM*Z+AZI)/(EEL-PCM*Z+AZJ) | |
550 | 252 CONTINUE | |
551 | END IF | |
552 | 251 CONTINUE | |
553 | SIG=4*PCM/128./PI/SH/EB*PHIZ | |
554 | SIG=FLEM*FLEP*SIG*UNITS/2. | |
555 | IF (IBREM.AND..NOT.IBEAM) THEN | |
556 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
557 | ELSE IF (IBEAM) THEN | |
558 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
559 | END IF | |
560 | CALL SIGFIL(SIG,0,0,47,36) | |
561 | ENDIF | |
562 | IF(GOQ(37,1).AND.GOQ(46,2)) THEN | |
563 | PCM=SQRT(SSXLAM(SH,MEL**2,MER**2))/4./EB | |
564 | EER=SQRT(PCM**2+MER**2) | |
565 | Z=-CTH(1) | |
566 | PHIZ=0. | |
567 | DO 253 IZ1=1,4 | |
568 | BEZS=BEZ(IZ1)*CONJG(BEZ(IZ1)) | |
569 | AEZS=AEZ(IZ1)*CONJG(AEZ(IZ1)) | |
570 | AZI=(AMZISS(IZ1)**2-MER**2)/2./EB | |
571 | PHIZ=PHIZ+AEZS*BEZS*AMZISS(IZ1)**2/(EER-PCM*Z+AZI)**2 | |
572 | IF (IZ1.LE.3) THEN | |
573 | DO 254 IP=IZ1+1,4 | |
574 | AZJ=(AMZISS(IP)**2-MER**2)/2./EB | |
575 | PHIZ=PHIZ+2*ABS(AMZISS(IZ1)*AMZISS(IP))* | |
576 | $ REAL(AEZ(IZ1)*CONJG(AEZ(IP))*CONJG(BEZ(IZ1))*BEZ(IP))/ | |
577 | $ (EER-PCM*Z+AZI)/(EER-PCM*Z+AZJ) | |
578 | 254 CONTINUE | |
579 | END IF | |
580 | 253 CONTINUE | |
581 | SIG=4*PCM/128./PI/SH/EB*PHIZ | |
582 | SIG=FREM*FREP*SIG*UNITS/2. | |
583 | IF (IBREM.AND..NOT.IBEAM) THEN | |
584 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
585 | ELSE IF (IBEAM) THEN | |
586 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
587 | END IF | |
588 | CALL SIGFIL(SIG,0,0,37,46) | |
589 | ENDIF | |
590 | 270 CONTINUE | |
591 | C | |
592 | C Chargino pair production | |
593 | C | |
594 | DO 300 IW1=1,4 | |
595 | JW1=(IW1+1)/2 | |
596 | AMWI=ABS(AMWISS(JW1)) | |
597 | JTW1=IW2JS(IW1) | |
598 | JTW2=IW2JS(MATCHL(IW1)) | |
599 | IF (.NOT.(GOQ(JTW1,1).AND.GOQ(JTW2,2))) GO TO 300 | |
600 | IF((2*AMWI).GE.RSH) GO TO 300 | |
601 | PCM=SQRT(SSXLAM(SH,AMWI**2,AMWI**2))/4./EB | |
602 | Z=CTH(1) | |
603 | IF (IW1.EQ.1.OR.IW1.EQ.3) Z=-CTH(1) | |
604 | SINGR=SIN(GAMMAR) | |
605 | COSGR=COS(GAMMAR) | |
606 | SINGL=SIN(GAMMAL) | |
607 | COSGL=COS(GAMMAL) | |
608 | XC=1.-(COSGL**2+COSGR**2)/4./CS2THW | |
609 | YC=(COSGR**2-COSGL**2)/4./CS2THW | |
610 | XS=1.-(SINGL**2+SINGR**2)/4./CS2THW | |
611 | YS=(SINGR**2-SINGL**2)/4./CS2THW | |
612 | IF (IW1.GE.3) THEN | |
613 | XC=XS | |
614 | YC=YS | |
615 | SINGR=COSGR | |
616 | END IF | |
617 | TGG=16*E**4/SH*(EB**2*(1.+Z**2)+AMWI**2*(1.-Z**2)) | |
618 | TZZLR=16*E**4*CTTHW**2*SH/PROPZ*((XC**2+YC**2)*(AE-BE)**2* | |
619 | $ (EB**2*(1.+Z**2)+AMWI**2*(1.-Z**2))- | |
620 | $ 2*YC**2*(AE-BE)**2*AMWI**2+4*XC*YC*(AE-BE)**2*EB*PCM*Z) | |
621 | TZZRL=16*E**4*CTTHW**2*SH/PROPZ*((XC**2+YC**2)*(AE+BE)**2* | |
622 | $ (EB**2*(1.+Z**2)+AMWI**2*(1.-Z**2))- | |
623 | $ 2*YC**2*(AE+BE)**2*AMWI**2-4*XC*YC*(AE+BE)**2*EB*PCM*Z) | |
624 | TGZLR=-32*E**4*CTTHW*(SH-AMZ**2)/PROPZ*((AE-BE)*XC* | |
625 | $ (EB**2*(1.+Z**2)+AMWI**2*(1.-Z**2))-2*(BE-AE)*YC*EB*PCM*Z) | |
626 | TGZRL=-32*E**4*CTTHW*(SH-AMZ**2)/PROPZ*((AE+BE)*XC* | |
627 | $ (EB**2*(1.+Z**2)+AMWI**2*(1.-Z**2))-2*(BE+AE)*YC*EB*PCM*Z) | |
628 | TNN=2*E**4*SINGR**4*SH*(EB-PCM*Z)**2/SN2THW**2/ | |
629 | $ (EB**2+PCM**2-2*EB*PCM*Z+MSNE**2)**2 | |
630 | TGN=-8*E**4*SINGR**2*((EB-PCM*Z)**2+AMWI**2)/SN2THW/ | |
631 | $ (EB**2+PCM**2-2*EB*PCM*Z+MSNE**2) | |
632 | TZN=8*E**4*CTTHW*SINGR**2*(SH-AMZ**2)*(AE-BE)*SH/ | |
633 | $ SN2THW/PROPZ*((XC-YC)*((EB-PCM*Z)**2+AMWI**2)+2*YC*AMWI**2)/ | |
634 | $ (EB**2+PCM**2-2*EB*PCM*Z+MSNE**2) | |
635 | SIGLR=2*PCM/128./PI/SH/EB*(TGG+TZZLR+TGZLR+TNN+TGN+TZN) | |
636 | SIGRL=2*PCM/128./PI/SH/EB*(TGG+TZZRL+TGZRL) | |
637 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
638 | IF (IBREM.AND..NOT.IBEAM) THEN | |
639 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
640 | ELSE IF (IBEAM) THEN | |
641 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
642 | END IF | |
643 | CALL SIGFIL(SIG,0,0,JTW1,JTW2) | |
644 | 300 CONTINUE | |
645 | C | |
646 | C Chargino_1 + chargino_2 pair production | |
647 | IF((ABS(AMW1SS)+ABS(AMW2SS)).GE.RSH) GO TO 340 | |
648 | PCM=SQRT(SSXLAM(SH,AMW1SS**2,AMW2SS**2))/4./EB | |
649 | XC=(THX*SIN(GAMMAL)*COS(GAMMAL)-THY*SIN(GAMMAR)*COS(GAMMAR))/2. | |
650 | YC=(THX*SIN(GAMMAL)*COS(GAMMAL)+THY*SIN(GAMMAR)*COS(GAMMAR))/2. | |
651 | DEL=(AMW2SS**2-AMW1SS**2)/4./EB | |
652 | XI=-1.*SIGN(1.,AMWISS(1))*SIGN(1.,AMWISS(2)) | |
653 | IF (.NOT.(GOQ(27,1).AND.GOQ(28,2))) GO TO 310 | |
654 | Z=CTH(1) | |
655 | TZZLR=4*(CTTHW+TNTHW)**2/PROPZ*((XC**2+YC**2)*(AE-BE)**2* | |
656 | $ (EB**2+PCM**2*Z**2-DEL**2-XI*ABS(AMW1SS*AMW2SS))+ | |
657 | $ 2*XC**2*XI*(AE-BE)**2*ABS(AMW1SS*AMW2SS)+ | |
658 | $ 4*XC*YC*(AE-BE)**2*EB*PCM*Z) | |
659 | TZZRL=4*(CTTHW+TNTHW)**2/PROPZ*((XC**2+YC**2)*(AE+BE)**2* | |
660 | $ (EB**2+PCM**2*Z**2-DEL**2-XI*ABS(AMW1SS*AMW2SS))+ | |
661 | $ 2*XC**2*XI*(AE+BE)**2*ABS(AMW1SS*AMW2SS)- | |
662 | $ 4*XC*YC*(AE+BE)**2*EB*PCM*Z) | |
663 | TNN=2*SIN(GAMMAR)**2*COS(GAMMAR)**2*((EB-PCM*Z)**2-DEL**2)/ | |
664 | $ SN2THW**2/(2*EB*(EB-DEL)-2*EB*PCM*Z+MSNE**2-AMW1SS**2)**2 | |
665 | TZN=-4*THY*(CTTHW+TNTHW)*SIN(GAMMAR)*COS(GAMMAR)*(SH-AMZ**2) | |
666 | $ *(AE-BE)/SN2THW/PROPZ*((XC-YC)*((EB-PCM*Z)**2-DEL**2- | |
667 | $ XI*ABS(AMW1SS*AMW2SS))+2*XC*XI*ABS(AMW1SS*AMW2SS))/ | |
668 | $ (2*EB*(EB-DEL)-2*EB*PCM*Z+MSNE**2-AMW1SS**2) | |
669 | SIGLR=2*E**4*PCM/128./PI/EB*(TZZLR+TNN+TZN) | |
670 | SIGRL=2*E**4*PCM/128./PI/EB*TZZRL | |
671 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
672 | IF (IBREM.AND..NOT.IBEAM) THEN | |
673 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
674 | ELSE IF (IBEAM) THEN | |
675 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
676 | END IF | |
677 | CALL SIGFIL(SIG,0,0,27,28) | |
678 | 310 CONTINUE | |
679 | IF (.NOT.(GOQ(28,1).AND.GOQ(27,2))) GO TO 320 | |
680 | Z=-CTH(1) | |
681 | TZZLR=4*(CTTHW+TNTHW)**2/PROPZ*((XC**2+YC**2)*(AE-BE)**2* | |
682 | $ (EB**2+PCM**2*Z**2-DEL**2-XI*ABS(AMW1SS*AMW2SS))+ | |
683 | $ 2*XC**2*XI*(AE-BE)**2*ABS(AMW1SS*AMW2SS)+ | |
684 | $ 4*XC*YC*(AE-BE)**2*EB*PCM*Z) | |
685 | TZZRL=4*(CTTHW+TNTHW)**2/PROPZ*((XC**2+YC**2)*(AE+BE)**2* | |
686 | $ (EB**2+PCM**2*Z**2-DEL**2-XI*ABS(AMW1SS*AMW2SS))+ | |
687 | $ 2*XC**2*XI*(AE+BE)**2*ABS(AMW1SS*AMW2SS)- | |
688 | $ 4*XC*YC*(AE+BE)**2*EB*PCM*Z) | |
689 | TNN=2*SIN(GAMMAR)**2*COS(GAMMAR)**2*((EB-PCM*Z)**2-DEL**2)/ | |
690 | $ SN2THW**2/(2*EB*(EB-DEL)-2*EB*PCM*Z+MSNE**2-AMW1SS**2)**2 | |
691 | TZN=-4*THY*(CTTHW+TNTHW)*SIN(GAMMAR)*COS(GAMMAR)*(SH-AMZ**2) | |
692 | $ *(AE-BE)/SN2THW/PROPZ*((XC-YC)*((EB-PCM*Z)**2-DEL**2- | |
693 | $ XI*ABS(AMW1SS*AMW2SS))+2*XC*XI*ABS(AMW1SS*AMW2SS))/ | |
694 | $ (2*EB*(EB-DEL)-2*EB*PCM*Z+MSNE**2-AMW1SS**2) | |
695 | SIGLR=2*E**4*PCM/128./PI/EB*(TZZLR+TNN+TZN) | |
696 | SIGRL=2*E**4*PCM/128./PI/EB*TZZRL | |
697 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
698 | IF (IBREM.AND..NOT.IBEAM) THEN | |
699 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
700 | ELSE IF (IBEAM) THEN | |
701 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
702 | END IF | |
703 | CALL SIGFIL(SIG,0,0,28,27) | |
704 | 320 CONTINUE | |
705 | IF (.NOT.(GOQ(29,1).AND.GOQ(26,2))) GO TO 330 | |
706 | Z=CTH(1) | |
707 | TZZLR=4*(CTTHW+TNTHW)**2/PROPZ*((XC**2+YC**2)*(AE-BE)**2* | |
708 | $ (EB**2+PCM**2*Z**2-DEL**2-XI*ABS(AMW1SS*AMW2SS))+ | |
709 | $ 2*XC**2*XI*(AE-BE)**2*ABS(AMW1SS*AMW2SS)+ | |
710 | $ 4*XC*YC*(AE-BE)**2*EB*PCM*Z) | |
711 | TZZRL=4*(CTTHW+TNTHW)**2/PROPZ*((XC**2+YC**2)*(AE+BE)**2* | |
712 | $ (EB**2+PCM**2*Z**2-DEL**2-XI*ABS(AMW1SS*AMW2SS))+ | |
713 | $ 2*XC**2*XI*(AE+BE)**2*ABS(AMW1SS*AMW2SS)- | |
714 | $ 4*XC*YC*(AE+BE)**2*EB*PCM*Z) | |
715 | TNN=2*SIN(GAMMAR)**2*COS(GAMMAR)**2*((EB-PCM*Z)**2-DEL**2)/ | |
716 | $ SN2THW**2/(2*EB*(EB-DEL)-2*EB*PCM*Z+MSNE**2-AMW1SS**2)**2 | |
717 | TZN=-4*THY*(CTTHW+TNTHW)*SIN(GAMMAR)*COS(GAMMAR)*(SH-AMZ**2) | |
718 | $ *(AE-BE)/SN2THW/PROPZ*((XC-YC)*((EB-PCM*Z)**2-DEL**2- | |
719 | $ XI*ABS(AMW1SS*AMW2SS))+2*XC*XI*ABS(AMW1SS*AMW2SS))/ | |
720 | $ (2*EB*(EB-DEL)-2*EB*PCM*Z+MSNE**2-AMW1SS**2) | |
721 | SIGLR=2*E**4*PCM/128./PI/EB*(TZZLR+TNN+TZN) | |
722 | SIGRL=2*E**4*PCM/128./PI/EB*TZZRL | |
723 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
724 | IF (IBREM.AND..NOT.IBEAM) THEN | |
725 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
726 | ELSE IF (IBEAM) THEN | |
727 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
728 | END IF | |
729 | CALL SIGFIL(SIG,0,0,29,26) | |
730 | 330 CONTINUE | |
731 | IF (.NOT.(GOQ(26,1).AND.GOQ(29,2))) GO TO 340 | |
732 | Z=-CTH(1) | |
733 | TZZLR=4*(CTTHW+TNTHW)**2/PROPZ*((XC**2+YC**2)*(AE-BE)**2* | |
734 | $ (EB**2+PCM**2*Z**2-DEL**2-XI*ABS(AMW1SS*AMW2SS))+ | |
735 | $ 2*XC**2*XI*(AE-BE)**2*ABS(AMW1SS*AMW2SS)+ | |
736 | $ 4*XC*YC*(AE-BE)**2*EB*PCM*Z) | |
737 | TZZRL=4*(CTTHW+TNTHW)**2/PROPZ*((XC**2+YC**2)*(AE+BE)**2* | |
738 | $ (EB**2+PCM**2*Z**2-DEL**2-XI*ABS(AMW1SS*AMW2SS))+ | |
739 | $ 2*XC**2*XI*(AE+BE)**2*ABS(AMW1SS*AMW2SS)- | |
740 | $ 4*XC*YC*(AE+BE)**2*EB*PCM*Z) | |
741 | TNN=2*SIN(GAMMAR)**2*COS(GAMMAR)**2*((EB-PCM*Z)**2-DEL**2)/ | |
742 | $ SN2THW**2/(2*EB*(EB-DEL)-2*EB*PCM*Z+MSNE**2-AMW1SS**2)**2 | |
743 | TZN=-4*THY*(CTTHW+TNTHW)*SIN(GAMMAR)*COS(GAMMAR)*(SH-AMZ**2) | |
744 | $ *(AE-BE)/SN2THW/PROPZ*((XC-YC)*((EB-PCM*Z)**2-DEL**2- | |
745 | $ XI*ABS(AMW1SS*AMW2SS))+2*XC*XI*ABS(AMW1SS*AMW2SS))/ | |
746 | $ (2*EB*(EB-DEL)-2*EB*PCM*Z+MSNE**2-AMW1SS**2) | |
747 | SIGLR=2*E**4*PCM/128./PI/EB*(TZZLR+TNN+TZN) | |
748 | SIGRL=2*E**4*PCM/128./PI/EB*TZZRL | |
749 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
750 | IF (IBREM.AND..NOT.IBEAM) THEN | |
751 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
752 | ELSE IF (IBEAM) THEN | |
753 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
754 | END IF | |
755 | CALL SIGFIL(SIG,0,0,26,29) | |
756 | 340 CONTINUE | |
757 | C | |
758 | C Neutralino pair production | |
759 | C | |
760 | DO 400 IZ1=1,4 | |
761 | AMZIZ1=ABS(AMZISS(IZ1)) | |
762 | JTYPZ1=IZ2JS(IZ1) | |
763 | DO 410 IZ2=1,4 | |
764 | AMZIZ2=ABS(AMZISS(IZ2)) | |
765 | JTYPZ2=IZ2JS(IZ2) | |
766 | IF(.NOT.(GOQ(JTYPZ1,1).AND.GOQ(JTYPZ2,2))) GO TO 410 | |
767 | IF((AMZIZ1+AMZIZ2).GE.RSH) GO TO 410 | |
768 | WIJ=SQRT(G**2+GP**2)*ZI**(ITHZ(IZ2))*(-ZI)**(ITHZ(IZ1))* | |
769 | $ (ZMIXSS(1,IZ1)*ZMIXSS(1,IZ2)-ZMIXSS(2,IZ1)* | |
770 | $ ZMIXSS(2,IZ2))/4. | |
771 | KK=SQRT(SH*SH+(AMZIZ1**2-AMZIZ2**2)**2-2*SH* | |
772 | $ (AMZIZ1**2+AMZIZ2**2))/4./EB | |
773 | Z=CTH(1) | |
774 | SIGLL=2*AEZ(IZ1)*CONJG(AEZ(IZ1))*AEZ(IZ2)*CONJG(AEZ(IZ2))* | |
775 | $ SSGT(SH,MEL,Z,IZ1,IZ2) | |
776 | SIGRR=2*BEZ(IZ1)*CONJG(BEZ(IZ1))*BEZ(IZ2)*CONJG(BEZ(IZ2))* | |
777 | $ SSGT(SH,MER,Z,IZ1,IZ2) | |
778 | SIGZZL=4*E**2*WIJ*CONJG(WIJ)*(AE-BE)**2* | |
779 | $ (SH*SH-(AMZIZ1**2-AMZIZ2**2)**2+4*(-1.)**(ITHZ(IZ1)+ | |
780 | $ ITHZ(IZ2)+1)*SH*AMZIZ1*AMZIZ2+4*SH*KK*KK*Z*Z)/PROPZ | |
781 | SIGZZR=4*E**2*WIJ*CONJG(WIJ)*(AE+BE)**2* | |
782 | $ (SH*SH-(AMZIZ1**2-AMZIZ2**2)**2+4*(-1.)**(ITHZ(IZ1)+ | |
783 | $ ITHZ(IZ2)+1)*SH*AMZIZ1*AMZIZ2+4*SH*KK*KK*Z*Z)/PROPZ | |
784 | SIGLZ=-E*(AE-BE)*(SH-AMZ**2)/2./PROPZ* | |
785 | $ (REAL(WIJ*CONJG(AEZ(IZ1))*AEZ(IZ2))* | |
786 | $ SSGST(SH,MEL,Z,IZ1,IZ2)+(-1.)**(ITHZ(IZ1)+ITHZ(IZ2))* | |
787 | $ REAL(WIJ*AEZ(IZ1)*CONJG(AEZ(IZ2)))* | |
788 | $ SSGST(SH,MEL,-Z,IZ1,IZ2)) | |
789 | SIGRZ=-E*(-1.)**(ITHZ(IZ1)+ITHZ(IZ2)+1)* | |
790 | $ (AE+BE)*(SH-AMZ**2)/2./PROPZ* | |
791 | $ (REAL(WIJ*CONJG(BEZ(IZ1))*BEZ(IZ2))* | |
792 | $ SSGST(SH,MER,Z,IZ1,IZ2)+(-1.)**(ITHZ(IZ1)+ITHZ(IZ2))* | |
793 | $ REAL(WIJ*BEZ(IZ1)*CONJG(BEZ(IZ2)))* | |
794 | $ SSGST(SH,MER,-Z,IZ1,IZ2)) | |
795 | SIGLR=2*KK/16./PI/SH/SQRT(SH)*(SIGLL+SIGZZL+SIGLZ) | |
796 | SIGRL=2*KK/16./PI/SH/SQRT(SH)*(SIGRR+SIGZZR+SIGRZ) | |
797 | C BELOW FACTOR OF 2 FOR ID PARTICLES AND JETTYP SWITCH | |
798 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
799 | IF (IBREM.AND..NOT.IBEAM) THEN | |
800 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
801 | ELSE IF (IBEAM) THEN | |
802 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
803 | END IF | |
804 | CALL SIGFIL(SIG,0,0,JTYPZ1,JTYPZ2) | |
805 | 410 CONTINUE | |
806 | 400 CONTINUE | |
807 | C | |
808 | C Higgs boson mechanisms | |
809 | C | |
810 | C E+ E- --> Z H_L; symmetric in cos(theta) | |
811 | IF((AMZ+AMHL).LT.RSH) THEN | |
812 | FACLR=E**2*G**2*(SIN(ALFAH+BETA))**2*(AE-BE)**2/CS2THW | |
813 | FACRL=E**2*G**2*(SIN(ALFAH+BETA))**2*(AE+BE)**2/CS2THW | |
814 | Z=CTH(1) | |
815 | PCM=SQRT(SSXLAM(SH,AMZ**2,AMHL**2))/4./EB | |
816 | EZ0=SQRT(PCM**2+AMZ**2) | |
817 | FAC1=AMZ**2+EZ0**2-PCM**2*Z**2 | |
818 | SIGLR=2*FACLR/32./PI/PROPZ/SQRT(SH)*PCM*FAC1 | |
819 | SIGRL=2*FACRL/32./PI/PROPZ/SQRT(SH)*PCM*FAC1 | |
820 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
821 | IF (IBREM.AND..NOT.IBEAM) THEN | |
822 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
823 | ELSE IF (IBEAM) THEN | |
824 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
825 | END IF | |
826 | IF(GOQ(80,1).AND.GOQ(81,2)) CALL SIGFIL(SIG,0,0,80,81) | |
827 | IF(GOQ(81,1).AND.GOQ(80,2)) CALL SIGFIL(SIG,0,0,81,80) | |
828 | ENDIF | |
829 | C E+ E- --> Z H_H; symmetric in cos(theta) | |
830 | IF((AMZ+AMHH).LT.RSH) THEN | |
831 | FACLR=E**2*G**2*(COS(ALFAH+BETA))**2*(AE-BE)**2/CS2THW | |
832 | FACRL=E**2*G**2*(COS(ALFAH+BETA))**2*(AE+BE)**2/CS2THW | |
833 | Z=CTH(1) | |
834 | PCM=SQRT(SSXLAM(SH,AMZ**2,AMHH**2))/4./EB | |
835 | EZ0=SQRT(PCM**2+AMZ**2) | |
836 | FAC1=AMZ**2+EZ0**2-PCM**2*Z**2 | |
837 | SIGLR=2*FACLR/32./PI/PROPZ/SQRT(SH)*PCM*FAC1 | |
838 | SIGRL=2*FACRL/32./PI/PROPZ/SQRT(SH)*PCM*FAC1 | |
839 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
840 | IF (IBREM.AND..NOT.IBEAM) THEN | |
841 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
842 | ELSE IF (IBEAM) THEN | |
843 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
844 | END IF | |
845 | IF(GOQ(80,1).AND.GOQ(82,2)) CALL SIGFIL(SIG,0,0,80,82) | |
846 | IF(GOQ(82,1).AND.GOQ(80,2)) CALL SIGFIL(SIG,0,0,82,80) | |
847 | ENDIF | |
848 | C E+ E- --> H_P H_L; symmetric in cos(theta) | |
849 | IF((AMHA+AMHL).LT.RSH) THEN | |
850 | PCM=SQRT(SSXLAM(SH,AMHA**2,AMHL**2))/4./EB | |
851 | Z=CTH(1) | |
852 | FAC1=PCM**3*(1.-Z**2) | |
853 | FACLR=E**4*(COS(ALFAH+BETA))**2*(AE-BE)**2*FAC1 | |
854 | FACRL=E**4*(COS(ALFAH+BETA))**2*(AE+BE)**2*FAC1 | |
855 | SIGLR=2*FACLR/32./PI/SQRT(SH)/SN2THW/CS2THW/PROPZ | |
856 | SIGRL=2*FACRL/32./PI/SQRT(SH)/SN2THW/CS2THW/PROPZ | |
857 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
858 | IF (IBREM.AND..NOT.IBEAM) THEN | |
859 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
860 | ELSE IF (IBEAM) THEN | |
861 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
862 | END IF | |
863 | IF(GOQ(81,1).AND.GOQ(83,2)) CALL SIGFIL(SIG,0,0,81,83) | |
864 | IF(GOQ(83,1).AND.GOQ(81,2)) CALL SIGFIL(SIG,0,0,83,81) | |
865 | ENDIF | |
866 | C E+ E- --> H_P H_H; SYMMETRIC IN COS(THETA) | |
867 | IF((AMHA+AMHH).LT.RSH) THEN | |
868 | PCM=SQRT(SSXLAM(SH,AMHA**2,AMHH**2))/4./EB | |
869 | Z=CTH(1) | |
870 | FAC1=PCM**3*(1.-Z**2) | |
871 | FACLR=E**4*(SIN(ALFAH+BETA))**2*(AE-BE)**2*FAC1 | |
872 | FACRL=E**4*(SIN(ALFAH+BETA))**2*(AE+BE)**2*FAC1 | |
873 | SIGLR=2*FACLR/32./PI/SQRT(SH)/SN2THW/CS2THW/PROPZ | |
874 | SIGRL=2*FACRL/32./PI/SQRT(SH)/SN2THW/CS2THW/PROPZ | |
875 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
876 | IF (IBREM.AND..NOT.IBEAM) THEN | |
877 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
878 | ELSE IF (IBEAM) THEN | |
879 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
880 | END IF | |
881 | IF(GOQ(82,1).AND.GOQ(83,2)) CALL SIGFIL(SIG,0,0,82,83) | |
882 | IF(GOQ(83,1).AND.GOQ(82,2)) CALL SIGFIL(SIG,0,0,83,82) | |
883 | ENDIF | |
884 | C E+ E- --> H^+ H^-; symmetric in cos(theta) | |
885 | IF((2*AMHC).LT.RSH) THEN | |
886 | PCM=SQRT(SSXLAM(SH,AMHC**2,AMHC**2))/4./EB | |
887 | Z=CTH(1) | |
888 | FAC1=PCM**3*(1.-Z**2) | |
889 | FACLR=FAC1*(1./SH**2+(2*SN2THW-1.)**2/SN2THW/CS2THW* | |
890 | $(AE-BE)**2/4./PROPZ+(2*SN2THW-1.)*(AE-BE)*(SH-AMZ**2)/SH/ | |
891 | $SQRT(SN2THW*CS2THW)/PROPZ) | |
892 | FACRL=FAC1*(1./SH**2+(2*SN2THW-1.)**2/SN2THW/CS2THW* | |
893 | $(AE+BE)**2/4./PROPZ+(2*SN2THW-1.)*(AE+BE)*(SH-AMZ**2)/SH/ | |
894 | $SQRT(SN2THW*CS2THW)/PROPZ) | |
895 | SIGLR=2*E**4*FACLR/8./PI/SQRT(SH) | |
896 | SIGRL=2*E**4*FACRL/8./PI/SQRT(SH) | |
897 | SIG=(FLEM*FREP*SIGLR+FREM*FLEP*SIGRL)*UNITS/2. | |
898 | IF (IBREM.AND..NOT.IBEAM) THEN | |
899 | SIG=SIG*ESTRUC(X1,QSQ)*ESTRUC(X2,QSQ)*JAC | |
900 | ELSE IF (IBEAM) THEN | |
901 | SIG=SIG*SSFEL(X1,0)*SSFEL(X2,0)*JAC | |
902 | END IF | |
903 | IF(GOQ(84,1).AND.GOQ(85,2)) CALL SIGFIL(SIG,0,0,84,85) | |
904 | IF(GOQ(85,1).AND.GOQ(84,2)) CALL SIGFIL(SIG,0,0,85,84) | |
905 | ENDIF | |
906 | C | |
907 | RETURN | |
908 | END |