IAP GITLAB

Commit d7b09ba8 authored by Colin Baus's avatar Colin Baus

copy full hepevt block for dpmjet

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@3899 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent b0abe2b0
......@@ -74,18 +74,19 @@ C general initialization
ITRY = 0
IREJ = 0
IRET=1
KKMAT = -1
NEVHKK = K
ITRY = ITRY+1
CALL dt_kkinc(NPMASS,NPCHAR,NTMASS,NTCHAR,IDP,EPROJ, KKMAT,IREJ)
IF(IREJ.NE.0) print *,"uff"
cGOTO 1001
IF(IREJ.NE.0) THEN
IF(ish.ge.2) WRITE(ifch,'(a)')
&'Bad return from DPMJET. Skipping event'
GOTO 1001
ENDIF
c PLEASE FILL!!
ncol=1
nevt=IREJ
nevt=1
kolevt=ncol
npjevt=maproj
ntgevt=matarg
......@@ -96,30 +97,45 @@ c PLEASE FILL!!
phievt=0.
phi=0.
IF(NHKK.GE.NMXHEP) THEN
IF(ish.ge.2) WRITE(ifch,'(a)')
&'NMAXHEP too small for event. Skipping'
GOTO 1001
ENDIF
NEVHEP = NEVHKK
NHEP = 0
nhep = nhkk
C particle copy loop
DO I=1,NHKK
isthep(I) = isthkk(I)
idhep(I) = idhkk(I)
phep(1,I) = phkk(1,I)
phep(2,I) = phkk(2,I)
phep(3,I) = phkk(3,I)
phep(4,I) = phkk(4,I)
phep(5,I) = phkk(5,I)
ENDIF
print *,'asd ',I,NHKK, phkk(4:5,I),ISTHKK(I)
IF(ISTHKK(I).NE.1 .AND. ISTHKK(I).NE.4) CYCLE
NHEP = NHEP + 1
ISTHEP(I) = ISTHKK(I)
IDHEP(I) = IDHKK(I)
JMOHEP(1,I) = JMOHKK(1,I)
JMOHEP(2,I) = JMOHKK(2,I)
JDAHEP(1,I) = JDAHKK(1,I)
JDAHEP(2,I) = JDAHKK(2,I)
PHEP(1,I) = PHKK(1,I)
PHEP(2,I) = PHKK(2,I)
PHEP(3,I) = PHKK(3,I)
PHEP(4,I) = PHKK(4,I)
PHEP(5,I) = PHKK(5,I)
VHEP(1,I) = VHKK(1,I)
VHEP(2,I) = VHKK(2,I)
VHEP(3,I) = VHKK(3,I)
VHEP(4,I) = VHKK(4,I)
ENDDO
print *,' !!!! HERE !!!! ',nevhep,phep(1,1),phep(2,1),phep(3,1)
& ,phep(4,1),phep(5,1)
print *,' !!!! HERE !!!! ',nevhep,phkk(1,1),phkk(2,1),phkk(3,1)
& ,phkk(4,1),phkk(5,1)
print *,' !!!! HERE !!!! ',nhep,nevhep,phep(4,i),phep(5,i),
&vhep(1,i),jmohep(1,i),jmohep(1,i)
print *,' !!!! HERE !!!! ',nhkk,nevhkk,phkk(1,1),phkk(2,1),
&phkk(3,1),phkk(4,1),phkk(5,1)
print *,itry,' !!!! NOT FINISHED !!!!!!!!!!!!!!!!!!!!!!!!!! !!!!',
& 'dtevt2 needs to be copied to hepvt here'
c nhkk = number of particles in event
c nevhkk = number of event
iret=0
1000 return
......
......@@ -1941,7 +1941,7 @@ C...Commonblocks.
COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
SAVE /PYJETS/,/PYDAT1/,/PYDAT2/
C...HEPEVT commonblock.
PARAMETER (NMXHEP=4000)
PARAMETER (NMXHEP=9990)
COMMON/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
&JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
DOUBLE PRECISION PHEP,VHEP
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