Review of technologies and procedures of clinical dosimetry for scanned ion beam radiotherapy

Study of a plastic scintillating plate-based quality assurance system for pencil beam scanning proton beams

INTRODUCTION

The purpose of this study was to evaluate a plastic scintillating plate-based beam monitoring system to perform quality assurance (QA) measurements in pencil beam scanning proton beam.

METHODS

Single spots and scanned fields were measured with the high-resolution dosimetry system, consisting of a plastic scintillation plate coupled to a camera in a dark box at the isocenter. The measurements were taken at 110-190 MeV beam energies with 30° gantry angle intervals at each energy. Spot positions were determined using the plastic scintillating plate-based dosimetry system at the isocenter for 70-230 MeV beam energies with 30° gantry angle intervals. The effect of gantry angle on dose distribution was also assessed by determining the scanning pattern for daily QA and 25 fields treated with intensity-modulated proton therapy.

RESULTS

Spot size, field flatness, and field symmetry of plastic scintillating plate-based dosimetry system were consistent with EBT3 at all investigated energies and angles. In all investigated energies and angles, the spot size measured was ±10% of the average size of each energy, the spot position measured was within ±2 mm, field flatness was within ±2%, and field symmetry was within ±1%. The mean gamma passing rates with the 3%/3 mm gamma criterion of the scanning pattern and 25 fields were 99.2% and 99.8%, respectively.

CONCLUSIONS

This system can be effective for QA determinations of spot size, spot position, field flatness, and field symmetry over 360° of gantry rotation in a time- and cost-effective manner, with spatial resolution comparable to that of EBT3 film.

Calibration method and performance of a time-of-flight detector to measure absolute beam energy in proton therapy

BACKGROUND

The beam energy is one of the most significant parameters in particle therapy since it is directly correlated to the particles' penetration depth inside the patient. Nowadays, the range accuracy is guaranteed by offline routine quality control checks mainly performed with water phantoms, 2D detectors with PMMA wedges, or multi-layer ionization chambers. The latter feature low sensitivity, slow collection time, and response dependent on external parameters, which represent limiting factors for the quality controls of beams delivered with fast energy switching modalities, as foreseen in future treatments. In this context, a device based on solid-state detectors technology, able to perform a direct and absolute beam energy measurement, is proposed as a viable alternative for quality assurance measurements and beam commissioning, paving the way for online range monitoring and treatment verification.

PURPOSE

This work follows the proof of concept of an energy monitoring system for clinical proton beams, based on Ultra Fast Silicon Detectors (featuring tenths of ps time resolution in 50 μm active thickness, and single particle detection capability) and time-of-flight techniques. An upgrade of such a system is presented here, together with the description of a dedicated self-calibration method, proving that this second prototype is able to assess the mean particles energy of a monoenergetic beam without any constraint on the beam temporal structure, neither any a priori knowledge of the beam energy for the calibration of the system.

METHODS

A new detector geometry, consisting of sensors segmented in strips, has been designed and implemented in order to enhance the statistics of coincident protons, thus improving the accuracy of the measured time differences. The prototype was tested on the cyclotron proton beam of the Trento Protontherapy Center (TPC). In addition, a dedicated self-calibration method, exploiting the measurement of monoenergetic beams crossing the two telescope sensors for different flight distances, was introduced to remove the systematic uncertainties independently from any external reference.

RESULTS

The novel calibration strategy was applied to the experimental data collected at TPC (Trento) and CNAO (Pavia). Deviations between measured and reference beam energies in the order of a few hundreds of keV with a maximum uncertainty of 0.5 MeV were found, in compliance with the clinically required water range accuracy of 1 mm.

CONCLUSIONS

The presented version of the telescope system, minimally perturbative of the beam, relies on a few seconds of acquisition time to achieve the required clinical accuracy and therefore represents a feasible solution for beam commission, quality assurance checks, and online beam energy monitoring.

Evaluation of International Atomic Energy Agency Technical Report Series-483 Detector-specific Output Correction Factor for Various Collimator Systems

