IAP GITLAB

Commit 5afb9d34 authored by Tanguy Pierog's avatar Tanguy Pierog

add glauber table calculation in DPMJET to get proper impact parameter


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@5166 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent b4c6935d
......@@ -47,6 +47,10 @@ c fresh common block
COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)
COMMON/PYDAT3/MDCY(500,3),MDME(4000,2),BRAT(4000),KFDP(4000,5)
COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
* properties of interacting particles
COMMON /DTPRTA/ IT,ITZ,IP,IPZ,IJPROJ,IBPROJ,IJTARG,IBTARG
dimension WHAT(6)
integer istatus
real taugm
......@@ -70,6 +74,11 @@ c will be treated as nucleus in DT_INIT and NPMASS,... will be used
CALL DT_INIT(NCASES,EPROJ,NPMASS,NPCHAR,NTMASS,NTCHAR,IDPDG
+,IGLAU)
* pre-initialization of profile function for nuclear collisions*
ECM=dble(ecms)
CALL DT_XSGLAU(NPMASS,NTMASS,IJPROJ,0D0,0D0,ECM,1,1,-1)
c set decay flag in Pythia for DPMJET (after DT_INIT otherwise default is used)
call IniDkyJetset
......@@ -132,6 +141,8 @@ c-----------------------------------------------------------------------
CHARACTER*8 ANAME
COMMON /DTPART/ ANAME(210),AAM(210),GA(210),TAU(210),
& IICH(210),IIBAR(210),K1(210),K2(210)
* central particle production, impact parameter biasing
COMMON /DTIMPA/ BIMIN,BIMAX,XSFRAC,ICENTR
* Glauber formalism: collision properties
COMMON /DTGLCP/ RPROJ,RTARG,BIMPAC,
& NWTSAM,NWASAM,NWBSAM,NWTACC,NWAACC,NWBACC
......
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