|Message: Cross sections of the reaction Al(p,x)26Al and Si(n,x)26Al calculated by Geant4 are different with experimental data.||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 calculated the cross sections of the reaction Al(p,x)26Al and Si(n,x)26Al based on the Example/Hadr03 with Geant4.10.1. In comparison with the experimental data in the published works (Nishiizumi et al., 2009; Reedy 2013) we can find there is a great difference. The comparisons are shown in the Figs attached at bottom.
Q1: Please have a look whether it is right to calculate the cross section in the following way. I modified the example/Hadr03. I set the material as Aluminum, chose the G4HadronPhysicsQGSP_BIC_HP as the Hadron Inelastic physics model and used proton as the particle source. I change the energy of incident protons each time. For example, when the energy of incident proton is 100MeV, we get the following output file.
>>Proton+Al >>>>>>>>Energy 100 MeV Number:10000 hadElastic= 4851 protonInelastic= 5149 MeanFreePath: 21.13 cm +- 20.94 cm massic: 57.052 g/cm2 CrossSection: 0.047325 cm^-1 massic: 1.7528 mm2/g crossSection per atom: 785.32 mbarn Verification: crossSections from G4HadronicProcessStore: hadElastic= 846.82 um2/mg 379.41 mbarn protonInelastic= 910.64 um2/mg 408 mbarn total= 1.7575 mm2/g 787.41 mbarn Al26(produced in the reactions): 942
We can see there are 5149 protonInelastic reactions, the cross section of protonInelastic is 408 mbarn and 942 Al26 were produced. Then I get the cross section of Al(p,x)26Al, CS=942/5149*408mbarn.
Q2. If my calculation is reasonable, what can I do with the difference between Geant4 and the experimental data? I need the accurate production of Al26. Is there any way to simulate the production of Al26 using the experimental data?
Sorry for my poor English and my long question, I hope you can get my mean.
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2015/08/20/19.45-35756-Al-px-Al26.png http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2015/08/20/19.45-90232-Si-nx-Al26.png
|Inline Depth:||Outline Depth:||Add message:|