IAP GITLAB

Commit e115d43f authored by Colin Baus's avatar Colin Baus

dpmejet: some fixes concerning beam particles and vertex of children

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4132 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 699bc52f
......@@ -93,7 +93,9 @@ c-----------------------------------------------------------------------
COMMON /DPMEVTINI/ EPROJ,NPMASS,NPCHAR,NTMASS,NTCHAR,IDP
integer nptlhep(NMXHKK)
call conre !! adds beam particles that are knonw from intial call of crmc
anintine=anintine+1.
call conre !! adds beam particles that are known from intial call of crmc
call conwr
IREJ = 0
......@@ -155,10 +157,6 @@ c PLEASE FILL!!
GOTO 1001
ENDIF
c set projectile and target as non final particles
istptl(1)=1
istptl(2)=1
do k=1,NHKK
c LIST is the code of final particle, P - its 4-momentum and mass.
......@@ -172,7 +170,7 @@ c LIST is the code of final particle, P - its 4-momentum and mass.
c final particle
nptl=nptl+1 !! add 1 particle to stack
nptlhep(k)=nptl
istptl(nptl)=ISTHKK(k)-1 !!somehow 0 is final?
istptl(nptl)=ISTHKK(k)-1 !!0 means last generation other codes are e.g. for pomerons, remnants...
else
nptlhep(k)=0
cycle
......@@ -189,10 +187,10 @@ c final particle
pptl(4,nptl)=sngl(PHKK(4,k)) !E
pptl(5,nptl)=sngl(PHKK(5,k)) !mass
ityptl(nptl)=0
iorptl(nptl)=0
if(JMOHKK(1,k).ne.0)iorptl(nptl)=nptlhep(JMOHKK(1,k))
jorptl(nptl)=0
if(JMOHKK(2,k).ne.0)jorptl(nptl)=nptlhep(JMOHKK(2,k))
iorptl(nptl)=1
c if(JMOHKK(1,k).ne.0)iorptl(nptl)=nptlhep(JMOHKK(1,k))
jorptl(nptl)=maproj+matarg
c if(JMOHKK(2,k).ne.0)jorptl(nptl)=nptlhep(JMOHKK(2,k))
ifrptl(1,nptl)=0
ifrptl(2,nptl)=0
xorptl(1,nptl)=sngl(VHKK(1,k))*1e12
......
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