Re: Particle Yields from mgdraw.f

From: Lindley Winslow (LAWinslow@lbl.gov)
Date: Tue Feb 27 2007 - 19:39:36 CET

  • Next message: Ercan Pilicer: "energy loss in air"

    Stefan,

    Thank you for the advice and the confirmation of my setup.

    I wanted to verify that for iCode=300 reactions that the number of
    created neutrons in one reaction is the weight of the neutrons going
    out minus the weight of the neutrons going in.

    Thank you,
    Lindley

    On Feb 27, 2007, at 6:37 AM, Stefan Roesler wrote:

    > 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


  • Next message: Ercan Pilicer: "energy loss in air"

    This archive was generated by hypermail 2.1.6 : Wed Feb 28 2007 - 09:19:48 CET