IAP GITLAB

Commit 15dd6d74 authored by Colin Baus's avatar Colin Baus

first working dpmjet version for pp

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4047 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent c4a75707
c 13.04.2010 Link routines between DPMJET and CRMC.
c authors C. Baus, A. Feydnitch
c-----------------------------------------------------------------------
subroutine IniDPMJET
c-----------------------------------------------------------------------
......@@ -16,6 +16,14 @@ c-----------------------------------------------------------------------
INTEGER NCASES,NPMASS,NPCHAR,NTMASS,NTCHAR,IDP,IGLAU
DOUBLE PRECISION EPROJ
c change to center of mass system
c common block var iframe renamed to dframe since it is used by epos
LOGICAL LEMCCK,LHADRO,LSEADI,LEVAPO
COMMON /DTFLG1/ IFRAG(2),IRESCO,IMSHL,IRESRJ,IOULEV(6),
& LEMCCK,LHADRO(0:9),LSEADI,LEVAPO,DFRAME,ITRSPT
DFRAME=2
C general initialization
NCASES = -1
NPMASS = maproj
......@@ -83,6 +91,35 @@ C general initialization
GOTO 1001
ENDIF
c nevt .......... error code. 1=valid event, 0=invalid event
c bimevt ........ absolute value of impact parameter
c phievt ........ angle of impact parameter
c kolevt ........ number of collisions
c koievt ........ number of inelastic collisions
c kohevt ........ number of hard collisions
c pmxevt ........ reference momentum
c egyevt ........ pp cm energy (hadron) or string energy (lepton)
c npjevt ........ number of primary projectile participants
c ntgevt ........ number of primary target participants
c npnevt ........ number of primary projectile neutron spectators
c nppevt ........ number of primary projectile proton spectators
c ntnevt ........ number of primary target neutron spectators
c ntpevt ........ number of primary target proton spectators
c jpnevt ........ number of absolute projectile neutron spectators
c jppevt ........ number of absolute projectile proton spectators
c jtnevt ........ number of absolute target neutron spectators
c jtpevt ........ number of absolute target proton spectators
c xbjevt ........ bjorken x for dis
c qsqevt ........ q**2 for dis
c sigtot ........ total cross section
c nglevt ........ number of collisions acc to Glauber
c zppevt ........ average Z-parton-proj
c zptevt ........ average Z-parton-targ
c ng1evt ........ number of Glauber participants with at least one IAs
c ng2evt ........ number of Glauber participants with at least two IAs
c ikoevt ........ number of elementary parton-parton scatterings
c typevt ........ type of event (1=Non Diff, 2=Double Diff, 3=Central Diff, 4=AB->XB, -4=AB->AX)
c PLEASE FILL!!
ncol=1
nevt=1
......@@ -106,26 +143,36 @@ c PLEASE FILL!!
NHEP = 0
C particle copy loop
DO I=1,NHKK
DO I=3,NHKK
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)
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)
ISTHEP(I-2) = ISTHKK(I)
IDHEP(I-2) = IDHKK(I)
JMOHEP(1,I-2) = JMOHKK(1,I)-2
JMOHEP(2,I-2) = JMOHKK(2,I)-2
JDAHEP(1,I-2) = JDAHKK(1,I)-2
JDAHEP(2,I-2) = JDAHKK(2,I)-2
IF(JMOHKK(1,I).EQ.0) JMOHEP(1,I-2) = 1
IF(JMOHKK(2,I).EQ.0) JMOHEP(2,I-2) = 1
IF(JDAHKK(1,I).EQ.0) JDAHEP(1,I-2) = 0
IF(JDAHKK(2,I).EQ.0) JDAHEP(2,I-2) = 0
PHEP(1,I-2) = PHKK(1,I)
PHEP(2,I-2) = PHKK(2,I)
PHEP(3,I-2) = PHKK(3,I)
PHEP(4,I-2) = PHKK(4,I)
PHEP(5,I-2) = PHKK(5,I)
VHEP(1,I-2) = VHKK(1,I)
VHEP(2,I-2) = VHKK(2,I)
VHEP(3,I-2) = VHKK(3,I)
VHEP(4,I-2) = VHKK(4,I)
ENDDO
ISTHEP(1) = 4
ISTHEP(2) = 4
JMOHEP(1,1) = -1
JMOHEP(2,1) = -1
JMOHEP(1,2) = -1
JMOHEP(2,2) = -1
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