Aim

In this study, a 6MV flattening filter (FF) and 6MV FF Free (FFF) photon beam small-field output factors (OF) were measured with various collimators using different detectors. The corrected OFs were compared with the treatment planning system (TPS) calculated OFs.

Materials and Methods

OF measurements were performed with four different types of collimators: Varian Millennium multi-leaf collimator (MLC), Elekta Agility MLC, Apex micro-MLC (mMLC) and a stereotactic cone. Ten detectors (four ionization chambers and six diodes) were used to perform the OF measurements at a depth of 10 cm with a source-to-surface distance of 90 cm. The corrected OF was calculated from the measurements. The corrected OFs were compared with existing TPS-generated OFs.

Results

The use of detector-specific output correction factor (OCF) in the PTW diode P detector reduced the OF uncertainty by <4.1% for 1 cm × 1 cm Sclin. The corrected OF was compared with TPS calculated OF; the maximum variation with the IBA CC01 chamber was 3.75%, 3.72%, 1.16%, and 0.90% for 5 mm stereotactic cone, 0.49 cm × 0.49 cm Apex mMLC, 1 cm × 1 cm Agility MLC, and 1 cm × 1 cm Millennium MLC, respectively.

Conclusion

The technical report series-483 protocol recommends that detector-specific OCF should be used to calculate the corrected OF from the measured OF. The implementation of OCF in the TPS commissioning will reduce the small-field OF variation by <3% for any type of detector.

Radio-fluorogenic nanoclay gel dosimeters with reduced linear energy transfer dependence for carbon-ion beam radiotherapy

PURPOSE

The precise assessment of the dose distribution of high linear energy transfer (LET) radiation remains a challenge, because the signal of most dosimeters will be saturated due to the high ionization density. Such measurements are particularly important for heavy-ion beam cancer therapy. On this basis, the present work examined the high LET effect associated with three-dimensional gel dosimetry based on radiation-induced chemical reactions. The purpose of this study was to create an ion beam radio-fluorogenic gel dosimeter with a reduced effect of LET.

METHODS

Nanoclay radio-fluorogenic gel (NC-RFG) dosimeters were prepared, typically containing 100 μM dihydrorhodamine 123 (DHR123) and 2.0 wt% nanoclay together with catalytic additives promoting Fenton or Fenton-like reactions. The radiological properties of NC-RFG dosimeters having different compositions in response to a carbon-ion beam were investigated using a fluorescence gel scanner.

RESULTS

An NC-RFG dosimeter capable of generating a fluorescence intensity distribution reflecting the carbon-ion beam dose profile was obtained. It was clarified that the reduction of the unfavorable LET dependence results from an acceleration of the reactions between DHR123 and H₂ O₂ , which is a molecular radiolysis product. The effects of varying the preparation conditions on the radiological properties of these gels were also examined. The optimum H₂ O₂ catalyst was determined to include 1 mM Fe³⁺ ions, and the addition of 100 mM pyridine was also found to increase the sensitivity.

CONCLUSIONS

This technique allows the first-ever evaluation of the depth-dose profile of a carbon-ion beam at typical therapeutic levels of several Gy without LET effect.

Implementation of Sphinx/Lynx as daily QA equipment for scanned proton and carbon ion beams

PURPOSE

Reporting on the first implementation of a proton dedicated commercial device (IBA Sphinx/Lynx) for daily Quality Assurance (QA) of scanned proton and carbon ion beams.

METHODS

Daily QA trendlines over more than 3 years for protons and more than 2 years for carbon ions have been acquired. Key daily QA parameters were reviewed, namely the spot size and position, beam range, Bragg peak width, coincidence (between beam and imaging system isocenters), homogeneity and dose.

RESULTS

The performance of the QA equipment for protons and carbon ions was evaluated. Daily QA trendlines allowed us to detect machine performance drifts and changes. The definition of tolerances and action levels is provided and compared with levels used in the literature.

CONCLUSION

