Extracting high resolution gamma decay spectrum from FISPACT-II output file

[zach-hartwig]

Hi,

I'm exploring pypact for my FISPACT-II workflows, and it looks like a great tool. Thanks for developing.

At present, I'm primarily interested in extracting, plotting, and analyzing a highly resolved gamma decay spectrum from a given inventory time step. Unfortunately, the default gamma spectrum in the FISPACT-II output file only provides the 22 or 24 group structure, which I assume is structured for gamma dose calculations but is useless for detailed gamma line and isotopic intensity analysis.

• It looks like I can use pypact to manually extract the dominant nuclides and intensities from an output file and then in a following loop use a PRINTLIB 5 output to reconstruct the highly resolved gamma spectrum. Is there a better way with pypact that might be implemented?

• The plot on the pypact homepage ('gammaspectrumanimation.py') shows a highly resolved gamma spectrum, which is just what I want! However, when I run that exact same example locally, I get the very low resolution 24 bin gamma structure. See the plot below. Can you let me know what I might be missing?

Thanks for any insight!

[thomasms]

Hello, the example file provided indeed only has 22 or 24 bins. I don't think there is an example file on this repository with more bins, but I did try it offline with an output file with many more bins and it works fine.

If you have F-II you need to use READGG keyword with a custom group structure in order to get something other than the 22 or 24 group.

Do you have a file with custom bins already to hand?

[zach-hartwig]

Hello! Thanks for the great advice. READGG was exactly what I needed to utilize; I had glanced over it previously without recognizing that it would indeed enable me to extract what I want. I've got everything working now.

[thomasms]

You may also be interested in using the other FISPACT-II python package: actigamma if you're just interested in extracting, plotting, and analyzing a highly resolved gamma decay spectrum from a given inventory time step.