IAP GITLAB

Commit a0162df4 authored by Colin Baus's avatar Colin Baus

fixed major bug in our dpmjet interface. impact parameter was not updated for...

fixed major bug in our dpmjet interface. impact parameter was not updated for each event. there exists a special common block (DTEVNO) that internally updates nevhkk. with this, glauber initialisation gets called again for each event.

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4284 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 4f002df4
......@@ -7,6 +7,8 @@ c Primary initialization for DPMJET
c-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION(A-H,O-Z)
include 'epos.inc'
c event number
COMMON /DTEVNO/ NDPMEVENT,ICASCA
c common block var iframe renamed to idframe since it is used by epos
LOGICAL LEMCCK,LHADRO,LSEADI,LEVAPO
......@@ -14,12 +16,13 @@ c common block var iframe renamed to idframe since it is used by epos
& LEMCCK,LHADRO(0:9),LSEADI,LEVAPO,IDFRAME,ITRSPT
IDFRAME=2 !dpmjet iframe variable. nucleon-nucleon frame
IDFRAME = 2 !dpmjet iframe variable. nucleon-nucleon frame
NDPMEVENT = 0 !event number needed in dpmjet
egymin=10. !min energy for model
egymax=5e4 !max energy for model
egymin = 10. !min energy for model
egymax = 5e4 !max energy for model
irescl=0 !don't rescale/skip events with wrong energy
irescl = 0 !don't rescale/skip events with wrong energy
END
......@@ -47,7 +50,7 @@ c fresh common block
C general initialization
NCASES = -1 !skip reading steering cards
EPROJ = dble(elab)
EPROJ = dble(elab)
NPMASS = maproj
NPCHAR = laproj
NTMASS = matarg
......@@ -59,7 +62,7 @@ C general initialization
c will be treated as nucleus in DT_INIT and NPMASS,... will be used
IDPDG = 0
ENDIF
IGLAU = 0
CALL DT_INIT(NCASES,EPROJ,NPMASS,NPCHAR,NTMASS,NTCHAR,IDPDG
+,IGLAU)
......@@ -70,7 +73,7 @@ c set decay flag in Pythia for DPMJET (after DT_INIT otherwise default is used)
c particle without breit-wigner decay will be decayed in EPOS to get full history
MDCY(13,1)=1 !force muon decay in Jetset
do i=100,500
if(PMAS(i,2).LT.PARP(41))MDCY(i,1)=0
if(PMAS(i,2).LT.PARP(41))MDCY(i,1)=0
c print *,i,PMAS(i,1),PMAS(i,2),KCHG(I,4),PMAS(i,2).GE.PARP(41)
enddo
......@@ -106,6 +109,8 @@ c-----------------------------------------------------------------------
& NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC
* properties of interacting particles
COMMON /DTPRTA/ IT,ITZ,IP,IPZ,IJPROJ,IBPROJ,IJTARG,IBTARG
c event number
COMMON /DTEVNO/ NDPMEVENT,ICASCA
INTEGER NPMASS,NPCHAR,NTMASS,NTCHAR,IDPDG
INTEGER IREJ, KKMAT
......@@ -122,6 +127,7 @@ c call conwr
KKMAT = -1
C call dpmjet event simulation
NDPMEVENT = NDPMEVENT+1 !needs to be updated for internal init. sets NEVHKK
IF (IDPDG.EQ.0) THEN
IDP = 1
ELSE
......@@ -193,11 +199,11 @@ c LIST is the code of final particle, P - its 4-momentum and mass.
if(ish.ge.7)write(ifch,'(a,i5,a,i5,2a,i3,1x,5(e10.4,1x),4i4)')
$ ' DPMJET particle ',k,' id :',ic,' before conversion'
$ , ' momentum :',ISTHKK(k),(sngl(PHKK(i,k)),i=1,5)
$ ,JMOHKK(1,k),JMOHKK(2,k),JDAHKK(1,k),JDAHKK(2,k)
$ ,JMOHKK(1,k),JMOHKK(2,k),JDAHKK(1,k),JDAHKK(2,k)
nptl=nptl+1 !! add 1 particle to stack
nptlhep(k)=nptl
IF(ISTHKK(k).GE.1 .AND. ISTHKK(k).LE.2 .AND. ic.NE.99999 )THEN !! if final particle
ist=ISTHKK(k)-1 !!0 means last generation other codes are e.g. for pomerons, remnants...
elseif(ic.eq.99999)then
......@@ -228,7 +234,7 @@ c LIST is the code of final particle, P - its 4-momentum and mass.
idptl(nptl)=id
c treatment of mothers and daughters
if(ist.ge.11.and.ist.le.18)then
c fix beam particles which are registered at rest in DPMJET
if(mod(ist,2).eq.0)then !target remnant (ist=12,14,16,18)
......@@ -283,7 +289,7 @@ c fix beam particles which are registered at rest in DPMJET
iorptl(nptl)=0
jorptl(nptl)=0
endif
if(JMOHKK(1,k).ne.0.AND.nptlhep(JMOHKK(1,k)).NE.0) then !mother
iorptl(nptl)=nptlhep(JMOHKK(1,k))
endif
......
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