I am trying to get a flux spectrum with very fine bins I have approximately 48,000 bins between 0-15 MeV. As I have custom bin width as per ICRP, I am doing this by command-based scoring in a Macro, Later I use MATLAB to parse the data. I am getting a very weird spectrum (attached). Here my source energy is only 14.1 MeV But I get flux above that energies, shouldn’t there be a cutoff at 14.1? Is this a bug?
My physics list is below.
G4VModularPhysicsList* phys = new QGSP_BIC_HP(1);
I also get lot of zeros for many bins.
My Macro is also available at:
My output file is available at:
Also, this run took 21 days to finish.
Please let me know your suggestions and advice on this result.
Thanks very much in advance.