Re: Particle Yields from mgdraw.f

From: Stefan Roesler (sroesler@mail.cern.ch)
Date: Tue Feb 27 2007 - 15:37:16 CET

  • Next message: Lindley Winslow: "Re: Particle Yields from mgdraw.f"

    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: Lindley Winslow: "Re: Particle Yields from mgdraw.f"

    This archive was generated by hypermail 2.1.6 : Tue Feb 27 2007 - 16:40:10 CET