The device has been successfully implemented for routine daily QA activities in a dual particle therapy facility for more than 2 years. It improved the efficiency of daily QA and provides a comprehensive QA process.

Investigating the lateral dose response functions of point detectors in proton beams

Objective. Point detector measurements in proton fields are perturbed by the volume effect originating from geometrical volume-averaging within the extended detector’s sensitive volume and density perturbations by non-water equivalent detector components. Detector specific lateral dose response functions K(x) can be used to characterize the volume effect within the framework of a mathematical convolution model, where K(x) is the convolution kernel transforming the true dose profile D(x) into the measured signal profile of a detector M(x). The aim of this work is to investigate K(x) for detectors in proton beams. Approach. The K(x) for five detectors were determined by iterative deconvolution of measurements of D(x) and M(x) profiles at 2 cm water equivalent depth of a narrow 150 MeV proton beam. Monte Carlo simulations were carried out for two selected detectors to investigate a potential energy dependence, and to study the contribution of volume-averaging and density perturbation to the volume effect. Main results. The Monte Carlo simulated and experimentally determined K(x) agree within 2.1% of the maximum value. Further simulations demonstrate that the main contribution to the volume effect is volume-averaging. The results indicate that an energy or depth dependence of K(x) is almost negligible in proton beams. While the signal reduction from a Semiflex 3D ionization chamber in the center of a gaussian shaped field with 2 mm sigma is 32% for photons, it is 15% for protons. When measuring the field with a microDiamond the trend is less pronounced and reversed with a signal reduction for protons of 3.9% and photons of 1.9%. Significance. The determined K(x) can be applied to characterize the influence of the volume effect on detectors measured signal profiles at all clinical proton energies and measurement depths. The functions can be used to derive the actual dose distribution from point detector measurements.

Resistive electrode array (REA) for radiotherapy beam monitoring and quality assurance

We have developed a new type of detector array for monitoring of radiation beams in radiotherapy. The detector has parallel-plane architecture with multiple large-area uniform thin-film electrodes. At least one of the electrodes is resistive and has multiple signal readouts spread out along its perimeter. The integral dose deposited in the detector gives rise to multiple signals that depend on the distribution of radiation with respect to resistive electrode array (REA) geometry. The purpose of the present study was to experimentally determine basic detector response to MLC collimated x-ray fields. Two detector arrays have been characterized: circular and rectangular. The current and electrostatic potential distribution within the resistive electrode are governed by the Laplace and continuity equations with boundary conditions at the border with the readouts. Measurements for pencil beams showed that signal strength depends primarily on the distances between the location of the pencil beam and the readouts. Measurements for larger irregular MLC showed that signals as a function of time are quasi-linear with respect to MLC position and are proportional to the MLC area. Derivation of clinically relevant radiation beam parameters from REA signals, such as MLC position, MLC gap size and monitor unit per MLC segment relies on the detector response model with empirical model parameters. An approximate analytical detector response model was proposed and used to fit experiment data.

Monte Carlo computation of 3D distributions of stopping power ratios in light ion beam therapy using GATE-RTion

PURPOSE

This paper presents a novel method for the calculation of three-dimensional (3D) Bragg-Gray water-to-detector stopping power ratio (s_w,det ) distributions for proton and carbon ion beams.

METHODS

Contrary to previously published fluence-based calculations of the stopping power ratio, the s_w,det calculation method used in this work is based on the specific way GATE/Geant4 scores the energy deposition. It only requires the use of the so-called DoseActor, as available in GATE, for the calculation of the s_w,det at any point of a 3D dose distribution. The simulations are performed using GATE-RTion v1.0, a dedicated GATE release that was validated for the clinical use in light ion beam therapy.

RESULTS

