]> git.uio.no Git - u/mrichter/AliRoot.git/blame - ISAJET/code/sigssl.F
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Changed the OpenESD method to take into account the path in the name of the file
[u/mrichter/AliRoot.git] / ISAJET / code / sigssl.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE SIGSSL
3C
4C Calculate d(sigma)/d(pt**2)d(y1)d(y2) for supersymmetric
5C sleptons and sneutrinos in MSSM using cross
6C sections from Baer and Tata.
7C
8C SIGMA = cross section summed over types allowed by
9C JETTYPE cards.
10C SIGS(I) = partial cross section for I1 + I2 --> I3 + I4
11C INOUT(I) = IOPAK**3*I4 + IOPAK**2*I3 + IOPAK*I2 +I1
12C JETTYP -> IDENT mapping:
13C GLSS, UPSSL, UBSSL, ..., UPSSR, UBSSR, ...,
14C W1SS+, W1SS-, WS22+, W2SS-, Z1SS, Z2SS, Z3SS, Z4SS
15C NUEL, ANUEL, EL-, ..., TAUL+
16C
17C Extra factor of 1/2 needed for nonidentical final jets.
18C Y=-log(tan(theta/2)) gives jacobean P1*P2/E1*E2
19C
20C Called from SIGSSY and so does not reinitialize /JETSIG/.
21C
22C
23#if defined(CERNLIB_IMPNONE)
24 IMPLICIT NONE
25#endif
26#include "isajet/itapes.inc"
27#include "isajet/const.inc"
28#include "isajet/jetpar.inc"
29#include "isajet/jetsig.inc"
30#include "isajet/primar.inc"
31#include "isajet/q1q2.inc"
32#include "isajet/qcdpar.inc"
33#include "isajet/sspar.inc"
34#include "isajet/sssm.inc"
35#include "isajet/sstype.inc"
36#include "isajet/wcon.inc"
37C
38 REAL X(2)
39 EQUIVALENCE (X(1),X1)
40 EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT)
41 INTEGER JS2JT(25),IW2JS(4),IW2IM(4),IZ2JS(4),IS2UD(25)
42 SAVE JS2JT,IW2JS,IW2IM,IZ2JS,IS2UD
43 INTEGER IDLSS(18)
44 SAVE IDLSS
45 INTEGER IL2JS(18),IS2LN(18),II
46 SAVE IL2JS,IS2LN
47 REAL SIG,S,T,U,FAC,AM22,AM12,TT,GP,G,
48 $E1,E2
49 INTEGER IQ,IQ1,IQ2,IH
50 REAL QFCN,STRUC,PSIFCN,AMASS
51 REAL SR2,AML,AMN,SIGW,PROPZ
52 REAL CS2THW,TNTHW,CTTHW,AL(2),BE(2),ESQ,XWI(2),YWI(2)
53 REAL ALL(2),BEL(2),EL1
54 REAL EQ1,XMGG,XMZZ,XMGZ,XM,CTH2L
55 REAL SIGUT,SIGTU,EHAT,PHAT,EBM,TPP,AMWI,AMQ,PROPW
56 REAL A,B,ASPBS,ASMBS,TM1,TM2,TM3,COTB,TANB
57 INTEGER JTYP1,JTYP2,IFLQ,IUD(13)
58 INTEGER IFLL,IL,IN,IDL,IDN,IL1,IL2,JTYPL1,JTYPL2,IDL1,IDL2
59C
60C IDENT codes from /SSTYPE/. (Fortran 77 allows - signs in
61C parameter statements but not data statements.)
62 INTEGER MSUPL,MSDNL,MSSTL,MSCHL,MSBT1,MSTP1,
63 $MSUPR,MSDNR,MSSTR,MSCHR,MSBT2,MSTP2,MSW1,MSW2,
64 $MSNEL,MSEL,MSNML,MSMUL,MSNTL,MSTAU1,MSER,MSMUR,MSTAU2
65 PARAMETER (MSUPL=-ISUPL)
66 PARAMETER (MSDNL=-ISDNL)
67 PARAMETER (MSSTL=-ISSTL)
68 PARAMETER (MSCHL=-ISCHL)
69 PARAMETER (MSBT1=-ISBT1)
70 PARAMETER (MSTP1=-ISTP1)
71 PARAMETER (MSUPR=-ISUPR)
72 PARAMETER (MSDNR=-ISDNR)
73 PARAMETER (MSSTR=-ISSTR)
74 PARAMETER (MSCHR=-ISCHR)
75 PARAMETER (MSBT2=-ISBT2)
76 PARAMETER (MSTP2=-ISTP2)
77 PARAMETER (MSW1=-ISW1)
78 PARAMETER (MSW2=-ISW2)
79 PARAMETER (MSNEL=-ISNEL)
80 PARAMETER (MSEL=-ISEL)
81 PARAMETER (MSNML=-ISNML)
82 PARAMETER (MSMUL=-ISMUL)
83 PARAMETER (MSNTL=-ISNTL)
84 PARAMETER (MSTAU1=-ISTAU1)
85 PARAMETER (MSER=-ISER)
86 PARAMETER (MSMUR=-ISMUR)
87 PARAMETER (MSTAU2=-ISTAU2)
88 DATA IDLSS/ISNEL,MSNEL,ISEL,MSEL,ISNML,MSNML,ISMUL,MSMUL,
89 $ISNTL,MSNTL,ISTAU1,MSTAU1,ISER,MSER,ISMUR,MSMUR,
90 $ISTAU2,MSTAU2/
91 DATA IUD/0,1,-1,2,-2,2,-2,1,-1,2,-2,1,-1/
92C
93C JS2JT: Susy jettype -> normal jettype
94 DATA JS2JT/1,
95 $2,3,4,5,6,7,8,9,10,11,12,13,2,3,4,5,6,7,8,9,10,11,12,13/
96C IW2JS: Wino index -> susy jettype
97 DATA IW2JS/26,27,28,29/
98C IW2IM: Wino index -> match code
99 DATA IW2IM/2,3,2,3/
100C IZ2JS: Zino index -> susy jettype
101 DATA IZ2JS/30,31,32,33/
102C IS2UD: Susy jettype -> u/d code
103 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/
104 DATA IS2LN/1,1,2,2,1,1,2,2,1,1,2,2,2,2,2,2,2,2/
105 DATA IL2JS/34,35,36,37,38,39,40,41,42,43,44,45,46,47,
106 $48,49,50,51/
107C
108C Functions
109 QFCN(IQ,IH)=STRUC(X(IH),QSQ,IQ,IDIN(IH))/X(IH)
110 PSIFCN(AM12,AM22,TT)=((S+TT-AM12)/(2*S)
111 $-AM12*(AM22-TT)/(AM12-TT)**2
112 $+(TT*(AM22-AM12)+AM22*(S-AM22+AM12))/(S*(AM12-TT)))
113C
114C Constants from Baer and Tata,
115C
116 G=SQRT(4*PI*ALFAEM/SN2THW)
117 GP=G*SQRT(SN2THW/(1.-SN2THW))
118C Quark couplings to Z
119 CS2THW=1.-SN2THW
120 TNTHW=SQRT(SN2THW/CS2THW)
121 CTTHW=1./TNTHW
122 AL(1)=(CTTHW/4.-5*TNTHW/12.)
123 AL(2)=(TNTHW/12.-CTTHW/4.)
124 BE(1)=-(CTTHW+TNTHW)/4.
125 BE(2)=-BE(1)
126 ALL(1)=(CTTHW+TNTHW)/4.
127 ALL(2)=(-CTTHW+3*TNTHW)/4.
128 BEL(1)=-(CTTHW+TNTHW)/4.
129 BEL(2)=-BEL(1)
130 ESQ=4*PI*ALFAEM
131 SR2=SQRT(2.)
132 COTB=RV2V1
133 TANB=1./COTB
134C
135C qk qb --> slss slbss
136C
137C
138C Left-leftbar slepton pair production
139C
140 DO 200 IL=1,6
141 IL1=2*IL-1
142 IL2=IL1+1
143 AML=AMASS(IDLSS(IL1))
144 JTYPL1=IL2JS(IL1)
145 JTYPL2=IL2JS(IL2)
146 IDL1=IDLSS(IL1)
147 IDL2=IDLSS(IL2)
148 IF (.NOT.(GOQ(JTYPL1,1).AND.GOQ(JTYPL2,2))) GO TO 210
149 CALL TWOKIN(0.,0.,AML,AML)
150 IF (X1.GE.1..OR.X2.GE.1.) GO TO 210
151 E1=SQRT(P(1)**2+AML**2)
152 E2=SQRT(P(2)**2+AML**2)
153 FAC=1./(16.*PI*S**2)
154 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
155 DO 220 IQ1=2,11
156 IFLQ=IS2UD(IQ1)
157 IFLL=IS2LN(IL1)
158 IF (IFLQ.EQ.1) THEN
159 EQ1=2./3.
160 ELSE
161 EQ1=-1./3.
162 END IF
163 IF (IFLL.EQ.1) THEN
164 EL1=0.
165 ELSE
166 EL1=-1.
167 END IF
168 IQ2=MATCH(IQ1,4)
169 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
170 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 220
171 XMGG=EL1**2*EQ1**2/S/S
172 CTH2L=1.
173 IF (JTYPL1.EQ.44) CTH2L=COS(2*THETAL)
174 XMZZ=(AL(IFLQ)**2+BE(IFLQ)**2)*(ALL(IFLL)-BEL(IFLL)*
175 $ CTH2L)**2/PROPZ
176 XMGZ=2*EL1*EQ1*AL(IFLQ)*(ALL(IFLL)-BEL(IFLL)*CTH2L)*
177 $ (S-AMZ**2)/S/PROPZ
178 XM=2*ESQ*ESQ*(U*T-AML**4)/3.
179 SIG=XM*(XMGG+XMZZ+XMGZ)
180 SIG=SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
181 SIG=.5*SIG
182 CALL SIGFIL(SIG,IQ1,IQ2,JTYPL1,JTYPL2)
183220 CONTINUE
184210 CONTINUE
185200 CONTINUE
186C stau_1 + stau_2 bar
187 IF (GOQ(44,1).AND.GOQ(51,2)) THEN
188 CALL TWOKIN(0.,0.,AML1SS,AML2SS)
189 IF(X1.GE.1..OR.X2.GE.1.) GO TO 231
190 E1=SQRT(P(1)**2+AML1SS**2)
191 E2=SQRT(P(2)**2+AML2SS**2)
192 FAC=1./(16.*PI*S**2)
193 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
194 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
195 DO 230 IQ1=2,11
196 IFLQ=IS2UD(IQ1)
197 IQ2=MATCH(IQ1,4)
198 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 230
199 SIG=2*ESQ**2*(AL(IFLQ)**2+BE(IFLQ)**2)*BEL(2)**2*
200 $ SIN(2*THETAL)**2*(U*T-AML1SS**2*AML2SS**2)/3./PROPZ
201 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
202 CALL SIGFIL(SIG,IQ1,IQ2,44,51)
203230 CONTINUE
204231 CONTINUE
205 END IF
206C
207C
208C Right-rightbar slepton pair production
209C
210 DO 300 IL=1,3
211 IL1=11+2*IL
212 IL2=IL1+1
213 AML=AMASS(IDLSS(IL1))
214 JTYPL1=IL2JS(IL1)
215 JTYPL2=IL2JS(IL2)
216 IDL1=IDLSS(IL1)
217 IDL2=IDLSS(IL2)
218 IF (.NOT.(GOQ(JTYPL1,1).AND.GOQ(JTYPL2,2))) GO TO 310
219 CALL TWOKIN(0.,0.,AML,AML)
220 IF (X1.GE.1..OR.X2.GE.1.) GO TO 310
221 E1=SQRT(P(1)**2+AML**2)
222 E2=SQRT(P(2)**2+AML**2)
223 FAC=1./(16.*PI*S**2)
224 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
225 DO 320 IQ1=2,11
226 IFLQ=IS2UD(IQ1)
227 IFLL=IS2LN(IL1)
228 IF (IFLQ.EQ.1) THEN
229 EQ1=2./3.
230 ELSE
231 EQ1=-1./3.
232 END IF
233 IF (IFLL.EQ.1) THEN
234 EL1=0.
235 ELSE
236 EL1=-1.
237 END IF
238 IQ2=MATCH(IQ1,4)
239 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
240 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 320
241 XMGG=EL1**2*EQ1**2/S/S
242 CTH2L=1.
243 IF (JTYPL1.EQ.50) CTH2L=COS(2*THETAL)
244 XMZZ=(AL(IFLQ)**2+BE(IFLQ)**2)*(ALL(IFLL)+BEL(IFLL)*
245 $ CTH2L)**2/PROPZ
246 XMGZ=2*EL1*EQ1*AL(IFLQ)*(ALL(IFLL)+BEL(IFLL)*CTH2L)*
247 $ (S-AMZ**2)/S/PROPZ
248 XM=2*ESQ*ESQ*(U*T-AML**4)/3.
249 SIG=XM*(XMGG+XMZZ+XMGZ)
250 SIG=SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
251 SIG=.5*SIG
252 CALL SIGFIL(SIG,IQ1,IQ2,JTYPL1,JTYPL2)
253320 CONTINUE
254310 CONTINUE
255300 CONTINUE
256C stau_2 bar + stau_1
257 IF (GOQ(51,1).AND.GOQ(44,2)) THEN
258 CALL TWOKIN(0.,0.,AML2SS,AML1SS)
259 IF(X1.GE.1..OR.X2.GE.1.) GO TO 331
260 E1=SQRT(P(1)**2+AML2SS**2)
261 E2=SQRT(P(2)**2+AML1SS**2)
262 FAC=1./(16.*PI*S**2)
263 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
264 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
265 DO 330 IQ1=2,11
266 IFLQ=IS2UD(IQ1)
267 IQ2=MATCH(IQ1,4)
268 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 330
269 SIG=2*ESQ**2*(AL(IFLQ)**2+BE(IFLQ)**2)*BEL(2)**2*
270 $ SIN(2*THETAL)**2*(U*T-AML1SS**2*AML2SS**2)/3./PROPZ
271 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
272 CALL SIGFIL(SIG,IQ1,IQ2,51,44)
273330 CONTINUE
274331 CONTINUE
275 END IF
276C
277C
278C Leftbar-left slepton pair production
279C
280 DO 400 IL=1,6
281 IL1=2*IL
282 IL2=IL1-1
283 AML=AMASS(IDLSS(IL1))
284 JTYPL1=IL2JS(IL1)
285 JTYPL2=IL2JS(IL2)
286 IDL1=IDLSS(IL1)
287 IDL2=IDLSS(IL2)
288 IF (.NOT.(GOQ(JTYPL1,1).AND.GOQ(JTYPL2,2))) GO TO 410
289 CALL TWOKIN(0.,0.,AML,AML)
290 IF (X1.GE.1..OR.X2.GE.1.) GO TO 410
291 E1=SQRT(P(1)**2+AML**2)
292 E2=SQRT(P(2)**2+AML**2)
293 FAC=1./(16.*PI*S**2)
294 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
295 DO 420 IQ1=2,11
296 IFLQ=IS2UD(IQ1)
297 IFLL=IS2LN(IL1)
298 IF (IFLQ.EQ.1) THEN
299 EQ1=2./3.
300 ELSE
301 EQ1=-1./3.
302 END IF
303 IF (IFLL.EQ.1) THEN
304 EL1=0.
305 ELSE
306 EL1=-1.
307 END IF
308 IQ2=MATCH(IQ1,4)
309 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
310 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 420
311 XMGG=EL1**2*EQ1**2/S/S
312 CTH2L=1.
313 IF (JTYPL1.EQ.45) CTH2L=COS(2*THETAL)
314 XMZZ=(AL(IFLQ)**2+BE(IFLQ)**2)*(ALL(IFLL)-BEL(IFLL)*
315 $ CTH2L)**2/PROPZ
316 XMGZ=2*EL1*EQ1*AL(IFLQ)*(ALL(IFLL)-BEL(IFLL)*CTH2L)*
317 $ (S-AMZ**2)/S/PROPZ
318 XM=2*ESQ*ESQ*(U*T-AML**4)/3.
319 SIG=XM*(XMGG+XMZZ+XMGZ)
320 SIG=SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
321 SIG=.5*SIG
322 CALL SIGFIL(SIG,IQ1,IQ2,JTYPL1,JTYPL2)
323420 CONTINUE
324410 CONTINUE
325400 CONTINUE
326C stau_1 bar + stau_2
327 IF (GOQ(45,1).AND.GOQ(50,2)) THEN
328 CALL TWOKIN(0.,0.,AML1SS,AML2SS)
329 IF(X1.GE.1..OR.X2.GE.1.) GO TO 431
330 E1=SQRT(P(1)**2+AML1SS**2)
331 E2=SQRT(P(2)**2+AML2SS**2)
332 FAC=1./(16.*PI*S**2)
333 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
334 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
335 DO 430 IQ1=2,11
336 IFLQ=IS2UD(IQ1)
337 IQ2=MATCH(IQ1,4)
338 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 430
339 SIG=2*ESQ**2*(AL(IFLQ)**2+BE(IFLQ)**2)*BEL(2)**2*
340 $ SIN(2*THETAL)**2*(U*T-AML1SS**2*AML2SS**2)/3./PROPZ
341 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
342 CALL SIGFIL(SIG,IQ1,IQ2,45,50)
343430 CONTINUE
344431 CONTINUE
345 END IF
346C
347C
348C Rightbar-right slepton pair production
349C
350 DO 500 IL=1,3
351 IL1=12+2*IL
352 IL2=IL1-1
353 AML=AMASS(IDLSS(IL1))
354 JTYPL1=IL2JS(IL1)
355 JTYPL2=IL2JS(IL2)
356 IDL1=IDLSS(IL1)
357 IDL2=IDLSS(IL2)
358 IF (.NOT.(GOQ(JTYPL1,1).AND.GOQ(JTYPL2,2))) GO TO 510
359 CALL TWOKIN(0.,0.,AML,AML)
360 IF (X1.GE.1..OR.X2.GE.1.) GO TO 510
361 E1=SQRT(P(1)**2+AML**2)
362 E2=SQRT(P(2)**2+AML**2)
363 FAC=1./(16.*PI*S**2)
364 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
365 DO 520 IQ1=2,11
366 IFLQ=IS2UD(IQ1)
367 IFLL=IS2LN(IL1)
368 IF (IFLQ.EQ.1) THEN
369 EQ1=2./3.
370 ELSE
371 EQ1=-1./3.
372 END IF
373 IF (IFLL.EQ.1) THEN
374 EL1=0.
375 ELSE
376 EL1=-1.
377 END IF
378 IQ2=MATCH(IQ1,4)
379 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
380 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 520
381 XMGG=EL1**2*EQ1**2/S/S
382 CTH2L=1.
383 IF (JTYPL1.EQ.51) CTH2L=COS(2*THETAL)
384 XMZZ=(AL(IFLQ)**2+BE(IFLQ)**2)*(ALL(IFLL)+BEL(IFLL)*
385 $ CTH2L)**2/PROPZ
386 XMGZ=2*EL1*EQ1*AL(IFLQ)*(ALL(IFLL)+BEL(IFLL)*CTH2L)*
387 $ (S-AMZ**2)/S/PROPZ
388 XM=2*ESQ*ESQ*(U*T-AML**4)/3.
389 SIG=XM*(XMGG+XMZZ+XMGZ)
390 SIG=SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
391 SIG=.5*SIG
392 CALL SIGFIL(SIG,IQ1,IQ2,JTYPL1,JTYPL2)
393520 CONTINUE
394510 CONTINUE
395500 CONTINUE
396C stau_2 + stau_1 bar
397 IF (GOQ(50,1).AND.GOQ(45,2)) THEN
398 CALL TWOKIN(0.,0.,AML2SS,AML1SS)
399 IF(X1.GE.1..OR.X2.GE.1.) GO TO 531
400 E1=SQRT(P(1)**2+AML2SS**2)
401 E2=SQRT(P(2)**2+AML1SS**2)
402 FAC=1./(16.*PI*S**2)
403 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
404 PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
405 DO 530 IQ1=2,11
406 IFLQ=IS2UD(IQ1)
407 IQ2=MATCH(IQ1,4)
408 IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 530
409 SIG=2*ESQ**2*(AL(IFLQ)**2+BE(IFLQ)**2)*BEL(2)**2*
410 $ SIN(2*THETAL)**2*(U*T-AML1SS**2*AML2SS**2)/3./PROPZ
411 SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
412 CALL SIGFIL(SIG,IQ1,IQ2,50,45)
413530 CONTINUE
414531 CONTINUE
415 END IF
416C
417C slepton+sneutrino-bar via W-*
418C
419 DO 600 II=1,3
420 IL=4*II-1
421 IN=IL-1
422 IDL=IDLSS(IL)
423 IDN=IDLSS(IN)
424 AML=AMASS(IDL)
425 AMN=AMASS(IDN)
426 JTYP1=IL2JS(IL)
427 JTYP2=IL2JS(IN)
428 IF(.NOT.(GOQ(JTYP1,1).AND.GOQ(JTYP2,2))) GO TO 610
429 CALL TWOKIN(0.,0.,AML,AMN)
430 IF(X1.GE.1..OR.X2.GE.1.) GO TO 610
431 E1=SQRT(P(1)**2+AML**2)
432 E2=SQRT(P(2)**2+AMN**2)
433 FAC=1./(16.*PI*S**2)
434 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
435 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
436 SIGW=G**4*(U*T-AML**2*AMN**2)/12./PROPW
437 IF (JTYP1.EQ.44) SIGW=SIGW*COS(THETAL)**2
438 SIG=.5*SIGW*FAC*QFCN(3,1)*QFCN(4,2)
439 CALL SIGFIL(SIG,3,4,JTYP1,JTYP2)
440 SIG=.5*SIGW*FAC*QFCN(4,1)*QFCN(3,2)
441 CALL SIGFIL(SIG,4,3,JTYP1,JTYP2)
442 SIG=.5*SIGW*FAC*QFCN(9,1)*QFCN(6,2)
443 CALL SIGFIL(SIG,9,6,JTYP1,JTYP2)
444 SIG=.5*SIGW*FAC*QFCN(6,1)*QFCN(9,2)
445 CALL SIGFIL(SIG,6,9,JTYP1,JTYP2)
446610 CONTINUE
447600 CONTINUE
448C stau_2 +nu_tau bar
449 IF (GOQ(50,1).AND.GOQ(43,2)) THEN
450 CALL TWOKIN(0.,0.,AML2SS,AMN3SS)
451 IF(X1.GE.1..OR.X2.GE.1.) GO TO 620
452 E1=SQRT(P(1)**2+AML2SS**2)
453 E2=SQRT(P(2)**2+AMN3SS**2)
454 FAC=1./(16.*PI*S**2)
455 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
456 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
457 SIGW=G**4*(U*T-AML2SS**2*AMN3SS**2)/12./PROPW
458 SIGW=SIGW*SIN(THETAL)**2
459 SIG=.5*SIGW*FAC*QFCN(3,1)*QFCN(4,2)
460 CALL SIGFIL(SIG,3,4,50,43)
461 SIG=.5*SIGW*FAC*QFCN(4,1)*QFCN(3,2)
462 CALL SIGFIL(SIG,4,3,50,43)
463 SIG=.5*SIGW*FAC*QFCN(9,1)*QFCN(6,2)
464 CALL SIGFIL(SIG,9,6,50,43)
465 SIG=.5*SIGW*FAC*QFCN(6,1)*QFCN(9,2)
466 CALL SIGFIL(SIG,6,9,50,43)
467620 CONTINUE
468 END IF
469C
470C sneutrino-bar+slepton via W-*
471C
472 DO 700 II=1,3
473 IN=4*II-2
474 IL=IN+1
475 IDL=IDLSS(IL)
476 IDN=IDLSS(IN)
477 AML=AMASS(IDL)
478 AMN=AMASS(IDN)
479 JTYP1=IL2JS(IN)
480 JTYP2=IL2JS(IL)
481 IF(.NOT.(GOQ(JTYP1,1).AND.GOQ(JTYP2,2))) GO TO 710
482 CALL TWOKIN(0.,0.,AMN,AML)
483 IF(X1.GE.1..OR.X2.GE.1.) GO TO 710
484 E1=SQRT(P(1)**2+AMN**2)
485 E2=SQRT(P(2)**2+AML**2)
486 FAC=1./(16.*PI*S**2)
487 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
488 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
489 SIGW=G**4*(U*T-AML**2*AMN**2)/12./PROPW
490 IF (JTYP2.EQ.44) SIGW=SIGW*COS(THETAL)**2
491 SIG=.5*SIGW*FAC*QFCN(3,1)*QFCN(4,2)
492 CALL SIGFIL(SIG,3,4,JTYP1,JTYP2)
493 SIG=.5*SIGW*FAC*QFCN(4,1)*QFCN(3,2)
494 CALL SIGFIL(SIG,4,3,JTYP1,JTYP2)
495 SIG=.5*SIGW*FAC*QFCN(9,1)*QFCN(6,2)
496 CALL SIGFIL(SIG,9,6,JTYP1,JTYP2)
497 SIG=.5*SIGW*FAC*QFCN(6,1)*QFCN(9,2)
498 CALL SIGFIL(SIG,6,9,JTYP1,JTYP2)
499710 CONTINUE
500700 CONTINUE
501C nu_tau bar + STAU_2
502 IF (GOQ(43,1).AND.GOQ(50,2)) THEN
503 CALL TWOKIN(0.,0.,AMN3SS,AML2SS)
504 IF(X1.GE.1..OR.X2.GE.1.) GO TO 720
505 E1=SQRT(P(1)**2+AMN3SS**2)
506 E2=SQRT(P(2)**2+AML2SS**2)
507 FAC=1./(16.*PI*S**2)
508 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
509 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
510 SIGW=G**4*(U*T-AML2SS**2*AMN3SS**2)/12./PROPW
511 SIGW=SIGW*SIN(THETAL)**2
512 SIG=.5*SIGW*FAC*QFCN(3,1)*QFCN(4,2)
513 CALL SIGFIL(SIG,3,4,43,50)
514 SIG=.5*SIGW*FAC*QFCN(4,1)*QFCN(3,2)
515 CALL SIGFIL(SIG,4,3,43,50)
516 SIG=.5*SIGW*FAC*QFCN(9,1)*QFCN(6,2)
517 CALL SIGFIL(SIG,9,6,43,50)
518 SIG=.5*SIGW*FAC*QFCN(6,1)*QFCN(9,2)
519 CALL SIGFIL(SIG,6,9,43,50)
520720 CONTINUE
521 END IF
522C
523C slepton-bar+sneutrino via W+*
524C
525 DO 800 II=1,3
526 IL=4*II
527 IN=IL-3
528 IDL=IDLSS(IL)
529 IDN=IDLSS(IN)
530 AML=AMASS(IDL)
531 AMN=AMASS(IDN)
532 JTYP1=IL2JS(IL)
533 JTYP2=IL2JS(IN)
534 IF(.NOT.(GOQ(JTYP1,1).AND.GOQ(JTYP2,2))) GO TO 810
535 CALL TWOKIN(0.,0.,AML,AMN)
536 IF(X1.GE.1..OR.X2.GE.1.) GO TO 810
537 E1=SQRT(P(1)**2+AML**2)
538 E2=SQRT(P(2)**2+AMN**2)
539 FAC=1./(16.*PI*S**2)
540 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
541 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
542 SIGW=G**4*(U*T-AML**2*AMN**2)/12./PROPW
543 IF (JTYP1.EQ.45) SIGW=SIGW*COS(THETAL)**2
544 SIG=.5*SIGW*FAC*QFCN(2,1)*QFCN(5,2)
545 CALL SIGFIL(SIG,2,5,JTYP1,JTYP2)
546 SIG=.5*SIGW*FAC*QFCN(5,1)*QFCN(2,2)
547 CALL SIGFIL(SIG,5,2,JTYP1,JTYP2)
548 SIG=.5*SIGW*FAC*QFCN(8,1)*QFCN(7,2)
549 CALL SIGFIL(SIG,8,7,JTYP1,JTYP2)
550 SIG=.5*SIGW*FAC*QFCN(7,1)*QFCN(8,2)
551 CALL SIGFIL(SIG,7,8,JTYP1,JTYP2)
552810 CONTINUE
553800 CONTINUE
554C stau_2 bar+nu_tau
555 IF (GOQ(51,1).AND.GOQ(42,2)) THEN
556 CALL TWOKIN(0.,0.,AML2SS,AMN3SS)
557 IF(X1.GE.1..OR.X2.GE.1.) GO TO 820
558 E1=SQRT(P(1)**2+AML2SS**2)
559 E2=SQRT(P(2)**2+AMN3SS**2)
560 FAC=1./(16.*PI*S**2)
561 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
562 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
563 SIGW=G**4*(U*T-AML2SS**2*AMN3SS**2)/12./PROPW
564 SIGW=SIGW*SIN(THETAL)**2
565 SIG=.5*SIGW*FAC*QFCN(2,1)*QFCN(5,2)
566 CALL SIGFIL(SIG,2,5,51,42)
567 SIG=.5*SIGW*FAC*QFCN(5,1)*QFCN(2,2)
568 CALL SIGFIL(SIG,5,2,51,42)
569 SIG=.5*SIGW*FAC*QFCN(8,1)*QFCN(7,2)
570 CALL SIGFIL(SIG,8,7,51,42)
571 SIG=.5*SIGW*FAC*QFCN(7,1)*QFCN(8,2)
572 CALL SIGFIL(SIG,7,8,51,42)
573820 CONTINUE
574 END IF
575C
576C sneutrino+slepton-bar via W+*
577C
578 DO 900 II=1,3
579 IN=4*II-3
580 IL=IN+3
581 IDL=IDLSS(IL)
582 IDN=IDLSS(IN)
583 AML=AMASS(IDL)
584 AMN=AMASS(IDN)
585 JTYP1=IL2JS(IN)
586 JTYP2=IL2JS(IL)
587 IF(.NOT.(GOQ(JTYP1,1).AND.GOQ(JTYP2,2))) GO TO 910
588 CALL TWOKIN(0.,0.,AMN,AML)
589 IF(X1.GE.1..OR.X2.GE.1.) GO TO 910
590 E1=SQRT(P(1)**2+AMN**2)
591 E2=SQRT(P(2)**2+AML**2)
592 FAC=1./(16.*PI*S**2)
593 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
594 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
595 SIGW=G**4*(U*T-AML**2*AMN**2)/12./PROPW
596 IF (JTYP2.EQ.45) SIGW=SIGW*COS(THETAL)**2
597 SIG=.5*SIGW*FAC*QFCN(2,1)*QFCN(5,2)
598 CALL SIGFIL(SIG,2,5,JTYP1,JTYP2)
599 SIG=.5*SIGW*FAC*QFCN(5,1)*QFCN(2,2)
600 CALL SIGFIL(SIG,5,2,JTYP1,JTYP2)
601 SIG=.5*SIGW*FAC*QFCN(8,1)*QFCN(7,2)
602 CALL SIGFIL(SIG,8,7,JTYP1,JTYP2)
603 SIG=.5*SIGW*FAC*QFCN(7,1)*QFCN(8,2)
604 CALL SIGFIL(SIG,7,8,JTYP1,JTYP2)
605910 CONTINUE
606900 CONTINUE
607C nu_tau + stau_2 bar
608 IF (GOQ(42,1).AND.GOQ(51,2)) THEN
609 CALL TWOKIN(0.,0.,AMN3SS,AML2SS)
610 IF(X1.GE.1..OR.X2.GE.1.) GO TO 920
611 E1=SQRT(P(1)**2+AMN3SS**2)
612 E2=SQRT(P(2)**2+AML2SS**2)
613 FAC=1./(16.*PI*S**2)
614 FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
615 PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
616 SIGW=G**4*(U*T-AML2SS**2*AMN3SS**2)/12./PROPW
617 SIGW=SIGW*SIN(THETAL)**2
618 SIG=.5*SIGW*FAC*QFCN(2,1)*QFCN(5,2)
619 CALL SIGFIL(SIG,2,5,42,51)
620 SIG=.5*SIGW*FAC*QFCN(5,1)*QFCN(2,2)
621 CALL SIGFIL(SIG,5,2,42,51)
622 SIG=.5*SIGW*FAC*QFCN(8,1)*QFCN(7,2)
623 CALL SIGFIL(SIG,8,7,42,51)
624 SIG=.5*SIGW*FAC*QFCN(7,1)*QFCN(8,2)
625 CALL SIGFIL(SIG,7,8,42,51)
626920 CONTINUE
627 END IF
628C
629 RETURN
630 END