Boost also resonances into LAB frame
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 16 May 2013 16:06:26 +0000 (16:06 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 16 May 2013 16:06:26 +0000 (16:06 +0000)
DPMJET/dpmjet3.0-5.f
DPMJET/dpmjet3.0-5F.f

index 303dc02..f876d4c 100644 (file)
@@ -3359,8 +3359,11 @@ C     STOP
       IF (LBEAM) THEN
          IF ( (NPOINT(4).EQ.0).OR.(NHKK.LT.NPOINT(4)) ) RETURN
          DO 20 I=NPOINT(4),NHKK
-            IF ((ABS(ISTHKK(I)).EQ.1).OR.(ISTHKK(I).EQ.1000).OR.
-     &                                   (ISTHKK(I).EQ.1001)) THEN
+            IF ((ABS(ISTHKK(I)).EQ.1)  .OR.
+     &           (ABS(ISTHKK(I)).EQ.2) .OR.
+     &           (ISTHKK(I).EQ.1000)   .OR.
+     &           (ISTHKK(I).EQ.1001)) THEN
+               
                CALL DT_MYTRAN(1,PHKK(1,I),PHKK(2,I),PHKK(3,I),
      &                     COD,SID,COF,SIF,PXCMS,PYCMS,PZCMS)
                PECMS = PHKK(4,I)
@@ -26780,6 +26783,7 @@ C     DO 1 I=1,NEND
       DO 1 I=NPOINT(4),NHKK
          IF ((ABS(ISTHKK(I)).EQ.1).OR.(ISTHKK(I).EQ.1000).OR.
      &                                (ISTHKK(I).EQ.1001)) THEN
+            
             CALL DT_LTNUC(PHKK(3,I),PHKK(4,I),PZ,PE,-3)
             PHKK(3,I) = PZ
             PHKK(4,I) = PE
index 0386485..9061bf8 100644 (file)
@@ -3592,8 +3592,12 @@ C     STOP
       IF (LBEAM) THEN
          IF ( (NPOINT(4).EQ.0).OR.(NHKK.LT.NPOINT(4)) ) RETURN
          DO 20 I=NPOINT(4),NHKK
-            IF ((ABS(ISTHKK(I)).EQ.1).OR.(ISTHKK(I).EQ.1000).OR.
-     &                                   (ISTHKK(I).EQ.1001)) THEN
+
+            IF ((ABS(ISTHKK(I)).EQ.1)  .OR.
+     &           (ABS(ISTHKK(I)).EQ.2) .OR.
+     &           (ISTHKK(I).EQ.1000)   .OR.
+     &           (ISTHKK(I).EQ.1001)) THEN
+
                CALL DT_MYTRAN(1,PHKK(1,I),PHKK(2,I),PHKK(3,I),
      &                     COD,SID,COF,SIF,PXCMS,PYCMS,PZCMS)
                PECMS = PHKK(4,I)