IAP GITLAB

Commit abca1e00 authored by Colin Baus's avatar Colin Baus

working version with dpmjet at least for pp

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@3901 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 54a91227
......@@ -147,7 +147,7 @@ c Fix final particles and some event parameters
call afinal
c Fill HEP common
if(model.ne.4.and.model.ne.5)call hepmcstore
if(model.ne.4.and.model.ne.5.and.model.ne.12)call hepmcstore
c optional Statistic information (only with debug level ish=1)
call astati
......
......@@ -25,8 +25,7 @@ C general initialization
IDP = 1
IGLAU = 0
EPROJ = (ECMS**2-AAM(NPMASS)**2-AAM(NTMASS)**2)
& +(2.0D0*AAM(NTMASS))
print *,' !!!! BLA !!!! ',AAM(NPMASS),AAM(NTMASS),ECMS,EPROJ
& /(2.0D0*AAM(NTMASS))
CALL DT_INIT(NCASES,EPROJ,NPMASS,NPCHAR,NTMASS,NTCHAR,IDP,IGLAU)
END
......@@ -69,7 +68,7 @@ C general initialization
NTCHAR = latarg
IDP = 1
EPROJ = (ECMS**2-AAM(NPMASS)**2-AAM(NTMASS)**2)
& +(2.0D0*AAM(NTMASS))
& /(2.0D0*AAM(NTMASS))
ITRY = 0
......@@ -108,8 +107,8 @@ c PLEASE FILL!!
C particle copy loop
DO I=1,NHKK
print *,'asd ',I,NHKK, phkk(4:5,I),ISTHKK(I)
IF(ISTHKK(I).NE.1 .AND. ISTHKK(I).NE.4) CYCLE
c Destroys HEPEVT structure if only beam and final (status={1,4}) particles are copied
c IF(ISTHKK(I).NE.1 .AND. ISTHKK(I).NE.4) CYCLE
NHEP = NHEP + 1
ISTHEP(I) = ISTHKK(I)
IDHEP(I) = IDHKK(I)
......@@ -128,15 +127,6 @@ C particle copy loop
VHEP(4,I) = VHKK(4,I)
ENDDO
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)
c nhkk = number of particles in event
c nevhkk = number of event
iret=0
1000 return
......
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