Calculation of radiation therapy dose using all particle Monte Carlo transport

A computer implemented process for producing a 3-dimensional map of a radiation dose delivered to a patient, comprising:

• constructing patient-dependent information necessary for a Monte-Carlo transport calculation, wherein the step of constructing patient dependent-information comprises:
• determining user-specified options from an input of (i) Monte Carlo parameters, (ii) physics options and (iii) output options, to set switches for code control;
• determining the number of energy groups for transport of each particle type from nuclear/atomic/electron data to provide the number of energy groups for each particle type;
• using computed tomography (CT) information to define dimensions and material composition for each CT voxel from an input of (i) user options comprising user-specified thresholds for processing CT scans and (ii) a CT scan array to contribute to the production of a material specification array;
• reading user-drawn contours that describe patient structures and modifying said material specification arrays from an input of a second set of user options comprising user-drawn contours, to complete the production of said material specification array and to provide a standard deviation zone identification array;
• reading user input specifying each radiation beam source from an input of radiation source specifications comprising external beam characteristics and modifiers to provide radiation source angular and energy distributions and arrays describing beam delivery components;
• reading user input specifying each internal (brachytherapy) source from an input of an internal radiation source specification to provide radiation source angular and energy distributions;
• completing the final setup for material arrays from an input of material composition data which is defined internally within the code to provide problem dependent material and isotope specification arrays;
• reading nuclear and atomic data and constructing transport arrays from an input of (i) said nuclear/atomic/electron data, (ii) said number of energy groups for each particle type and (iii) said problem-dependent material and isotope specification arrays to provide (i) nuclear and atomic transport data arrays, (ii) heavy charged particle transport data arrays and (iii) energy group structure for each particle type;
• and
• reading electron data and constructing transport arrays from an input of (i) said nuclear/atomic/electron data and (ii) said problem-dependent material and isotope specification arrays to produce electron transport data arrays;
• executing said Monte-Carlo transport calculation, wherein the step of executing said Monte Carlo transport calculation comprises:
• selecting particle attributes for a primary particle arising from an external radiation beam from an input of (i) radiation source angular and energy distributions, (ii) arrays describing beam delivery components, (iii) material data, (iv) nuclear and atomic transport data arrays and (v) number and energy group structure for each particle type, and to contribute a first portion to the provision of attributes of one particle, wherein said attributes comprise energy, location, direction and type;
• selecting particle attributes for a primary particle arising from an internal radiation source from an input of radiation source angular and energy distributions to contribute a second portion to the provision of said attributes of one particle;
• selecting a particle that has been created by an interaction of another particle in a transport mesh from an input of secondary particle arrays to provide a third portion to and complete the provision of said attributes of one particle;
• (i) tracking a neutron through said transport mesh, (ii) recording energy deposited by said neutron and (iii) storing attributes of secondary particles produced in a secondary particle array from an input of (i) attributes of one particle (said neutron), (ii) switches set for code control, (iii) a material specification array, (iv) material data, (v) nuclear transport data arrays, and (vi) the number and energy group structure for said neutrons, to provide secondary particle arrays for neutrons and a 3-D energy deposit map for neutrons;
• (i) tracking a photon through said transport mesh, (ii) recording energy deposited by said photon and (iii) storing attributes of secondary particles produced in secondary particle arrays for photons from an input of (i) attributes of one photon particle, (ii) said switches set for code control, (iii) said material specification arrays, (iv) said material data, (v) said atomic transport data arrays, and (vi) number and energy group structure for photons, to provide secondary particle arrays for photons and a 3-D energy deposit map for photons;
• (i) tracking a heavy charged particle through said transport mesh, (ii) recording the energy deposited by said heavy charged particle and (iii) storing the attributes of secondary particles produced in secondary particle arrays for heavy charged particles from an input of (i) said switch settings for code control, (ii) said material specification array, (iii) said material data, (iv) said nuclear transport data arrays, (v) heavy charged particle transport data arrays and (vi) the number and energy group structure for heavy charged particles to provide secondary particle arrays for heavy charged particles and a 3-D deposit map for heavy charged particles;
• (i) tracking a primary electron through said transport mesh, (ii) recording the energy deposited by said electron and (iii) storing the attributes of secondary particles produced in secondary particle arrays for primary electrons from an input of (i) the attributes of one particle (primary electron), (ii) said switch settings for code control, (iii) said material specification array, (iv) said material data, (v) electron transport data arrays and (vi) said number and energy group structure for said electrons to provide secondary particle arrays for primary electrons and a 3-D energy deposit map for primary electrons;
• (i) tracking secondary electrons through said transport mesh, (ii) recording the energy deposited by said secondary electron and (iii) storing the attributes of secondary particles produced (secondary electrons) in secondary particle arrays for secondary electrons from an input of (i) the attributes of one particle (a secondary electron), (ii) said switch settings for code control, (iii) said material specification array, (iv) said material data, (v) said electron transport data arrays and (vi) said number and energy group structure for electrons to provide said secondary particle arrays for secondary electrons and a 3-D energy deposit map for secondary electrons;
• adding all 3-D energy deposit maps calculated over a batch to a 3-D energy deposit map for the problem from an input of a 3-D energy deposit map calculated for a single batch to provide an integral 3-D energy deposit map;
• updating the arrays necessary for a standard deviation calculation with energy deposit information determined from each said batch from an input of said integral 3-D energy deposit map and a standard deviation zone ID array to provide standard deviation precalculation arrays;
• and
• calculating standard deviation from an input of said standard deviation precalculation arrays to provide standard deviation arrays;
• and
• producing, from said patient-dependent information and said Monte-Carlo transport calculation, a 3-dimensional map of the dose delivered to said patient.

