From: Francesco Cerutti (Francesco.Cerutti@cern.ch)
Date: Fri Jan 26 2007 - 14:44:01 CET
Dear Hamide,
unfortunately in your input and source there are serious
mistakes.
TKEFLK (NPFLKA) = ENERGY is the TOTAL kinetic energy of your carbon
nucleus and you write
ENERGY = (ENEUP(K) - ENELOW(K))/2
instead of
ENERGY = (ENEUP(K) + ENELOW(K))/2 (better 2d0 or TWOTWO)
thus getting non-sense low values.
By the way your bins in fact range from 1.016 MeV (i.e. 85 keV/n) and
please note that the carbon nuclear interactions will be simulated only
for the last 5 bins (above 100 MeV/n -> 1.2 GeV).
ENELOW(13) should be 0.98629D+02 and not 0.98629D+03.
MOREOVER it is not enough to link DPMJET3 and RQMD event generators (by
$FLUPRO/flutil/ldpm3qmd) to have nuclear interactions. In addition you
have to put into your input an EVENTYPE card with WHAT(3)=2.0 and
SDUM=DPMJET (you can find this problem in some previous messages).
Finally (?), if you are using your own source routine, you cannot
get an isotropic distribution by WHAT(3) of the BEAM card because
the distribution is expected to be defined in your source itself.
Good luck
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 : Fri Jan 26 2007 - 14:50:04 CET