diff --git a/src/crmc-aaa.f b/src/crmc-aaa.f index c8b4b611382cf969b15cd17a6ee28316f1fed4e7..5c204a81138d0f94fc9ea0c12104fa8fd180606b 100644 --- a/src/crmc-aaa.f +++ b/src/crmc-aaa.f @@ -350,8 +350,8 @@ c Projecticle definitions if (abs(ipart) .eq. 1) then c proton idprojin = sign(1120,ipart) !proton - laproj = sign(1,ipart) !proj Z - maproj = 1 !proj A + laproj = -1 !<0 means not a nucleus + maproj = 1 !single particle elseif (ipart .eq. 12) then c carbon idprojin=1120 @@ -365,8 +365,8 @@ c lead elseif (abs(ipart) .eq. 120) then c pi+/- idprojin = ipart !pi+/- - laproj = -1 !proj Z - maproj = 1 !proj A + laproj = -1 !<0 means not a nucleus + maproj = 1 !single particle elseif (ipart.gt.10000)then c nuclei idprojin=1120 @@ -375,8 +375,8 @@ c nuclei c PDG else idprojin=idtrafo('pdg','nxs',ipart) - laproj = -1 !proj Z - maproj = 1 !proj A + laproj = -1 !<0 means not a nucleus + maproj = 1 !single particle endif if(idprojin.eq.99)then @@ -393,8 +393,8 @@ c idtarg = 1120 !proton if ( abs(itarg) .eq. 1 ) then c proton idtargin = sign(1120,itarg) - latarg = sign(1,itarg) !targ Z - matarg = 1 !targ A + latarg = -1 !<0 means not a nucleus + matarg = 1 !single particle elseif ( itarg .eq. 12 ) then c carbon idtargin=1120 @@ -413,8 +413,8 @@ c nuclei c PDG elseif (abs(itarg).eq.2112.or.abs(itarg).eq.2212)then idtargin=idtrafo('pdg','nxs',itarg) - latarg = 1 !targ Z - matarg = 1 !targ A + latarg = -1 !<0 means not a nucleus + matarg = 1 !single particle else print *,'Warning : target particle not known : ',itarg print *,' id particle must be +/-120(pi+/-)'