Beam quality correction factors provided in current codes of practice for proton beams are approximated using the water-to-air mass stopping power ratio and by assuming the proton beam quality related perturbation correction factors to be unity. The aim of this work is to use Monte Carlo simulations to calculate energy dependent beam quality and perturbation correction factors for a set of nine ionization chambers in proton beams. Methods: The Monte Carlo code EGSnrc was used to determine the ratio of the absorbed dose to water and the absorbed dose to the sensitive air volume of ionization chambers f Q 0 related to the reference photon beam quality (60 Co). For proton beams, the quantity f Q was simulated with GATE/ Geant4 for five monoenergetic beam energies between 70 MeV and 250 MeV. The perturbation correction factors for the air cavity, chamber wall, chamber stem, central electrode, and displacement effect in proton radiation were investigated separately. Additionally, the correction factors of cylindrical chambers were investigated with and without consideration of the effective point of measurement. Results: The perturbation factors p Q were shown to deviate from unity for the investigated chambers, contradicting the assumptions made in dosimetry protocols. The beam quality correction factors for both plane-parallel and cylindrical chambers positioned with the effective point of measurement at the measurement depth were constant within 0.8%. An increase of the beam quality correction factors determined for cylindrical ionization chambers placed with their reference point at the measurement depth with decreasing energy is attributed to the displacement perturbation correction factors p dis , which were up to 1.045 AE 0.1% for the lowest energy and 1.005 AE 0.1% for the highest energy investigated. Besides p dis , the largest perturbation was found for the chamber wall where the smallest p wall determined was 0.981 AE 0.3%. Conclusions: Beam quality correction factors applied in dosimetry with cylindrical chambers in monoenergetic proton beams strongly depend on the positioning method used. We found perturbation correction factors different from unity. Consequently, the approximation of ionization chamber perturbations in proton beams by the respective water-to-air mass stopping power ratio shall be revised.

The introduction of advanced treatment techniques in proton therapy, such as intensitymodulated proton therapy, leads to an increased need for patient-specific quality assurance, especially an accurate treatment plan verification becomes inevitable. In this study, signal theoretical analysis of dose distributions in scanned proton therapy is performed to investigate the feasibility and limits of two-dimensional (2D) detector arrays for treatment plan verification. Methods: 2D detector arrays are characterized by two main aspects: the distance between the single detectors on the array or the sampling frequency; and the lateral response functions of a single detector. The analysis is based on single spots, reference fields and on measured and calculated dose distributions of typical intensity-modulated proton therapy treatment plans with and without range shifter. Measurements were performed with Gafchromic EBT3 films (Ashland Speciality Ingredients G.P., Bridgewater, NJ, USA), the MatriXX PT detector array (IBA Dosimetry, Schwarzenbruck, Germany) and the OCTAVIUS detector array 1500XDR (PTW-Freiburg, Germany) at an IBA Proteus PLUS proton therapy system (Ion Beam Applications, Louvain-la-Neuve, Belgium). Dose calculations were performed with the treatment planning system RayStation 6 or 8 (RaySearch Laboratories, Sweden). Results: The Fourier analysis of the data of the treatment planning system and film measurements show maximum frequencies of 0.06/mm for the plan with range shifter and 0.083/mm for the plan without range shifter. According to the Nyquist theorem, this corresponds to minimum required sampling distances of 8.3 and 6 mm, respectively. By comparison, the sampling distances of the arrays of 7.6 mm (MatriXX PT) and 7.1 mm (OD1500XDR) are sufficient to reconstruct the dose distributions adequately from measurements if range shifters are used, whereas some fields of the plans without range shifter violated the Nyquist requirement. The lateral dose response functions of the single detectors within the arrays have clearly higher frequencies than the treatment plans and thus the volume effect only slightly influences the measurements. Consequently, the array measurements show high gamma passing rates with at least 96 % and a good agreement between the investigated line profiles. Conclusion: The results indicate that the detector dimensions and sampling distances of the arrays are in most studied cases adequate not to substantially influence the measurement process when they are used for analyzing typical intensity-modulated proton therapy treatment plans. Nevertheless, clinical conditions have been identified, for instance treatment plans without range shifter, under which the Nyquist theorem is violated such that a full representation of the dose distributions with the measurements is not feasible. In these cases, analysis of measurements is limited to pointwise comparisons.

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.

No abstract

