|Message: Re: Doubt about simulation type||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)|
To do your job first you have to get the number of times each reaction happens. This is not easy in Geant4, becuase you should use the neutron_hp physics and in this physics there is no information about the nucleus that suffered the collision. You can get it by summing up the secondary particles produced; this works always for neutron_hp except for a few exceptions, none of them Si or O. The exposure time has to be translated to number of events, changing the initial time of each event. The half lifes are also not easy to take into account, because Geant4 may decay the isotopes for you, but assigning them a time sampled from the decay time distribution, so you will have to count by yourself how many are inside the time you selected.
I have developed some tools to do what you want, which I have put in the GAMOS framework, and I honestly think you would save a lot of time by using it. Counting isotopes can be easily done by getting the number of times each reaction happens. This job will be done by a simple script like this:
/gamos/setParam GmGeometryFromText:FileName test.geom /gamos/geometry GmGeometryFromText /gamos/physicsList G4QGSP_BIC_HP /gamos/generator GmGenerator
/gamos/generator/addSingleParticleSource source neutron 10.*MeV
/gamos/userAction GmCountTracksAndStepsUA GmClassifierByHadronicReaction
and a geometry using Geant4 ASCII format:
:ELEM Silicon Si 14 28.0855 :ELEM Oxygen O 8 15.9994 :MIXT_BY_NATOMS SiO2 2.648 2 Silicon 1 Oxygen 2 :VOLU world BOX 1000. 1000. 1000. SiO2
(put your cylinder instead...) which will produce a table of this sort:
%%% COUNT_STEPS: Si28(n,p)Al28 = 123 %%% COUNT_STEPS: Si28(n,n)Si28 = 448 %%% COUNT_STEPS: O16(n,a)C13 = 180 %%% COUNT_STEPS: Si28(n,a)Mg25 = 72 ...
To take into account the time, look at the activation example, exercise 4a of tutorials/ShieldingTutorial