Adding macro to run the test pt task also on offline ESDs, plus corresponding documen...

[u/mrichter/AliRoot.git] / PYTHIA6 / pythia6214.f

diff --git a/PYTHIA6/pythia6214.f b/PYTHIA6/pythia6214.f
index 6a9b73c5c894e70ae604b899039d835d5287caa5..28e5d1e04ea89b3b4e45506ed6d72319f1c60b21 100644 (file)
--- a/PYTHIA6/pythia6214.f
+++ b/PYTHIA6/pythia6214.f
@@ -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
@@ -50112,26 +50113,26 @@ 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
@@ -56238,6 +56238,7 @@ C...Loop over all particles. Find cell that was hit by given particle.
      &    KC.EQ.18) GOTO 110
           IF(MSTU(41).GE.3.AND.KCHG(KC,2).EQ.0.AND.PYCHGE(K(I,2)).EQ.0)
      &    GOTO 110
+
         ENDIF
         NP=NP+1
         PT=SQRT(P(I,1)**2+P(I,2)**2)