WRITE(MSTU(11),5100)
STOP
ENDIF
-
C...Initial values for some counters.
N=0
MINT(5)=MINT(5)+1
C for the Lambda(1520).
C -------------------------------------------------------------------------
- IF (IABS(KF).EQ.3122) THEN
+C IF (IABS(KF).EQ.3122) THEN
C Converting a fraction (0.20) of Lambda0 to Lambda(1520) + c.c.
C This fraction is based on the experimental measurement at ISR
C Bobbink 83, NP B217,11 (1983)
C The region 0.5 < XF < 1.0 has been extrapolated to XF=0
- IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
- ENDIF
+C IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
+C ENDIF
- IF(IABS(KF).EQ.3212) THEN
+C IF(IABS(KF).EQ.3212) THEN
C Converting a fraction (0.20) of Sigma0 to Lambda(1520) + c.c.
C We suppose the same fraction as for Lambda0
- IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
- ENDIF
+C IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
+C ENDIF
- IF (IABS(KF).EQ.3214) THEN
+C IF (IABS(KF).EQ.3214) THEN
C Converting a fraction (0.30) of Sigma0(1385) to Lambda(1520) + c.c.
C This is conservative extimate supposing that the ratio
C scales as (M_Sigma1385/M_Lambda0)^2 ~ 1.5
- IF(PYR(0).LE.0.30) KF=ISIGN(3124,KF)
- ENDIF
+C IF(PYR(0).LE.0.30) KF=ISIGN(3124,KF)
+C ENDIF
RETURN
C...Use tabulated probabilities to select new flavour and hadron.
C...Integral of Altarelli-Parisi z kernel for QCD.
C...(Includes squark and gluino; with factor N_C/C_F extra for latter).
+
FMED = PARJ(200)
IF(MSTJ(49).EQ.0.AND.KFL(1).EQ.21) THEN
C Nestor