pythiaeRHIC/pythia/pysfdc.f

0001
0002 C*********************************************************************
0003
0004 C...PYSFDC
0005 C...Calculates decays of sfermions.
0006
0007       SUBROUTINE PYSFDC(KFIN,XLAM,IDLAM,IKNT)
0008
0009 C...Double precision and integer declarations.
0010       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0011       IMPLICIT INTEGER(I-N)
0012       INTEGER PYK,PYCHGE,PYCOMP
0013 C...Parameter statement to help give large particle numbers.
0014       PARAMETER (KSUSY1=1000000,KSUSY2=2000000,KTECHN=3000000,
0015      &KEXCIT=4000000,KDIMEN=5000000)
0016 C...Commonblocks.
0017       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0018       COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
0019       COMMON/PYMSSM/IMSS(0:99),RMSS(0:99)
0020       COMMON/PYSSMT/ZMIX(4,4),UMIX(2,2),VMIX(2,2),SMZ(4),SMW(2),
0021      &SFMIX(16,4),ZMIXI(4,4),UMIXI(2,2),VMIXI(2,2)
0022       SAVE /PYDAT1/,/PYDAT2/,/PYMSSM/,/PYSSMT/
0023
0024 C...Local variables.
0025       COMPLEX*16 ZMIXC(4,4),VMIXC(2,2),UMIXC(2,2)
0026       COMPLEX*16 CAL,CAR,CBL,CBR,CALP,CARP,CBLP,CBRP,CA,CB
0027       INTEGER KFIN,KCIN
0028       DOUBLE PRECISION XMI,XMJ,XMF,XMSF1,XMSF2,XMW,XMW2,XMZ,AXMJ
0029       DOUBLE PRECISION XMI2,XMI3,XMA2,XMB2,XMFP
0030       DOUBLE PRECISION PYLAMF,XL
0031       DOUBLE PRECISION TANW,XW,AEM,C1,AS
0032       DOUBLE PRECISION AL,AR,BL,BR
0033       DOUBLE PRECISION CH1,CH2,CH3,CH4
0034       DOUBLE PRECISION XMBOT,XMTOP
0035       DOUBLE PRECISION XLAM(0:400)
0036       INTEGER IDLAM(400,3)
0037       INTEGER LKNT,IX,ILR,IDU,J,I,IKNT,IFL,II
0038       DOUBLE PRECISION SR2
0039       DOUBLE PRECISION CBETA,SBETA
0040       DOUBLE PRECISION CW
0041       DOUBLE PRECISION BETA,ALFA,XMU,AT,AB,ATRIT,ATRIB,ATRIL
0042       DOUBLE PRECISION COSA,SINA,TANB
0043       DOUBLE PRECISION PYALEM,PI,PYALPS,EI
0044       DOUBLE PRECISION GHRR,GHLL,GHLR,XMB,BLR
0045       INTEGER IG,KF1,KF2
0046       INTEGER IGG(4),KFNCHI(4),KFCCHI(2)
0047       DATA IGG/23,25,35,36/
0048       DATA PI/3.141592654D0/
0049       DATA SR2/1.4142136D0/
0050       DATA KFNCHI/1000022,1000023,1000025,1000035/
0051       DATA KFCCHI/1000024,1000037/
0052
0053 C...COUNT THE NUMBER OF DECAY MODES
0054       LKNT=0
0055
0056 C...NO NU_R DECAYS
0057       IF(KFIN.EQ.KSUSY2+12.OR.KFIN.EQ.KSUSY2+14.OR.
0058      &KFIN.EQ.KSUSY2+16) RETURN
0059
0060       XMW=PMAS(24,1)
0061       XMW2=XMW**2
0062       XMZ=PMAS(23,1)
0063       XW=PARU(102)
0064       TANW = SQRT(XW/(1D0-XW))
0065       CW=SQRT(1D0-XW)
0066
0067       DO 110 I=1,4
0068         DO 100 J=1,4
0069           ZMIXC(J,I)=DCMPLX(ZMIX(J,I),ZMIXI(J,I))
0070   100   CONTINUE
0071   110 CONTINUE
0072       DO 130 I=1,2
0073         DO 120 J=1,2
0074            VMIXC(J,I)=DCMPLX(VMIX(J,I),VMIXI(J,I))
0075            UMIXC(J,I)=DCMPLX(UMIX(J,I),UMIXI(J,I))
0076   120   CONTINUE
0077   130 CONTINUE
0078
0079 C...KCIN
0080       KCIN=PYCOMP(KFIN)
0081 C...ILR is 1 for left and 2 for right.
0082       ILR=KFIN/KSUSY1
0083 C...IFL is matching non-SUSY flavour.
0084       IFL=MOD(KFIN,KSUSY1)
0085 C...IDU is weak isospin, 1 for down and 2 for up.
0086       IDU=2-MOD(IFL,2)
0087
0088       XMI=PMAS(KCIN,1)
0089       XMI2=XMI**2
0090       AEM=PYALEM(XMI2)
0091       AS =PYALPS(XMI2)
0092       C1=AEM/XW
0093       XMI3=XMI**3
0094       EI=KCHG(IFL,1)/3D0
0095
0096       XMBOT=PYMRUN(5,XMI2)
0097       XMTOP=PYMRUN(6,XMI2)
0098
0099       TANB=RMSS(5)
0100       BETA=ATAN(TANB)
0101       ALFA=RMSS(18)
0102       CBETA=COS(BETA)
0103       SBETA=TANB*CBETA
0104       SINA=SIN(ALFA)
0105       COSA=COS(ALFA)
0106       XMU=-RMSS(4)
0107       ATRIT=RMSS(16)
0108       ATRIB=RMSS(15)
0109       ATRIL=RMSS(17)
0110
0111 C...2-BODY DECAYS OF SFERMION -> GRAVITINO + FERMION
0112
0113       IF(IMSS(11).EQ.1) THEN
0114         XMP=RMSS(29)
0115         IDG=39+KSUSY1
0116         XMGR=PMAS(PYCOMP(IDG),1)
0117         XFAC=(XMI2/(XMP*XMGR))**2*XMI/48D0/PI
0118         IF(IFL.EQ.5) THEN
0119           XMF=XMBOT
0120         ELSEIF(IFL.EQ.6) THEN
0121           XMF=XMTOP
0122         ELSE
0123           XMF=PMAS(IFL,1)
0124         ENDIF
0125         IF(XMI.GT.XMGR+XMF) THEN
0126           LKNT=LKNT+1
0127           IDLAM(LKNT,1)=IDG
0128           IDLAM(LKNT,2)=IFL
0129           IDLAM(LKNT,3)=0
0130           XLAM(LKNT)=XFAC*(1D0-XMF**2/XMI2)**4
0131         ENDIF
0132       ENDIF
0133
0134 C...2-BODY DECAYS OF SFERMION -> FERMION + GAUGE/GAUGINO
0135
0136 C...CHARGED DECAYS:
0137       DO 140 IX=1,2
0138 C...DI -> U CHI1-,CHI2-
0139         IF(IDU.EQ.1) THEN
0140           XMFP=PMAS(IFL+1,1)
0141           XMF =PMAS(IFL,1)
0142 C...UI -> D CHI1+,CHI2+
0143         ELSE
0144           XMFP=PMAS(IFL-1,1)
0145           XMF =PMAS(IFL,1)
0146         ENDIF
0147         XMJ=SMW(IX)
0148         AXMJ=ABS(XMJ)
0149         IF(XMI.GE.AXMJ+XMFP) THEN
0150           XMA2=XMJ**2
0151           XMB2=XMFP**2
0152           IF(IDU.EQ.2) THEN
0153             IF(IFL.EQ.6) THEN
0154               XMFP=XMBOT
0155               XMF =XMTOP
0156             ELSEIF(IFL.LT.6) THEN
0157               XMF=0D0
0158               XMFP=0D0
0159             ENDIF
0160             CBL=VMIXC(IX,1)
0161             CAL=-XMFP*UMIXC(IX,2)/SR2/XMW/CBETA
0162             CBR=-XMF*VMIXC(IX,2)/SR2/XMW/SBETA
0163             CAR=0D0
0164           ELSE
0165             IF(IFL.EQ.5) THEN
0166               XMF =XMBOT
0167               XMFP=XMTOP
0168             ELSEIF(IFL.LT.5) THEN
0169               XMF=0D0
0170               XMFP=0D0
0171             ENDIF
0172             CBL=UMIXC(IX,1)
0173             CAL=-XMFP*VMIXC(IX,2)/SR2/XMW/SBETA
0174             CBR=-XMF*UMIXC(IX,2)/SR2/XMW/CBETA
0175             CAR=0D0
0176           ENDIF
0177
0178           CALP=SFMIX(IFL,1)*CAL + SFMIX(IFL,2)*CAR
0179           CBLP=SFMIX(IFL,1)*CBL + SFMIX(IFL,2)*CBR
0180           CARP=SFMIX(IFL,4)*CAR + SFMIX(IFL,3)*CAL
0181           CBRP=SFMIX(IFL,4)*CBR + SFMIX(IFL,3)*CBL
0182           CAL=CALP
0183           CBL=CBLP
0184           CAR=CARP
0185           CBR=CBRP
0186
0187 C...F1 -> F` CHI
0188           IF(ILR.EQ.1) THEN
0189             CA=CAL
0190             CB=CBL
0191 C...F2 -> F` CHI
0192           ELSE
0193             CA=CAR
0194             CB=CBR
0195           ENDIF
0196           LKNT=LKNT+1
0197           XL=PYLAMF(XMI2,XMA2,XMB2)
0198 C...SPIN AVERAGE = 1/1 NOT 1/2....NO COLOR ENHANCEMENT
0199           XLAM(LKNT)=2D0*C1/8D0/XMI3*SQRT(XL)*((XMI2-XMB2-XMA2)*
0200      &    (ABS(CA)**2+ABS(CB)**2)-4D0*DBLE(CA*DCONJG(CB))*XMJ*XMFP)
0201           IDLAM(LKNT,3)=0
0202           IF(IDU.EQ.1) THEN
0203             IDLAM(LKNT,1)=-KFCCHI(IX)
0204             IDLAM(LKNT,2)=IFL+1
0205           ELSE
0206             IDLAM(LKNT,1)=KFCCHI(IX)
0207             IDLAM(LKNT,2)=IFL-1
0208           ENDIF
0209         ENDIF
0210   140 CONTINUE
0211
0212 C...NEUTRAL DECAYS
0213       DO 150 IX=1,4
0214 C...DI -> D CHI10
0215         XMF=PMAS(IFL,1)
0216         XMJ=SMZ(IX)
0217         AXMJ=ABS(XMJ)
0218         IF(XMI.GE.AXMJ+XMF) THEN
0219           XMA2=XMJ**2
0220           XMB2=XMF**2
0221           IF(IDU.EQ.1) THEN
0222             IF(IFL.EQ.5) THEN
0223               XMF=XMBOT
0224             ELSEIF(IFL.LT.5) THEN
0225               XMF=0D0
0226             ENDIF
0227             CBL=-ZMIXC(IX,2)+TANW*ZMIXC(IX,1)*(2D0*EI+1)
0228             CAL=XMF*ZMIXC(IX,3)/XMW/CBETA
0229             CAR=-2D0*EI*TANW*ZMIXC(IX,1)
0230             CBR=CAL
0231           ELSE
0232             IF(IFL.EQ.6) THEN
0233               XMF=XMTOP
0234             ELSEIF(IFL.LT.5) THEN
0235               XMF=0D0
0236             ENDIF
0237             CBL=ZMIXC(IX,2)+TANW*ZMIXC(IX,1)*(2D0*EI-1)
0238             CAL=XMF*ZMIXC(IX,4)/XMW/SBETA
0239             CAR=-2D0*EI*TANW*ZMIXC(IX,1)
0240             CBR=CAL
0241           ENDIF
0242
0243           CALP=SFMIX(IFL,1)*CAL + SFMIX(IFL,2)*CAR
0244           CBLP=SFMIX(IFL,1)*CBL + SFMIX(IFL,2)*CBR
0245           CARP=SFMIX(IFL,4)*CAR + SFMIX(IFL,3)*CAL
0246           CBRP=SFMIX(IFL,4)*CBR + SFMIX(IFL,3)*CBL
0247           CAL=CALP
0248           CBL=CBLP
0249           CAR=CARP
0250           CBR=CBRP
0251
0252 C...F1 -> F CHI
0253           IF(ILR.EQ.1) THEN
0254             CA=CAL
0255             CB=CBL
0256 C...F2 -> F CHI
0257           ELSE
0258             CA=CAR
0259             CB=CBR
0260           ENDIF
0261           LKNT=LKNT+1
0262           XL=PYLAMF(XMI2,XMA2,XMB2)
0263 C...SPIN AVERAGE = 1/1 NOT 1/2....NO COLOR ENHANCEMENT
0264           XLAM(LKNT)=C1/8D0/XMI3*SQRT(XL)*((XMI2-XMB2-XMA2)*
0265      &    (ABS(CA)**2+ABS(CB)**2)-4D0*DBLE(CA*DCONJG(CB))*XMJ*XMF)
0266           IDLAM(LKNT,1)=KFNCHI(IX)
0267           IDLAM(LKNT,2)=IFL
0268           IDLAM(LKNT,3)=0
0269         ENDIF
0270   150 CONTINUE
0271
0272 C...2-BODY DECAYS TO SM GAUGE AND HIGGS BOSONS
0273 C...IG=23,25,35,36
0274       DO 160 II=1,4
0275         IG=IGG(II)
0276         IF(ILR.EQ.1) GOTO 160
0277         XMB=PMAS(IG,1)
0278         XMSF1=PMAS(PYCOMP(KFIN-KSUSY1),1)
0279         IF(XMI.LT.XMSF1+XMB) GOTO 160
0280         IF(IG.EQ.23) THEN
0281           BL=-SIGN(.5D0,EI)/CW+EI*XW/CW
0282           BR=EI*XW/CW
0283           BLR=0D0
0284         ELSEIF(IG.EQ.25) THEN
0285           IF(IFL.EQ.5) THEN
0286             XMF=XMBOT
0287           ELSEIF(IFL.EQ.6) THEN
0288             XMF=XMTOP
0289           ELSEIF(IFL.LT.5) THEN
0290             XMF=0D0
0291           ELSE
0292             XMF=PMAS(IFL,1)
0293           ENDIF
0294           IF(IDU.EQ.2) THEN
0295             GHLL=XMZ/CW*(0.5D0-EI*XW)*(-SIN(ALFA+BETA))+
0296      &      XMF**2/XMW*COSA/SBETA
0297             GHRR=XMZ/CW*(EI*XW)*(-SIN(ALFA+BETA))+
0298      &      XMF**2/XMW*COSA/SBETA
0299           ELSE
0300             GHLL=XMZ/CW*(0.5D0-EI*XW)*(-SIN(ALFA+BETA))+
0301      &      XMF**2/XMW*(-SINA)/CBETA
0302             GHRR=XMZ/CW*(EI*XW)*(-SIN(ALFA+BETA))+
0303      &      XMF**2/XMW*(-SINA)/CBETA
0304           ENDIF
0305           IF(IFL.EQ.5) THEN
0306             AT=ATRIB
0307           ELSEIF(IFL.EQ.6) THEN
0308             AT=ATRIT
0309           ELSEIF(IFL.EQ.15) THEN
0310             AT=ATRIL
0311           ELSE
0312             AT=0D0
0313           ENDIF
0314 C.........need to complexify
0315           IF(IDU.EQ.2) THEN
0316             GHLR=XMF/2D0/XMW/SBETA*(-XMU*SINA+
0317      &      AT*COSA)
0318           ELSE
0319             GHLR=XMF/2D0/XMW/CBETA*(XMU*COSA-
0320      &      AT*SINA)
0321           ENDIF
0322           BL=GHLL
0323           BR=GHRR
0324           BLR=-GHLR
0325         ELSEIF(IG.EQ.35) THEN
0326           IF(IFL.EQ.5) THEN
0327             XMF=XMBOT
0328           ELSEIF(IFL.EQ.6) THEN
0329             XMF=XMTOP
0330           ELSEIF(IFL.LT.5) THEN
0331             XMF=0D0
0332           ELSE
0333             XMF=PMAS(IFL,1)
0334           ENDIF
0335           IF(IDU.EQ.2) THEN
0336             GHLL=XMZ/CW*(0.5D0-EI*XW)*COS(ALFA+BETA)+
0337      &      XMF**2/XMW*SINA/SBETA
0338             GHRR=XMZ/CW*(EI*XW)*COS(ALFA+BETA)+
0339      &      XMF**2/XMW*SINA/SBETA
0340           ELSE
0341             GHLL=XMZ/CW*(0.5D0-EI*XW)*COS(ALFA+BETA)+
0342      &      XMF**2/XMW*COSA/CBETA
0343             GHRR=XMZ/CW*(EI*XW)*COS(ALFA+BETA)+
0344      &      XMF**2/XMW*COSA/CBETA
0345           ENDIF
0346           IF(IFL.EQ.5) THEN
0347             AT=ATRIB
0348           ELSEIF(IFL.EQ.6) THEN
0349             AT=ATRIT
0350           ELSEIF(IFL.EQ.15) THEN
0351             AT=ATRIL
0352           ELSE
0353             AT=0D0
0354           ENDIF
0355 C.........Need to complexify
0356           IF(IDU.EQ.2) THEN
0357             GHLR=XMF/2D0/XMW/SBETA*(XMU*COSA+
0358      &      AT*SINA)
0359           ELSE
0360             GHLR=XMF/2D0/XMW/CBETA*(XMU*SINA+
0361      &      AT*COSA)
0362           ENDIF
0363           BL=GHLL
0364           BR=GHRR
0365           BLR=GHLR
0366         ELSEIF(IG.EQ.36) THEN
0367           GHLL=0D0
0368           GHRR=0D0
0369           IF(IFL.EQ.5) THEN
0370             XMF=XMBOT
0371           ELSEIF(IFL.EQ.6) THEN
0372             XMF=XMTOP
0373           ELSEIF(IFL.LT.5) THEN
0374             XMF=0D0
0375           ELSE
0376             XMF=PMAS(IFL,1)
0377           ENDIF
0378           IF(IFL.EQ.5) THEN
0379             AT=ATRIB
0380           ELSEIF(IFL.EQ.6) THEN
0381             AT=ATRIT
0382           ELSEIF(IFL.EQ.15) THEN
0383             AT=ATRIL
0384           ELSE
0385             AT=0D0
0386           ENDIF
0387 C.........Need to complexify
0388           IF(IDU.EQ.2) THEN
0389             GHLR=XMF/2D0/XMW*(-XMU+AT/TANB)
0390           ELSE
0391             GHLR=XMF/2D0/XMW/(-XMU+AT*TANB)
0392           ENDIF
0393           BL=GHLL
0394           BR=GHRR
0395           BLR=GHLR
0396         ENDIF
0397         AL=SFMIX(IFL,1)*SFMIX(IFL,3)*BL+
0398      &  SFMIX(IFL,2)*SFMIX(IFL,4)*BR+
0399      &  (SFMIX(IFL,1)*SFMIX(IFL,4)+SFMIX(IFL,3)*SFMIX(IFL,2))*BLR
0400         XL=PYLAMF(XMI2,XMSF1**2,XMB**2)
0401         LKNT=LKNT+1
0402         IF(IG.EQ.23) THEN
0403           XLAM(LKNT)=C1/4D0/XMI3*XL**1.5D0/XMB**2*AL**2
0404         ELSE
0405           XLAM(LKNT)=C1/4D0/XMI3*SQRT(XL)*AL**2
0406         ENDIF
0407         IDLAM(LKNT,3)=0
0408         IDLAM(LKNT,1)=KFIN-KSUSY1
0409         IDLAM(LKNT,2)=IG
0410   160 CONTINUE
0411
0412 C...SF -> SF' + W
0413       XMB=PMAS(24,1)
0414       IF(MOD(IFL,2).EQ.0) THEN
0415         KF1=KSUSY1+IFL-1
0416       ELSE
0417         KF1=KSUSY1+IFL+1
0418       ENDIF
0419       KF2=KF1+KSUSY1
0420       XMSF1=PMAS(PYCOMP(KF1),1)
0421       XMSF2=PMAS(PYCOMP(KF2),1)
0422       IF(XMI.GT.XMB+XMSF1) THEN
0423         IF(MOD(IFL,2).EQ.0) THEN
0424           IF(ILR.EQ.1) THEN
0425             AL=1D0/SR2*SFMIX(IFL,1)*SFMIX(IFL-1,1)
0426           ELSE
0427             AL=1D0/SR2*SFMIX(IFL,3)*SFMIX(IFL-1,1)
0428           ENDIF
0429         ELSE
0430           IF(ILR.EQ.1) THEN
0431             AL=1D0/SR2*SFMIX(IFL,1)*SFMIX(IFL+1,1)
0432           ELSE
0433             AL=1D0/SR2*SFMIX(IFL,3)*SFMIX(IFL+1,1)
0434           ENDIF
0435         ENDIF
0436         XL=PYLAMF(XMI2,XMSF1**2,XMB**2)
0437         LKNT=LKNT+1
0438         XLAM(LKNT)=C1/4D0/XMI3*XL**1.5D0/XMB**2*AL**2
0439         IDLAM(LKNT,3)=0
0440         IDLAM(LKNT,1)=KF1
0441         IDLAM(LKNT,2)=SIGN(24,KCHG(IFL,1))
0442       ENDIF
0443       IF(XMI.GT.XMB+XMSF2) THEN
0444         IF(MOD(IFL,2).EQ.0) THEN
0445           IF(ILR.EQ.1) THEN
0446             AL=1D0/SR2*SFMIX(IFL,1)*SFMIX(IFL-1,3)
0447           ELSE
0448             AL=1D0/SR2*SFMIX(IFL,3)*SFMIX(IFL-1,3)
0449           ENDIF
0450         ELSE
0451           IF(ILR.EQ.1) THEN
0452             AL=1D0/SR2*SFMIX(IFL,1)*SFMIX(IFL+1,3)
0453           ELSE
0454             AL=1D0/SR2*SFMIX(IFL,3)*SFMIX(IFL+1,3)
0455           ENDIF
0456         ENDIF
0457         XL=PYLAMF(XMI2,XMSF2**2,XMB**2)
0458         LKNT=LKNT+1
0459         XLAM(LKNT)=C1/4D0/XMI3*XL**1.5D0/XMB**2*AL**2
0460         IDLAM(LKNT,3)=0
0461         IDLAM(LKNT,1)=KF2
0462         IDLAM(LKNT,2)=SIGN(24,KCHG(IFL,1))
0463       ENDIF
0464
0465 C...SF -> SF' + HC
0466       XMB=PMAS(37,1)
0467       IF(MOD(IFL,2).EQ.0) THEN
0468         KF1=KSUSY1+IFL-1
0469       ELSE
0470         KF1=KSUSY1+IFL+1
0471       ENDIF
0472       KF2=KF1+KSUSY1
0473       XMSF1=PMAS(PYCOMP(KF1),1)
0474       XMSF2=PMAS(PYCOMP(KF2),1)
0475       IF(XMI.GT.XMB+XMSF1) THEN
0476         XMF=0D0
0477         XMFP=0D0
0478         AT=0D0
0479         AB=0D0
0480         IF(MOD(IFL,2).EQ.0) THEN
0481 C...T1-> B1 HC
0482           IF(ILR.EQ.1) THEN
0483             CH1=-SFMIX(IFL,1)*SFMIX(IFL-1,1)
0484             CH2= SFMIX(IFL,2)*SFMIX(IFL-1,2)
0485             CH3=-SFMIX(IFL,1)*SFMIX(IFL-1,2)
0486             CH4=-SFMIX(IFL,2)*SFMIX(IFL-1,1)
0487 C...T2-> B1 HC
0488           ELSE
0489             CH1= SFMIX(IFL,3)*SFMIX(IFL-1,1)
0490             CH2=-SFMIX(IFL,4)*SFMIX(IFL-1,2)
0491             CH3= SFMIX(IFL,3)*SFMIX(IFL-1,2)
0492             CH4= SFMIX(IFL,4)*SFMIX(IFL-1,1)
0493           ENDIF
0494           IF(IFL.EQ.6) THEN
0495             XMF=XMTOP
0496             XMFP=XMBOT
0497             AT=ATRIT
0498             AB=ATRIB
0499           ENDIF
0500         ELSE
0501 C...B1 -> T1 HC
0502           IF(ILR.EQ.1) THEN
0503             CH1=-SFMIX(IFL+1,1)*SFMIX(IFL,1)
0504             CH2= SFMIX(IFL+1,2)*SFMIX(IFL,2)
0505             CH3=-SFMIX(IFL+1,1)*SFMIX(IFL,2)
0506             CH4=-SFMIX(IFL+1,2)*SFMIX(IFL,1)
0507 C...B2-> T1 HC
0508           ELSE
0509             CH1= SFMIX(IFL,3)*SFMIX(IFL+1,1)
0510             CH2=-SFMIX(IFL,4)*SFMIX(IFL+1,2)
0511             CH3= SFMIX(IFL,4)*SFMIX(IFL+1,1)
0512             CH4= SFMIX(IFL,3)*SFMIX(IFL+1,2)
0513           ENDIF
0514           IF(IFL.EQ.5) THEN
0515             XMF=XMTOP
0516             XMFP=XMBOT
0517             AT=ATRIT
0518             AB=ATRIB
0519           ENDIF
0520         ENDIF
0521         XL=PYLAMF(XMI2,XMSF1**2,XMB**2)
0522         LKNT=LKNT+1
0523 C.......Need to complexify
0524         AL=CH1*(XMW2*2D0*CBETA*SBETA-XMFP**2*TANB-XMF**2/TANB)+
0525      &  CH2*2D0*XMF*XMFP/(2D0*CBETA*SBETA)+
0526      &  CH3*XMFP*(-XMU+AB*TANB)+CH4*XMF*(-XMU+AT/TANB)
0527         XLAM(LKNT)=C1/8D0/XMI3*SQRT(XL)/XMW2*AL**2
0528         IDLAM(LKNT,3)=0
0529         IDLAM(LKNT,1)=KF1
0530         IDLAM(LKNT,2)=SIGN(37,KCHG(IFL,1))
0531       ENDIF
0532       IF(XMI.GT.XMB+XMSF2) THEN
0533         XMF=0D0
0534         XMFP=0D0
0535         AT=0D0
0536         AB=0D0
0537         IF(MOD(IFL,2).EQ.0) THEN
0538 C...T1-> B2 HC
0539           IF(ILR.EQ.1) THEN
0540             CH1= SFMIX(IFL-1,3)*SFMIX(IFL,1)
0541             CH2=-SFMIX(IFL-1,4)*SFMIX(IFL,2)
0542             CH3= SFMIX(IFL-1,4)*SFMIX(IFL,1)
0543             CH4= SFMIX(IFL-1,3)*SFMIX(IFL,2)
0544 C...T2-> B2 HC
0545           ELSE
0546             CH1= -SFMIX(IFL,3)*SFMIX(IFL-1,3)
0547             CH2= SFMIX(IFL,4)*SFMIX(IFL-1,4)
0548             CH3= -SFMIX(IFL,3)*SFMIX(IFL-1,4)
0549             CH4= -SFMIX(IFL,4)*SFMIX(IFL-1,3)
0550           ENDIF
0551           IF(IFL.EQ.6) THEN
0552             XMF=XMTOP
0553             XMFP=XMBOT
0554             AT=ATRIT
0555             AB=ATRIB
0556           ENDIF
0557         ELSE
0558 C...B1 -> T2 HC
0559           IF(ILR.EQ.1) THEN
0560             CH1= SFMIX(IFL+1,3)*SFMIX(IFL,1)
0561             CH2=-SFMIX(IFL+1,4)*SFMIX(IFL,2)
0562             CH3= SFMIX(IFL+1,3)*SFMIX(IFL,2)
0563             CH4= SFMIX(IFL+1,4)*SFMIX(IFL,1)
0564 C...B2-> T2 HC
0565           ELSE
0566             CH1= -SFMIX(IFL+1,3)*SFMIX(IFL,3)
0567             CH2= SFMIX(IFL+1,4)*SFMIX(IFL,4)
0568             CH3= -SFMIX(IFL+1,3)*SFMIX(IFL,4)
0569             CH4= -SFMIX(IFL+1,4)*SFMIX(IFL,3)
0570           ENDIF
0571           IF(IFL.EQ.5) THEN
0572             XMF=XMTOP
0573             XMFP=XMBOT
0574             AT=ATRIT
0575             AB=ATRIB
0576           ENDIF
0577         ENDIF
0578         XL=PYLAMF(XMI2,XMSF1**2,XMB**2)
0579         LKNT=LKNT+1
0580 C.......Need to complexify
0581         AL=CH1*(XMW2*2D0*CBETA*SBETA-XMFP**2*TANB-XMF**2/TANB)+
0582      &  CH2*2D0*XMF*XMFP/(2D0*CBETA*SBETA)+
0583      &  CH3*XMFP*(-XMU+AB*TANB)+CH4*XMF*(-XMU+AT/TANB)
0584         XLAM(LKNT)=C1/8D0/XMI3*SQRT(XL)/XMW2*AL**2
0585         IDLAM(LKNT,3)=0
0586         IDLAM(LKNT,1)=KF2
0587         IDLAM(LKNT,2)=SIGN(37,KCHG(IFL,1))
0588       ENDIF
0589
0590 C...2-BODY DECAYS OF SQUARK -> QUARK GLUINO
0591
0592       IF(IFL.LE.6) THEN
0593         XMFP=0D0
0594         XMF=0D0
0595         IF(IFL.EQ.6) XMF=PMAS(6,1)
0596         IF(IFL.EQ.5) XMF=PMAS(5,1)
0597         XMJ=PMAS(PYCOMP(KSUSY1+21),1)
0598         AXMJ=ABS(XMJ)
0599         IF(XMI.GE.AXMJ+XMF) THEN
0600           AL=-SFMIX(IFL,3)
0601           BL=SFMIX(IFL,1)
0602           AR=-SFMIX(IFL,4)
0603           BR=SFMIX(IFL,2)
0604 C...F1 -> F CHI
0605           IF(ILR.EQ.1) THEN
0606             XCA=AL
0607             XCB=BL
0608 C...F2 -> F CHI
0609           ELSE
0610             XCA=AR
0611             XCB=BR
0612           ENDIF
0613           LKNT=LKNT+1
0614           XMA2=XMJ**2
0615           XMB2=XMF**2
0616           XL=PYLAMF(XMI2,XMA2,XMB2)
0617           XLAM(LKNT)=4D0/3D0*AS/2D0/XMI3*SQRT(XL)*((XMI2-XMB2-XMA2)*
0618      &    (XCA**2+XCB**2)+4D0*XCA*XCB*XMJ*XMF)
0619           IDLAM(LKNT,1)=KSUSY1+21
0620           IDLAM(LKNT,2)=IFL
0621           IDLAM(LKNT,3)=0
0622         ENDIF
0623       ENDIF
0624
0625 C...IF NOTHING ELSE FOR T1, THEN T1* -> C+CHI0
0626       IF(KFIN.EQ.KSUSY1+6.AND.PMAS(KCIN,1).GT.
0627      &PMAS(PYCOMP(KSUSY1+22),1)+PMAS(4,1)) THEN
0628 C...THIS IS A BACK-OF-THE-ENVELOPE ESTIMATE
0629 C...M = 1/(16PI**2)G**3 = G*2/(4PI) G/(4PI) = C1 * G/(4PI)
0630 C...M*M = C1**2 * G**2/(16PI**2)
0631 C...G = 1/(8PI)P/MI**2 * M*M = C1**3/(32PI**2)*LAM/(2*MI**3)
0632         LKNT=LKNT+1
0633         XL=PYLAMF(XMI2,0D0,PMAS(PYCOMP(KSUSY1+22),1)**2)
0634         XLAM(LKNT)=C1**3/64D0/PI**2/XMI3*SQRT(XL)
0635         IF(XLAM(LKNT).EQ.0) XLAM(LKNT)=1D-3
0636         IDLAM(LKNT,1)=KSUSY1+22
0637         IDLAM(LKNT,2)=4
0638         IDLAM(LKNT,3)=0
0639       ENDIF
0640
0641 C...R-violating sfermion decays (SKANDS).
0642       CALL PYRVSF(KFIN,XLAM,IDLAM,LKNT)
0643
0644       IKNT=LKNT
0645       XLAM(0)=0D0
0646       DO 170 I=1,IKNT
0647         IF(XLAM(I).LT.0D0) XLAM(I)=0D0
0648         XLAM(0)=XLAM(0)+XLAM(I)
0649   170 CONTINUE
0650       IF(XLAM(0).EQ.0D0) XLAM(0)=1D-3
0651
0652       RETURN
0653       END