IAP GITLAB

Commit 673a015c authored by Tanguy Pierog's avatar Tanguy Pierog

correct istmax definition in dpmjet-crmc.f and try to copy properly mother...

correct istmax definition in dpmjet-crmc.f and try to copy properly mother particles if EPOS is not used.


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4165 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 8952c0ad
......@@ -146,7 +146,8 @@ c will be treated as nucleus in DT_INIT and NPMASS,... will be used
c initialize cross-sections by calling epos-bas.f function -> models.F -> dpmjet-crmc.f
call xsigma
istmax=0 !don't copy decayed particles atm
c istmax=0 !uncomment only if you want to be sure to have only final particles
END
c-----------------------------------------------------------------------
......@@ -262,7 +263,7 @@ c LIST is the code of final particle, P - its 4-momentum and mass.
$ ' DPMJET particle ',k,' id :',ic,' before conversion'
$ , ' momentum :',(sngl(PHKK(i,k)),i=1,5)
imaxhepstatus = MAX(1,MIN(ISTMAX,2))
imaxhepstatus = MAX(1,MIN(ISTMAX+1,2))
IF(ISTHKK(k).GE.1 .AND. ISTHKK(k).LE.imaxhepstatus)THEN !! if final particle
nptl=nptl+1 !! add 1 particle to stack
nptlhep(k)=nptl
......
......@@ -1375,7 +1375,7 @@ c copy first list in final list to define daughters of beam particles
nhep=0
nhepio=0
if(istmaxhep.ne.0)nhep=maproj+matarg
if(istmaxhep.ne.0.and.model.eq.1)nhep=maproj+matarg
lrcor0=.true. !link spectator remnants to core only once
lrcore=.false.
......@@ -1383,7 +1383,7 @@ c start with beam particles (except spectators producing fragments)
do j=1,maproj+matarg
if(istmaxhep.eq.0)then
if(istmaxhep.eq.0.or.model.ne.1)then
c when no daughter/mother informations, simply copy beam particles
nhep=nhep+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