IAP GITLAB

Commit 4e1f1419 authored by Colin Baus's avatar Colin Baus

copy decayed particles disabled

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4138 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 3eb5b9c3
......@@ -146,6 +146,7 @@ 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
END
c-----------------------------------------------------------------------
......@@ -202,7 +203,7 @@ C call dpmjet event simulation
GOTO 1001
ENDIF
call dt_evtout(1)
c call dt_evtout(1)
c nevt .......... error code. 1=valid event, 0=invalid event
c bimevt ........ absolute value of impact parameter
......@@ -261,7 +262,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 = MIN(ISTMAX,2)
imaxhepstatus = MAX(1,MIN(ISTMAX,2))
print*,imaxhepstatus
IF(ISTHKK(k).GE.1 .AND. ISTHKK(k).LE.imaxhepstatus)THEN !! if final particle
nptl=nptl+1 !! add 1 particle to stack
......@@ -298,12 +299,12 @@ c treatment of mothers and daughters
iorptl(nptl)=1 !by default set to first and last beam particle
jorptl(nptl)=maproj+matarg
c if(JMOHKK(1,k).ne.0.AND.nptlhep(JMOHKK(1,k)).NE.0) then !mother
c iorptl(nptl)=nptlhep(JMOHKK(1,k))
c endif
c if(JMOHKK(2,k).ne.0.AND.nptlhep(JMOHKK(2,k)).NE.0) then !father
c jorptl(nptl)=nptlhep(JMOHKK(2,k))
c endif
if(JMOHKK(1,k).ne.0.AND.nptlhep(JMOHKK(1,k)).NE.0) then !mother
iorptl(nptl)=nptlhep(JMOHKK(1,k))
endif
if(JMOHKK(2,k).ne.0.AND.nptlhep(JMOHKK(2,k)).NE.0) then !father
jorptl(nptl)=nptlhep(JMOHKK(2,k))
endif
print*,istmax,k,JMOHKK(1,k),nptlhep(k),iorptl(nptl),
$ nptlhep(JMOHKK(1,k))
......
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