IAP GITLAB

Commit 6a10cbd5 authored by Colin Baus's avatar Colin Baus

moved common block to right subroutine and enabled output only in debug mode

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4127 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent c6ba9801
......@@ -8,10 +8,6 @@ c-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION(A-H,O-Z)
include 'epos.inc'
* particle properties (BAMJET index convention)
CHARACTER*8 ANAME
COMMON /DTPART/ ANAME(210),AAM(210),GA(210),TAU(210),
& IICH(210),IIBAR(210),K1(210),K2(210)
c change to center of mass system
c common block var iframe renamed to dframe since it is used by epos
......@@ -35,11 +31,16 @@ C set some variables for decays handled by epos
c-----------------------------------------------------------------------
subroutine IniEvtDpm
c-----------------------------------------------------------------------
c Setting energy,... for each event. Usefull for e.g. conex
c Setting energy, primaries,... for each event. Usefull for e.g. conex
c-----------------------------------------------------------------------
IMPLICIT DOUBLE PRECISION(A-H,O-Z)
include 'epos.inc'
* particle properties (BAMJET index convention)
CHARACTER*8 ANAME
COMMON /DTPART/ ANAME(210),AAM(210),GA(210),TAU(210),
& IICH(210),IIBAR(210),K1(210),K2(210)
INTEGER NPMASS,NPCHAR,NTMASS,NTCHAR,IDP
DOUBLE PRECISION EPROJ
COMMON /DPMEVTINI/ EPROJ,NPMASS,NPCHAR,NTMASS,NTCHAR,IDP
......@@ -146,7 +147,7 @@ c PLEASE FILL!!
phi=0.
IF(NHKK.GE.NMXHKK) THEN
WRITE(ifch,'(a)')
if(ish.ge.2)WRITE(ifch,'(a)')
&'NMAXHKK too small for event. Skipping'
GOTO 1001
ENDIF
......@@ -171,8 +172,6 @@ c final particle
istptl(nptl)=ISTHKK(k)-1 !!somehow 0 is final?
else
nptlhep(k)=0
write(ifch,'(a)')
& 'Not final particle: skipped ...'
cycle
endif !! all non final particles skipped
......@@ -202,7 +201,7 @@ c final particle
idptl(nptl)=id
write(ifch,'(a,i5,a,i5,a,i5,a,4(e10.4,1x),f6.3)')
if(ish.ge.5)write(ifch,'(a,i5,a,i5,a,i5,a,4(e10.4,1x),f6.3)')
$ ' particle from DPMJET ',nptl,' id :',idptl(nptl)
$ ,' st :',istptl(nptl)
$ , ' momentum :',(pptl(i,nptl),i=1,5)
......
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