A computer implemented process for producing a 3-dimensional map of a radiation dose delivered to a patient, comprising:

• constructing patient-dependent information necessary for a Monte-Carlo transport calculation, wherein the step of constructing patient dependent-information comprises:
• determining user-specified options from an input of (i) Monte Carlo parameters, (ii) physics options and (iii) output options, to set switches for code control;
• determining the number of energy groups for transport of each particle type from nuclear/atomic/electron data to provide the number of energy groups for each particle type;
• using computed tomography (CT) information to define dimensions and material composition for each CT voxel from an input of (i) user options comprising user-specified thresholds for processing CT scans and (ii) a CT scan array to contribute to the production of a material specification array;
• reading user-drawn contours that describe patient structures and modifying said material specification arrays from an input of a second set of user options comprising user-drawn contours, to complete the production of said material specification array and to provide a standard deviation zone identification array;
• reading user input specifying each radiation beam source from an input of radiation source specifications comprising external beam characteristics and modifiers to provide radiation source angular and energy distributions and arrays describing beam delivery components;
• reading user input specifying each internal (brachytherapy) source from an input of an internal radiation source specification to provide radiation source angular and energy distributions;
• completing the final setup for material arrays from an input of material composition data which is defined internally within the code to provide problem dependent material and isotope specification arrays;
• reading nuclear and atomic data and constructing transport arrays from an input of (i) said nuclear/atomic/electron data, (ii) said number of energy groups for each particle type and (iii) said problem-dependent material and isotope specification arrays to provide (i) nuclear and atomic transport data arrays, (ii) heavy charged particle transport data arrays and (iii) energy group structure for each particle type;
• and
• reading electron data and constructing transport arrays from an input of (i) said nuclear/atomic/electron data and (ii) said problem-dependent material and isotope specification arrays to produce electron transport data arrays;
• executing said Monte-Carlo transport calculation;
• and
• producing, from said patient-dependent information and said Monte-Carlo transport calculation, a 3-dimensional map of the dose delivered to said patient.

A computer implemented process for producing a 3-dimensional map of a radiation dose delivered to a patient, comprising:

