Small mod to run with boundary checking

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

diff --git a/PYTHIA6/pythia6214.f b/PYTHIA6/pythia6214.f
index 6f78fec60c00d4154c5845300e74acf82e22f16c..02a2088d5fcd9f703eee440c1e0109c5fd6e1755 100644 (file)
--- a/PYTHIA6/pythia6214.f
+++ b/PYTHIA6/pythia6214.f
@@ -54963,24 +54963,26 @@ C...Identify particle code and whether already defined  (for MUPDA=3).
   150         CONTINUE
             ENDIF
 C...Remove duplicate old decay data.
-            IF(KCREP.NE.0.AND.MDCY(KCREP,3).GT.0) THEN
-              IDCREP=MDCY(KCREP,2)
-              NDCREP=MDCY(KCREP,3)
-              DO 160 I=1,KCC
-                IF(MDCY(I,2).GT.IDCREP) MDCY(I,2)=MDCY(I,2)-NDCREP
-  160         CONTINUE
-              DO 180 I=IDCREP,NDC-NDCREP
-                MDME(I,1)=MDME(I+NDCREP,1)
-                MDME(I,2)=MDME(I+NDCREP,2)
-                BRAT(I)=BRAT(I+NDCREP)
-                DO 170 J=1,5
-                  KFDP(I,J)=KFDP(I+NDCREP,J)
-  170           CONTINUE
-  180         CONTINUE
-              NDC=NDC-NDCREP
-              KC=KCREP
-            ELSEIF(KCREP.NE.0) THEN
-              KC=KCREP
+            IF(KCREP.NE.0) THEN
+               IF(MDCY(KCREP,3).GT.0) THEN
+                  IDCREP=MDCY(KCREP,2)
+                  NDCREP=MDCY(KCREP,3)
+                  DO 160 I=1,KCC
+                     IF(MDCY(I,2).GT.IDCREP) MDCY(I,2)=MDCY(I,2)-NDCREP
+ 160              CONTINUE
+                  DO 180 I=IDCREP,NDC-NDCREP
+                     MDME(I,1)=MDME(I+NDCREP,1)
+                     MDME(I,2)=MDME(I+NDCREP,2)
+                     BRAT(I)=BRAT(I+NDCREP)
+                     DO 170 J=1,5
+                        KFDP(I,J)=KFDP(I+NDCREP,J)
+ 170                 CONTINUE
+ 180              CONTINUE
+                  NDC=NDC-NDCREP
+                  KC=KCREP
+               ELSE
+                  KC=KCREP
+               ENDIF
             ELSE
               KCC=KCC+1
               KC=KCC