IAP GITLAB

Commit 6d1e1260 authored by Tanguy Pierog's avatar Tanguy Pierog

update QII04 to now if collision are ND, DD or SD


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@3167 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 35b78a0e
......@@ -54,7 +54,7 @@ C C
C last modification: 26.06.2012 C
C Version qgsjet-II-04 (for CONEX) C
C C
C small corrections to adapt to CRMC : 23.07.2012 by T.Pierog C
C small corrections to adapt to CORSIKA : 29.06.2012 by T.Pierog C
C=======================================================================
......@@ -5398,6 +5398,12 @@ c-----------------------------------------------------------------------------
*,bqgs,bmaxqgs,bmaxnex,bminnex
common/jdiff/bdiff,jdiff !for external use: impact parameter
!for diffraction, diffraction type
ctp from epos
integer ng1evt,ng2evt,ikoevt
real rglevt,sglevt,eglevt,fglevt,typevt
common/c2evt/ng1evt,ng2evt,rglevt,sglevt,eglevt,fglevt,ikoevt
*,typevt !in epos.inc
external qgran
if(debug.ge.1)write (moniou,201)
......@@ -5990,6 +5996,25 @@ c check diffractive cross sections (hadron-proton only)
endif !end of diffr. check
12 bdiff=b
ctp define collision type
typevt=1
if(ia(1).eq.1.and.ia(2).eq.1)then !only for h-h
if(lqa(1).gt.0.and.lqb(1).gt.0)then !ND (at least one cut P attached)
typevt=1
elseif(nwp.eq.0.and.nwt.eq.0)then !no interaction
typevt=0
elseif((lqa(1).gt.0.and.iwt(1).eq.2)
* .or.(lqb(1).gt.0.and.iwp(1).eq.2)
* .or.(iwp(1).eq.2.and.iwt(1).eq.2))then !DD
typevt=2
elseif(lqa(1).gt.0.or.lqb(1).gt.0
* .or.iwp(1).eq.2.or.iwt(1).eq.2)then !SD
typevt=3
else !CD
typevt=1
endif
endif
c form projectile spectator part
if(debug.ge.1)write (moniou,227)
nspec=0
......
......@@ -27,7 +27,7 @@ c-----------------------------------------------------------------------
c model parameter setting
call qgset
c common initialization procedure
DATDIR="qgsjetII"
DATDIR="tabs/"
call qgaini(DATDIR)
BMAXNEX=dble(bmaxim)
BMINNEX=dble(bminim)
......
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