IAP GITLAB

Commit fecfc4c6 authored by Tanguy Pierog's avatar Tanguy Pierog

correct bug in momentum beam definition for DPMJET (just for the output)


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4192 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 27994e23
......@@ -76,7 +76,7 @@ c call DPMJET to simulate interaction
c-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION(A-H,O-Z)
include 'epos.inc'
real yprj,ytrg
* event history
PARAMETER (NMXHKK=200000)
COMMON /DTEVT1/ NHKK,NEVHKK,ISTHKK(NMXHKK),IDHKK(NMXHKK),
......@@ -220,7 +220,9 @@ c treatment of mothers and daughters
if(ist.ge.11.and.ist.le.18)then
c fix beam particles which are registered at rest in DPMJET
if(mod(ist,2).eq.0)then !target remnant (ist=12,14,16,18)
call utlob5(yhaha,pptl(1,nptl)
ytrg=log((sqrt(pnullx**2+pptl(5,nptl)**2)+pnullx)
. /pptl(5,nptl))
call utlob5(ytrg,pptl(1,nptl)
. ,pptl(2,nptl),pptl(3,nptl),pptl(4,nptl),pptl(5,nptl)) !boost in cms frame
if(ist.eq.12)then !inelastic nucleon
istptl(nptl)=1
......@@ -241,7 +243,9 @@ c fix beam particles which are registered at rest in DPMJET
jorptl(nptl)=0 !defined later
endif
else !projectile remnant (ist=11,13,15,17)
call utlob5(-yhaha,pptl(1,nptl)
yprj=log((sqrt(pnullx**2+pptl(5,nptl)**2)+pnullx)
. /pptl(5,nptl))
call utlob5(-yprj,pptl(1,nptl)
. ,pptl(2,nptl),pptl(3,nptl),pptl(4,nptl),pptl(5,nptl)) !boost in cms frame
if(ist.eq.11)then !inelastic nucleon
istptl(nptl)=1
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment