IAP GITLAB

Commit 953d4188 authored by Felix Riehn's avatar Felix Riehn

updated event type definition for sibyll

git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@5282 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent af327c4b
......@@ -197,12 +197,40 @@ C SIBYLL
if(ish.ge.5)write(ifch,'(a,i5)')
$ ' number of particles from Sibyll :',NP
c save interaction type
if(JDIF(1).eq.0)then
typevt=1 !ND
elseif(JDIF(1).eq.3)then
typevt=2 !DD
if(nwd.eq.1)then
c single nucleon interaction
if(JDIF(1).eq.0)then
typevt=1 !ND
elseif(JDIF(1).eq.3)then
typevt=2 !DD
elseif(JDIF(1).eq.2)then
typevt=-4 !SD tar
else
typevt=4 !SD pro
endif
else
typevt=4 !SD
c multiple nucleon interaction
indif=0
idif_b=0
idif_t=0
idif_d=0
do jj=1,nwd
if(jdif(jj).eq.0) indif=1
if(jdif(jj).eq.1) idif_b=1
if(jdif(jj).eq.2) idif_t=1
if(jdif(jj).eq.3) idif_d=1
enddo
if(indif.eq.1)then
typevt=1
else
if(idif_d.eq.1.or.(idif_t.eq.1.and.idif_b.eq.1))then
typevt=1
elseif(idif_t.eq.1)then
typevt=-4
else
typevt=4
endif
endif
endif
do k=1,NP
......
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