fluxesfile.py: Writing fluxes files using setValue method doesn't set correct groups when small rounding errors present about pypact HOT 4 CLOSED

Hi Tom, the values attribute works fine.

Thanks,

Ross

from pypact.

Congratulations for submitting your first issue on pypact. Thanks for contributing!

from pypact.

We could add a tolerance, but I think that would defeat the purpose of that API. The setValue is intended as a convenience method for simple fluxes i.e setting a monoenergetic flux as:

# will make a 14 MeV beam in 3 lines for 175 group
ff = pp.FluxesFile()
ff.setGroup(175)
ff.setValue(14.0e6, 1.0)

# to file
pp.to_file(ff, "fluxes")

or the other use case of setValue is to permute one or two bin values after setting a base spectrum i.e. doing total MC.

ff = pp.FluxesFile()
ff.setGroup(175)
ff.values = [.....]
# change the 14 MeV bin to 0.8
ff.setValue(14.0e6, 0.8)

# to file
pp.to_file(ff, "fluxes")

If you simply just want to set a complete spectrum (all the bins) just use the values attribute.
For example:

# an array of 175 values set or retrieved from somewhere else.
myvalues = [......]

ff = pp.FluxesFile()
ff.setGroup(175)
# values are not protected - can read & write
ff.values = myvalues

# to file
pp.to_file(ff, "fluxes")

Does this help with your use case?

from pypact.

But if you still want to add the tolerance, then please make a pull request :)

from pypact.