The Bragg-Gray water-to-air stopping power ratio (s_w,air ) was calculated for monoenergetic proton and carbon ion beams with the default stopping power data in GATE-RTion v1.0 and the new ICRU90 recommendation. The s_w,air differences between the use of the default and the ICRU90 configuration were 0.6% and 5.4% at the physical range (R₈₀ - 80% dose level in the distal dose fall-off) for a 70 MeV proton beam and a 120 MeV/u carbon ion beam, respectively. For protons, the s_w,det results for lithium fluoride, silicon, gadolinium oxysulfide, and the active layer material of EBT2 (radiochromic film) were compared with the literature and a reasonable agreement was found. For a real patient treatment plan, the 3D distributions of s_w,det in proton beams were calculated.

CONCLUSIONS

Our method was validated by comparison with available literature data. Its equivalence with Bragg-Gray cavity theory was demonstrated mathematically. The capability of GATE-RTion v1.0 for the s_w,det calculation at any point of a 3D dose distribution for simple and complex proton and carbon ion plans was presented.

A novel hybrid 3D dose reconstruction approach for pre-treatment verification of intensity modulated proton therapy plans

AIM

A novel hybrid three-dimensional (3D) dose reconstruction method, based on planar dose measured at a single shallower depth, was developed for use as patient-specific quality assurance (PSQA) of intensity modulated proton therapy (IMPT) plans. The accuracy, robustness and sensitivity of the presented method were validated for multiple IMPT plans of varying complexities.

METHODS AND MATERIALS

An in-house MATLAB program was developed to reconstruct 3D dose distribution from the planar dose (GyRBE) measured at 3 g cm^-2 depth in water or solid phantom using a MatriXX PT ion chamber array. The presented method was validated extensively for 11 single-field optimization (SFO) and multi-field optimization (MFO) plans on Proteus Plus. A total of 47 reconstructed planar doses at different depths were compared against the corresponding RayStation treatment planning system (TPS) and MatriXX PT measurement using a gamma passing rate (γ%) evaluated for 3%/3 mm. The robustness of the reconstruction method with respect to depth, energy layers, field dimensions and complexities in the spot intensity map (SIM) were analysed and compared against the standard PSQA. The sensitivity of the reconstruction method was tested for plans with intentional errors.

RESULTS

The presented reconstruction method showed excellent agreement (mean γ% > 98%) and robustness with both TPS-calculated and measured dose planes at all depths (2.97-30 g cm^-2), energy layers (82.1-225.5 MeV), field dimensions, target volume (17.7-1000 cm³) and SIMs from both SFO and MFO plans. In comparison to the overall mean ± SD γ% from standard PSQA, the reconstruction method showed reductions in mean γ% within 1% for both standard cubes and clinical plans. The reconstruction method was sensitive enough to detect intentional spot positional errors in a selected energy layer of a plan.

CONCLUSION

The presented hybrid reconstruction method is sufficiently accurate, robust and sensitive to estimate planar dose at any user-defined depth. It simplifies the measurement setup and eliminates multiple depth measurements.

Reconstruction of dose distributions for fine carbon-ion beams using iterative approximation toward carbon-knife

For the practical application of carbon-knife with fine carbon-ion beams, the quantification of the dose distribution is essential and requires a high spatial resolution. We propose a novel method to quantify dose distributions with a spatial resolution smaller than the dosimeter size. The proposed method innovates the iterative reconstruction technique. Using a diode dosimeter with a sensitive area of 1 mm², two-dimensional dose-area-product (DAP) distributions were measured at a 0.1 mm step at the surface and near the Bragg peak depths for fine carbon-ion beams of ∼1 mm size at the full width at half maximum (FWHM). Then, the dose distributions were reconstructed with a spatial resolution of 0.1 × 0.1 mm² from the measured DAP distributions. However, an unnaturally high noise was observed in the reconstructed dose distributions, which were considered to originate from the measurement reproducibility errors of the DAP distributions estimated to be 2.5%-3%. Therefore, a low-pass filtering process was implemented to reduce the errors on the reconstructed dose distributions. The optimum cut-off frequencies of the low-pass filter were estimated depending on the amplitude of the induced noise. Using the filtering process with the obtained optimum cut-off frequency, the dose distribution was quantified with an average error of approximately 3% or less with respect to the peak value, when the actual measurement had an error of 3%. In the reconstructed dose rate distributions, a steep penumbra P_80-20 ∼ 0.2 mm was observed at the surface, and a dose rate at the center axis of ∼90 Gy s^-1 and a beam size of ∼1.1 mm at FWHM near the Bragg peak were obtained. The proposed method is expected to be useful for the measurement-based determination of microbeam models for commissioning and dose distribution calculations toward carbon-knife applications.

