From: Francesco Cerutti (Francesco.Cerutti@cern.ch)
Date: Mon Jan 29 2007 - 18:37:25 CET
Dear Lindley,
> Comparing my results for events with iCode=101 and RESNUC scoring with
> What(1)=1, I get the same results for all nuclei except 6Li, 7Li, 6He. These
> nuclei always show higher numbers in the RESNUC output.
my guess is that you are missing the heavy fragments produced besides the
residual nucleus. You identify the residual nucleus through the statement
WRITE( 41.0, * ) NPHEAV, ICRES, IBRES, ERES, EKRES
then you loop over the heavy secondaries created in the nuclear
evaporation
do 20 ip = 1, NPHEAV
WRITE( 41.0, *) ip, KHEAVY(ip), TKHEAV(ip), WHEAVY(ip)
20 continue
but, if kheavy(ip) > 6 (it is the case of 6Li, 7Li, 6He, ...), you
do not record charge (icheav(KHEAVY(ip))) and mass number
(ibheav(KHEAVY(ip))) of the fragment.
Let me know if this is not the reason
Cheers
Francesco
**************************************************
Francesco Cerutti
CERN-AB
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854
This archive was generated by hypermail 2.1.6 : Tue Jan 30 2007 - 09:12:05 CET