IAP GITLAB

Commit 67586912 authored by Tanguy Pierog's avatar Tanguy Pierog

correct beam particle definition to have same events using 1 or 2212. No bug,...

correct beam particle definition to have same events using 1 or 2212. No bug, just different seed used before.


git-svn-id: https://devel-ik.fzk.de/svn/mc/crmc/trunk@4721 c7a5e08c-de06-0410-9364-b41cf42a0b17
parent 35a14372
...@@ -350,8 +350,8 @@ c Projecticle definitions ...@@ -350,8 +350,8 @@ c Projecticle definitions
if (abs(ipart) .eq. 1) then if (abs(ipart) .eq. 1) then
c proton c proton
idprojin = sign(1120,ipart) !proton idprojin = sign(1120,ipart) !proton
laproj = sign(1,ipart) !proj Z laproj = -1 !<0 means not a nucleus
maproj = 1 !proj A maproj = 1 !single particle
elseif (ipart .eq. 12) then elseif (ipart .eq. 12) then
c carbon c carbon
idprojin=1120 idprojin=1120
...@@ -365,8 +365,8 @@ c lead ...@@ -365,8 +365,8 @@ c lead
elseif (abs(ipart) .eq. 120) then elseif (abs(ipart) .eq. 120) then
c pi+/- c pi+/-
idprojin = ipart !pi+/- idprojin = ipart !pi+/-
laproj = -1 !proj Z laproj = -1 !<0 means not a nucleus
maproj = 1 !proj A maproj = 1 !single particle
elseif (ipart.gt.10000)then elseif (ipart.gt.10000)then
c nuclei c nuclei
idprojin=1120 idprojin=1120
...@@ -375,8 +375,8 @@ c nuclei ...@@ -375,8 +375,8 @@ c nuclei
c PDG c PDG
else else
idprojin=idtrafo('pdg','nxs',ipart) idprojin=idtrafo('pdg','nxs',ipart)
laproj = -1 !proj Z laproj = -1 !<0 means not a nucleus
maproj = 1 !proj A maproj = 1 !single particle
endif endif
if(idprojin.eq.99)then if(idprojin.eq.99)then
...@@ -393,8 +393,8 @@ c idtarg = 1120 !proton ...@@ -393,8 +393,8 @@ c idtarg = 1120 !proton
if ( abs(itarg) .eq. 1 ) then if ( abs(itarg) .eq. 1 ) then
c proton c proton
idtargin = sign(1120,itarg) idtargin = sign(1120,itarg)
latarg = sign(1,itarg) !targ Z latarg = -1 !<0 means not a nucleus
matarg = 1 !targ A matarg = 1 !single particle
elseif ( itarg .eq. 12 ) then elseif ( itarg .eq. 12 ) then
c carbon c carbon
idtargin=1120 idtargin=1120
...@@ -413,8 +413,8 @@ c nuclei ...@@ -413,8 +413,8 @@ c nuclei
c PDG c PDG
elseif (abs(itarg).eq.2112.or.abs(itarg).eq.2212)then elseif (abs(itarg).eq.2112.or.abs(itarg).eq.2212)then
idtargin=idtrafo('pdg','nxs',itarg) idtargin=idtrafo('pdg','nxs',itarg)
latarg = 1 !targ Z latarg = -1 !<0 means not a nucleus
matarg = 1 !targ A matarg = 1 !single particle
else else
print *,'Warning : target particle not known : ',itarg print *,'Warning : target particle not known : ',itarg
print *,' id particle must be +/-120(pi+/-)' print *,' id particle must be +/-120(pi+/-)'
......
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