Adding material decay gamma dose
Description
This PR adds a simplified gamma dose for a material. This allows users to find the dose at a distance from a gamma emitting material.
In this implementation I sample the gammas but I would be interested knowing if there is a way to do this without sampling the gammas.
The user can provide a distance from the material to get dose in pSv or the distance is assumed to be at the surface of a sphere of the material.
@tokasamwin and @samha-pf I think this feature will make comparing activated materials easier.
Fixes # (issue)
Checklist
- [x] I have performed a self-review of my own code
- [x] I have followed the style guidelines for Python source files (if applicable)
- [x] I have made corresponding changes to the documentation (if applicable)
- [x] I have added tests that prove my fix is effective or that my feature works (if applicable)
You can avoid sampling in the special (but very common) case of a discrete gamma distribution.
re: I would be interested knowing if there is a way to do this without sampling the gammas.
Can't we use Gamma-Ray Dose Constants to do that? (been thinking about it myself, and I think they'd be super useful for rapid evals) https://www.nrc.gov/docs/ML1121/ML11210B521.pdf page 10 http://www.iem-inc.com/information/tools/gamma-ray-dose-constants https://www.osti.gov/servlets/purl/1607209 <--- these guys use ICRP116 and AP https://www.osti.gov/biblio/6246345
This way it depends solely on the nuclide inventories.
re: I would be interested knowing if there is a way to do this without sampling the gammas.
Can't we use Gamma-Ray Dose Constants to do that? (been thinking about it myself, and I think they'd be super useful for rapid evals) https://www.nrc.gov/docs/ML1121/ML11210B521.pdf page 10 http://www.iem-inc.com/information/tools/gamma-ray-dose-constants https://www.osti.gov/servlets/purl/1607209 <--- these guys use ICRP116 and AP https://www.osti.gov/biblio/6246345
This way it depends solely on the nuclide inventories.
This is super interesting, shall we put this in as a separate alternative PR. Could also be nice to compare answers for each option. I do really like your proposed idea, I guess it would be super quick to calculate. It looks like these Gamma-Ray Dose Constants use a distance of 1m.
You can avoid sampling in the special (but very common) case of a discrete gamma distribution.
Thanks for this suggestion, I have had a go at this.
I've create an 'issue' feature request for this I think you can change the distance and mass, just need data on air attenuation.
re: I would be interested knowing if there is a way to do this without sampling the gammas. Can't we use Gamma-Ray Dose Constants to do that? (been thinking about it myself, and I think they'd be super useful for rapid evals) https://www.nrc.gov/docs/ML1121/ML11210B521.pdf page 10 http://www.iem-inc.com/information/tools/gamma-ray-dose-constants https://www.osti.gov/servlets/purl/1607209 <--- these guys use ICRP116 and AP https://www.osti.gov/biblio/6246345 This way it depends solely on the nuclide inventories.
This is super interesting, shall we put this in as a separate alternative PR. Could also be nice to compare answers for each option. I do really like your proposed idea, I guess it would be super quick to calculate. It looks like these Gamma-Ray Dose Constants use a distance of 1m.