A new detector for the beam energy measurement in proton therapy: a feasibility study

The proof of concept of a new device, capable of determining in a few seconds the energy of clinical proton beams by measuring the time of flight (ToF) of protons, is presented. The prototype consists of two thin ultra fast silicon detector (UFSD) pads, aligned along the beam direction in a telescope configuration and readout by a digitizer. The method developed for extracting the energy at the isocenter from the measured ToF, validated by Monte Carlo simulations, and the procedure used to calibrate the system are also presented and discussed in detail. The prototype was tested at the Centro Nazionale di Adroterapia Oncologica (CNAO, Pavia, Italy), at several beam energies, covering the entire clinical range, and using different distances between the sensors. The measured beam energies were benchmarked against the nominal CNAO energy values, obtained during the commissioning of the centre from the measured ranges in water. Deviations of few hundreds of keV have been achieved for all considered proton beam energies for distances between the two sensors larger than 60 cm, indicating a sensitivity to the corresponding beam range in water smaller than the clinical tolerance of 1 mm. Moreover, few seconds of irradiation were necessary to collect the required statistics. These preliminary results indicate that a telescope of UFSDs could achieve in a short time the accuracy required for the clinical application and therefore encourage further investigations towards the improvement and the optimization of the present prototype.

A GEMPix-based integrated system for measurements of 3D dose distributions in water for carbon ion scanning beam radiotherapy

Purpose

Commercially available systems for ion beam reference dosimetry in water are mainly based on ionization chambers. In those systems, a large number of small detectors are typically arranged in a two‐dimensional (2D) array or matrix to achieve high spatial resolution (order of several millimeters) and large field coverage at the same time. The goal of this work was to investigate the reliability of a detector of superior spatial resolution to perform three‐dimensional (3D) ionization measurements in carbon ion pencil beams.

Methods

The GEMPix is a small gaseous detector with a highly pixelated readout, consisting of a drift region (with 2.8 cm³ × 2.8 cm³ × 0.3 cm³ volume), three gas electron multipliers (GEMs) for signal amplification and four Timepix ASICs with 55 µm pixel pitch and a total of 262,144 pixels. An integrated system was designed and built, which consists of a commercial water phantom with a three‐axis motorized arm, a reference large‐area ionization chamber for signal normalization to the beam output and the GEMPix itself. Measurements at different depths in water have been performed at the Italian National Centre for Oncological Hadrontherapy (CNAO) with three carbon ion beam energies. Lateral beam profiles measured with the GEMPix at the shallowest depth were compared to those measured with radiochromic EBT3 films in air in the position of the reference ionization chamber. The Timepix readout was calibrated in energy by using one independent depth scan with carbon ions of 150 mm range. Bragg peak curves were also simulated using the Monte Carlo FLUKA code as a reference.

Results

Beam profiles measured with the GEMPix were smooth and showed similar shape and full width at half maximum when compared to those measured with radiochromic EBT3 films. Smooth, reproducible Bragg curves were obtained with statistical uncertainties of about 2%, matching FLUKA simulations of the Bragg curves within 15% for most data points. This difference is partially explained for the measurement with carbon ions of 150 mm range by a saturation effect in the GEMs. The high granularity of the readout allowed to produce 2D images of the deposited dose at different depths, as well as 3D data distributions.

Conclusions

This paper demonstrates the capability of the GEMPix detector to measure the 3D dose distribution of carbon ions in water for a clinical pencil beam reliably. In the future, the detector area will be increased to cover fields of scanned beams. Measurements at higher beam intensities and with protons are planned.

