ISAJET/code/sigssl.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE SIGSSL
   3 C
   4 C          Calculate d(sigma)/d(pt**2)d(y1)d(y2) for supersymmetric
   5 C          sleptons and sneutrinos in MSSM using cross
   6 C          sections from Baer and Tata.
   7 C
   8 C          SIGMA    = cross section summed over types allowed by
   9 C                     JETTYPE cards.
  10 C          SIGS(I)  = partial cross section for I1 + I2 --> I3 + I4
  11 C          INOUT(I) = IOPAK**3*I4 + IOPAK**2*I3 + IOPAK*I2 +I1
  12 C          JETTYP -> IDENT mapping:
  13 C          GLSS, UPSSL, UBSSL, ..., UPSSR, UBSSR, ...,
  14 C          W1SS+, W1SS-, WS22+, W2SS-, Z1SS, Z2SS, Z3SS, Z4SS
  15 C          NUEL, ANUEL, EL-, ..., TAUL+
  16 C
  17 C          Extra factor of 1/2 needed for nonidentical final jets.
  18 C          Y=-log(tan(theta/2)) gives jacobean P1*P2/E1*E2
  19 C
  20 C          Called from SIGSSY and so does not reinitialize /JETSIG/.
  21 C
  22 C
  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"
  37 C
  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
  59 C
  60 C          IDENT codes from /SSTYPE/. (Fortran 77 allows - signs in
  61 C          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/
  92 C
  93 C          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/
  96 C          IW2JS: Wino index -> susy jettype
  97       DATA IW2JS/26,27,28,29/
  98 C          IW2IM: Wino index -> match code
  99       DATA IW2IM/2,3,2,3/
 100 C          IZ2JS: Zino index -> susy jettype
 101       DATA IZ2JS/30,31,32,33/
 102 C          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/
 107 C
 108 C          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)))
 113 C
 114 C          Constants from Baer and Tata,
 115 C
 116       G=SQRT(4*PI*ALFAEM/SN2THW)
 117       GP=G*SQRT(SN2THW/(1.-SN2THW))
 118 C          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
 134 C
 135 C         qk qb --> slss slbss
 136 C
 137 C
 138 C     Left-leftbar slepton pair production
 139 C
 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)
 183 220       CONTINUE
 184 210     CONTINUE
 185 200   CONTINUE
 186 C          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)
 203 230         CONTINUE
 204 231       CONTINUE
 205       END IF
 206 C
 207 C
 208 C     Right-rightbar slepton pair production
 209 C
 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)
 253 320       CONTINUE
 254 310     CONTINUE
 255 300   CONTINUE
 256 C          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)
 273 330         CONTINUE
 274 331       CONTINUE
 275       END IF
 276 C
 277 C
 278 C     Leftbar-left slepton pair production
 279 C
 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)
 323 420       CONTINUE
 324 410     CONTINUE
 325 400   CONTINUE
 326 C          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)
 343 430         CONTINUE
 344 431       CONTINUE
 345       END IF
 346 C
 347 C
 348 C     Rightbar-right slepton pair production
 349 C
 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)
 393 520       CONTINUE
 394 510     CONTINUE
 395 500   CONTINUE
 396 C          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)
 413 530         CONTINUE
 414 531       CONTINUE
 415       END IF
 416 C
 417 C          slepton+sneutrino-bar via W-*
 418 C
 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)
 446 610     CONTINUE
 447 600   CONTINUE
 448 C          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)
 467 620       CONTINUE
 468       END IF
 469 C
 470 C          sneutrino-bar+slepton via W-*
 471 C
 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)
 499 710     CONTINUE
 500 700   CONTINUE
 501 C          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)
 520 720       CONTINUE
 521       END IF
 522 C
 523 C          slepton-bar+sneutrino via W+*
 524 C
 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)
 552 810     CONTINUE
 553 800   CONTINUE
 554 C          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)
 573 820       CONTINUE
 574       END IF
 575 C
 576 C          sneutrino+slepton-bar via W+*
 577 C
 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)
 605 910     CONTINUE
 606 900   CONTINUE
 607 C          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)
 626 920       CONTINUE
 627       END IF
 628 C
 629       RETURN
 630       END