Hello, I’m beginner at Geant4. I’m trying to calculate the neutron flux through the 70 cm of graphite moderator using Hadr04 example. The result has a data of number of collisions, track length, time of flight etc. Can i see neutron flux from result or should use another example?
You have to implement flux calculation in the example. Or in another one, eg. Hadr07
Thanks for your response.
Can you be more specific ? My input data is 5x10^7 neutrons go through the 70 cm thickness of graphite and i need to know that how many neutrons will transmit to the graphite (neutron flux in a n/cm^2 s ) and how many of them will have a thermal energy.
As a first approach, I send macros for examples Hadr06, Hadr07, Activation (in radioactivedecay subdirectory). And plots from Hadr06 and Activation.
Processing: munk06.mac.txt…
Uploading: munk07.mac.txt…
Uploading: munkActiv.mac.txt…
Uploading: munk06.gif…
Uploading: munActiv.gif…
For time being, my system is unable to upload files to Forum
if you give me an email, I can send them directly to you.
Macro file for “Hadr06.cc”
/control/verbose 1
/run/verbose 0
/testhadr/det/setMat graphite
/testhadr/det/setRadius 70 cm
/testhadr/phys/thermalScattering true
/process/em/verbose 0
/process/had/verbose 1
/run/initialize
/gun/particle neutron
/gun/energy 5 MeV
/analysis/setFileName munk06
/analysis/h1/set 16 200 1.e-7 100. keV log10 #emerging neutrons
/run/printProgress 1000
/run/beamOn 10000
Macro file for “Hadr07.cc”
/control/verbose 1
/run/verbose 0
/testhadr/det/setNbOfAbsor 1
/testhadr/det/setAbsor 1 graphite 70 cm
/testhadr/det/setSizeYZ 70 cm
/testhadr/phys/thermalScattering true
/process/em/verbose 0
/process/had/verbose 1
/run/initialize
/testhadr/gun/setDefault
/gun/particle neutron
/gun/energy 5 MeV
/run/beamOn 10000
Macro file for “Activation.cc”
/control/verbose 1
/run/verbose 1
/testhadr/det/setMat graphite
/testhadr/det/setThickness 70 cm
/testhadr/det/setSizeYZ 70 cm
/testhadr/phys/thermalScattering true
/process/em/verbose 0
/process/had/verbose 1
/run/initialize
/gun/particle neutron
/gun/energy 5 MeV
/testhadr/gun/beamSize 10 cm
/testhadr/gun/beamTime 0 s
/analysis/setFileName munkActiv
/analysis/h1/set 6 200 1.e-7 10000. keV log10 #emerging neutron
/analysis/h1/set 16 100 1.e-4 10. ms log10 #energing neutron
/run/printProgress 10000
/run/beamOn 100000
I’m just trying out uploading a file here to see if I can reproduce @maire’s problem:
test.mac.txt (409 Bytes)
This topic was automatically closed 90 days after the last reply. New replies are no longer allowed.