pythiaeRHIC/pythia/pysgsu.f

0001
0002 C*********************************************************************
0003
0004 C...PYSGSU
0005 C...Subprocess cross sections for SUSY processes,
0006 C...including Higgs pair production.
0007 C...Auxiliary to PYSIGH.
0008
0009       SUBROUTINE PYSGSU(NCHN,SIGS)
0010
0011 C...Double precision and integer declarations
0012       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0013       IMPLICIT INTEGER(I-N)
0014       INTEGER PYK,PYCHGE,PYCOMP
0015 C...Parameter statement to help give large particle numbers.
0016       PARAMETER (KSUSY1=1000000,KSUSY2=2000000,KTECHN=3000000,
0017      &KEXCIT=4000000,KDIMEN=5000000)
0018 C...Commonblocks
0019       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0020       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
0021       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0022       COMMON/PYINT1/MINT(400),VINT(400)
0023       COMMON/PYINT2/ISET(500),KFPR(500,2),COEF(500,20),ICOL(40,4,2)
0024       COMMON/PYINT3/XSFX(2,-40:40),ISIG(1000,3),SIGH(1000)
0025       COMMON/PYINT4/MWID(500),WIDS(500,5)
0026       COMMON/PYMSSM/IMSS(0:99),RMSS(0:99)
0027       COMMON/PYSSMT/ZMIX(4,4),UMIX(2,2),VMIX(2,2),SMZ(4),SMW(2),
0028      &SFMIX(16,4),ZMIXI(4,4),UMIXI(2,2),VMIXI(2,2)
0029       COMMON/PYSGCM/ISUB,ISUBSV,MMIN1,MMAX1,MMIN2,MMAX2,MMINA,MMAXA,
0030      &KFAC(2,-40:40),COMFAC,FACK,FACA,SH,TH,UH,SH2,TH2,UH2,SQM3,SQM4,
0031      &SHR,SQPTH,TAUP,BE34,CTH,X(2),SQMZ,SQMW,GMMZ,GMMW,
0032      &AEM,AS,XW,XW1,XWC,XWV,POLL,POLR,POLLL,POLRR
0033       SAVE /PYDAT1/,/PYDAT2/,/PYPARS/,/PYINT1/,/PYINT2/,/PYINT3/,
0034      &/PYINT4/,/PYMSSM/,/PYSSMT/,/PYSGCM/
0035 C...Local arrays and complex variables
0036       DIMENSION WDTP(0:400),WDTE(0:400,0:5)
0037       COMPLEX*16 OLPP,ORPP,OLP,ORP,OL,OR,QLL,QLR
0038       COMPLEX*16 QRR,QRL,GLIJ,GRIJ,PROPW,PROPZ
0039       COMPLEX*16 ZMIXC(4,4),UMIXC(2,2),VMIXC(2,2)
0040
0041 CMRENNA++
0042 C...Z and W width, combinations of weak mixing angle
0043       ZWID=PMAS(23,2)
0044       WWID=PMAS(24,2)
0045       TANW=SQRT(XW/XW1)
0046       CT2W=(1D0-2D0*XW)/(2D0*XW/TANW)
0047
0048 C...Convert almost equivalent SUSY processes into each other
0049 C...Extract differences in flavours and couplings
0050
0051 C...Sleptons and sneutrinos
0052       IF(ISUB.EQ.201.OR.ISUB.EQ.204.OR.ISUB.EQ.207) THEN
0053         KFID=MOD(KFPR(ISUB,1),KSUSY1)
0054         ISUB=201
0055         ILR=0
0056       ELSEIF(ISUB.EQ.202.OR.ISUB.EQ.205.OR.ISUB.EQ.208) THEN
0057         KFID=MOD(KFPR(ISUB,1),KSUSY1)
0058         ISUB=201
0059         ILR=1
0060       ELSEIF(ISUB.EQ.203.OR.ISUB.EQ.206.OR.ISUB.EQ.209) THEN
0061         KFID=MOD(KFPR(ISUB,1),KSUSY1)
0062         ISUB=203
0063       ELSEIF(ISUB.GE.210.AND.ISUB.LE.212) THEN
0064         IF(ISUB.EQ.210) THEN
0065           RKF=2.0D0
0066         ELSEIF(ISUB.EQ.211) THEN
0067           RKF=SFMIX(15,1)**2
0068         ELSEIF(ISUB.EQ.212) THEN
0069           RKF=SFMIX(15,2)**2
0070         ENDIF
0071           ISUB=210
0072       ELSEIF(ISUB.EQ.213.OR.ISUB.EQ.214) THEN
0073         IF(ISUB.EQ.213) THEN
0074           KFID=MOD(KFPR(ISUB,1),KSUSY1)
0075           RKF=2.0D0
0076         ELSEIF(ISUB.EQ.214) THEN
0077           KFID=16
0078           RKF=1.0D0
0079         ENDIF
0080         ISUB=213
0081
0082 C...Neutralinos
0083       ELSEIF(ISUB.GE.216.AND.ISUB.LE.225) THEN
0084         IF(ISUB.EQ.216) THEN
0085           IZID1=1
0086           IZID2=1
0087         ELSEIF(ISUB.EQ.217) THEN
0088           IZID1=2
0089           IZID2=2
0090         ELSEIF(ISUB.EQ.218) THEN
0091           IZID1=3
0092           IZID2=3
0093         ELSEIF(ISUB.EQ.219) THEN
0094           IZID1=4
0095           IZID2=4
0096         ELSEIF(ISUB.EQ.220) THEN
0097           IZID1=1
0098           IZID2=2
0099         ELSEIF(ISUB.EQ.221) THEN
0100           IZID1=1
0101           IZID2=3
0102         ELSEIF(ISUB.EQ.222) THEN
0103           IZID1=1
0104           IZID2=4
0105         ELSEIF(ISUB.EQ.223) THEN
0106           IZID1=2
0107           IZID2=3
0108         ELSEIF(ISUB.EQ.224) THEN
0109           IZID1=2
0110           IZID2=4
0111         ELSEIF(ISUB.EQ.225) THEN
0112           IZID1=3
0113           IZID2=4
0114         ENDIF
0115         ISUB=216
0116
0117 C...Charginos
0118       ELSEIF(ISUB.GE.226.AND.ISUB.LE.228) THEN
0119         IF(ISUB.EQ.226) THEN
0120           IZID1=1
0121           IZID2=1
0122         ELSEIF(ISUB.EQ.227) THEN
0123           IZID1=2
0124           IZID2=2
0125         ELSEIF(ISUB.EQ.228) THEN
0126           IZID1=1
0127           IZID2=2
0128         ENDIF
0129         ISUB=226
0130
0131 C...Neutralino + chargino
0132       ELSEIF(ISUB.GE.229.AND.ISUB.LE.236) THEN
0133         IF(ISUB.EQ.229) THEN
0134           IZID1=1
0135           IZID2=1
0136         ELSEIF(ISUB.EQ.230) THEN
0137           IZID1=1
0138           IZID2=2
0139         ELSEIF(ISUB.EQ.231) THEN
0140           IZID1=1
0141           IZID2=3
0142         ELSEIF(ISUB.EQ.232) THEN
0143           IZID1=1
0144           IZID2=4
0145         ELSEIF(ISUB.EQ.233) THEN
0146           IZID1=2
0147           IZID2=1
0148         ELSEIF(ISUB.EQ.234) THEN
0149           IZID1=2
0150           IZID2=2
0151         ELSEIF(ISUB.EQ.235) THEN
0152           IZID1=2
0153           IZID2=3
0154         ELSEIF(ISUB.EQ.236) THEN
0155           IZID1=2
0156           IZID2=4
0157         ENDIF
0158         ISUB=229
0159
0160 C...Gluino + neutralino
0161       ELSEIF(ISUB.GE.237.AND.ISUB.LE.240) THEN
0162         IF(ISUB.EQ.237) THEN
0163           IZID=1
0164         ELSEIF(ISUB.EQ.238) THEN
0165           IZID=2
0166         ELSEIF(ISUB.EQ.239) THEN
0167           IZID=3
0168         ELSEIF(ISUB.EQ.240) THEN
0169           IZID=4
0170         ENDIF
0171         ISUB=237
0172
0173 C...Gluino + chargino
0174       ELSEIF(ISUB.GE.241.AND.ISUB.LE.242) THEN
0175         IF(ISUB.EQ.241) THEN
0176           IZID=1
0177         ELSEIF(ISUB.EQ.242) THEN
0178           IZID=2
0179         ENDIF
0180         ISUB=241
0181
0182 C...Squark + neutralino
0183       ELSEIF(ISUB.GE.246.AND.ISUB.LE.253) THEN
0184         ILR=0
0185         IF(MOD(ISUB,2).NE.0) ILR=1
0186         IF(ISUB.LE.247) THEN
0187           IZID=1
0188         ELSEIF(ISUB.LE.249) THEN
0189           IZID=2
0190         ELSEIF(ISUB.LE.251) THEN
0191           IZID=3
0192         ELSEIF(ISUB.LE.253) THEN
0193           IZID=4
0194         ENDIF
0195         ISUB=246
0196         RKF=5D0
0197
0198 C...Squark + chargino
0199       ELSEIF(ISUB.GE.254.AND.ISUB.LE.257) THEN
0200         IF(ISUB.LE.255) THEN
0201           IZID=1
0202         ELSEIF(ISUB.LE.257) THEN
0203           IZID=2
0204         ENDIF
0205         IF(MOD(ISUB,2).EQ.0) THEN
0206           ILR=0
0207         ELSE
0208           ILR=1
0209         ENDIF
0210         ISUB=254
0211         RKF=5D0
0212
0213 C...Squark + gluino
0214       ELSEIF(ISUB.EQ.258.OR.ISUB.EQ.259) THEN
0215         ISUB=258
0216         RKF=4D0
0217
0218 C...Stops
0219       ELSEIF(ISUB.EQ.261.OR.ISUB.EQ.262) THEN
0220         ILR=0
0221         IF(ISUB.EQ.262) ILR=1
0222         ISUB=261
0223       ELSEIF(ISUB.EQ.265) THEN
0224         ISUB=264
0225
0226 C...Squarks
0227       ELSEIF(ISUB.GE.271.AND.ISUB.LE.280) THEN
0228         ILR=0
0229         IF(ISUB.LE.273) THEN
0230           IF(ISUB.EQ.273) ILR=1
0231           ISUB=271
0232           RKF=16D0
0233         ELSEIF(ISUB.LE.276) THEN
0234           IF(ISUB.EQ.276) ILR=1
0235           ISUB=274
0236           RKF=16D0
0237         ELSEIF(ISUB.LE.278) THEN
0238           IF(ISUB.EQ.278) ILR=1
0239           ISUB=277
0240           RKF=4D0
0241         ELSE
0242           IF(ISUB.EQ.280) ILR=1
0243           ISUB=279
0244           RKF=4D0
0245         ENDIF
0246 C...Sbottoms
0247       ELSEIF(ISUB.GE.281.AND.ISUB.LE.296) THEN
0248         ILR=0
0249         IF(ISUB.LE.283) THEN
0250           IF(ISUB.EQ.283) ILR=1
0251           ISUB=271
0252           RKF=4D0
0253         ELSEIF(ISUB.LE.286) THEN
0254           IF(ISUB.EQ.286) ILR=1
0255           ISUB=274
0256           RKF=4D0
0257         ELSEIF(ISUB.LE.288) THEN
0258           IF(ISUB.EQ.288) ILR=1
0259           ISUB=277
0260           RKF=1D0
0261         ELSEIF(ISUB.LE.290) THEN
0262           IF(ISUB.EQ.290) ILR=1
0263           ISUB=279
0264           RKF=1D0
0265         ELSEIF(ISUB.LE.293) THEN
0266           IF(ISUB.EQ.293) ILR=1
0267           ISUB=271
0268           RKF=1D0
0269         ELSEIF(ISUB.EQ.296) THEN
0270           ILR=1
0271           ISUB=274
0272           RKF=1D0
0273 C...Squark + gluino
0274         ELSEIF(ISUB.EQ.294.OR.ISUB.EQ.295) THEN
0275           ISUB=258
0276           RKF=1D0
0277         ENDIF
0278 C...H+/- + H0
0279       ELSEIF(ISUB.EQ.297.OR.ISUB.EQ.298) THEN
0280         IF(ISUB.EQ.297) THEN
0281           RKF=.5D0*PARU(195)**2
0282         ELSEIF(ISUB.EQ.298) THEN
0283           RKF=.5D0*(1D0-PARU(195)**2)
0284         ENDIF
0285         ISUB=210
0286 C...A0 + H0
0287       ELSEIF(ISUB.EQ.299.OR.ISUB.EQ.300) THEN
0288         IF(ISUB.EQ.299) THEN
0289           RKF=PARU(186)**2
0290           KFID=25
0291         ELSEIF(ISUB.EQ.300) THEN
0292           RKF=PARU(187)**2
0293           KFID=35
0294         ENDIF
0295         ISUB=213
0296 C...H+ + H-
0297       ELSEIF(ISUB.EQ.301) THEN
0298         KFID=37
0299         RKF=1D0
0300         ISUB=201
0301       ENDIF
0302
0303 C...Supersymmetric processes - all of type 2 -> 2 :
0304 C...correct final-state Breit-Wigners from fixed to running width.
0305       IF(MSTP(42).GT.0) THEN
0306         DO 100 I=1,2
0307         KFLW=KFPR(ISUBSV,I)
0308         KCW=PYCOMP(KFLW)
0309         IF(PMAS(KCW,2).LT.PARP(41)) GOTO 100
0310         IF(I.EQ.1) SQMI=SQM3
0311         IF(I.EQ.2) SQMI=SQM4
0312         SQMS=PMAS(KCW,1)**2
0313         GMMS=PMAS(KCW,1)*PMAS(KCW,2)
0314         HBWS=GMMS/((SQMI-SQMS)**2+GMMS**2)
0315         CALL PYWIDT(KFLW,SQMI,WDTP,WDTE)
0316         GMMI=SQRT(SQMI)*WDTP(0)
0317         HBWI=GMMI/((SQMI-SQMS)**2+GMMI**2)
0318         COMFAC=COMFAC*(HBWI/HBWS)
0319   100   CONTINUE
0320       ENDIF
0321
0322 C...Differential cross section expressions.
0323
0324       IF(ISUB.LE.210) THEN
0325         IF(ISUB.EQ.201) THEN
0326 C...f + fbar -> e_L + e_Lbar
0327           COMFAC=COMFAC*WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
0328           DO 130 I=MMIN1,MMAX1
0329             IA=IABS(I)
0330             IF(I.EQ.0.OR.KFAC(1,I)*KFAC(2,-I).EQ.0) GOTO 130
0331             EI=KCHG(IA,1)/3D0
0332             TT3I=SIGN(1D0,EI+1D-6)/2D0
0333             EJ=-1D0
0334             TT3J=-1D0/2D0
0335             FCOL=1D0
0336 C...Color factor for e+ e-
0337             IF(IA.GE.11) FCOL=3D0
0338             IF(ISUBSV.EQ.301) THEN
0339               A1=1D0
0340               A2=0D0
0341             ELSEIF(ILR.EQ.1) THEN
0342               A1=SFMIX(KFID,3)**2
0343               A2=SFMIX(KFID,4)**2
0344             ELSEIF(ILR.EQ.0) THEN
0345               A1=SFMIX(KFID,1)**2
0346               A2=SFMIX(KFID,2)**2
0347             ENDIF
0348             XLQ=(TT3J-EJ*XW)*A1
0349             XRQ=(-EJ*XW)*A2
0350             XLF=(TT3I-EI*XW)
0351             XRF=(-EI*XW)
0352             TAA=(EI*EJ)**2*(POLL+POLR)
0353             TZZ=(XLF**2*POLL+XRF**2*POLR)*(XLQ+XRQ)**2/XW**2/XW1**2
0354             TZZ=TZZ/((1D0-SQMZ/SH)**2+SQMZ*ZWID/SH**2)
0355             TAZ=2D0*EI*EJ*(XLQ+XRQ)*(XLF*POLL+XRF*POLR)/XW/XW1
0356             TAZ=TAZ/((1D0-SQMZ/SH)**2+SQMZ*(ZWID/SH)**2)*(1D0-SQMZ/SH)
0357             TNN=0.0D0
0358             TAN=0.0D0
0359             TZN=0.0D0
0360             IF(IA.GE.11.AND.IA.LE.18.AND.KFID.EQ.IA) THEN
0361               FAC2=SQRT(2D0)
0362               TNN1=0D0
0363               TNN2=0D0
0364               TNN3=0D0
0365               DO 120 II=1,4
0366                 DK=1D0/(TH-SMZ(II)**2)
0367                 FLEK=-FAC2*(TT3I*ZMIX(II,2)-TANW*(TT3I-EI)*
0368      &          ZMIX(II,1))
0369                 FREK=FAC2*TANW*EI*ZMIX(II,1)
0370                 TNN1=TNN1+FLEK**2*DK
0371                 TNN2=TNN2+FREK**2*DK
0372                 DO 110 JJ=1,4
0373                   DL=1D0/(TH-SMZ(JJ)**2)
0374                   FLEL=-FAC2*(TT3J*ZMIX(JJ,2)-TANW*(TT3J-EJ)*
0375      &            ZMIX(JJ,1))
0376                   FREL=FAC2*TANW*EJ*ZMIX(JJ,1)
0377                   TNN3=TNN3+FLEK*FREK*FLEL*FREL*DK*DL*SMZ(II)*SMZ(JJ)
0378   110           CONTINUE
0379   120         CONTINUE
0380               TNN=(UH*TH-SQM3*SQM4)*(A1**2*TNN1**2*POLL+
0381      &        A2**2*TNN2**2*POLR)
0382               TNN=(TNN+SH*A1*A2*TNN3*((1D0-PARJ(131))*(1D0-PARJ(132))+
0383      &        (1D0+PARJ(131))*(1D0+PARJ(132))))/4D0/XW**2
0384               TZN=(UH*TH-SQM3*SQM4)*(XLQ+XRQ)*
0385      &        (TNN1*XLF*A1*POLL+TNN2*XRF*A2*POLR)
0386               TZN=TZN/((1D0-SQMZ/SH)**2+SQMZ*(ZWID/SH)**2)*
0387      &        (1D0-SQMZ/SH)/SH
0388               TZN=TZN/XW**2/XW1
0389               TAN=EI*EJ*(UH*TH-SQM3*SQM4)/SH*(A1*TNN1*POLL+
0390      &        A2*TNN2*POLR)/XW
0391             ENDIF
0392             FACQQ1=COMFAC*AEM**2*(TAA+TZZ+TAZ)*FCOL/3D0
0393             FACQQ1=FACQQ1*( UH*TH-SQM3*SQM4 )/SH**2
0394             FACQQ2=COMFAC*AEM**2*(TNN+TZN+TAN)*FCOL/3D0
0395             NCHN=NCHN+1
0396             ISIG(NCHN,1)=I
0397             ISIG(NCHN,2)=-I
0398             ISIG(NCHN,3)=1
0399             SIGH(NCHN)=FACQQ1+FACQQ2
0400   130     CONTINUE
0401
0402         ELSEIF(ISUB.EQ.203) THEN
0403 C...f + fbar -> e_L + e_Rbar
0404           DO 160 I=MMIN1,MMAX1
0405             IA=IABS(I)
0406             IF(I.EQ.0.OR.KFAC(1,I)*KFAC(2,-I).EQ.0) GOTO 160
0407             EI=KCHG(IABS(I),1)/3D0
0408             TT3I=SIGN(1D0,EI)/2D0
0409             EJ=-1
0410             TT3J=-1D0/2D0
0411             FCOL=1D0
0412 C...Color factor for e+ e-
0413             IF(IA.GE.11) FCOL=3D0
0414             A1=SFMIX(KFID,1)**2
0415             A2=SFMIX(KFID,2)**2
0416             XLQ=(TT3J-EJ*XW)
0417             XRQ=(-EJ*XW)
0418             XLF=(TT3I-EI*XW)
0419             XRF=(-EI*XW)
0420             TZZ=(XLF**2*POLL+XRF**2*POLR)*(XLQ-XRQ)**2
0421      &      /XW**2/XW1**2*A1*A2
0422             TZZ=TZZ/((1D0-SQMZ/SH)**2+SQMZ*(ZWID/SH)**2)
0423             TNN=0.0D0
0424             TZN=0.0D0
0425             TNNA=0D0
0426             TNNB=0D0
0427             IF(IA.GE.11.AND.IA.LE.18.AND.KFID.EQ.IA) THEN
0428               FAC2=SQRT(2D0)
0429               TNN1=0D0
0430               TNN2=0D0
0431               TNN3=0D0
0432               DO 150 II=1,4
0433                 DK=1D0/(TH-SMZ(II)**2)
0434                 FLEK=-FAC2*(TT3I*ZMIX(II,2)-TANW*(TT3I-EI)*
0435      &          ZMIX(II,1))
0436                 FREK=FAC2*TANW*EI*ZMIX(II,1)
0437                 TNN1=TNN1+FLEK**2*DK
0438                 TNN2=TNN2+FREK**2*DK
0439                 DO 140 JJ=1,4
0440                   DL=1D0/(TH-SMZ(JJ)**2)
0441                   FLEL=-FAC2*(TT3J*ZMIX(JJ,2)-TANW*(TT3J-EJ)*
0442      &            ZMIX(JJ,1))
0443                   FREL=FAC2*TANW*EJ*ZMIX(JJ,1)
0444                   TNN3=TNN3+FLEK*FREK*FLEL*FREL*DK*DL*SMZ(II)*SMZ(JJ)
0445   140           CONTINUE
0446   150         CONTINUE
0447               TNN=(UH*TH-SQM3*SQM4)*A1*A2*(TNN2**2*POLR+TNN1**2*POLL)
0448               TNNA=(TNN+SH*(A1**2*POLLL+A2**2*POLRR)*TNN3)/4D0
0449               TNNB=(TNN+SH*(A1**2*POLRR+A2**2*POLLL)*TNN3)/4D0
0450               TZN=(UH*TH-SQM3*SQM4)*A1*A2
0451               TZN=TZN*(XLQ-XRQ)*(XLF*TNN1*POLL-XRF*TNN2*POLR)/XW1
0452               TZN=TZN/((1D0-SQMZ/SH)**2+SQMZ*(ZWID/SH)**2)*
0453      &        (1D0-SQMZ/SH)/SH
0454             ENDIF
0455             FACQQ0=COMFAC*AEM**2*TZZ*FCOL/3D0*(UH*TH-SQM3*SQM4)/SH2
0456             FACQQ2=COMFAC*AEM**2/XW**2*(TNNA+TZN)*FCOL/3D0
0457             FACQQ1=COMFAC*AEM**2/XW**2*(TNNB+TZN)*FCOL/3D0
0458 C%%%%%%%%%%%
0459             NCHN=NCHN+1
0460             ISIG(NCHN,1)=I
0461             ISIG(NCHN,2)=-I
0462             ISIG(NCHN,3)=1
0463             SIGH(NCHN)=(FACQQ0+FACQQ1)*WIDS(PYCOMP(KFPR(ISUBSV,1)),2)*
0464      &      WIDS(PYCOMP(KFPR(ISUBSV,2)),3)
0465             NCHN=NCHN+1
0466             ISIG(NCHN,1)=I
0467             ISIG(NCHN,2)=-I
0468             ISIG(NCHN,3)=2
0469             SIGH(NCHN)=(FACQQ0+FACQQ2)*WIDS(PYCOMP(KFPR(ISUBSV,1)),3)*
0470      &      WIDS(PYCOMP(KFPR(ISUBSV,2)),2)
0471   160     CONTINUE
0472
0473         ELSEIF(ISUB.EQ.210) THEN
0474 C...q + qbar' -> W*- > ~l_L + ~nu_L
0475           FAC0=RKF*COMFAC*AEM**2/XW**2/12D0
0476           FAC1=(TH*UH-SQM3*SQM4)/((SH-SQMW)**2+WWID**2*SQMW)
0477           DO 180 I=MMIN1,MMAX1
0478             IA=IABS(I)
0479             IF(I.EQ.0.OR.IA.GT.10.OR.KFAC(1,I).EQ.0) GOTO 180
0480             DO 170 J=MMIN2,MMAX2
0481               JA=IABS(J)
0482               IF(J.EQ.0.OR.JA.GT.10.OR.KFAC(2,J).EQ.0) GOTO 170
0483               IF(I*J.GT.0.OR.MOD(IA+JA,2).EQ.0) GOTO 170
0484               FCKM=3D0
0485               IF(IA.LE.10) FCKM=VCKM((IA+1)/2,(JA+1)/2)
0486               KCHSUM=KCHG(IA,1)*ISIGN(1,I)+KCHG(JA,1)*ISIGN(1,J)
0487               KCHW=2
0488               IF(KCHSUM.LT.0) KCHW=3
0489               NCHN=NCHN+1
0490               ISIG(NCHN,1)=I
0491               ISIG(NCHN,2)=J
0492               ISIG(NCHN,3)=1
0493               IF(ISUBSV.EQ.297.OR.ISUBSV.EQ.298) THEN
0494                 FACR=WIDS(PYCOMP(KFPR(ISUBSV,1)),5-KCHW)*
0495      &          WIDS(PYCOMP(KFPR(ISUBSV,2)),2)
0496               ELSE
0497                 FACR=WIDS(PYCOMP(KFPR(ISUBSV,1)),5-KCHW)*
0498      &          WIDS(PYCOMP(KFPR(ISUBSV,2)),KCHW)
0499               ENDIF
0500               SIGH(NCHN)=FAC0*FAC1*FCKM*FACR
0501   170       CONTINUE
0502   180     CONTINUE
0503         ENDIF
0504
0505       ELSEIF(ISUB.LE.220) THEN
0506         IF(ISUB.EQ.213) THEN
0507 C...f + fbar -> ~nu_L + ~nu_Lbar
0508           IF(ISUBSV.EQ.299.OR.ISUBSV.EQ.300) THEN
0509             FACR=WIDS(PYCOMP(KFPR(ISUBSV,1)),2)*
0510      &      WIDS(PYCOMP(KFPR(ISUBSV,2)),2)
0511           ELSE
0512             FACR=WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
0513           ENDIF
0514           COMFAC=COMFAC*FACR
0515           PROPZ2=(SH-SQMZ)**2+ZWID**2*SQMZ
0516           XLL=0.5D0
0517           XLR=0.0D0
0518           DO 190 I=MMIN1,MMAX1
0519             IA=IABS(I)
0520             IF(I.EQ.0.OR.KFAC(1,I)*KFAC(2,-I).EQ.0) GOTO 190
0521             EI=KCHG(IA,1)/3D0
0522             FCOL=1D0
0523 C...Color factor for e+ e-
0524             IF(IA.GE.11) FCOL=3D0
0525             XLQ=(SIGN(1D0,EI)-2D0*EI*XW)/2D0
0526             XRQ=-EI*XW
0527             TZC=0.0D0
0528             TCC=0.0D0
0529             IF(IA.GE.11.AND.KFID.EQ.IA+1) THEN
0530               TZC=VMIX(1,1)**2/(TH-SMW(1)**2)+VMIX(2,1)**2/
0531      &        (TH-SMW(2)**2)
0532               TCC=TZC**2
0533               TZC=TZC/XW1*(SH-SQMZ)/PROPZ2*XLQ*XLL
0534             ENDIF
0535             FACQQ1=(XLQ**2+XRQ**2)*(XLL+XLR)**2/XW1**2/PROPZ2
0536             FACQQ2=TZC+TCC/4D0
0537             NCHN=NCHN+1
0538             ISIG(NCHN,1)=I
0539             ISIG(NCHN,2)=-I
0540             ISIG(NCHN,3)=1
0541             SIGH(NCHN)=(FACQQ1+FACQQ2)*RKF*(UH*TH-SQM3*SQM4)*COMFAC
0542      &      *AEM**2*FCOL/3D0/XW**2
0543   190     CONTINUE
0544
0545         ELSEIF(ISUB.EQ.216) THEN
0546 C...q + qbar -> ~chi0_1 + ~chi0_1
0547           IF(IZID1.EQ.IZID2) THEN
0548             COMFAC=COMFAC*WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
0549           ELSE
0550             COMFAC=COMFAC*WIDS(PYCOMP(KFPR(ISUBSV,1)),2)*
0551      &      WIDS(PYCOMP(KFPR(ISUBSV,2)),2)
0552           ENDIF
0553           FACXX=COMFAC*AEM**2/3D0/XW**2
0554           IF(IZID1.EQ.IZID2) FACXX=FACXX/2D0
0555           ZM12=SQM3
0556           ZM22=SQM4
0557           WU2 = (UH-ZM12)*(UH-ZM22)
0558           WT2 = (TH-ZM12)*(TH-ZM22)
0559           WS2 = SMZ(IZID1)*SMZ(IZID2)*SH
0560           PROPZ2 = (SH-SQMZ)**2 + SQMZ*ZWID**2
0561           PROPZ=DCMPLX(SH-SQMZ,-ZWID*PMAS(23,1))/DCMPLX(PROPZ2)
0562           DO 200 I=1,4
0563             ZMIXC(IZID1,I)=DCMPLX(ZMIX(IZID1,I),ZMIXI(IZID1,I))
0564             IF(IZID2.NE.IZID1) THEN
0565               ZMIXC(IZID2,I)=DCMPLX(ZMIX(IZID2,I),ZMIXI(IZID2,I))
0566             ENDIF
0567   200     CONTINUE
0568           OLPP=(ZMIXC(IZID1,3)*DCONJG(ZMIXC(IZID2,3))-
0569      &    ZMIXC(IZID1,4)*DCONJG(ZMIXC(IZID2,4)))/2D0
0570           ORPP=DCONJG(OLPP)
0571           DO 210 I=MMINA,MMAXA
0572             IF(I.EQ.0.OR.KFAC(1,I)*KFAC(2,-I).EQ.0) GOTO 210
0573             EI=KCHG(IABS(I),1)/3D0
0574             T3I=SIGN(1D0,EI+1D-6)/2D0
0575             XML2=PMAS(PYCOMP(KSUSY1+IABS(I)),1)**2
0576             XMR2=PMAS(PYCOMP(KSUSY2+IABS(I)),1)**2
0577             GLIJ=(T3I*ZMIXC(IZID1,2)-TANW*(T3I-EI)*ZMIXC(IZID1,1))*
0578      &      DCONJG(T3I*ZMIXC(IZID2,2)-TANW*(T3I-EI)*ZMIXC(IZID2,1))
0579             GRIJ=ZMIXC(IZID1,1)*DCONJG(ZMIXC(IZID2,1))*(EI*TANW)**2
0580             QLL=DCMPLX((T3I-EI*XW)/XW1)*OLPP*PROPZ-GLIJ/DCMPLX(UH-XML2)
0581             QLR=-DCMPLX((T3I-EI*XW)/XW1)*ORPP*PROPZ+DCONJG(GLIJ)
0582      &      /DCMPLX(TH-XML2)
0583             QRL=-DCMPLX((EI*XW)/XW1)*OLPP*PROPZ+GRIJ/DCMPLX(TH-XMR2)
0584             QRR=DCMPLX((EI*XW)/XW1)*ORPP*PROPZ
0585      &      -DCONJG(GRIJ)/DCMPLX(UH-XMR2)
0586             FCOL=1D0
0587             IF(IABS(I).GE.11) FCOL=3D0
0588             FACGG1=(ABS(QLL)**2*POLL+ABS(QRR)**2*POLR)*WU2+
0589      &      (ABS(QRL)**2*POLR+ABS(QLR)**2*POLL)*WT2+
0590      &      2D0*DBLE(QLR*DCONJG(QLL)*POLL+
0591      &      QRL*DCONJG(QRR)*POLR)*WS2
0592             NCHN=NCHN+1
0593             ISIG(NCHN,1)=I
0594             ISIG(NCHN,2)=-I
0595             ISIG(NCHN,3)=1
0596             SIGH(NCHN)=FACXX*FACGG1*FCOL
0597   210     CONTINUE
0598         ENDIF
0599
0600       ELSEIF(ISUB.LE.230) THEN
0601         IF(ISUB.EQ.226) THEN
0602 C...f + fbar -> ~chi+_1 + ~chi-_1
0603           FACXX=COMFAC*AEM**2/3D0
0604           ZM12=SQM3
0605           ZM22=SQM4
0606           WU2 = (UH-ZM12)*(UH-ZM22)
0607           WT2 = (TH-ZM12)*(TH-ZM22)
0608           WS2 = SMW(IZID1)*SMW(IZID2)*SH
0609           PROPZ2 = (SH-SQMZ)**2 + SQMZ*ZWID**2
0610           PROPZ=DCMPLX(SH-SQMZ,-ZWID*PMAS(23,1))/DCMPLX(PROPZ2)
0611           DIFF=0D0
0612           IF(IZID1.EQ.IZID2) DIFF=1D0
0613           DO 220 I=1,2
0614             VMIXC(IZID1,I)=DCMPLX(VMIX(IZID1,I),VMIXI(IZID1,I))
0615             UMIXC(IZID1,I)=DCMPLX(UMIX(IZID1,I),UMIXI(IZID1,I))
0616             IF(IZID2.NE.IZID1) THEN
0617               VMIXC(IZID2,I)=DCMPLX(VMIX(IZID2,I),VMIXI(IZID2,I))
0618               UMIXC(IZID2,I)=DCMPLX(UMIX(IZID2,I),UMIXI(IZID2,I))
0619             ENDIF
0620   220     CONTINUE
0621           OLP=-VMIXC(IZID2,1)*DCONJG(VMIXC(IZID1,1))-
0622      &    VMIXC(IZID2,2)*DCONJG(VMIXC(IZID1,2))/2D0+DCMPLX(XW*DIFF)
0623           ORP=-UMIXC(IZID1,1)*DCONJG(UMIXC(IZID2,1))-
0624      &    UMIXC(IZID1,2)*DCONJG(UMIXC(IZID2,2))/2D0+DCMPLX(XW*DIFF)
0625           DO 230 I=MMINA,MMAXA
0626             IF(I.EQ.0.OR.KFAC(1,I)*KFAC(2,-I).EQ.0) GOTO 230
0627             EI=KCHG(IABS(I),1)/3D0
0628             T3I=SIGN(1D0,EI+1D-6)/2D0
0629             QRL=DCMPLX(-EI/SH*DIFF)-DCMPLX(EI/XW1)*PROPZ*ORP
0630             QLL=DCMPLX(-EI/SH*DIFF)+DCMPLX((T3I-XW*EI)/XW/XW1)*PROPZ*ORP
0631             QRR=DCMPLX(-EI/SH*DIFF)-DCMPLX(EI/XW1)*PROPZ*OLP
0632             IF(MOD(I,2).EQ.0) THEN
0633               XML2=PMAS(PYCOMP(KSUSY1+IABS(I)-1),1)**2
0634               QLR=DCMPLX(-EI/SH*DIFF)+DCMPLX((T3I-XW*EI)/XW/XW1)*
0635      &        PROPZ*OLP-UMIXC(IZID2,1)*DCONJG(UMIXC(IZID1,1))*
0636      &        DCMPLX(T3I/XW/(TH-XML2))
0637             ELSE
0638               XML2=PMAS(PYCOMP(KSUSY1+IABS(I)+1),1)**2
0639               QLR=DCMPLX(-EI/SH*DIFF)+DCMPLX((T3I-XW*EI)/XW/XW1)*
0640      &        PROPZ*OLP-VMIXC(IZID2,1)*DCONJG(VMIXC(IZID1,1))*
0641      &        DCMPLX(T3I/XW/(TH-XML2))
0642             ENDIF
0643             FCOL=1D0
0644             IF(IABS(I).GE.11) FCOL=3D0
0645             FACSUM=((ABS(QLL)**2*POLL+ABS(QRR)**2*POLR)*WU2+
0646      &      (ABS(QRL)**2*POLR+ABS(QLR)**2*POLL)*WT2+
0647      &      2D0*DBLE(QLR*DCONJG(QLL)*POLL+
0648      &      QRL*DCONJG(QRR)*POLR)*WS2)*FACXX*FCOL
0649             NCHN=NCHN+1
0650             ISIG(NCHN,1)=I
0651             ISIG(NCHN,2)=-I
0652             ISIG(NCHN,3)=1
0653             IF(IZID1.EQ.IZID2) THEN
0654               SIGH(NCHN)=FACSUM*WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
0655             ELSE
0656               SIGH(NCHN)=FACSUM*WIDS(PYCOMP(KFPR(ISUBSV,1)),3)*
0657      &        WIDS(PYCOMP(KFPR(ISUBSV,2)),2)
0658               NCHN=NCHN+1
0659               ISIG(NCHN,1)=I
0660               ISIG(NCHN,2)=-I
0661               ISIG(NCHN,3)=2
0662               SIGH(NCHN)=FACSUM*WIDS(PYCOMP(KFPR(ISUBSV,1)),2)*
0663      &        WIDS(PYCOMP(KFPR(ISUBSV,2)),3)
0664             ENDIF
0665   230     CONTINUE
0666
0667         ELSEIF(ISUB.EQ.229) THEN
0668 C...q + qbar' -> ~chi0_1 + ~chi+-_1
0669           FACXX=COMFAC*AEM**2/6D0/XW**2
0670           ZM12=SQM3
0671           ZM22=SQM4
0672           WU2 = (UH-ZM12)*(UH-ZM22)
0673           WT2 = (TH-ZM12)*(TH-ZM22)
0674           WS2 = SMW(IZID1)*SMZ(IZID2)*SH
0675           RT2I = 1D0/SQRT(2D0)
0676           PROPW = DCMPLX(SH-SQMW,-WWID*PMAS(24,1))/
0677      &    DCMPLX((SH-SQMW)**2+WWID**2*SQMW,0D0)
0678           DO 240 I=1,2
0679             VMIXC(IZID1,I)=DCMPLX(VMIX(IZID1,I),VMIXI(IZID1,I))
0680             UMIXC(IZID1,I)=DCMPLX(UMIX(IZID1,I),UMIXI(IZID1,I))
0681   240     CONTINUE
0682           DO 250 I=1,4
0683             ZMIXC(IZID2,I)=DCMPLX(ZMIX(IZID2,I),ZMIXI(IZID2,I))
0684   250     CONTINUE
0685           OL=(DCONJG(ZMIXC(IZID2,2))*VMIXC(IZID1,1)-
0686      &    DCONJG(ZMIXC(IZID2,4))*VMIXC(IZID1,2)*RT2I)*PROPW
0687           OR=(ZMIXC(IZID2,2)*DCONJG(UMIXC(IZID1,1))+
0688      &    ZMIXC(IZID2,3)*DCONJG(UMIXC(IZID1,2))*RT2I)*PROPW
0689
0690           DO 270 I=MMIN1,MMAX1
0691             IA=IABS(I)
0692             IF(I.EQ.0.OR.IA.GT.20.OR.KFAC(1,I).EQ.0) GOTO 270
0693             EI=KCHG(IA,1)/3D0
0694             T3I=SIGN(1D0,EI+1D-6)/2D0
0695             DO 260 J=MMIN2,MMAX2
0696               JA=IABS(J)
0697               IF(J.EQ.0.OR.JA.GT.20.OR.KFAC(2,J).EQ.0) GOTO 260
0698               IF(I*J.GT.0.OR.MOD(IA+JA,2).EQ.0) GOTO 260
0699               EJ=KCHG(JA,1)/3D0
0700               T3J=SIGN(1D0,EJ+1D-6)/2D0
0701               FCKM=3D0
0702               IF(IA.LE.10) FCKM=VCKM((IA+1)/2,(JA+1)/2)
0703               KCHSUM=KCHG(IA,1)*ISIGN(1,I)+KCHG(JA,1)*ISIGN(1,J)
0704               KCHW=2
0705               IF(KCHSUM.LT.0) KCHW=3
0706               IF(MOD(IA,2).EQ.0) THEN
0707                 ZMI2  = PMAS(PYCOMP(KSUSY1+IA),1)**2
0708                 ZMJ2  = PMAS(PYCOMP(KSUSY1+JA),1)**2
0709                 QLL=OL+VMIXC(IZID1,1)*DCONJG(ZMIXC(IZID2,1)*(EI-T3I)*
0710      &          TANW+ZMIXC(IZID2,2)*T3I)/DCMPLX(UH-ZMI2)
0711                 QLR=OR-DCONJG(UMIXC(IZID1,1))*(
0712      &          ZMIXC(IZID2,1)*(EJ-T3J)*TANW+ZMIXC(IZID2,2)*T3J)
0713      &          /DCMPLX(TH-ZMJ2)
0714               ELSE
0715                 ZMI2  = PMAS(PYCOMP(KSUSY1+JA),1)**2
0716                 ZMJ2  = PMAS(PYCOMP(KSUSY1+IA),1)**2
0717                 QLL=OL+VMIXC(IZID1,1)*DCONJG(ZMIXC(IZID2,1)*(EJ-T3J)*
0718      &          TANW+ZMIXC(IZID2,2)*T3J)/DCMPLX(UH-ZMJ2)
0719                 QLR=OR-DCONJG(UMIXC(IZID1,1))*(
0720      &          ZMIXC(IZID2,1)*(EI-T3I)*TANW+ZMIXC(IZID2,2)*T3I)
0721      &          /DCMPLX(TH-ZMI2)
0722               ENDIF
0723               ZINTR=DBLE(QLR*DCONJG(QLL))
0724               FACGG1=FACXX*(ABS(QLL)**2*WU2+ABS(QLR)**2*WT2+
0725      &        2D0*ZINTR*WS2)
0726               NCHN=NCHN+1
0727               ISIG(NCHN,1)=I
0728               ISIG(NCHN,2)=J
0729               ISIG(NCHN,3)=1
0730               SIGH(NCHN)=FACGG1*FCKM*WIDS(PYCOMP(KFPR(ISUBSV,1)),2)*
0731      &        WIDS(PYCOMP(KFPR(ISUBSV,2)),KCHW)
0732   260       CONTINUE
0733   270     CONTINUE
0734         ENDIF
0735
0736       ELSEIF(ISUB.LE.240) THEN
0737         IF(ISUB.EQ.237) THEN
0738 C...q + qbar -> gluino + ~chi0_1
0739           COMFAC=COMFAC*WIDS(PYCOMP(KFPR(ISUBSV,1)),2)*
0740      &    WIDS(PYCOMP(KFPR(ISUBSV,2)),2)
0741           ASYUK=RMSS(42)*AS
0742           FAC0=COMFAC*ASYUK*AEM*4D0/9D0/XW
0743           GM2=SQM3
0744           ZM2=SQM4
0745           DO 280 I=MMINA,MMAXA
0746             IF(I.EQ.0.OR.IABS(I).GT.MSTP(58)) GOTO 280
0747             EI=KCHG(IABS(I),1)/3D0
0748             IA=IABS(I)
0749             XLQC = -TANW*EI*ZMIX(IZID,1)
0750             XRQC =(SIGN(1D0,EI)*ZMIX(IZID,2)-TANW*
0751      &      (SIGN(1D0,EI)-2D0*EI)*ZMIX(IZID,1))/2D0
0752             XLQ2=XLQC**2
0753             XRQ2=XRQC**2
0754             XML2=PMAS(PYCOMP(KSUSY1+IA),1)**2
0755             XMR2=PMAS(PYCOMP(KSUSY2+IA),1)**2
0756             ATKIN=(TH-GM2)*(TH-ZM2)/(TH-XML2)**2
0757             AUKIN=(UH-GM2)*(UH-ZM2)/(UH-XML2)**2
0758             ATUKIN=SMZ(IZID)*SQRT(GM2)*SH/(TH-XML2)/(UH-XML2)
0759             SGCHIL=XLQ2*(ATKIN+AUKIN-2D0*ATUKIN)
0760             ATKIN=(TH-GM2)*(TH-ZM2)/(TH-XMR2)**2
0761             AUKIN=(UH-GM2)*(UH-ZM2)/(UH-XMR2)**2
0762             ATUKIN=SMZ(IZID)*SQRT(GM2)*SH/(TH-XMR2)/(UH-XMR2)
0763             SGCHIR=XRQ2*(ATKIN+AUKIN-2D0*ATUKIN)
0764             NCHN=NCHN+1
0765             ISIG(NCHN,1)=I
0766             ISIG(NCHN,2)=-I
0767             ISIG(NCHN,3)=1
0768             SIGH(NCHN)=FAC0*(SGCHIL+SGCHIR)
0769   280     CONTINUE
0770         ENDIF
0771
0772       ELSEIF(ISUB.LE.250) THEN
0773         IF(ISUB.EQ.241) THEN
0774 C...q + qbar' -> ~chi+-_1 + gluino
0775           FACWG=COMFAC*AS*AEM/XW*2D0/9D0
0776           GM2=SQM3
0777           ZM2=SQM4
0778           FAC01=2D0*UMIX(IZID,1)*VMIX(IZID,1)
0779           FAC0=UMIX(IZID,1)**2
0780           FAC1=VMIX(IZID,1)**2
0781           DO 300 I=MMIN1,MMAX1
0782             IA=IABS(I)
0783             IF(I.EQ.0.OR.IA.GT.10.OR.KFAC(1,I).EQ.0) GOTO 300
0784             DO 290 J=MMIN2,MMAX2
0785               JA=IABS(J)
0786               IF(J.EQ.0.OR.JA.GT.10.OR.KFAC(2,J).EQ.0) GOTO 290
0787               IF(I*J.GT.0.OR.MOD(IA+JA,2).EQ.0) GOTO 290
0788               FCKM=1D0
0789               IF(IA.LE.10) FCKM=VCKM((IA+1)/2,(JA+1)/2)
0790               KCHSUM=KCHG(IA,1)*ISIGN(1,I)+KCHG(JA,1)*ISIGN(1,J)
0791               KCHW=2
0792               IF(KCHSUM.LT.0) KCHW=3
0793               XMU2=PMAS(PYCOMP(KSUSY1+2),1)**2
0794               XMD2=PMAS(PYCOMP(KSUSY1+1),1)**2
0795               ATKIN=(TH-GM2)*(TH-ZM2)/(TH-XMU2)**2
0796               AUKIN=(UH-GM2)*(UH-ZM2)/(UH-XMD2)**2
0797               ATUKIN=SMW(IZID)*SQRT(GM2)*SH/(TH-XMU2)/(UH-XMD2)
0798               XMU2=PMAS(PYCOMP(KSUSY2+2),1)**2
0799               XMD2=PMAS(PYCOMP(KSUSY2+1),1)**2
0800               ATKIN=(ATKIN+(TH-GM2)*(TH-ZM2)/(TH-XMU2)**2)/2D0
0801               AUKIN=(AUKIN+(UH-GM2)*(UH-ZM2)/(UH-XMD2)**2)/2D0
0802               ATUKIN=(ATUKIN+SMW(IZID)*SQRT(GM2)*
0803      &        SH/(TH-XMU2)/(UH-XMD2))/2D0
0804               NCHN=NCHN+1
0805               ISIG(NCHN,1)=I
0806               ISIG(NCHN,2)=J
0807               ISIG(NCHN,3)=1
0808               SIGH(NCHN)=FACWG*FCKM*(FAC0*ATKIN+FAC1*AUKIN-
0809      &        FAC01*ATUKIN)*WIDS(PYCOMP(KFPR(ISUBSV,1)),2)*
0810      &        WIDS(PYCOMP(KFPR(ISUBSV,2)),KCHW)
0811   290       CONTINUE
0812   300     CONTINUE
0813
0814         ELSEIF(ISUB.EQ.243) THEN
0815 C...q + qbar -> gluino + gluino
0816           COMFAC=COMFAC*WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
0817           XMT=SQM3-TH
0818           XMU=SQM3-UH
0819           DO 310 I=MMINA,MMAXA
0820             IF(I.EQ.0.OR.IABS(I).GT.MSTP(58).OR.
0821      &      KFAC(1,I)*KFAC(2,-I).EQ.0) GOTO 310
0822             NCHN=NCHN+1
0823             XSU=PMAS(PYCOMP(KSUSY1+IABS(I)),1)**2-UH
0824             XST=PMAS(PYCOMP(KSUSY1+IABS(I)),1)**2-TH
0825             FACGG1=COMFAC*AS**2*8D0/3D0*( (XMT**2+XMU**2+
0826      &      2D0*SQM3*SH)/SH2 + RMSS(42)**2*(4D0/9D0*(XMT**2/XST**2+
0827      &      XMU**2/XSU**2) + SQM3*SH/XST/XSU/9D0) - RMSS(42)*(
0828      &      (XMT**2+SH*SQM3)/SH/XST + (XMU**2+SH*SQM3)/SH/XSU ))
0829             XSU=PMAS(PYCOMP(KSUSY2+IABS(I)),1)**2-UH
0830             XST=PMAS(PYCOMP(KSUSY2+IABS(I)),1)**2-TH
0831             FACGG2=COMFAC*AS**2*8D0/3D0*( (XMT**2+XMU**2+
0832      &      2D0*SQM3*SH)/SH2 + RMSS(42)**2*(4D0/9D0*(XMT**2/XST**2+
0833      &      XMU**2/XSU**2) + SQM3*SH/XST/XSU/9D0) - RMSS(42)*(
0834      &      (XMT**2+SH*SQM3)/SH/XST + (XMU**2+SH*SQM3)/SH/XSU ))
0835             ISIG(NCHN,1)=I
0836             ISIG(NCHN,2)=-I
0837             ISIG(NCHN,3)=1
0838 C...1/2 for identical particles
0839             SIGH(NCHN)=0.25D0*(FACGG1+FACGG2)
0840   310     CONTINUE
0841
0842         ELSEIF(ISUB.EQ.244) THEN
0843 C...g + g -> gluino + gluino
0844           COMFAC=COMFAC*WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
0845           XMT=SQM3-TH
0846           XMU=SQM3-UH
0847           FACQQ1=COMFAC*AS**2*9D0/4D0*(
0848      &    (XMT*XMU-2D0*SQM3*(TH+SQM3))/XMT**2 -
0849      &    (XMT*XMU+SQM3*(UH-TH))/SH/XMT )
0850           FACQQ2=COMFAC*AS**2*9D0/4D0*(
0851      &    (XMU*XMT-2D0*SQM3*(UH+SQM3))/XMU**2 -
0852      &    (XMU*XMT+SQM3*(TH-UH))/SH/XMU )
0853           FACQQ3=COMFAC*AS**2*9D0/4D0*(2D0*XMT*XMU/SH2 +
0854      &    SQM3*(SH-4D0*SQM3)/XMT/XMU)
0855           IF(KFAC(1,21)*KFAC(2,21).EQ.0) GOTO 320
0856           NCHN=NCHN+1
0857           ISIG(NCHN,1)=21
0858           ISIG(NCHN,2)=21
0859           ISIG(NCHN,3)=1
0860           SIGH(NCHN)=FACQQ1/2D0
0861           NCHN=NCHN+1
0862           ISIG(NCHN,1)=21
0863           ISIG(NCHN,2)=21
0864           ISIG(NCHN,3)=2
0865           SIGH(NCHN)=FACQQ2/2D0
0866           NCHN=NCHN+1
0867           ISIG(NCHN,1)=21
0868           ISIG(NCHN,2)=21
0869           ISIG(NCHN,3)=3
0870           SIGH(NCHN)=FACQQ3/2D0
0871   320     CONTINUE
0872
0873         ELSEIF(ISUB.EQ.246) THEN
0874 C...g + q_j -> ~chi0_1 + ~q_j
0875           FAC0=COMFAC*AS*AEM/6D0/XW
0876           ZM2=SQM4
0877           QM2=SQM3
0878           FACZQ0=FAC0*( (ZM2-TH)/SH +
0879      &    (UH-ZM2)*(UH+QM2)/(UH-QM2)**2 -
0880      &    (SH*(UH+ZM2)+2D0*(QM2-ZM2)*(ZM2-UH))/SH/(UH-QM2) )
0881           KFNSQ=MOD(KFPR(ISUBSV,1),KSUSY1)
0882           DO 340 I=-KFNSQ,KFNSQ,2*KFNSQ
0883             IF(I.LT.MMINA.OR.I.GT.MMAXA) GOTO 340
0884             IF(I.EQ.0.OR.IABS(I).GT.MSTP(58)) GOTO 340
0885             EI=KCHG(IABS(I),1)/3D0
0886             IA=IABS(I)
0887             XRQZ = -TANW*EI*ZMIX(IZID,1)
0888             XLQZ =(SIGN(1D0,EI)*ZMIX(IZID,2)-TANW*
0889      &      (SIGN(1D0,EI)-2D0*EI)*ZMIX(IZID,1))/2D0
0890             IF(ILR.EQ.0) THEN
0891               BS=XLQZ**2*SFMIX(IA,1)**2+XRQZ**2*SFMIX(IA,2)**2
0892             ELSE
0893               BS=XLQZ**2*SFMIX(IA,3)**2+XRQZ**2*SFMIX(IA,4)**2
0894             ENDIF
0895             FACZQ=FACZQ0*BS
0896             KCHQ=2
0897             IF(I.LT.0) KCHQ=3
0898             DO 330 ISDE=1,2
0899               IF(ISDE.EQ.1.AND.KFAC(1,I)*KFAC(2,21).EQ.0) GOTO 330
0900               IF(ISDE.EQ.2.AND.KFAC(1,21)*KFAC(2,I).EQ.0) GOTO 330
0901               NCHN=NCHN+1
0902               ISIG(NCHN,ISDE)=I
0903               ISIG(NCHN,3-ISDE)=21
0904               ISIG(NCHN,3)=1
0905               SIGH(NCHN)=FACZQ*RKF*WIDS(PYCOMP(KFPR(ISUBSV,1)),KCHQ)*
0906      &        WIDS(PYCOMP(KFPR(ISUBSV,2)),2)
0907   330       CONTINUE
0908   340     CONTINUE
0909         ENDIF
0910
0911       ELSEIF(ISUB.LE.260) THEN
0912         IF(ISUB.EQ.254) THEN
0913 C...g + q_j -> ~chi1_1 + ~q_i
0914           FAC0=COMFAC*AS*AEM/12D0/XW
0915           ZM2=SQM4
0916           QM2=SQM3
0917           AU=UMIX(IZID,1)**2
0918           AD=VMIX(IZID,1)**2
0919           FACZQ0=FAC0*( (ZM2-TH)/SH +
0920      &    (UH-ZM2)*(UH+QM2)/(UH-QM2)**2 -
0921      &    (SH*(UH+ZM2)+2D0*(QM2-ZM2)*(ZM2-UH))/SH/(UH-QM2) )
0922           KFNSQ1=MOD(KFPR(ISUBSV,1),KSUSY1)
0923           IF(MOD(KFNSQ1,2).EQ.0) THEN
0924             KFNSQ=KFNSQ1-1
0925             KCHW=2
0926           ELSE
0927             KFNSQ=KFNSQ1+1
0928             KCHW=3
0929           ENDIF
0930           DO 360 I=-KFNSQ,KFNSQ,2*KFNSQ
0931             IF(I.LT.MMINA.OR.I.GT.MMAXA) GOTO 360
0932             IF(I.EQ.0.OR.IABS(I).GT.MSTP(58)) GOTO 360
0933             IA=IABS(I)
0934             IF(MOD(IA,2).EQ.0) THEN
0935               FACZQ=FACZQ0*AU
0936             ELSE
0937               FACZQ=FACZQ0*AD
0938             ENDIF
0939             FACZQ=FACZQ*SFMIX(KFNSQ1,1+2*ILR)**2
0940             KCHQ=2
0941             IF(I.LT.0) KCHQ=3
0942             KCHWQ=KCHW
0943             IF(I.LT.0) KCHWQ=5-KCHW
0944             DO 350 ISDE=1,2
0945               IF(ISDE.EQ.1.AND.KFAC(1,I)*KFAC(2,21).EQ.0) GOTO 350
0946               IF(ISDE.EQ.2.AND.KFAC(1,21)*KFAC(2,I).EQ.0) GOTO 350
0947               NCHN=NCHN+1
0948               ISIG(NCHN,ISDE)=I
0949               ISIG(NCHN,3-ISDE)=21
0950               ISIG(NCHN,3)=1
0951               SIGH(NCHN)=FACZQ*RKF*WIDS(PYCOMP(KFPR(ISUBSV,1)),KCHQ)*
0952      &        WIDS(PYCOMP(KFPR(ISUBSV,2)),KCHWQ)
0953   350       CONTINUE
0954   360     CONTINUE
0955
0956         ELSEIF(ISUB.EQ.258) THEN
0957 C...g + q_j -> gluino + ~q_i
0958           XG2=SQM4
0959           XQ2=SQM3
0960           XMT=XG2-TH
0961           XMU=XG2-UH
0962           XST=XQ2-TH
0963           XSU=XQ2-UH
0964           FACQG1=0.5D0*4D0/9D0*XMT/SH + (XMT*SH+2D0*XG2*XST)/XMT**2 -
0965      &    ( (SH-XQ2+XG2)*(-XST)-SH*XG2 )/SH/(-XMT) +
0966      &    0.5D0*1D0/2D0*( XST*(TH+2D0*UH+XG2)-XMT*(SH-2D0*XST) +
0967      &    (-XMU)*(TH+XG2+2D0*XQ2) )/2D0/XMT/XSU
0968           FACQG2= 4D0/9D0*(-XMU)*(UH+XQ2)/XSU**2 + 1D0/18D0*
0969      &    (SH*(UH+XG2)
0970      &    +2D0*(XQ2-XG2)*XMU)/SH/(-XSU) + 0.5D0*4D0/9D0*XMT/SH +
0971      &    0.5D0*1D0/2D0*(XST*(TH+2D0*UH+XG2)-XMT*(SH-2D0*XST)+
0972      &    (-XMU)*(TH+XG2+2D0*XQ2))/2D0/XMT/XSU
0973           ASYUK=RMSS(42)*AS
0974           FACQG1=COMFAC*AS*ASYUK*FACQG1/2D0
0975           FACQG2=COMFAC*AS*ASYUK*FACQG2/2D0
0976           KFNSQ=MOD(KFPR(ISUBSV,1),KSUSY1)
0977           DO 380 I=-KFNSQ,KFNSQ,2*KFNSQ
0978             IF(I.LT.MMINA.OR.I.GT.MMAXA) GOTO 380
0979             IF(I.EQ.0.OR.IABS(I).GT.10) GOTO 380
0980             KCHQ=2
0981             IF(I.LT.0) KCHQ=3
0982             FACSEL=RKF*WIDS(PYCOMP(KFPR(ISUBSV,1)),KCHQ)*
0983      &      WIDS(PYCOMP(KFPR(ISUBSV,2)),2)
0984             DO 370 ISDE=1,2
0985               IF(ISDE.EQ.1.AND.KFAC(1,I)*KFAC(2,21).EQ.0) GOTO 370
0986               IF(ISDE.EQ.2.AND.KFAC(1,21)*KFAC(2,I).EQ.0) GOTO 370
0987               NCHN=NCHN+1
0988               ISIG(NCHN,ISDE)=I
0989               ISIG(NCHN,3-ISDE)=21
0990               ISIG(NCHN,3)=1
0991               SIGH(NCHN)=FACQG1*FACSEL
0992               NCHN=NCHN+1
0993               ISIG(NCHN,ISDE)=I
0994               ISIG(NCHN,3-ISDE)=21
0995               ISIG(NCHN,3)=2
0996               SIGH(NCHN)=FACQG2*FACSEL
0997   370       CONTINUE
0998   380     CONTINUE
0999         ENDIF
1000
1001       ELSEIF(ISUB.LE.270) THEN
1002         IF(ISUB.EQ.261) THEN
1003 C...q_i + q_ibar -> ~t_1 + ~t_1bar
1004           FACQQ1=COMFAC*( (UH*TH-SQM3*SQM4)/ SH**2 )*
1005      &    WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
1006           KFNSQ=MOD(KFPR(ISUBSV,1),KSUSY1)
1007           FAC0=AS**2*4D0/9D0
1008           DO 390 I=MMIN1,MMAX1
1009             IA=IABS(I)
1010             IF(I.EQ.0.OR.KFAC(1,I)*KFAC(2,-I).EQ.0) GOTO 390
1011             IF(IA.GE.11.AND.IA.LE.18) THEN
1012               EI=KCHG(IA,1)/3D0
1013               EJ=KCHG(KFNSQ,1)/3D0
1014               T3I=SIGN(1D0,EI)/2D0
1015               T3J=SIGN(1D0,EJ)/2D0
1016               XLQ=2D0*(T3J-EJ*XW)*SFMIX(KFNSQ,2*ILR+1)**2
1017               XRQ=2D0*(-EJ*XW)*SFMIX(KFNSQ,2*ILR+2)**2
1018               XLF=2D0*(T3I-EI*XW)
1019               XRF=2D0*(-EI*XW)
1020               TAA=0.5D0*(EI*EJ)**2
1021               TZZ=(XLF**2+XRF**2)*(XLQ+XRQ)**2/64D0/XW**2/XW1**2
1022               TZZ=TZZ/((1D0-SQMZ/SH)**2+SQMZ*(ZWID/SH)**2)
1023               TAZ=EI*EJ*(XLQ+XRQ)*(XLF+XRF)/8D0/XW/XW1
1024               TAZ=TAZ/((1D0-SQMZ/SH)**2+SQMZ*(ZWID/SH)**2)*(1D0-SQMZ/SH)
1025               FAC0=AEM**2*12D0*(TAA+TZZ+TAZ)
1026             ENDIF
1027             NCHN=NCHN+1
1028             ISIG(NCHN,1)=I
1029             ISIG(NCHN,2)=-I
1030             ISIG(NCHN,3)=1
1031             SIGH(NCHN)=FACQQ1*FAC0
1032   390     CONTINUE
1033
1034         ELSEIF(ISUB.EQ.263) THEN
1035 C...f + fbar -> ~t1 + ~t2bar
1036           DO 400 I=MMIN1,MMAX1
1037             IA=IABS(I)
1038             IF(I.EQ.0.OR.KFAC(1,I)*KFAC(2,-I).EQ.0) GOTO 400
1039             EI=KCHG(IABS(I),1)/3D0
1040             TT3I=SIGN(1D0,EI)/2D0
1041             EJ=2D0/3D0
1042             TT3J=1D0/2D0
1043             FCOL=1D0
1044 C...Color factor for e+ e-
1045             IF(IA.GE.11) FCOL=3D0
1046             XLQ=2D0*(TT3J-EJ*XW)
1047             XRQ=2D0*(-EJ*XW)
1048             XLF=2D0*(TT3I-EI*XW)
1049             XRF=2D0*(-EI*XW)
1050             TZZ=(XLF**2+XRF**2)*(XLQ-XRQ)**2/64D0/XW**2/XW1**2
1051             TZZ=TZZ*(SFMIX(6,1)*SFMIX(6,2))**2
1052             TZZ=TZZ/((1D0-SQMZ/SH)**2+SQMZ*(ZWID/SH)**2)
1053 C...Factor of 2 for t1 t2bar + t2 t1bar
1054             FACQQ1=2D0*COMFAC*AEM**2*TZZ*FCOL*4D0
1055             FACQQ1=FACQQ1*( UH*TH-SQM3*SQM4 )/SH2
1056             NCHN=NCHN+1
1057             ISIG(NCHN,1)=I
1058             ISIG(NCHN,2)=-I
1059             ISIG(NCHN,3)=1
1060             SIGH(NCHN)=FACQQ1*WIDS(PYCOMP(KFPR(ISUBSV,1)),2)*
1061      &      WIDS(PYCOMP(KFPR(ISUBSV,2)),3)
1062             NCHN=NCHN+1
1063             ISIG(NCHN,1)=I
1064             ISIG(NCHN,2)=-I
1065             ISIG(NCHN,3)=2
1066             SIGH(NCHN)=FACQQ1*WIDS(PYCOMP(KFPR(ISUBSV,1)),3)*
1067      &      WIDS(PYCOMP(KFPR(ISUBSV,2)),2)
1068   400     CONTINUE
1069
1070         ELSEIF(ISUB.EQ.264) THEN
1071 C...g + g -> ~t_1 + ~t_1bar
1072           XSU=SQM3-UH
1073           XST=SQM3-TH
1074           FAC0=COMFAC*AS**2*(7D0/48D0+3D0*(UH-TH)**2/16D0/SH2 )*0.5D0*
1075      &    WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
1076           FACQQ1=FAC0*(0.5D0+2D0*SQM3*TH/XST**2 + 2D0*SQM3**2/XSU/XST)
1077           FACQQ2=FAC0*(0.5D0+2D0*SQM3*UH/XSU**2 + 2D0*SQM3**2/XSU/XST)
1078           IF(KFAC(1,21)*KFAC(2,21).EQ.0) GOTO 410
1079           NCHN=NCHN+1
1080           ISIG(NCHN,1)=21
1081           ISIG(NCHN,2)=21
1082           ISIG(NCHN,3)=1
1083           SIGH(NCHN)=FACQQ1
1084           NCHN=NCHN+1
1085           ISIG(NCHN,1)=21
1086           ISIG(NCHN,2)=21
1087           ISIG(NCHN,3)=2
1088           SIGH(NCHN)=FACQQ2
1089   410     CONTINUE
1090         ENDIF
1091
1092       ELSEIF(ISUB.LE.280) THEN
1093         IF(ISUB.EQ.271) THEN
1094 C...q + q' -> ~q + ~q' (~g exchange)
1095           XMG2=PMAS(PYCOMP(KSUSY1+21),1)**2
1096           XMT=XMG2-TH
1097           XMU=XMG2-UH
1098           XSU1=SQM3-UH
1099           XSU2=SQM4-UH
1100           XST1=SQM3-TH
1101           XST2=SQM4-TH
1102           ASYUK=RMSS(42)*AS
1103           IF(ILR.EQ.1) THEN
1104             FACQQ1=COMFAC*ASYUK**2*4D0/9D0*( -(XST1*XST2+SH*TH)/XMT**2 )
1105             FACQQ2=COMFAC*ASYUK**2*4D0/9D0*( -(XSU1*XSU2+SH*UH)/XMU**2 )
1106             FACQQB=0.0D0
1107           ELSE
1108             FACQQ1=0.5D0*COMFAC*ASYUK**2*4D0/9D0*( SH*XMG2/XMT**2 )
1109             FACQQ2=0.5D0*COMFAC*ASYUK**2*4D0/9D0*( SH*XMG2/XMU**2 )
1110             FACQQB=0.5D0*COMFAC*ASYUK**2*4D0/9D0*( -2D0*SH*XMG2/3D0/
1111      &      XMT/XMU )
1112           ENDIF
1113           KFNSQI=MOD(KFPR(ISUBSV,1),KSUSY1)
1114           KFNSQJ=MOD(KFPR(ISUBSV,2),KSUSY1)
1115           DO 430 I=-KFNSQI,KFNSQI,2*KFNSQI
1116             IF(I.LT.MMIN1.OR.I.GT.MMAX1) GOTO 430
1117             IA=IABS(I)
1118             IF(I.EQ.0.OR.IA.GT.MSTP(58).OR.KFAC(1,I).EQ.0) GOTO 430
1119             KCHQ=2
1120             IF(I.LT.0) KCHQ=3
1121             DO 420 J=-KFNSQJ,KFNSQJ,2*KFNSQJ
1122               IF(J.LT.MMIN2.OR.J.GT.MMAX2) GOTO 420
1123               JA=IABS(J)
1124               IF(J.EQ.0.OR.JA.GT.MSTP(58).OR.KFAC(2,J).EQ.0) GOTO 420
1125               IF(I*J.LT.0) GOTO 420
1126               NCHN=NCHN+1
1127               ISIG(NCHN,1)=I
1128               ISIG(NCHN,2)=J
1129               ISIG(NCHN,3)=1
1130               SIGH(NCHN)=FACQQ1*RKF*WIDS(PYCOMP(KFPR(ISUBSV,1)),KCHQ)*
1131      &        WIDS(PYCOMP(KFPR(ISUBSV,2)),KCHQ)
1132               IF(I.EQ.J) THEN
1133                 IF(ILR.EQ.0) THEN
1134                   SIGH(NCHN)=0.5D0*(FACQQ1+0.5D0*FACQQB)*RKF*
1135      &            WIDS(PYCOMP(KFPR(ISUBSV,1)),KCHQ+2)
1136                 ELSE
1137                   SIGH(NCHN)=0.5D0*FACQQ1*RKF*
1138      &            WIDS(PYCOMP(KFPR(ISUBSV,1)),KCHQ)*
1139      &            WIDS(PYCOMP(KFPR(ISUBSV,2)),KCHQ)
1140                 ENDIF
1141                 NCHN=NCHN+1
1142                 ISIG(NCHN,1)=I
1143                 ISIG(NCHN,2)=J
1144                 ISIG(NCHN,3)=2
1145                 IF(ILR.EQ.0) THEN
1146                   SIGH(NCHN)=0.5D0*(FACQQ2+0.5D0*FACQQB)*RKF*
1147      &            WIDS(PYCOMP(KFPR(ISUBSV,1)),KCHQ+2)
1148                 ELSE
1149                   SIGH(NCHN)=0.5D0*FACQQ2*RKF*
1150      &            WIDS(PYCOMP(KFPR(ISUBSV,1)),KCHQ)*
1151      &            WIDS(PYCOMP(KFPR(ISUBSV,2)),KCHQ)
1152                 ENDIF
1153               ENDIF
1154   420       CONTINUE
1155   430     CONTINUE
1156
1157         ELSEIF(ISUB.EQ.274) THEN
1158 C...q + qbar' -> ~q + ~qbar'
1159           XMG2=PMAS(PYCOMP(KSUSY1+21),1)**2
1160           XMT=XMG2-TH
1161           XMU=XMG2-UH
1162           IF(ILR.EQ.0) THEN
1163 C...Mrenna...Normalization.and.1/XMT
1164             FACQQ1=COMFAC*AS**2*2D0/9D0*(
1165      &      (UH*TH-SQM3*SQM4)/XMT**2 )*RMSS(42)**2
1166             FACQQB=COMFAC*AS**2*4D0/9D0*(
1167      &      (UH*TH-SQM3*SQM4)/SH2 )
1168             FACQQI=-COMFAC*AS**2*4D0/27D0*(
1169      &      (UH*TH-SQM3*SQM4)/SH/XMT )*RMSS(42)
1170             FACQQB=FACQQB+FACQQ1+FACQQI
1171           ELSE
1172             FACQQ1=COMFAC*AS**2*4D0/9D0*( XMG2*SH/XMT**2 )*RMSS(42)**2
1173             FACQQB=FACQQ1
1174           ENDIF
1175           KFNSQI=MOD(KFPR(ISUBSV,1),KSUSY1)
1176           KFNSQJ=MOD(KFPR(ISUBSV,2),KSUSY1)
1177           DO 450 I=-KFNSQI,KFNSQI,2*KFNSQI
1178             IF(I.LT.MMIN1.OR.I.GT.MMAX1) GOTO 450
1179             IA=IABS(I)
1180             IF(I.EQ.0.OR.IA.GT.MSTP(58).OR.KFAC(1,I).EQ.0) GOTO 450
1181             KCHQ=2
1182             IF(I.LT.0) KCHQ=3
1183             DO 440 J=-KFNSQJ,KFNSQJ,2*KFNSQJ
1184               IF(J.LT.MMIN2.OR.J.GT.MMAX2) GOTO 440
1185               JA=IABS(J)
1186               IF(J.EQ.0.OR.JA.GT.MSTP(58).OR.KFAC(2,J).EQ.0) GOTO 440
1187               IF(I*J.GT.0) GOTO 440
1188               NCHN=NCHN+1
1189               ISIG(NCHN,1)=I
1190               ISIG(NCHN,2)=J
1191               ISIG(NCHN,3)=1
1192               SIGH(NCHN)=FACQQ1*RKF*WIDS(PYCOMP(KFPR(ISUBSV,1)),KCHQ)*
1193      &        WIDS(PYCOMP(KFPR(ISUBSV,2)),5-KCHQ)
1194               IF(ILR.EQ.0.AND.I.EQ.-J) SIGH(NCHN)=FACQQB*RKF*
1195      &        WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
1196   440       CONTINUE
1197   450     CONTINUE
1198
1199         ELSEIF(ISUB.EQ.277) THEN
1200 C...q_i + q_ibar -> ~q_j + ~q_jbar ,i .ne. j
1201 C...if i .eq. j covered in 274
1202           FACQQ1=COMFAC*( (UH*TH-SQM3*SQM4)/ SH**2 )
1203           KFNSQ=MOD(KFPR(ISUBSV,1),KSUSY1)
1204           FAC0=0D0
1205           DO 460 I=MMIN1,MMAX1
1206             IA=IABS(I)
1207             IF(I.EQ.0.OR.(IA.GT.MSTP(58).AND.IA.LE.10).OR.
1208      &      KFAC(1,I)*KFAC(2,-I).EQ.0) GOTO 460
1209             IF(IA.EQ.KFNSQ) GOTO 460
1210             IF(IA.EQ.11.OR.IA.EQ.13.OR.IA.EQ.15) THEN
1211               EI=KCHG(IA,1)/3D0
1212               EJ=KCHG(KFNSQ,1)/3D0
1213               T3J=SIGN(0.5D0,EJ)
1214               T3I=SIGN(1D0,EI)/2D0
1215               IF(ILR.EQ.0) THEN
1216                 XLQ=2D0*(T3J-EJ*XW)*SFMIX(KFNSQ,1)
1217                 XRQ=2D0*(-EJ*XW)*SFMIX(KFNSQ,2)
1218               ELSE
1219                 XLQ=2D0*(T3J-EJ*XW)*SFMIX(KFNSQ,3)
1220                 XRQ=2D0*(-EJ*XW)*SFMIX(KFNSQ,4)
1221               ENDIF
1222               XLF=2D0*(T3I-EI*XW)
1223               XRF=2D0*(-EI*XW)
1224               IF(ILR.EQ.0) THEN
1225                 XRQ=0D0
1226               ELSE
1227                 XLQ=0D0
1228               ENDIF
1229               TAA=0.5D0*(EI*EJ)**2
1230               TZZ=(XLF**2+XRF**2)*(XLQ+XRQ)**2/64D0/XW**2/XW1**2
1231               TZZ=TZZ/((1D0-SQMZ/SH)**2+SQMZ*(ZWID/SH)**2)
1232               TAZ=EI*EJ*(XLQ+XRQ)*(XLF+XRF)/8D0/XW/XW1
1233               TAZ=TAZ/((1D0-SQMZ/SH)**2+SQMZ*(ZWID/SH)**2)*(1D0-SQMZ/SH)
1234               FAC0=AEM**2*12D0*(TAA+TZZ+TAZ)
1235             ELSEIF(IA.LE.6) THEN
1236               FAC0=AS**2*8D0/9D0/2D0
1237             ENDIF
1238             NCHN=NCHN+1
1239             ISIG(NCHN,1)=I
1240             ISIG(NCHN,2)=-I
1241             ISIG(NCHN,3)=1
1242             SIGH(NCHN)=FACQQ1*FAC0*RKF*WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
1243   460     CONTINUE
1244
1245         ELSEIF(ISUB.EQ.279) THEN
1246 C...g + g -> ~q_j + ~q_jbar
1247           XSU=SQM3-UH
1248           XST=SQM3-TH
1249 C...5=RKF because ~t ~tbar treated separately
1250           FAC0=RKF*COMFAC*AS**2*( 7D0/48D0+3D0*(UH-TH)**2/16D0/SH2 )
1251           FACQQ1=FAC0*(0.5D0+2D0*SQM3*TH/XST**2 + 2D0*SQM3**2/XSU/XST)
1252           FACQQ2=FAC0*(0.5D0+2D0*SQM3*UH/XSU**2 + 2D0*SQM3**2/XSU/XST)
1253           IF(KFAC(1,21)*KFAC(2,21).EQ.0) GOTO 470
1254           NCHN=NCHN+1
1255           ISIG(NCHN,1)=21
1256           ISIG(NCHN,2)=21
1257           ISIG(NCHN,3)=1
1258           SIGH(NCHN)=FACQQ1/2D0*WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
1259           NCHN=NCHN+1
1260           ISIG(NCHN,1)=21
1261           ISIG(NCHN,2)=21
1262           ISIG(NCHN,3)=2
1263           SIGH(NCHN)=FACQQ2/2D0*WIDS(PYCOMP(KFPR(ISUBSV,1)),1)
1264   470     CONTINUE
1265
1266         ENDIF
1267       ENDIF
1268 CMRENNA--
1269
1270       RETURN
1271       END