Purpose: The collimator monitoring fill factor (CM-FF) introduced by Stelljes et al. (2017) and the FWHM fill factor (FWHM-FF) introduced by Gago-Arias et al. (2012) were determined using the measured photon fluence response functions of various 2D-arrays. The error detection capabilities of 2D-arrays were studied by comparing detector signal changes and local gamma index passing rates in different field setups with introduced collimation errors. Methods: The fill factor is defined as the ratio of the sensitive detector area and the cell area of a detector, defined by the detector arrangement on a 2D-array. Gago-Arias et al. calculated the FWHM-FF, using the FWHM² of a detector's fluence response function K M (x) as the sensitive detector area. For the CM-FF a sensitive detector width w(D mm, d%) is calculated. The sensitive detector width is the lateral extent of K M (x), lying inside the detector cell area, along which a collimator error of D mm yields a signal change exceeding a detection threshold of d%. The sensitive area for a single detector is calculated using w(D mm, d%)². The CM-FF is then calculated as the ratio of the sensitive area of a detector within its cell area and the detector cell area. The fluence response functions of the central detector of the OCTAVIUS 729, 1500, and 1000 SRS array (all PTW-Freiburg, Freiburg, Germany) and the Map-CHECK 2 array (Sun Nuclear, Melbourne, US) were measured using a photon slit beam. The FWHM-FF and the CM-FF were calculated and compared for all 2D-arrays under investigation. The error detection capabilities of 2D-arrays in quadratic fields were studied by investigating the signal changes in the detectors adjacent to the collimator edge when changing the collimator position. The change in local gamma index passing rate with respect to the introduced collimator error was investigated for an ionization chamber and a diode array in quadratic and two intensity modulated fields. Results: Values for the CM-FF and FWHM-FF were 1.0 and 0.35, respectively for the area of the liquid-filled 1000 SRS ionization chamber array with a detector to detector distance of 5 mm and 0.32 and 0.04, respectively, for the MapCHECK 2 diode array. For the vented ionization chamber array OCTAVIUS 729 fill factors were calculated as CM-FF = 0.59 and FWHM-FF = 0.53, while the OCTAVIUS 1500 array yielded fill factors of CM-FF = 0.77 and FWHM-FF = 0.72. Signal changes in vented ionization chambers for collimator errors of 1 mm surpassed those of diodes by a factor of 2 in quadratic fields. The gamma index passing rates in quadratic fields reflect those findings. In intensity modulated fields, the decline of the gamma index passing rate is bigger for the ionization chamber array compared to the diode array when introducing collimator errors. Conclusions:The calculated values of the CM-FF correlate with the signal changes in quadratic field setups with introduced collimator position errors of 1 mm, while the FWHM-FF underestimates the error detection capabilities of 2D-arrays. An increased er...

Purpose This study evaluates the clinical use of the RUBY modular QA phantom for linac QA to validate the integrity of IGRT workflows including the congruence of machine isocenter, imaging isocenter, and room lasers. The results have been benchmarked against those obtained with widely used systems. Additionally, the RUBY phantom has been implemented to perform system QA (End‐to‐End testing) from imaging to radiation for IGRT‐based VMAT and stereotactic radiations at an Elekta Synergy linac. Material and Methods The daily check of IGRT workflow was performed using the RUBY phantom, the Penta‐Guide, and the STEEV phantom. Furthermore, Winston–Lutz tests was carried out with the RUBY phantom and a ball‐bearing phantom to determine the offsets and the diameters of the isospheres of gantry, collimator, and couch rotations, with respect to the room lasers and kV‐imaging isocenter. System QA was performed with the RUBY phantom and STEEV phantom for eight VMAT treatment plans. Additionally, the visibility of the embedded objects within these phantoms in the images and the results of CT and MR image fusions were evaluated. Results All systems used for daily QA of IGRT workflows show comparable results. Calculated shifts based on CBCT imaging agree within 1 mm to the expected values. The results of the Winston–Lutz test based on kV imaging (2D planar and CBCT) or room lasers are consistent regardless of the system tested. The point dose values in the RUBY phantom agree to the expected values calculated using algorithms in Masterplan and Monte Carlo engine in Monaco within 3% of the clinical acceptance criteria. Conclusion All the systems evaluated in this study yielded comparable results in terms of linac QA and system QA procedures. A system QA protocol has been derived using the RUBY phantom to check the IGRT‐based VMAT and stereotactic radiations workflow at an Elekta Synergy linac.

No abstract

Background and purpose As a part of the commissioning and quality assurance in proton beam therapy, lateral dose profiles and output factors have to be acquired. Such measurements can be performed with point detectors and are especially challenging in small fields or steep lateral penumbra regions as the detector's volume effect may lead to perturbations. To address this issue, this work aims to quantify and correct for such perturbations of six point detectors in small proton fields created via three different delivery techniques. Methods Lateral dose profile and output measurements of three proton beam delivery techniques (pencil beam scanning, pencil beam scanning combined with collimators, passive scattering with collimators) were performed using high‐resolution EBT3 films, a PinPoint 3D 31022 ionization chamber, a microSilicon diode 60023 and a microDiamond detector 60019 (all PTW Freiburg, Germany). Detector specific lateral dose response functions K(x,y) acting as the convolution kernel transforming the undisturbed dose distribution D(x,y) into the measured signal profiles M(x,y) were applied to quantify perturbations of the six investigated detectors in the proton fields and correct the measurements. A signal theoretical analysis in Fourier space of the dose distributions and detector's K(x,y) was performed to aid the understanding of the measurement process with regard to the combination of detector choice and delivery technique. Results Quantification of the lateral penumbra broadening and signal reduction at the fields center revealed that measurements in the pencil beam scanning fields are only compromised slightly even by large volume ionization chambers with maximum differences in the lateral penumbra of 0.25 mm and 4% signal reduction at the field center. In contrast, radiation techniques with collimation are not accurately represented by the investigated detectors as indicated by a penumbra broadening up to 1.6 mm for passive scattering with collimators and 2.2 mm for pencil beam scanning with collimators. For a 3 mm diameter collimator field, a signal reduction at field center between 7.6% and 60.7% was asserted. Lateral dose profile measurements have been corrected via deconvolution with the corresponding K(x,y) to obtain the undisturbed D(x,y). Corrected output ratios of the passively scattered collimated fields obtained for the microDiamond, microSilicon and PinPoint 3D show agreement better than 0.9% (one standard deviation) for the smallest field size of 3 mm. Conclusion Point detector perturbations in small proton fields created with three delivery techniques were quantified and found to be especially pronounced for collimated small proton fields with steep dose gradients. Among all investigated detectors, the microSilicon diode showed the smallest perturbations. The correction strategies based on detector's K(x,y) were found suitable for obtaining unperturbed lateral dose profiles and output factors. Approximation of K(x,y) by considering only the geometrical averaging effect has ...