ISAJET/code/sigsse.F

   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