|Message: Troubles with hadronic simulations at energies above 20 MeV||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)|
Hi there I'm simulating a Voxel-Phantom and the simulations work pretty well for Photons, Electrons and Muons. But I have serial problems with Protons and Neutrons. Below about 20 - 30 MeV the simulation runs perfect. But having neutrons as primary particles with primary energies above say 50 MeV the simulations start and keep running but most of the time suddenly the calculations seem to get stuck. They don't go on. I use the EventManager to print a message every 10,000 events. That's why I can see that the calculations seem to get stuck at the simulation at 1 particle because for days (!) no additional message from the EventManager is printed in the run.out-file. I made my own hadronic physics list and here's what I use: * Neutron Elastic Processes - G4NeutronHPThermalScattering till 4 eV - G4NeutronHPElastic from 4 eV to 20 MeV - G4LElastic above 20 MeV * Neutron Inelastic Processes - G4NeutronHPInelastic till 20 MeV - G4PreCompoundModel from 20 MeV to 70 MeV as used in the LISA-Simulation - G4CascadeInterface from 70 MeV to 5 GeV (Bertini-Cascade) - G4HENeutronInelastic above 5 GeV
I really have no clue what's wrong here and why the calculations get stuck. From the outside everything looks okay. CPU running and no errors occur. Sometimes an error occurs writing: *** Error in G4HEInelastic::GenerateNBodyEvent total mass (1.87654) >= total energy (1.85209) Could this be the reason?
|Inline Depth:||Outline Depth:||Add message:|