Getting the uncertainty on total heating with pypact

[shimwell]

There is a section in the fispact ouptus file that has uncertainty estimates for each timestep, is it possible to access those values with pypact?

Snippet of fispact output file

                                                              UNCERTAINTY ESTIMATES (cross sections only)
                                                              -------------------------------------------

    Uncertainty estimates are based on pathway analysis for the irradiation phase
    Total Activity is          2.23849E+07  +/-   1.72E+06  Bq.        Error is   7.66E+00 % of the total.

    Total Heat Production is   2.00178E-09  +/-   1.58E-10  kW.        Error is   7.90E+00 % of the total.

    Total Gamma Dose Rate is   4.37457E-01  +/-   4.71E-02  Sv/hr.     Error is   1.08E+01 % of the total.

    Total Ingestion Dose is    1.05633E-04  +/-   1.18E-05  Sv.        Error is   1.12E+01 % of the total.

Is this +/- value accessible with pypact, I'm particularly interested in the +/- value for total heating uncertainty

I am able to read in the total heat like this, but can appear to find the uncertainty on this value

import pypact as pp
filename='fns/Ag/TENDL-2017_2000exp_5min.out'
with pp.Reader(filename) as output:
    for id in output.inventory_data:
        print(id.total_heat)

[thomasms]

Sadly this functionality is not in pypact, no, but this would be a great addition. I remember that we never implemented the pathways/uncertainty stuff. We intended to add it later but never found the time.

I can have a look at how one could achieve this but in all honesty I don't think I have the time.

I've recently started to wonder about rewriting pypact using a tool like Lark instead, and to include the pathways from the start, as this is one of the most important features of F-II.

But feel free to open a PR for this if you have something already.