From 430525dd1ef3a15848f073518387accfbbeeab8d Mon Sep 17 00:00:00 2001 From: morsch Date: Thu, 16 May 2013 16:06:26 +0000 Subject: [PATCH] Boost also resonances into LAB frame --- DPMJET/dpmjet3.0-5.f | 8 ++++++-- DPMJET/dpmjet3.0-5F.f | 8 ++++++-- 2 files changed, 12 insertions(+), 4 deletions(-) diff --git a/DPMJET/dpmjet3.0-5.f b/DPMJET/dpmjet3.0-5.f index 303dc02cf32..f876d4c14fd 100644 --- a/DPMJET/dpmjet3.0-5.f +++ b/DPMJET/dpmjet3.0-5.f @@ -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 diff --git a/DPMJET/dpmjet3.0-5F.f b/DPMJET/dpmjet3.0-5F.f index 03864854124..9061bf8c937 100644 --- a/DPMJET/dpmjet3.0-5F.f +++ b/DPMJET/dpmjet3.0-5F.f @@ -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) -- 2.43.0