I’m trying to score neutrons due to spontaneous fission from uranium fuel. I am able to score alphas successfully but I’m not getting any data for neutrons. I tried running example rdecay02 with macro u238c.mac which should model the decay of stationary u238 nuclei but I’m having the same problem where there are no neutrons scored in the secondary particle output file. In my own model I have implemented the physics list and biasing physics from rdecay02 and also used /process/had/rdm/BRbias true in my macro to switch on branching ratio biasing but no luck so far. Am I missing something?
I recently updated my code from version 10.7.4 where it was working previously. I suspect it might be a physics list issue so I have attached the older version that showed the proper neutron/gamma multiplicity and the newer version (G4_11) wherein I am only able to see properly the gamma ray multiplicity and no neutron emission. I have separately verified that neutrons can be simulated correctly and set:
I have tried to set every single verbose flag that I can, but am looking for further suggestions. It compiles and no flags are raised. Even in UI mode, I can see the tracks of ions, gammas, etc. but no neutrons.
I have a hunch that there is a conflict or flag that must be set that I am skipping over in reading the release notes.
Thanks,
Justin