pythiaeRHIC/pythia/pynjdc.f

0001
0002 C*********************************************************************
0003
0004 C...PYNJDC
0005 C...Calculates decay widths for the neutralinos (admixtures of
0006 C...Bino, W3-ino, Higgs1-ino, Higgs2-ino)
0007
0008 C...Input:  KCIN = KF code for particle
0009 C...Output: XLAM = widths
0010 C...        IDLAM = KF codes for decay particles
0011 C...        IKNT = number of decay channels defined
0012 C...AUTHOR: STEPHEN MRENNA
0013 C...Last change:
0014 C...10-15-95:  force decay chi^0_2 -> chi^0_1 + gamma
0015 C...when CHIGAMMA .NE. 0
0016 C...10 FEB 96:  Calculate this decay for small tan(beta)
0017
0018       SUBROUTINE PYNJDC(KFIN,XLAM,IDLAM,IKNT)
0019
0020 C...Double precision and integer declarations.
0021       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0022       IMPLICIT INTEGER(I-N)
0023       INTEGER PYK,PYCHGE,PYCOMP
0024 C...Parameter statement to help give large particle numbers.
0025       PARAMETER (KSUSY1=1000000,KSUSY2=2000000,KTECHN=3000000,
0026      &KEXCIT=4000000,KDIMEN=5000000)
0027 C...Commonblocks.
0028       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0029       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
0030       COMMON/PYMSSM/IMSS(0:99),RMSS(0:99)
0031 c      COMMON/PYSSMT/ZMIX(4,4),UMIX(2,2),VMIX(2,2),SMZ(4),SMW(2),
0032 c     &SFMIX(16,4)
0033       COMMON/PYSSMT/ZMIX(4,4),UMIX(2,2),VMIX(2,2),SMZ(4),SMW(2),
0034      &SFMIX(16,4),ZMIXI(4,4),UMIXI(2,2),VMIXI(2,2)
0035 C      COMMON/PYINTS/XXM(20)
0036       COMPLEX*16 CXC
0037       COMMON/PYINTC/XXC(10),CXC(8)
0038       SAVE /PYDAT1/,/PYDAT2/,/PYMSSM/,/PYSSMT/,/PYINTC/
0039
0040 C...Local variables.
0041       COMPLEX*16 ZMIXC(4,4),VMIXC(2,2),UMIXC(2,2),OLPP,ORPP,GLIJ,GRIJ
0042       COMPLEX*16 QIJ,RIJ,F21K,F12K,CAL,CAR,CBL,CBR,CA,CB
0043       INTEGER KFIN
0044       DOUBLE PRECISION XMI,XMJ,XMF,XMSF1,XMSF2,XMW,XMW2,
0045      &XMZ,XMZ2,AXMJ,AXMI
0046       DOUBLE PRECISION S12MIN,S12MAX
0047       DOUBLE PRECISION XMI2,XMI3,XMJ2,XMH,XMH2,XMHP,XMA2,XMB2
0048       DOUBLE PRECISION PYLAMF,XL
0049       DOUBLE PRECISION TANW,XW,AEM,C1,AS,EI,T3I
0050       DOUBLE PRECISION PYX2XH,PYX2XG
0051       DOUBLE PRECISION XLAM(0:400)
0052       INTEGER IDLAM(400,3)
0053       INTEGER LKNT,IX,IH,J,IJ,I,IKNT,FID
0054       INTEGER ITH(3),KF1,KF2
0055       INTEGER ITHC
0056       DOUBLE PRECISION DH(3),EH(3)
0057       DOUBLE PRECISION SR2
0058       DOUBLE PRECISION CBETA,SBETA
0059       DOUBLE PRECISION GAMCON,XMT1,XMT2
0060       DOUBLE PRECISION PYALEM,PI,PYALPS
0061       DOUBLE PRECISION RAT1,RAT2
0062       DOUBLE PRECISION T3T,FCOL
0063       DOUBLE PRECISION ALFA,BETA,TANB
0064       DOUBLE PRECISION PYXXGA
0065       EXTERNAL PYGAUS,PYXXZ6
0066       DOUBLE PRECISION PYGAUS,PYXXZ6
0067       DOUBLE PRECISION PREC
0068       INTEGER KFNCHI(4),KFCCHI(2)
0069       DATA ITH/25,35,36/
0070       DATA ITHC/37/
0071       DATA PREC/1D-2/
0072       DATA PI/3.141592654D0/
0073       DATA SR2/1.4142136D0/
0074       DATA KFNCHI/1000022,1000023,1000025,1000035/
0075       DATA KFCCHI/1000024,1000037/
0076
0077 C...COUNT THE NUMBER OF DECAY MODES
0078       LKNT=0
0079
0080       XMW=PMAS(24,1)
0081       XMW2=XMW**2
0082       XMZ=PMAS(23,1)
0083       XMZ2=XMZ**2
0084       XW=1D0-XMW2/XMZ2
0085       XW1=1D0-XW
0086       TANW = SQRT(XW/XW1)
0087
0088 C...IX IS 1 - 4 DEPENDING ON SEQUENCE NUMBER
0089       IX=1
0090       IF(KFIN.EQ.KFNCHI(2)) IX=2
0091       IF(KFIN.EQ.KFNCHI(3)) IX=3
0092       IF(KFIN.EQ.KFNCHI(4)) IX=4
0093
0094       XMI=SMZ(IX)
0095       XMI2=XMI**2
0096       AXMI=ABS(XMI)
0097       AEM=PYALEM(XMI2)
0098       AS =PYALPS(XMI2)
0099       C1=AEM/XW
0100       XMI3=ABS(XMI**3)
0101
0102       TANB=RMSS(5)
0103       BETA=ATAN(TANB)
0104       ALFA=RMSS(18)
0105       CBETA=COS(BETA)
0106       SBETA=TANB*CBETA
0107       CALFA=COS(ALFA)
0108       SALFA=SIN(ALFA)
0109
0110       DO 110 I=1,4
0111         DO 100 J=1,4
0112           ZMIXC(J,I)=DCMPLX(ZMIX(J,I),ZMIXI(J,I))
0113   100   CONTINUE
0114   110 CONTINUE
0115       DO 130 I=1,2
0116         DO 120 J=1,2
0117            VMIXC(J,I)=DCMPLX(VMIX(J,I),VMIXI(J,I))
0118            UMIXC(J,I)=DCMPLX(UMIX(J,I),UMIXI(J,I))
0119   120   CONTINUE
0120   130 CONTINUE
0121
0122 C...CHECK ALL 2-BODY DECAYS TO GAUGE AND HIGGS BOSONS
0123       IF(IX.EQ.1.AND.IMSS(11).EQ.0) GOTO 300
0124
0125 C...FORCE CHI0_2 -> CHI0_1 + GAMMA
0126       IF(IX.EQ.2 .AND. IMSS(10).NE.0 ) THEN
0127         XMJ=SMZ(1)
0128         AXMJ=ABS(XMJ)
0129         LKNT=LKNT+1
0130         GAMCON=AEM**3/8D0/PI/XMW2/XW
0131         XMT1=(PMAS(PYCOMP(KSUSY1+6),1)/PMAS(6,1))**2
0132         XMT2=(PMAS(PYCOMP(KSUSY2+6),1)/PMAS(6,1))**2
0133         XLAM(LKNT)=PYXXGA(GAMCON,AXMI,AXMJ,XMT1,XMT2)
0134         IDLAM(LKNT,1)=KSUSY1+22
0135         IDLAM(LKNT,2)=22
0136         IDLAM(LKNT,3)=0
0137         WRITE(MSTU(11),*) 'FORCED N2 -> N1 + GAMMA ',XLAM(LKNT)
0138         GOTO 340
0139       ENDIF
0140
0141 C...GRAVITINO DECAY MODES
0142
0143       IF(IMSS(11).EQ.1) THEN
0144         XMP=RMSS(29)
0145         IDG=39+KSUSY1
0146         XMGR=PMAS(PYCOMP(IDG),1)
0147         SINW=SQRT(XW)
0148         COSW=SQRT(1D0-XW)
0149         XFAC=(XMI2/(XMP*XMGR))**2*AXMI/48D0/PI
0150         IF(AXMI.GT.XMGR+PMAS(22,1)) THEN
0151           LKNT=LKNT+1
0152           IDLAM(LKNT,1)=IDG
0153           IDLAM(LKNT,2)=22
0154           IDLAM(LKNT,3)=0
0155           XLAM(LKNT)=XFAC*ABS(ZMIXC(IX,1)*COSW+ZMIXC(IX,2)*SINW)**2
0156         ENDIF
0157         IF(AXMI.GT.XMGR+XMZ) THEN
0158           LKNT=LKNT+1
0159           IDLAM(LKNT,1)=IDG
0160           IDLAM(LKNT,2)=23
0161           IDLAM(LKNT,3)=0
0162           XLAM(LKNT)=XFAC*(ABS(ZMIXC(IX,1)*SINW-ZMIXC(IX,2)*COSW)**2 +
0163      $  .5D0*ABS(ZMIXC(IX,3)*CBETA-ZMIXC(IX,4)*SBETA)**2)*
0164      &  (1D0-XMZ2/XMI2)**4
0165         ENDIF
0166         IF(AXMI.GT.XMGR+PMAS(25,1)) THEN
0167           LKNT=LKNT+1
0168           IDLAM(LKNT,1)=IDG
0169           IDLAM(LKNT,2)=25
0170           IDLAM(LKNT,3)=0
0171           XLAM(LKNT)=XFAC*(ABS(ZMIXC(IX,3)*SALFA-ZMIXC(IX,4)*CALFA)**2)*
0172      $  .5D0*(1D0-PMAS(25,1)**2/XMI2)**4
0173         ENDIF
0174         IF(AXMI.GT.XMGR+PMAS(35,1)) THEN
0175           LKNT=LKNT+1
0176           IDLAM(LKNT,1)=IDG
0177           IDLAM(LKNT,2)=35
0178           IDLAM(LKNT,3)=0
0179           XLAM(LKNT)=XFAC*(ABS(ZMIXC(IX,3)*CALFA+ZMIXC(IX,4)*SALFA)**2)*
0180      $  .5D0*(1D0-PMAS(35,1)**2/XMI2)**4
0181         ENDIF
0182         IF(AXMI.GT.XMGR+PMAS(36,1)) THEN
0183           LKNT=LKNT+1
0184           IDLAM(LKNT,1)=IDG
0185           IDLAM(LKNT,2)=36
0186           IDLAM(LKNT,3)=0
0187           XLAM(LKNT)=XFAC*(ABS(ZMIXC(IX,3)*SBETA+ZMIXC(IX,4)*CBETA)**2)*
0188      $  .5D0*(1D0-PMAS(36,1)**2/XMI2)**4
0189         ENDIF
0190         IF(IX.EQ.1) GOTO 300
0191       ENDIF
0192
0193       DO 220 IJ=1,IX-1
0194         XMJ=SMZ(IJ)
0195         AXMJ=ABS(XMJ)
0196         XMJ2=XMJ**2
0197
0198 C...CHI0_I -> CHI0_J + GAMMA
0199         IF(AXMI.GE.AXMJ.AND.SBETA/CBETA.LE.2D0) THEN
0200           RAT1=ABS(ZMIXC(IJ,1))**2+ABS(ZMIXC(IJ,2))**2
0201           RAT1=RAT1/( 1D-6+ABS(ZMIXC(IX,3))**2+ABS(ZMIXC(IX,4))**2 )
0202           RAT2=ABS(ZMIXC(IX,1))**2+ABS(ZMIXC(IX,2))**2
0203           RAT2=RAT2/( 1D-6+ABS(ZMIXC(IJ,3))**2+ABS(ZMIXC(IJ,4))**2 )
0204           IF((RAT1.GT. 0.90D0 .AND. RAT1.LT. 1.10D0) .OR.
0205      &    (RAT2.GT. 0.90D0 .AND. RAT2.LT. 1.10D0)) THEN
0206             LKNT=LKNT+1
0207             IDLAM(LKNT,1)=KFNCHI(IJ)
0208             IDLAM(LKNT,2)=22
0209             IDLAM(LKNT,3)=0
0210             GAMCON=AEM**3/8D0/PI/XMW2/XW
0211             XMT1=(PMAS(PYCOMP(KSUSY1+6),1)/PMAS(6,1))**2
0212             XMT2=(PMAS(PYCOMP(KSUSY2+6),1)/PMAS(6,1))**2
0213             XLAM(LKNT)=PYXXGA(GAMCON,AXMI,AXMJ,XMT1,XMT2)
0214           ENDIF
0215         ENDIF
0216
0217 C...CHI0_I -> CHI0_J + Z0
0218         IF(AXMI.GE.AXMJ+XMZ) THEN
0219           LKNT=LKNT+1
0220           OLPP=(ZMIXC(IX,3)*DCONJG(ZMIXC(IJ,3))-
0221      &    ZMIXC(IX,4)*DCONJG(ZMIXC(IJ,4)))/2D0
0222           ORPP=-DCONJG(OLPP)
0223           GX2=ABS(OLPP)**2+ABS(ORPP)**2
0224           GLR=DBLE(OLPP*DCONJG(ORPP))
0225           XLAM(LKNT)=PYX2XG(C1/XMW2,XMI,XMJ,XMZ,GX2,GLR)
0226           IDLAM(LKNT,1)=KFNCHI(IJ)
0227           IDLAM(LKNT,2)=23
0228           IDLAM(LKNT,3)=0
0229         ELSEIF(AXMI.GE.AXMJ) THEN
0230           XXC(1)=0D0
0231           XXC(2)=XMJ
0232           XXC(3)=0D0
0233           XXC(4)=XMI
0234           XXC(9)=XMZ
0235           XXC(10)=PMAS(23,2)
0236           OLPP=(ZMIXC(IX,3)*DCONJG(ZMIXC(IJ,3))-
0237      &    ZMIXC(IX,4)*DCONJG(ZMIXC(IJ,4)))/2D0
0238           ORPP=DCONJG(OLPP)
0239 C...CHARGED LEPTONS
0240           FID=11
0241           XXC(5)=PMAS(PYCOMP(KSUSY1+FID),1)
0242           XXC(6)=PMAS(PYCOMP(KSUSY2+FID),1)
0243           EI=KCHG(FID,1)/3D0
0244           T3I=SIGN(1D0,EI+1D-6)/2D0
0245           GLIJ=(T3I*ZMIXC(IX,2)-TANW*(T3I-EI)*ZMIXC(IX,1))*
0246      &    DCONJG(T3I*ZMIXC(IJ,2)-TANW*(T3I-EI)*ZMIXC(IJ,1))
0247           GRIJ=ZMIXC(IX,1)*DCONJG(ZMIXC(IJ,1))*(EI*TANW)**2
0248           CXC(1)=DCMPLX((T3I-EI*XW)/XW1)*OLPP
0249           CXC(2)=-GLIJ
0250           CXC(3)=-DCMPLX((T3I-EI*XW)/XW1)*ORPP
0251           CXC(4)=DCONJG(GLIJ)
0252           CXC(5)=-DCMPLX((EI*XW)/XW1)*OLPP
0253           CXC(6)=GRIJ
0254           CXC(7)=DCMPLX((EI*XW)/XW1)*ORPP
0255           CXC(8)=-DCONJG(GRIJ)
0256           S12MIN=0D0
0257           S12MAX=(AXMI-AXMJ)**2
0258           IF( XXC(5).LT.AXMI ) THEN
0259             XXC(5)=1D6
0260           ENDIF
0261           IF(XXC(6).LT.AXMI ) THEN
0262             XXC(6)=1D6
0263           ENDIF
0264           XXC(7)=XXC(5)
0265           XXC(8)=XXC(6)
0266
0267           IF(AXMI.GE.AXMJ+2D0*PMAS(11,1)) THEN
0268             LKNT=LKNT+1
0269             XLAM(LKNT)=C1**2/XMI3/(16D0*PI)*
0270      &      PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)
0271             IDLAM(LKNT,1)=KFNCHI(IJ)
0272             IDLAM(LKNT,2)=FID
0273             IDLAM(LKNT,3)=-FID
0274             IF(AXMI.GE.AXMJ+2D0*PMAS(13,1)) THEN
0275               LKNT=LKNT+1
0276               XLAM(LKNT)=XLAM(LKNT-1)
0277               IDLAM(LKNT,1)=KFNCHI(IJ)
0278               IDLAM(LKNT,2)=13
0279               IDLAM(LKNT,3)=-13
0280             ENDIF
0281           ENDIF
0282   140     CONTINUE
0283           IF(ABS(SFMIX(15,1)).GT.ABS(SFMIX(15,2))) THEN
0284             XXC(5)=PMAS(PYCOMP(KSUSY1+15),1)
0285             XXC(6)=PMAS(PYCOMP(KSUSY2+15),1)
0286           ELSE
0287             XXC(6)=PMAS(PYCOMP(KSUSY1+15),1)
0288             XXC(5)=PMAS(PYCOMP(KSUSY2+15),1)
0289           ENDIF
0290           IF( XXC(5).LT.AXMI ) THEN
0291             XXC(5)=1D6
0292           ENDIF
0293           IF(XXC(6).LT.AXMI ) THEN
0294             XXC(6)=1D6
0295           ENDIF
0296           XXC(7)=XXC(5)
0297           XXC(8)=XXC(6)
0298
0299           IF(AXMI.GE.AXMJ+2D0*PMAS(15,1)) THEN
0300             LKNT=LKNT+1
0301             XLAM(LKNT)=C1**2/XMI3/(16D0*PI)*
0302      &      PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)
0303             IDLAM(LKNT,1)=KFNCHI(IJ)
0304             IDLAM(LKNT,2)=15
0305             IDLAM(LKNT,3)=-15
0306           ENDIF
0307
0308 C...NEUTRINOS
0309   150     CONTINUE
0310           FID=12
0311           XXC(5)=PMAS(PYCOMP(KSUSY1+FID),1)
0312           XXC(6)=PMAS(PYCOMP(KSUSY2+FID),1)
0313           EI=KCHG(FID,1)/3D0
0314           T3I=SIGN(1D0,EI+1D-6)/2D0
0315           GLIJ=(T3I*ZMIXC(IX,2)-TANW*(T3I-EI)*ZMIXC(IX,1))*
0316      &    DCONJG(T3I*ZMIXC(IJ,2)-TANW*(T3I-EI)*ZMIXC(IJ,1))
0317           GRIJ=ZMIXC(IX,1)*DCONJG(ZMIXC(IJ,1))*(EI*TANW)**2
0318           CXC(1)=DCMPLX((T3I-EI*XW)/XW1)*OLPP
0319           CXC(2)=-GLIJ
0320           CXC(3)=-DCMPLX((T3I-EI*XW)/XW1)*ORPP
0321           CXC(4)=DCONJG(GLIJ)
0322           CXC(5)=-DCMPLX((EI*XW)/XW1)*OLPP
0323           CXC(6)=GRIJ
0324           CXC(7)=DCMPLX((EI*XW)/XW1)*ORPP
0325           CXC(8)=-DCONJG(GRIJ)
0326           S12MIN=0D0
0327           S12MAX=(AXMI-AXMJ)**2
0328           IF( XXC(5).LT.AXMI ) THEN
0329             XXC(5)=1D6
0330           ENDIF
0331           IF( XXC(6).LT.AXMI ) THEN
0332             XXC(6)=1D6
0333           ENDIF
0334           XXC(7)=XXC(5)
0335           XXC(8)=XXC(6)
0336
0337           LKNT=LKNT+1
0338           XLAM(LKNT)=C1**2/XMI3/(16D0*PI)*
0339      &    PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)
0340           IDLAM(LKNT,1)=KFNCHI(IJ)
0341           IDLAM(LKNT,2)=12
0342           IDLAM(LKNT,3)=-12
0343           LKNT=LKNT+1
0344           XLAM(LKNT)=XLAM(LKNT-1)
0345           IDLAM(LKNT,1)=KFNCHI(IJ)
0346           IDLAM(LKNT,2)=14
0347           IDLAM(LKNT,3)=-14
0348   160     CONTINUE
0349
0350           IF(PMAS(PYCOMP(KSUSY1+16),1).NE.PMAS(PYCOMP(KSUSY1+12),1))
0351      &    THEN
0352             XXC(5)=PMAS(PYCOMP(KSUSY1+16),1)
0353             IF( XXC(5).LT.AXMI ) THEN
0354               XXC(5)=1D6
0355             ENDIF
0356             XXC(7)=XXC(5)
0357             LKNT=LKNT+1
0358             XLAM(LKNT)=C1**2/XMI3/(16D0*PI)*
0359      &      PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)
0360           ELSE
0361             LKNT=LKNT+1
0362             XLAM(LKNT)=XLAM(LKNT-1)
0363           ENDIF
0364           IDLAM(LKNT,1)=KFNCHI(IJ)
0365           IDLAM(LKNT,2)=16
0366           IDLAM(LKNT,3)=-16
0367 C...D-TYPE QUARKS
0368   170     CONTINUE
0369           FID=1
0370           XXC(5)=PMAS(PYCOMP(KSUSY1+FID),1)
0371           XXC(6)=PMAS(PYCOMP(KSUSY2+FID),1)
0372           EI=KCHG(FID,1)/3D0
0373           T3I=SIGN(1D0,EI+1D-6)/2D0
0374           GLIJ=(T3I*ZMIXC(IX,2)-TANW*(T3I-EI)*ZMIXC(IX,1))*
0375      &    DCONJG(T3I*ZMIXC(IJ,2)-TANW*(T3I-EI)*ZMIXC(IJ,1))
0376           GRIJ=ZMIXC(IX,1)*DCONJG(ZMIXC(IJ,1))*(EI*TANW)**2
0377           CXC(1)=DCMPLX((T3I-EI*XW)/XW1)*OLPP
0378           CXC(2)=-GLIJ
0379           CXC(3)=-DCMPLX((T3I-EI*XW)/XW1)*ORPP
0380           CXC(4)=DCONJG(GLIJ)
0381           CXC(5)=-DCMPLX((EI*XW)/XW1)*OLPP
0382           CXC(6)=GRIJ
0383           CXC(7)=DCMPLX((EI*XW)/XW1)*ORPP
0384           CXC(8)=-DCONJG(GRIJ)
0385           S12MIN=0D0
0386           S12MAX=(AXMI-AXMJ)**2
0387           IF( XXC(5).LT.AXMI ) THEN
0388             XXC(5)=1D6
0389           ENDIF
0390           IF( XXC(6).LT.AXMI ) THEN
0391             XXC(6)=1D6
0392           ENDIF
0393           XXC(7)=XXC(5)
0394           XXC(8)=XXC(6)
0395
0396           IF(AXMI.GE.AXMJ+2D0*PMAS(1,1)) THEN
0397             LKNT=LKNT+1
0398             XLAM(LKNT)=C1**2/XMI3/(16D0*PI)*
0399      &      PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)*3D0
0400             IDLAM(LKNT,1)=KFNCHI(IJ)
0401             IDLAM(LKNT,2)=1
0402             IDLAM(LKNT,3)=-1
0403             IF(AXMI.GE.AXMJ+2D0*PMAS(3,1)) THEN
0404               LKNT=LKNT+1
0405               XLAM(LKNT)=XLAM(LKNT-1)
0406               IDLAM(LKNT,1)=KFNCHI(IJ)
0407               IDLAM(LKNT,2)=3
0408               IDLAM(LKNT,3)=-3
0409             ENDIF
0410           ENDIF
0411   180     CONTINUE
0412           IF(ABS(SFMIX(5,1)).GT.ABS(SFMIX(5,2))) THEN
0413             XXC(5)=PMAS(PYCOMP(KSUSY1+5),1)
0414             XXC(6)=PMAS(PYCOMP(KSUSY2+5),1)
0415           ELSE
0416             XXC(6)=PMAS(PYCOMP(KSUSY1+5),1)
0417             XXC(5)=PMAS(PYCOMP(KSUSY2+5),1)
0418           ENDIF
0419           IF( XXC(5).LT.AXMI .AND. XXC(6).LT.AXMI ) GOTO 190
0420           IF(XXC(5).LT.AXMI) THEN
0421             XXC(5)=1D6
0422           ELSEIF(XXC(6).LT.AXMI) THEN
0423             XXC(6)=1D6
0424           ENDIF
0425           XXC(7)=XXC(5)
0426           XXC(8)=XXC(6)
0427           IF(AXMI.GE.AXMJ+2D0*PMAS(5,1)) THEN
0428             LKNT=LKNT+1
0429             XLAM(LKNT)=C1**2/XMI3/(16D0*PI)*
0430      &      PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)*3D0
0431             IDLAM(LKNT,1)=KFNCHI(IJ)
0432             IDLAM(LKNT,2)=5
0433             IDLAM(LKNT,3)=-5
0434           ENDIF
0435
0436 C...U-TYPE QUARKS
0437   190     CONTINUE
0438           FID=2
0439           XXC(5)=PMAS(PYCOMP(KSUSY1+FID),1)
0440           XXC(6)=PMAS(PYCOMP(KSUSY2+FID),1)
0441           EI=KCHG(FID,1)/3D0
0442           T3I=SIGN(1D0,EI+1D-6)/2D0
0443           GLIJ=(T3I*ZMIXC(IX,2)-TANW*(T3I-EI)*ZMIXC(IX,1))*
0444      &    DCONJG(T3I*ZMIXC(IJ,2)-TANW*(T3I-EI)*ZMIXC(IJ,1))
0445           GRIJ=ZMIXC(IX,1)*DCONJG(ZMIXC(IJ,1))*(EI*TANW)**2
0446           CXC(1)=DCMPLX((T3I-EI*XW)/XW1)*OLPP
0447           CXC(2)=-GLIJ
0448           CXC(3)=-DCMPLX((T3I-EI*XW)/XW1)*ORPP
0449           CXC(4)=DCONJG(GLIJ)
0450           CXC(5)=-DCMPLX((EI*XW)/XW1)*OLPP
0451           CXC(6)=GRIJ
0452           CXC(7)=DCMPLX((EI*XW)/XW1)*ORPP
0453           CXC(8)=-DCONJG(GRIJ)
0454
0455           IF( XXC(5).LT.AXMI .AND. XXC(6).LT.AXMI ) GOTO 200
0456           IF(XXC(5).LT.AXMI) THEN
0457             XXC(5)=1D6
0458           ELSEIF(XXC(6).LT.AXMI) THEN
0459             XXC(6)=1D6
0460           ENDIF
0461           XXC(7)=XXC(5)
0462           XXC(8)=XXC(6)
0463
0464           IF(AXMI.GE.AXMJ+2D0*PMAS(2,1)) THEN
0465             LKNT=LKNT+1
0466             XLAM(LKNT)=C1**2/XMI3/(16D0*PI)*
0467      &      PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)*3D0
0468             IDLAM(LKNT,1)=KFNCHI(IJ)
0469             IDLAM(LKNT,2)=2
0470             IDLAM(LKNT,3)=-2
0471             IF(AXMI.GE.AXMJ+2D0*PMAS(4,1)) THEN
0472               LKNT=LKNT+1
0473               XLAM(LKNT)=XLAM(LKNT-1)
0474               IDLAM(LKNT,1)=KFNCHI(IJ)
0475               IDLAM(LKNT,2)=4
0476               IDLAM(LKNT,3)=-4
0477             ENDIF
0478           ENDIF
0479   200     CONTINUE
0480         ENDIF
0481
0482 C...CHI0_I -> CHI0_J + H0_K
0483         EH(1)=SIN(ALFA)
0484         EH(2)=COS(ALFA)
0485         EH(3)=-SIN(BETA)
0486         DH(1)=COS(ALFA)
0487         DH(2)=-SIN(ALFA)
0488         DH(3)=COS(BETA)
0489         QIJ=ZMIXC(IX,3)*DCONJG(ZMIXC(IJ,2))+
0490      &  DCONJG(ZMIXC(IJ,3))*ZMIXC(IX,2)-
0491      &  TANW*(ZMIXC(IX,3)*DCONJG(ZMIXC(IJ,1))+
0492      &  DCONJG(ZMIXC(IJ,3))*ZMIXC(IX,1))
0493         RIJ=DCONJG(ZMIXC(IX,4))*ZMIXC(IJ,2)+
0494      &  ZMIXC(IJ,4)*DCONJG(ZMIXC(IX,2))-
0495      &  TANW*(DCONJG(ZMIXC(IX,4))*ZMIXC(IJ,1)+
0496      &  ZMIXC(IJ,4)*DCONJG(ZMIXC(IX,1)))
0497         DO 210 IH=1,3
0498           XMH=PMAS(ITH(IH),1)
0499           XMH2=XMH**2
0500           IF(AXMI.GE.AXMJ+XMH) THEN
0501             LKNT=LKNT+1
0502             XL=PYLAMF(XMI2,XMJ2,XMH2)
0503             F21K=0.5D0*(QIJ*EH(IH)+RIJ*DH(IH))
0504             F12K=F21K
0505 C...SIGN OF MASSES I,J
0506             XMK=XMJ
0507             IF(IH.EQ.3) XMK=-XMK
0508             GX2=ABS(F21K)**2+ABS(F12K)**2
0509             GLR=DBLE(F21K*DCONJG(F12K))
0510             XLAM(LKNT)=PYX2XH(C1,XMI,XMK,XMH,GX2,GLR)
0511             IDLAM(LKNT,1)=KFNCHI(IJ)
0512             IDLAM(LKNT,2)=ITH(IH)
0513             IDLAM(LKNT,3)=0
0514           ENDIF
0515   210   CONTINUE
0516   220 CONTINUE
0517
0518 C...CHI0_I -> CHI+_J + W-
0519       DO 260 IJ=1,2
0520         XMJ=SMW(IJ)
0521         AXMJ=ABS(XMJ)
0522         XMJ2=XMJ**2
0523         IF(AXMI.GE.AXMJ+XMW) THEN
0524           LKNT=LKNT+1
0525           CXC(1)=(DCONJG(ZMIXC(IX,2))*VMIXC(IJ,1)-
0526      &    DCONJG(ZMIXC(IX,4))*VMIXC(IJ,2)/SR2)
0527           CXC(3)=(ZMIXC(IX,2)*DCONJG(UMIXC(IJ,1))+
0528      &    ZMIXC(IX,3)*DCONJG(UMIXC(IJ,2))/SR2)
0529           GX2=ABS(CXC(1))**2+ABS(CXC(3))**2
0530           GLR=DBLE(CXC(1)*DCONJG(CXC(3)))
0531           XLAM(LKNT)=PYX2XG(C1/XMW2,XMI,XMJ,XMW,GX2,GLR)
0532           IDLAM(LKNT,1)=KFCCHI(IJ)
0533           IDLAM(LKNT,2)=-24
0534           IDLAM(LKNT,3)=0
0535           LKNT=LKNT+1
0536           XLAM(LKNT)=XLAM(LKNT-1)
0537           IDLAM(LKNT,1)=-KFCCHI(IJ)
0538           IDLAM(LKNT,2)=24
0539           IDLAM(LKNT,3)=0
0540         ELSEIF(AXMI.GE.AXMJ) THEN
0541           S12MIN=0D0
0542           S12MAX=(AXMI-AXMJ)**2
0543           RT2I = 1D0/SQRT(2D0)
0544           CXC(1)=(DCONJG(ZMIXC(IX,2))*VMIXC(IJ,1)-
0545      &    DCONJG(ZMIXC(IX,4))*VMIXC(IJ,2)*RT2I)*RT2I
0546           CXC(3)=(ZMIXC(IX,2)*DCONJG(UMIXC(IJ,1))+
0547      &    ZMIXC(IX,3)*DCONJG(UMIXC(IJ,2))*RT2I)*RT2I
0548           CXC(5)=DCMPLX(0D0,0D0)
0549           CXC(7)=DCMPLX(0D0,0D0)
0550           IA=11
0551           JA=12
0552           EI=KCHG(IA,1)/3D0
0553           T3I=SIGN(1D0,EI+1D-6)/2D0
0554           EJ=KCHG(JA,1)/3D0
0555           T3J=SIGN(1D0,EJ+1D-6)/2D0
0556           CXC(2)=VMIXC(IJ,1)*DCONJG(ZMIXC(IX,1)*(EJ-T3J)*
0557      &    TANW+ZMIXC(IX,2)*T3J)*RT2I
0558           CXC(4)=-DCONJG(UMIXC(IJ,1))*(
0559      &    ZMIXC(IX,1)*(EI-T3I)*TANW+ZMIXC(IX,2)*T3I)*RT2I
0560           CXC(6)=DCMPLX(0D0,0D0)
0561           CXC(8)=DCMPLX(0D0,0D0)
0562           XXC(1)=0D0
0563           XXC(2)=XMJ
0564           XXC(3)=0D0
0565           XXC(4)=XMI
0566           XXC(5)=PMAS(PYCOMP(KSUSY1+JA),1)
0567           XXC(6)=PMAS(PYCOMP(KSUSY1+IA),1)
0568           XXC(9)=PMAS(24,1)
0569           XXC(10)=PMAS(24,2)
0570           IF( XXC(5).LT.AXMI .AND. XXC(6).LT.AXMI ) GOTO 230
0571           IF(XXC(5).LT.AXMI) THEN
0572             XXC(5)=1D6
0573           ELSEIF(XXC(6).LT.AXMI) THEN
0574             XXC(6)=1D6
0575           ENDIF
0576           XXC(7)=XXC(6)
0577           XXC(8)=XXC(5)
0578           IF(AXMI.GE.AXMJ+PMAS(11,1)+PMAS(12,1)) THEN
0579             LKNT=LKNT+1
0580             XLAM(LKNT)=C1**2/XMI3/(16D0*PI)*
0581      &      PYGAUS(PYXXZ6,S12MIN,S12MAX,PREC)
0582             IDLAM(LKNT,1)=KFCCHI(IJ)
0583             IDLAM(LKNT,2)=11
0584             IDLAM(LKNT,3)=-12
0585             LKNT=LKNT+1
0586             XLAM(LKNT)=XLAM(LKNT-1)
0587             IDLAM(LKNT,1)=-IDLAM(LKNT-1,1)
0588             IDLAM(LKNT,2)=-IDLAM(LKNT-1,2)
0589             IDLAM(LKNT,3)=-IDLAM(LKNT-1,3)
0590             IF(AXMI.GE.AXMJ+PMAS(13,1)+PMAS(14,1)) THEN
0591               LKNT=LKNT+1
0592               XLAM(LKNT)=XLAM(LKNT-1)
0593               IDLAM(LKNT,1)=KFCCHI(IJ)
0594               IDLAM(LKNT,2)=13
0595               IDLAM(LKNT,3)=-14
0596               LKNT=LKNT+1
0597               XLAM(LKNT)=XLAM(LKNT-1)
0598               IDLAM(LKNT,1)=-IDLAM(LKNT-1,1)
0599               IDLAM(LKNT,2)=-IDLAM(LKNT-1,2)
0600               IDLAM(LKNT,3)=-IDLAM(LKNT-1,3)
0601             ENDIF
0602           ENDIF
0603   230     CONTINUE
0604           IF(ABS(SFMIX(15,1)).GT.ABS(SFMIX(15,2))) THEN
0605             XXC(5)=PMAS(PYCOMP(KSUSY1+15),1)
0606             XXC(6)=PMAS(PYCOMP(KSUSY1+16),1)
0607           ELSE
0608             XXC(5)=PMAS(PYCOMP(KSUSY2+15),1)
0609             XXC(6)=PMAS(PYCOMP(KSUSY1+16),1)
0610           ENDIF
0611           IF(XXC(5).LT.AXMI) THEN
0612             XXC(5)=1D6
0613           ENDIF
0614           IF(XXC(6).LT.AXMI) THEN
0615             XXC(6)=1D6
0616           ENDIF
0617           XXC(7)=XXC(6)
0618           XXC(8)=XXC(5)
0619           IF(AXMI.GE.AXMJ+PMAS(15,1)+PMAS(16,1)) THEN
0620             LKNT=LKNT+1
0621             XLAM(LKNT)=C1**2/XMI3/(16D0*PI)*
0622      &      PYGAUS(PYXXZ6,S12MIN,S12MAX,PREC)
0623             XLAM(LKNT)=XLAM(LKNT-1)
0624             IDLAM(LKNT,1)=KFCCHI(IJ)
0625             IDLAM(LKNT,2)=15
0626             IDLAM(LKNT,3)=-16
0627             LKNT=LKNT+1
0628             XLAM(LKNT)=XLAM(LKNT-1)
0629             IDLAM(LKNT,1)=-IDLAM(LKNT-1,1)
0630             IDLAM(LKNT,2)=-IDLAM(LKNT-1,2)
0631             IDLAM(LKNT,3)=-IDLAM(LKNT-1,3)
0632           ENDIF
0633
0634 C...NOW, DO THE QUARKS
0635   240     CONTINUE
0636           IA=1
0637           JA=2
0638           EI=KCHG(IA,1)/3D0
0639           T3I=SIGN(1D0,EI+1D-6)/2D0
0640           EJ=KCHG(JA,1)/3D0
0641           T3J=SIGN(1D0,EJ+1D-6)/2D0
0642           CXC(2)=VMIXC(IJ,1)*DCONJG(ZMIXC(IX,1)*(EJ-T3J)*
0643      &    TANW+ZMIXC(IX,2)*T3J)
0644           CXC(4)=-DCONJG(UMIXC(IJ,1))*(
0645      &    ZMIXC(IX,1)*(EI-T3I)*TANW+ZMIXC(IX,2)*T3I)
0646           XXC(5)=PMAS(PYCOMP(KSUSY1+IA),1)
0647           XXC(6)=PMAS(PYCOMP(KSUSY1+JA),1)
0648           IF(XXC(5).LT.AXMI) THEN
0649             XXC(5)=1D6
0650           ENDIF
0651           IF(XXC(6).LT.AXMI) THEN
0652             XXC(6)=1D6
0653           ENDIF
0654           XXC(7)=XXC(6)
0655           XXC(8)=XXC(5)
0656           IF(AXMI.GE.AXMJ+PMAS(2,1)+PMAS(1,1)) THEN
0657             LKNT=LKNT+1
0658             XLAM(LKNT)=3D0*C1**2/XMI3/(16D0*PI)*
0659      &      PYGAUS(PYXXZ6,S12MIN,S12MAX,PREC)
0660             IDLAM(LKNT,1)=KFCCHI(IJ)
0661             IDLAM(LKNT,2)=1
0662             IDLAM(LKNT,3)=-2
0663             LKNT=LKNT+1
0664             XLAM(LKNT)=XLAM(LKNT-1)
0665             IDLAM(LKNT,1)=-IDLAM(LKNT-1,1)
0666             IDLAM(LKNT,2)=-IDLAM(LKNT-1,2)
0667             IDLAM(LKNT,3)=-IDLAM(LKNT-1,3)
0668             IF(AXMI.GE.AXMJ+PMAS(3,1)+PMAS(4,1)) THEN
0669               LKNT=LKNT+1
0670               XLAM(LKNT)=XLAM(LKNT-1)
0671               IDLAM(LKNT,1)=KFCCHI(IJ)
0672               IDLAM(LKNT,2)=3
0673               IDLAM(LKNT,3)=-4
0674               LKNT=LKNT+1
0675               XLAM(LKNT)=XLAM(LKNT-1)
0676               IDLAM(LKNT,1)=-IDLAM(LKNT-1,1)
0677               IDLAM(LKNT,2)=-IDLAM(LKNT-1,2)
0678               IDLAM(LKNT,3)=-IDLAM(LKNT-1,3)
0679             ENDIF
0680           ENDIF
0681   250     CONTINUE
0682         ENDIF
0683   260 CONTINUE
0684   270 CONTINUE
0685
0686 C...CHI0_I -> CHI+_I + H-
0687       DO 280 IJ=1,2
0688         XMJ=SMW(IJ)
0689         AXMJ=ABS(XMJ)
0690         XMJ2=XMJ**2
0691         XMHP=PMAS(ITHC,1)
0692         IF(AXMI.GE.AXMJ+XMHP) THEN
0693           LKNT=LKNT+1
0694           OLPP=CBETA*(ZMIXC(IX,4)*DCONJG(VMIXC(IJ,1))+(ZMIXC(IX,2)+
0695      &    ZMIXC(IX,1)*TANW)*DCONJG(VMIXC(IJ,2))/SR2)
0696           ORPP=SBETA*(DCONJG(ZMIXC(IX,3))*UMIXC(IJ,1)-
0697      &    (DCONJG(ZMIXC(IX,2))+DCONJG(ZMIXC(IX,1))*TANW)*
0698      &    UMIXC(IJ,2)/SR2)
0699           GX2=ABS(OLPP)**2+ABS(ORPP)**2
0700           GLR=DBLE(OLPP*DCONJG(ORPP))
0701           XLAM(LKNT)=PYX2XH(C1,XMI,XMJ,XMHP,GX2,GLR)
0702           IDLAM(LKNT,1)=KFCCHI(IJ)
0703           IDLAM(LKNT,2)=-ITHC
0704           IDLAM(LKNT,3)=0
0705           LKNT=LKNT+1
0706           XLAM(LKNT)=XLAM(LKNT-1)
0707           IDLAM(LKNT,1)=-IDLAM(LKNT-1,1)
0708           IDLAM(LKNT,2)=-IDLAM(LKNT-1,2)
0709           IDLAM(LKNT,3)=-IDLAM(LKNT-1,3)
0710         ELSE
0711
0712         ENDIF
0713   280 CONTINUE
0714
0715 C...2-BODY DECAYS TO FERMION SFERMION
0716       DO 290 J=1,16
0717         IF(J.GE.7.AND.J.LE.10) GOTO 290
0718         KF1=KSUSY1+J
0719         KF2=KSUSY2+J
0720         XMSF1=PMAS(PYCOMP(KF1),1)
0721         XMSF2=PMAS(PYCOMP(KF2),1)
0722         XMF=PMAS(J,1)
0723         IF(J.LE.6) THEN
0724           FCOL=3D0
0725         ELSE
0726           FCOL=1D0
0727         ENDIF
0728
0729         EI=KCHG(J,1)/3D0
0730         T3T=SIGN(1D0,EI)
0731         IF(J.EQ.12.OR.J.EQ.14.OR.J.EQ.16) T3T=1D0
0732         IF(MOD(J,2).EQ.0) THEN
0733           CBL=T3T*ZMIXC(IX,2)+TANW*ZMIXC(IX,1)*(2D0*EI-T3T)
0734           CAL=XMF*ZMIXC(IX,4)/XMW/SBETA
0735           CAR=-2D0*EI*TANW*ZMIXC(IX,1)
0736           CBR=CAL
0737         ELSE
0738           CBL=T3T*ZMIXC(IX,2)+TANW*ZMIXC(IX,1)*(2D0*EI-T3T)
0739           CAL=XMF*ZMIXC(IX,3)/XMW/CBETA
0740           CAR=-2D0*EI*TANW*ZMIXC(IX,1)
0741           CBR=CAL
0742         ENDIF
0743
0744 C...D~ D_L
0745         IF(AXMI.GE.XMF+XMSF1) THEN
0746           LKNT=LKNT+1
0747           XMA2=XMSF1**2
0748           XMB2=XMF**2
0749           XL=PYLAMF(XMI2,XMA2,XMB2)
0750           CA=CAL*SFMIX(J,1)+CAR*SFMIX(J,2)
0751           CB=CBL*SFMIX(J,1)+CBR*SFMIX(J,2)
0752           XLAM(LKNT)=0.5D0*FCOL*C1/8D0/XMI3*SQRT(XL)*( (XMI2+XMB2-XMA2)*
0753      &    (ABS(CA)**2+ABS(CB)**2)+4D0*DBLE(CA*DCONJG(CB))*XMF*XMI)
0754           IDLAM(LKNT,1)=KF1
0755           IDLAM(LKNT,2)=-J
0756           IDLAM(LKNT,3)=0
0757           LKNT=LKNT+1
0758           XLAM(LKNT)=XLAM(LKNT-1)
0759           IDLAM(LKNT,1)=-IDLAM(LKNT-1,1)
0760           IDLAM(LKNT,2)=-IDLAM(LKNT-1,2)
0761           IDLAM(LKNT,3)=0
0762         ENDIF
0763
0764 C...D~ D_R
0765         IF(AXMI.GE.XMF+XMSF2) THEN
0766           LKNT=LKNT+1
0767           XMA2=XMSF2**2
0768           XMB2=XMF**2
0769           CA=CAL*SFMIX(J,3)+CAR*SFMIX(J,4)
0770           CB=CBL*SFMIX(J,3)+CBR*SFMIX(J,4)
0771           XL=PYLAMF(XMI2,XMA2,XMB2)
0772           XLAM(LKNT)=0.5D0*FCOL*C1/8D0/XMI3*SQRT(XL)*( (XMI2+XMB2-XMA2)*
0773      &    (ABS(CA)**2+ABS(CB)**2)+4D0*DBLE(CA*DCONJG(CB))*XMF*XMI)
0774           IDLAM(LKNT,1)=KF2
0775           IDLAM(LKNT,2)=-J
0776           IDLAM(LKNT,3)=0
0777           LKNT=LKNT+1
0778           XLAM(LKNT)=XLAM(LKNT-1)
0779           IDLAM(LKNT,1)=-IDLAM(LKNT-1,1)
0780           IDLAM(LKNT,2)=-IDLAM(LKNT-1,2)
0781           IDLAM(LKNT,3)=0
0782         ENDIF
0783   290 CONTINUE
0784   300 CONTINUE
0785 C...3-BODY DECAY TO Q Q~ GLUINO
0786       XMJ=PMAS(PYCOMP(KSUSY1+21),1)
0787       IF(AXMI.GE.XMJ) THEN
0788         RT2I = 1D0/SQRT(2D0)
0789         OLPP=DCMPLX(COS(RMSS(32)),SIN(RMSS(32)))*RT2I
0790         ORPP=DCONJG(OLPP)
0791         AXMJ=ABS(XMJ)
0792         XXC(1)=0D0
0793         XXC(2)=XMJ
0794         XXC(3)=0D0
0795         XXC(4)=XMI
0796         FID=1
0797         XXC(5)=PMAS(PYCOMP(KSUSY1+FID),1)
0798         XXC(6)=PMAS(PYCOMP(KSUSY2+FID),1)
0799         IF( XXC(5).LT.AXMI .OR. XXC(6).LT.AXMI ) GOTO 310
0800         XXC(7)=XXC(5)
0801         XXC(8)=XXC(6)
0802         XXC(9)=1D6
0803         XXC(10)=0D0
0804         EI=KCHG(FID,1)/3D0
0805         T3I=SIGN(1D0,EI+1D-6)/2D0
0806         GLIJ=(T3I*ZMIXC(IX,2)-TANW*(T3I-EI)*ZMIXC(IX,1))*OLPP
0807         GRIJ=ZMIXC(IX,1)*(EI*TANW)*ORPP
0808         CXC(1)=0D0
0809         CXC(2)=-GLIJ
0810         CXC(3)=0D0
0811         CXC(4)=DCONJG(GLIJ)
0812         CXC(5)=0D0
0813         CXC(6)=GRIJ
0814         CXC(7)=0D0
0815         CXC(8)=-DCONJG(GRIJ)
0816         S12MIN=0D0
0817         S12MAX=(AXMI-AXMJ)**2
0818 C...ALL QUARKS BUT T
0819         IF(AXMI.GE.AXMJ+2D0*PMAS(1,1)) THEN
0820           LKNT=LKNT+1
0821           XLAM(LKNT)=4D0*C1*AS/XMI3/(16D0*PI)*
0822      &    PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)
0823           IDLAM(LKNT,1)=KSUSY1+21
0824           IDLAM(LKNT,2)=1
0825           IDLAM(LKNT,3)=-1
0826           IF(AXMI.GE.AXMJ+2D0*PMAS(3,1)) THEN
0827             LKNT=LKNT+1
0828             XLAM(LKNT)=XLAM(LKNT-1)
0829             IDLAM(LKNT,1)=KSUSY1+21
0830             IDLAM(LKNT,2)=3
0831             IDLAM(LKNT,3)=-3
0832           ENDIF
0833         ENDIF
0834   310   CONTINUE
0835         IF(ABS(SFMIX(5,1)).GT.ABS(SFMIX(5,2))) THEN
0836           XXC(5)=PMAS(PYCOMP(KSUSY1+5),1)
0837           XXC(6)=PMAS(PYCOMP(KSUSY2+5),1)
0838         ELSE
0839           XXC(6)=PMAS(PYCOMP(KSUSY1+5),1)
0840           XXC(5)=PMAS(PYCOMP(KSUSY2+5),1)
0841         ENDIF
0842         IF( XXC(5).LT.AXMI .OR. XXC(6).LT.AXMI ) GOTO 320
0843         XXC(7)=XXC(5)
0844         XXC(8)=XXC(6)
0845         IF(AXMI.GE.AXMJ+2D0*PMAS(5,1)) THEN
0846           LKNT=LKNT+1
0847           XLAM(LKNT)=0.5D0*C1*AS/XMI3/(16D0*PI)*
0848      &    PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)
0849           IDLAM(LKNT,1)=KSUSY1+21
0850           IDLAM(LKNT,2)=5
0851           IDLAM(LKNT,3)=-5
0852         ENDIF
0853 C...U-TYPE QUARKS
0854   320   CONTINUE
0855         FID=2
0856         XXC(5)=PMAS(PYCOMP(KSUSY1+FID),1)
0857         XXC(6)=PMAS(PYCOMP(KSUSY2+FID),1)
0858         IF( XXC(5).LT.AXMI .OR. XXC(6).LT.AXMI ) GOTO 330
0859         XXC(7)=XXC(5)
0860         XXC(8)=XXC(6)
0861         EI=KCHG(FID,1)/3D0
0862         T3I=SIGN(1D0,EI+1D-6)/2D0
0863         GLIJ=(T3I*ZMIXC(IX,2)-TANW*(T3I-EI)*ZMIXC(IX,1))*OLPP
0864         GRIJ=ZMIXC(IX,1)*(EI*TANW)*ORPP
0865         CXC(2)=-GLIJ
0866         CXC(4)=DCONJG(GLIJ)
0867         CXC(6)=GRIJ
0868         CXC(8)=-DCONJG(GRIJ)
0869         IF(AXMI.GE.AXMJ+2D0*PMAS(2,1)) THEN
0870           LKNT=LKNT+1
0871           XLAM(LKNT)=0.5D0*C1*AS/XMI3/(16D0*PI)*
0872      &    PYGAUS(PYXXZ6,S12MIN,S12MAX,1D-3)
0873           IDLAM(LKNT,1)=KSUSY1+21
0874           IDLAM(LKNT,2)=2
0875           IDLAM(LKNT,3)=-2
0876           IF(AXMI.GE.AXMJ+2D0*PMAS(4,1)) THEN
0877             LKNT=LKNT+1
0878             XLAM(LKNT)=XLAM(LKNT-1)
0879             IDLAM(LKNT,1)=KSUSY1+21
0880             IDLAM(LKNT,2)=4
0881             IDLAM(LKNT,3)=-4
0882           ENDIF
0883         ENDIF
0884   330   CONTINUE
0885       ENDIF
0886
0887 C...R-violating decay modes (SKANDS).
0888       CALL PYRVNE(KFIN,XLAM,IDLAM,LKNT)
0889
0890   340 IKNT=LKNT
0891       XLAM(0)=0D0
0892       DO 350 I=1,IKNT
0893         IF(XLAM(I).LT.0D0) XLAM(I)=0D0
0894         XLAM(0)=XLAM(0)+XLAM(I)
0895   350 CONTINUE
0896       IF(XLAM(0).EQ.0D0) XLAM(0)=1D-6
0897
0898       RETURN
0899       END