|Message: Bug in dnaphysics example||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 would like to report a bug in
examples -> extended -> medical -> dna -> dnaphysics
example of Geant4 10.04.p01. I would appreciate if one of devs would confirm that it is a bug and reply to this message.
To reproduce this bug one should make these changes in "dnaphysics.in" control file:
/tracking/verbose 0 -> /tracking/verbose 1
/gun/particle proton -> /gun/particle e-
/gun/energy 100 keV -> /gun/energy 10 keV
Then run example with "dnaphysics.in" file as argument. For example:
You will see that electrons with kinetic energy of EXACTLY 10 keV can only suffer elastic scattering and cannot ionise other electrons or deposit any energy in water medium at all. The only process that occurs is "e-_G4DNAElastic". If one changes energy of electrons to 9.9 keV or 10.1 keV this problem dissapears.
Activating physics with
in PhysicsList instead of using mixed standard and DNA physics with G4EmDNAPhysicsActivator also solves the problem.
The source of the problem is likely in G4EmDNAPhysicsActivator who attaches two models for e-_G4DNAIonisation process: DNAEmfietzoglouIonisationModel (from 0 eV to 10 keV) and DNABornIonisationModel (from 10 keV to 1 MeV). So these models probably do not cover energy of exactly 10 keV so that e-_G4DNAIonisation process does not occur in simulation process at all.
This bug probably will not affect results of any simulations unless user sets energy of source to exactly 10 keV. But if he does he will obtain wrong results.
|Inline Depth:||Outline Depth:||Add message:|