IAP GITLAB

Commit 87cbab94 authored by Tanguy Pierog's avatar Tanguy Pierog

delete unnecessary table for QII-04 (compressed table is enough)

add correction for history of particle in DPMJet. EPOS decay is used for particle without Breit-Wigner width because DPMJet is not able to provide more than one decay step from a final particle.


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/tags/crmc.v1.5@4274 c7a5e08c-de06-0410-9364-b41cf42a0b17
parents 4f62d659 2a798f19
......@@ -41,6 +41,10 @@ c fresh common block
DOUBLE PRECISION EPROJ
COMMON /DPMEVTINI/ EPROJ,NPMASS,NPCHAR,NTMASS,NTCHAR,IDPDG
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)
C general initialization
NCASES = -1 !skip reading steering cards
EPROJ = dble(elab)
......@@ -63,6 +67,13 @@ c will be treated as nucleus in DT_INIT and NPMASS,... will be used
c set decay flag in Pythia for DPMJET (after DT_INIT otherwise default is used)
call IniDkyJetset
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
c print *,i,PMAS(i,1),PMAS(i,2),KCHG(I,4),PMAS(i,2).GE.PARP(41)
enddo
c initialize cross-sections by calling epos-bas.f function -> models.F -> dpmjet-crmc.f
call xsigma
......
This diff is collapsed.
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