A simple tool to visualize radiation damage in humans. - Toying around with pyQt.
This tools aims at visualizing the damage of different radiation types in the human body (or matter in general). It is by no means a way to reliably calculate a radiation hazard.
- Added specific pre-defined radiation types
- Energies can now be specified in ranges: "E1-E2". The energy distribution is always assumed to be logarithmically flat. I.e. E^{-1}.
In addition to the human torso and the cancer therapy example the geometry of the Radiation Assessment Detector (RAD) (see http://mslrad.boulder.swri.edu/) is also available. The simulation setup can be defined in LD50.py.
Set a geometry via importing HUMAN, CANCER, RPI, or RAD as geometry.
The RAD and RPI geometry the can be executed directly via rad.pyw, rpi.pyw.
LD50 is written in Python 2.7 and makes use of Py QT and numpy/scipy.
To get a complete and ready-to-use python setup, you can download https://code.google.com/p/pythonxy/
If you don't want to download/install the whole phythonxy package, you have to download/install
To Start the LD50, just execute LD50.py
at this stage LD50 uses a lot of simplifications.
- The model is only 2D
- The Bethe-Bloch equation is integrated via Euler method
- Neutron/gamma angular scattering are dummy values
- The magnitude of the radiation damage is given by the radius of the circles. This can give the impression that a large area of the Body is affected.
- The Torso is located in Vacuum, not in air. I.e. energy loss of particles in air is neglected.
This tool is free to use and can be freely modified and distributed. I'm grateful for any contribution to this project.
Jan Köhler