IAP GITLAB

Commit e7efdae3 authored by Tanguy Pierog's avatar Tanguy Pierog

according to Sergey, energy input for QGSJETII is elab and not ekin


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4923 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 8cbb1001
......@@ -62,7 +62,7 @@ c QGSJET-II Common
if(abs(icp).gt.5)
& call utstop('Projectile not allowed in QGSJET-II !&')
c Important to calculate e0 (kinetic energy) in double precision
e0=dble(elab)-dble(amproj)
e0=dble(elab)!-dble(amproj) !lab energy is used according to Sergey
call qgini(e0,icp,maproj,matarg)
call qgini(e0,icp,maproj,matarg) !again to set bm properly
bmax=BMAXQGS
......
......@@ -57,7 +57,7 @@ c QGSJET-II Common
call iclass(idtarg,icltar)
icp=idtrafo('nxs','qgs',idproj)
if(icp.eq.0)icp=1-2*int(rangen()+0.5) !pi0=pi+ or p-
e0=dble(elab)-dble(amproj)
e0=dble(elab)!-dble(amproj) !lab energy is used according to Sergey
call qgini(e0,icp,maproj,matarg)
call qgini(e0,icp,maproj,matarg) !again to set bm properly
bmax=BMAXQGS
......
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