From: Stefan Roesler (sroesler@mail.cern.ch)
Date: Tue Feb 27 2007 - 15:37:16 CET
Hi Lindley
I've had a look; the routine should be ok and give you the correct neutron
number.
If photoneutrons are important and biasing is compatible with you analysis
you may want to bias photonuclear interactions with
LAM-BIAS 0.0 0.02 0.0 7.0
It could also be better not to use coalescense since you have no model for
nucleus re-interaction. You may loose neutrons which would get liberated
by the latter.
Regards
Stefan
On Fri, 23 Feb 2007, Lindley Winslow wrote:
> Hello All,
>
> It seems from previous discussions that the correct way to get out particle
> yields is in fact from the user written routine mgdraw.f. I have managed to
> get fairly good agreement on pions and previous work with my current setup and
> my silly mgdraw.f routine that dumps things to a text file for analysis later.
> This of course leaves issues with low energy neutron reactions and the
> "multi-group" treatment.
>
> I have tried simple things with the incoming and outgoing particle weights for
> iCode=300 reactions but I am not confident in the results.
>
> I was hoping for some enlightenment as to the correct way to pull out the
> number of neutrons produced in a low energy neutron interaction.
>
> I have attached my mgdraw.f and .inp file as they are usually useful. I have
> been trying to explicitly turn on all the parts of the physics model to
> guarantee that I have photo-nuclear interactions, muon interactions and
> residual nuclei production as they are all important to this problem. If the
> experts or other on the discussion list have comments they would be
> appreciated as usual.
>
> Thank you,
> Lindley
>
>
-- ___________________________________ Stefan Roesler CERN, SC/RP CH-1211 Geneva 23 Switzerland Phone: +41-22-7679891 Fax: +41-22-7669639 E-mail: Stefan.Roesler@cern.ch
This archive was generated by hypermail 2.1.6 : Tue Feb 27 2007 - 16:40:10 CET