IAP GITLAB

Commit 2a798f19 authored by Tanguy Pierog's avatar Tanguy Pierog

Put decay of particles without Breit-Wigner width in EPOS when running DPMJet...

Put decay of particles without Breit-Wigner width in EPOS when running DPMJet to keep full history in particle list.
DPMJet register only the first mother of final particles so if a particle "A" decay and one of its daughter
decay too then this daughter will not be in the list of the daughters of "A". If the second daughter of "A" do not decay either, 
then "A" won't be in the list at all.


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4273 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent b1cf870a
......@@ -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
......
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