|Message: Re: Neutron energy deposit on hafnium dioxide||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)|
On Mon, 09 Oct 2017 16:54:32 GMT, michel maire wrote:
> 1- It is not directly related. Kinetic energy of secondaries which are not generated (e.g. below threshold) are counted in the dE/dx of the primary particle. I would suggest to keep default cut ( ~ 1 mm). Playing with cuts will be a later refinement. > > 2- Make sure that you have neutronHP in your physics list (compare with hafnium.out.txt). > > 3- Around 1 eV, neutron capture dominates, cross section increases, mean free path is "only" ~0.5 mm > Then ~10^7 events would allow to see few interactions. > > 4- At most energies, elastic collisions may be important. The recoil nucleus contributes to energy deposition. But, for various reasons, this recoil nucleus is not generated by default. > To create it, make sure that your physics list includes the statement setCutValue(0*mm, "proton") > Instead, you can use the command /run/setCutForAGivenParticle proton 0. mm > (if you do not understand well what I am speaking about, play with Hadr03 ...) > > 5- I tried your exercice with Hadr07. I attach macro, output and plot. As you will see, there are 94 interactions for 10^7 events. >
Wow thank you so much for your help. That is exactly what I wanted to see. I tried to replicate your simulation but I did not have hadr07 example because I am using 10.2 so I upgrade to 10.3 (Windows) but now I am having compiling error at ElectromagneticPhysics.cc where it says I am missing two members from G4EmParameters and indentifer "GetParticleIterator" being undefined.
I have updated all the data files as well and am downloading the G4TENDL1.3 just in case this data set is needed for Hadr07. I apologize if my questions are incredibly basic c++ language stuff. I am trying to google/search to find answers but running into walls.
Thank you again sincerely!
|Inline Depth:||Outline Depth:||Add message:|