|Message: Re: Photonuclear with Bertini in g4.10.00.p01||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 Thu, 31 Jul 2014 22:05:30 GMT, Dennis H. Wright wrote:
> The link you list describes part of the implemented code. The rest is > described in Chapter IV, section 38 (page 547) of the Physics Reference > Manual ( > http://geant4.web.cern.ch/geant4/UserDocumentation/UsersGuides/PhysicsReferenceManual/fo/PhysicsReferenceManual.pdf > ), which deals with the Fermi breakup model. That model is applied only > for A < 17. For Mg the evaporation model is used, which means alphas > will again be rare at low energies. The difference in thresholds you see > is at least somewhat due to the change in physics models in going from A > < 17 to A > 17. >
Thank you for the explanation. The reason for different thresholds is clear now. I have already read the Bertini part of Physics reference manual which has a different "A threshold" for final process selection (A < 12 && "not sure I fully get it" stuff) so that it also explains why you don't have 3alpha reaction with bertini standard handling of deexcitation (is it correct?).
I'm a bit confused because I have invoked Bertini model but then forced PreCompound deexcitation inside it and I thought Bertini limits were used for the choice of the model for the "ending" process, but clearly it is not so. Therefore, the PreCompound model imposes other selection mechanism for the physics models (A > 17 or not, is it the one described in the Binary cascade?) So in the end, forcing the use of Fermi break-up will probably give not physical results in terms of rates but will trigger proper alpha kinematics as output, is it?
Do you suggest to try forcing Fermi break-up? Finding the part of the code handling that doesn't look easy to me, but I will find it eventually.
Thank you again and sorry for the sparse questions.
|Inline Depth:||Outline Depth:||Add message:|