|Message: gamma emission from inelastic scattering of fast neutrons on C nucleus||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
I'd like to simulate the gamma emission resulting from the inelastic scattering of fast neutrons (14.1 MeV) on
carbon nucleus (12C). To do so, I've added a G4HPorLEInelasticModel instance in the inelastic process manager
attached to the neutrons particles.
As a result, I get the distribution of the emitted gamma, which, in my mind, should be more or less exclusively
made of 4.3 MeV carbon ray.
The problem is that, in addition to this ray, I also get some other peaks, corresponding to other carbon
nuclear excitation states, but with an intensity much higher than predicted (wrt nuclear database and MCNP simulations).
I tried to change the neutron inelastic model (NeutronHPorLEInelastic, Precompound, Bertini), as well as using
the LHEP_PRECO and LHEP_PRECO_HP physics list. I also verified that the G4NDL and the PhotonEvaporation
data are readable and are well pointed out by environment variables.
The results are very different from a model to another, but although the 4.3 MeV peak is always present, there is
always other peaks or continuous spectrum superimposed to this 4.3 MeV peak.
For exemple, if I use LHEP_PRECO_HP or my own physic list, where only NeutronHPorLEInelastic process is set (plus transportation),
I obtain the 4.3 MeV peak, but also 4 others peaks (namely, the 8.268, 9.637, 12.703 and 15.100 MeV peaks)
with intensities largely higher than expected.
So, my question is : does anybody already encountered this kind of problems or have a tip to help me with it ?
A possible reason would be that I don't plot only the gamma emitted from neutron inelastic interaction.
To get this spectrum, I've wrote a G4UserTrackingAction, where I record the kinetic energy of track if and only if
the track particle is a photon (PDG=22), the volume is the carbon target, the parent is track 1 (corresponding to
the 14.1 MeV neutron), and the process is "NeutronInelastic". Do you think I'm doing something wrong in this step ?
Thank in advance for any of your suggestions...
|Inline Depth:||Outline Depth:||Add message:|