Accuracy assessment of the CNAO dose delivery system in the initial period of clinical activity and impact of later improvements on delivered dose distributions

PURPOSE

A retrospective analysis of the dose delivery system (DDS) performances of the initial clinical operation at CNAO (Centro Nazionale di Adroterapia Oncologica) is reported, and compared with the dose delivery accuracy following the implementation of a position feedback control.

METHODS

Log files and raw data of the DDS were analyzed for every field of patients treated with protons and carbon ions between January 2012 and April 2013 (~3800 fields). To investigate the DDS accuracy, the spot positions and the number of particles per spot measured by the DDS and prescribed by the treatment planning system were compared for each field. The impact of deviations on dose distributions was studied by comparing, through the gamma-index method, 2 three-dimensional (3D) physical dose maps (one for prescribed, one for measured data), generated by a validated dose computation software. The maximum gamma and the percentage of points with gamma ≤ 1 (passing volume) were studied as a function of the treatment day, and correlated with the deviations from the prescription in the measured number of particles and spot positions. Finally, delivered dose distributions of same treatment plans were compared before and after the implementation of a feedback algorithm for the correction of small position deviations, to study the effect on the delivery quality. A double comparison of prescribed and measured 3D maps, before and after feedback implementation, is reported and studied for a representative treatment delivered in 2012, redelivered on a polymethyl methacrylate (PMMA) block in 2018.

RESULTS

Systematic deviations of spot positions, mainly due to beam lateral offsets, were always found within 1.5 mm, with the exception of the initial clinical period. The number of particles was very stable, as possible deviations are exclusively related to the quantization error in the conversion from monitor counts to particles. For the chosen representative patient treatment, the gamma-index evaluation of prescribed and measured dose maps, before and after feedback implementation, showed a higher variability of maximum gamma for the 2012 irradiation, with respect to the reirradiation of 2018. However, the 2012 passing volume is >99.8% for the sum of all fields, which is comparable to the value of 2018, with the exception of one day with 98.2% passing volume, probably related to an instability of the accelerating system.

CONCLUSIONS

A detailed retrospective analysis of the DDS performances in the initial period of CNAO clinical activity is reported. The spot position deviations are referable to beam lateral offset fluctuations, while almost no deviation was found in the number of particles. The impact of deviations on dose distributions showed that the position feedback implementation and the increased beam control capability acquired after the first years of clinical experience led to an evident improvement in the DDS stability, evaluated in terms of gamma-index as a measure of the impact on dose distributions. However, the clinical effect of the maximum gamma variability found in the 2012 representative irradiation is mitigated by averaging along the number of fractions, and the high percentage of passing volumes confirmed the accuracy of the delivery even before the feedback implementation.

Evaluation of the IAEA-TRS 483 protocol for the dosimetry of small fields (square and stereotactic cones) using multiple detectors

The IAEA TRS 483 protocol¹ for the dosimetry of small static fields in radiotherapy was used to calculate output factors for the Elekta Synergy linac at the Australian Radiation Protection and Nuclear Safety Agency (ARPANSA). Small field output factors for both square and circular fields were measured using nine different detectors. The “corrected” output factors (ratio of detector readings multiplied by the appropriate correction factor from the protocol) showed better consistency compared to the “uncorrected” output factors (ratio of detector readings only), with the relative standard deviation decreasing by approximately 1% after the application of the relevant correction factors. Comparisons relative to an arbitrarily chosen PTW 60019 microDiamond detector showed a reduction of maximal variation for the corrected values of approximately 3%. A full uncertainty budget was prepared to analyze the consistency of the output factors. Agreement within uncertainties between all detectors and field sizes was found, except for the 15 mm circular field. The results of this study show that the application of IAEA TRS 483¹ when measuring small fields will improve the consistency of small field measurements when using multiple detectors contained within the protocol.

Dynamic lung phantom commissioning for 4D dose assessment in proton therapy