• constructing patient-dependent information necessary for a Monte-Carlo transport calculation;
• executing said Monte-Carlo transport calculation;
• wherein the step of executing the Monte Carlo transport calculation comprises:
• selecting particle attributes for a primary particle arising from an external radiation beam from an input of (i) radiation source angular and energy distributions, (ii) arrays describing beam delivery components, (iii) material data, (iv) nuclear and atomic transport data arrays and (v) number and energy group structure for each particle type, and to contribute a first portion to the provision of attributes of one particle, wherein said attributes comprise energy, location, direction and type;
• selecting particle attributes for a primary particle arising from an internal radiation source from an input of radiation source angular and energy distributions to contribute a second portion to the provision of said attributes of one particle;
• selecting a particle that has been created by an interaction of another particle in a transport mesh from an input of secondary particle arrays to provide a third portion to and complete the provision of said attributes of one particle;
• (i) tracking a neutron through said transport mesh, (ii) recording energy deposited by said neutron and (iii) storing attributes of secondary particles produced in a secondary particle array from an input of (i) attributes of one particle (said neutron), (ii) switches set for code control, (iii) a material specification array, (iv) material data, (v) nuclear transport data arrays, and (vi) the number and energy group structure for said neutrons, to provide secondary particle arrays for neutrons and a 3-D energy deposit map for neutrons;
• (i) tracking a photon through said transport mesh, (ii) recording energy deposited by said photon and (iii) storing attributes of secondary particles produced in secondary particle arrays for photons from an input of (i) attributes of one photon particle, (ii) said switches set for code control, (iii) said material specification arrays, (iv) said material data, (v) said atomic transport data arrays, and (vi) number and energy group structure for photons, to provide secondary particle arrays for photons and a 3-D energy deposit map for photons;
• (i) tracking a heavy charged particle through said transport mesh, (ii) recording the energy deposited by said heavy charged particle and (iii) storing the attributes of secondary particles produced in secondary particle arrays for heavy charged particles from an input of (i) said switch settings for code control, (ii) said material specification array, (iii) said material data, (iv) said nuclear transport data arrays, (v) heavy charged particle transport data arrays and (vi) the number and energy group structure for heavy charged particles to provide secondary particle arrays for heavy charged particles and a 3-D deposit map for heavy charged particles;
• (i) tracking a primary electron through said transport mesh, (ii) recording the energy deposited by said electron and (iii) storing the attributes of secondary particles produced in secondary particle arrays for primary electrons from an input of (i) the attributes of one particle (primary electron), (ii) said switch settings for code control, (iii) said material specification array, (iv) said material data, (v) electron transport data arrays and (vi) said number and energy group structure for said electrons to provide secondary particle arrays for primary electrons and a 3-D energy deposit map for primary electrons;
• (i) tracking secondary electrons through said transport mesh, (ii) recording the energy deposited by said secondary electron and (iii) storing the attributes of secondary particles produced (secondary electrons) in secondary particle arrays for secondary electrons from an input of (i) the attributes of one particle (a secondary electron), (ii) said switch settings for code control, (iii) said material specification array, (iv) said material data, (v) said electron transport data arrays and (vi) said number and energy group structure for electrons to provide said secondary particle arrays for secondary electrons and a 3-D energy deposit map for secondary electrons;
• adding all 3-D energy deposit maps calculated over a batch to a 3-D energy deposit map for the problem from an input of a 3-D energy deposit map calculated for a single batch to provide an integral 3-D energy deposit map;
• updating the arrays necessary for a standard deviation calculation with energy deposit information determined from each said batch from an input of said integral 3-D energy deposit map and a standard deviation zone ID array to provide standard deviation precalculation arrays;
• and
• calculating standard deviation from an input of said standard deviation precalculation arrays to provide standard deviation arrays;
• and
• producing, from said patient-dependent information and said Monte-Carlo transport calculation, a 3-dimensional map of the dose delivered to said patient.