IAP GITLAB

Commit b0abe2b0 authored by Anatoli Fedynitch's avatar Anatoli Fedynitch

DPMJET should fill the HEPEVT common block after these changes.

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@3896 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent e2b7dc32
......@@ -24,8 +24,8 @@ C general initialization
NTCHAR = latarg
IDP = 1
IGLAU = 0
EPROJ = (ECMS**2-AAM(NPMASS)**2-AAM(NTMASS)**2)/
+(2.0D0*AAM(NTMASS))
EPROJ = (ECMS**2-AAM(NPMASS)**2-AAM(NTMASS)**2)
& +(2.0D0*AAM(NTMASS))
print *,' !!!! BLA !!!! ',AAM(NPMASS),AAM(NTMASS),ECMS,EPROJ
CALL DT_INIT(NCASES,EPROJ,NPMASS,NPCHAR,NTMASS,NTCHAR,IDP,IGLAU)
......@@ -68,8 +68,8 @@ C general initialization
NTMASS = matarg
NTCHAR = latarg
IDP = 1
EPROJ = (ECMS**2-AAM(NPMASS)**2-AAM(NTMASS)**2)/
+(2.0D0*AAM(NTMASS))
EPROJ = (ECMS**2-AAM(NPMASS)**2-AAM(NTMASS)**2)
& +(2.0D0*AAM(NTMASS))
ITRY = 0
......@@ -85,7 +85,7 @@ cGOTO 1001
c PLEASE FILL!!
ncol=1
nevt=1
nevt=IREJ
kolevt=ncol
npjevt=maproj
ntgevt=matarg
......@@ -96,17 +96,18 @@ c PLEASE FILL!!
phievt=0.
phi=0.
nhep = nhkk
C particle copy loop
DO I=1,NHKK
IF(ISTHKK(I).EQ.1) THEN
C write(6,*) nevhkk, idhkk(I),nobam(I), idbam(I) ,idch(I)
PX = phkk(1,I)
PY = phkk(2,I)
PZ = phkk(3,I)
EE = phkk(4,I)
AM = phkk(5,I)
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
ENDDO
......@@ -116,8 +117,9 @@ C write(6,*) nevhkk, idhkk(I),nobam(I), idbam(I) ,idch(I)
& ,phkk(4,1),phkk(5,1)
print *,itry,' !!!! NOT FINISHED !!!!!!!!!!!!!!!!!!!!!!!!!! !!!! '
$, 'dtevt2 needs to be copied to hepvt here'
print *,itry,' !!!! NOT FINISHED !!!!!!!!!!!!!!!!!!!!!!!!!! !!!!',
& 'dtevt2 needs to be copied to hepvt here'
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