System and method for modelling of scattering in ion radiotherapy treatment planning
A method of evaluating a radiotherapy treatment plan for ion based radiotherapy, is proposed, comprising the steps of determining multiple elastic scattering of ions for scattering angles in a first angular interval having an upper limit at a selected cut-off angle by means of a model for Coulomb scattering; determining multiple elastic scattering of ions for scattering angles in a second angular interval having a lower limit at the selected cut-off angle; determining the scattering for angles in a range comprising at least a part of the first angular interval and at least a part of the second angular interval, based on the results obtained for the first and second angular interval, respectively. The method avoids the double counting of particles at large scattering angles that occurs when using conventional methods.
1. A method of modelling ion multiple scattering for planning of ion based radiotherapy, comprising the following steps, performed in a computer:
determining multiple elastic scattering of ions for scattering angles in a first angular interval having an upper limit at a selected cut-off angle by means of a model for Coulomb scattering;
determining multiple elastic scattering of ions for scattering angles in a second angular interval having a lower limit at the selected cut-off angle to reflect the net effect of both a Coulomb force and a nuclear force in the second angular interval;
determining, as a model for ion multiple scattering, the scattering for angles in a range comprising at least a part of the first angular interval and at least a part of the second angular interval, based on the results obtained for the first and second angular intervals, respectively; and
using the model for ion multiple scattering to calculate dose for use in creating ion based radiotherapy treatment plans.
2. The method according to claim 1 , wherein the cut-off angle is selected such that the nuclear elastic scattering results in a higher scattering angle than the cut-off angle for a majority of particles scattered to angles where the Coulomb interaction does not dominate.
3. The method according to claim 1 , wherein the multiple elastic scattering in the first angular interval is performed by means of a model for Coulomb multiple scattering.
4. The method according to claim 1 , wherein the multiple elastic scattering in the first angular interval is determined according to the Goudsmit-Saunderson model.
5. The method according to claim 1 , wherein the scattering angular distribution in the second angular interval is determined based on angular differential cross sections from a database or calculated by a suitable model.
6. The method according to claim 1 , wherein the cut-off angle is selected in the interval between 1 and 15 degrees.
7. A computer program product comprising computer readable code means which, when executed in a computer, will cause the computer to perform the method according to claim 1 .
8. A non-transitory computer readable medium encoded with computer executable instructions which, when run in a first computer device will cause the device to perform the method according to claim 1 .
9. A computer system comprising a processor, a data memory and a program memory, wherein the program memory comprises a non-transitory computer readable medium according to claim 8 .
10. The method according to claim 6 , wherein the selected cut-off angle is 5 degrees.