author hristov Tue, 1 Mar 2011 08:46:26 +0000 (08:46 +0000) committer hristov Tue, 1 Mar 2011 08:46:26 +0000 (08:46 +0000)
 PYTHIA6/pythia6214.f patch | blob | blame | history

index 6a9b73c..6deccbf 100644 (file)
@@ -2856,6 +2856,7 @@ C...Stop if no subprocesses on.
WRITE(MSTU(11),5100)
STOP
ENDIF
C...Initial values for some counters.
N=0
MINT(5)=MINT(5)+1
C 2) Use the PYUPDA routine (see manual) to include particle and decay data
C for the Lambda(1520).
C -------------------------------------------------------------------------

-C      IF (IABS(KF).EQ.3122) THEN
+      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
-C         IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
-C      ENDIF
+         IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
+      ENDIF

-C      IF(IABS(KF).EQ.3212) THEN
+      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
-C         IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
-C      ENDIF
+         IF(PYR(0).LE.0.20) KF=ISIGN(3124,KF)
+      ENDIF

-C      IF (IABS(KF).EQ.3214) THEN
+      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
-C         IF(PYR(0).LE.0.30) KF=ISIGN(3124,KF)
-C      ENDIF
+         IF(PYR(0).LE.0.30) KF=ISIGN(3124,KF)
+      ENDIF
RETURN

C...Use tabulated probabilities to select new flavour and hadron.
@@ -51413,7 +51414,6 @@ C...Calculate allowed z range.

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