Anthropomorphic phantoms mimicking organ and tumor motion of patients are essential for end-to-end testing of motion mitigation techniques in ion beam therapy. In this work a commissioning procedure developed with the in-house designed respiratory phantom ARDOS (Advanced Radiation DOSimetry system) is presented. The phantom was tested and benchmarked for 4D dose verification in proton therapy, which included: characterization of the tissue equivalent materials from computed tomography (CT) imaging, assessment of dose calculation accuracy in critical structures of the phantom, and testing various detectors for proton dosimetry in the ARDOS phantom. To prove the validity of the CT calibration curve, measured relative stopping powers (RSP) of the ARDOS materials were compared with values from CTs: original and overwritten with known material parameters. Override of rib- and soft-tissue phantom components improved RSP accuracy while inhomogeneous lung tissue, represented by the balsa wood, was better modelled by the CT Hounsfield units. Monte Carlo (MC) dose calculations were benchmarked against measurements with a reference Farmer chamber embedded in ARDOS material samples showing less than 3% relative dose difference. Differences between MC calculated dose distributions and those calculated by analytical algorithms for the ARDOS geometry were higher than 20% of the prescribed dose, depending on the position in the phantom. Pinpoint ionization chambers and thermoluminescence dosimeters showed differences of up to 5.5% compared to MC dose calculations for all lung setups in the static phantom. They were also able to detect dose distortions due to motion. EBT3 film dosimetry was shown to be suitable for 2D relative dose characterization, which could provide extended information on dose distributions in the penumbra area. The presented methodology and results can be used for drafting general recommendations for dynamic phantom commissioning, which is an essential step towards end-to-end evaluation of motion mitigation techniques in ion beam therapy.

RIDOS: A new system for online computation of the delivered dose distributions in scanning ion beam therapy

PURPOSE

To describe a new system for scanned ion beam therapy, named RIDOS (Real-time Ion DOse planning and delivery System), which performs real time delivered dose verification integrating the information from a clinical beam monitoring system with a Graphic Processing Unit (GPU) based dose calculation in patient Computed Tomography.

METHODS

A benchmarked dose computation algorithm for scanned ion beams has been parallelized and adapted to run on a GPU architecture. A workstation equipped with a NVIDIA GPU has been interfaced through a National Instruments PXI-crate with the dose delivery system of the Italian National Center of Oncological Hadrontherapy (CNAO) to receive in real-time the measured beam parameters. Data from a patient monitoring system are also collected to associate the respiratory phases with each spot during the delivery of the dose. Using both measured and planned spot properties, RIDOS evaluates during the few seconds of inter-spill time the cumulative delivered and prescribed dose distributions and compares them through a fast γ-index algorithm.

RESULTS

The accuracy of the GPU-based algorithms was assessed against the CPU-based ones and the differences were found below 1‰. The cumulative planned and delivered doses are computed at the end of each spill in about 300 ms, while the dose comparison takes approximatively 400 ms. The whole operation provides the results before the next spill starts.

CONCLUSIONS

RIDOS system is able to provide a fast computation of the delivered dose in the inter-spill time of the CNAO facility and allows to monitor online the dose deposition accuracy all along the treatment.

Beam parameters optimization and characterization for a TUrning LInac for Protontherapy

TULIP (TUrning LInac for Protontherapy) is a novel compact accelerator system for protontherapy mounted on a rotating gantry (Amaldi et al., 2013, 2010, 2009). Its high-energy Linac has the unique property of being able to modulate the beam energy from one pulse to the next, in only a couple of milliseconds. The main purpose of this study is to optimize the properties of the beam exiting the Linac to make them compatible to medical therapy and to characterize their medical physics properties for later implementation in a Treatment Planning System. For this purpose, multi-particle tracking and Monte Carlo (MC) simulations are used to follow the particles through their path up to the treatment isocenter, following the so-called phase-space method. The data compiled includes particle fluences in air and depth-dose curves and provides the basis for a specific model of the TULIP beam.