Clarification on why only (n,g) reaction channel is identified and others reactions are missing

Scoring
8

Dear Zeenat,

I have a few further points that I suggest investigating:

  • With respect to the literature plot you provided, a large portion of the data appears to be missing below 1 keV.
    • In your code (line 134), the energy is formatted as a floating-point number. For very low energies, this results in printing only zeros, because the number of decimal places is insufficient (this can be verified by visually inspecting one of your output files). Please switch to exponential format (I/O formatting)
  • You can simplify the PKA selection logic (see here).
    • For low-energy neutrons (ICODE = 300), the residual nuclei of interest are stored in FHEAVY.
    • It is sufficient to inspect that stack only (e.g. PKA_E = TKHEAV(1), and so on).
    • You can therefore restrict the logic to the “priority 2” branch of your routine (line 98).
  • You are normalizing the results in your Python post-processing. Could you please clarify the assumptions underlying this normalization?

Thanks and best regards,
Tommaso