Three-voltage linear method to determine ion recombination in proton and light-ion beams.

A new practical method to determine the ion recombination correction factor (k s ) for plane-parallel and Farmer-type cylindrical chambers in particle beams is investigated. Experimental data were acquired in passively scattered and scanned particle beams and compared with theoretical models developed by Boag and/or Jaffé. The new method, named the three-voltage linear method (3VL-method), is simple and consists of determining the saturation current using the current measured at three voltages in a linear region and dividing it by the current at the operating voltage (V) (even if it is not in the linear region) to obtain k s . For plane-parallel chambers, comparing k s -values obtained by model fits to values obtained using the 3VL-method, an excellent agreement is found. For cylindrical chambers, recombination is due to volume recombination only. At low voltages, volume recombination is too large and Boag's models are not applicable. However, for Farmer-type chambers (NE2571), using a smaller voltage range, limited down to 100 V, we observe a linear variation of k s with 1/V 2 or 1/V for continuous or pulsed beams, respectively. This linearity trend allows applying the 3VL-method to determine k s at any polarizing voltage. For the particle beams used, the 3VL-method gives an accurate determination of k s at any polarizing voltage. The choice of the three voltages must to be done with care to ensure to be in a linear region. For Roos-type or Markus-type chambers (i.e. chambers with an electrode spacing of 2 mm) and NE2571 chambers, the use of the 3VL-method with 300 V, 200 V and 150 V is adequate. A difference with the 2V-method and some 3V-methods in the literature is that in the 3VL-method the operational voltage does not have to be one of the three voltages. An advantage over a 2V-method is that the 3VL-method can inherently verify if the linearity condition is fulfilled.

Experimental assessment of proton dose calculation accuracy in inhomogeneous media.

PURPOSE: Proton therapy with Pencil Beam Scanning (PBS) has the potential to improve radiotherapy treatments. Unfortunately, its promises are jeopardized by the sensitivity of the dose distributions to uncertainties, including dose calculation accuracy in inhomogeneous media. Monte Carlo dose engines (MC) are expected to handle heterogeneities better than analytical algorithms like the pencil-beam convolution algorithm (PBA). In this study, an experimental phantom has been devised to maximize the effect of heterogeneities and to quantify the capability of several dose engines (MC and PBA) to handle these. METHODS: An inhomogeneous phantom made of water surrounding a long insert of bone tissue substitute (1×10×10 cm3) was irradiated with a mono-energetic PBS field (10×10 cm2). A 2D ion chamber array (MatriXX, IBA Dosimetry GmbH) lied right behind the bone. The beam energy was such that the expected range of the protons exceeded the detector position in water and did not attain it in bone. The measurement was compared to the following engines: Geant4.9.5, PENH, MCsquare, as well as the MC and PBA algorithms of RayStation (RaySearch Laboratories AB). RESULTS: For a γ-index criteria of 2%/2mm, the passing rates are 93.8% for Geant4.9.5, 97.4% for PENH, 93.4% for MCsquare, 95.9% for RayStation MC, and 44.7% for PBA. The differences in γ-index passing rates between MC and RayStation PBA calculations can exceed 50%. CONCLUSION: The performance of dose calculation algorithms in highly inhomogeneous media was evaluated in a dedicated experiment. MC dose engines performed overall satisfactorily while large deviations were observed with PBA as expected.

Ion recombination correction factor in scanned light-ion beams for absolute dose measurement using plane-parallel ionisation chambers.

Based on international reference dosimetry protocols for light-ion beams, a correction factor (k s) has to be applied to the response of a plane-parallel ionisation chamber, to account for recombination of negative and positive charges in its air cavity before these charges can be collected on the electrodes. In this work, k s for IBA PPC40 Roos-type chambers is investigated in four scanned light-ion beams (proton, helium, carbon and oxygen). To take into account the high dose-rates used with scanned beams and LET-values, experimental results are compared to a model combining two theories. One theory, developed by Jaffé, describes the variation of k s with the ionization density within the ion track (initial recombination) and the other theory, developed by Boag, describes the variation of k s with the dose rate (volume recombination). Excellent agreement is found between experimental and theoretical k s-values. All results confirm that k s cannot be neglected. The solution to minimise k s is to use the ionisation chamber at high voltage. However, one must be aware that charge multiplication may complicate the interpretation of the measurement. For the chamber tested, it was found that a voltage of 300 V can be used without further complication. As the initial recombination has a logarithmic variation as a function of 1/V, the two-voltage method is not applicable to these scanned beams.

Development and application of a water calorimeter for the absolute dosimetry of short-range particle beams.

In this work, we describe a new design of water calorimeter built to measure absorbed dose in non-standard radiation fields with reference depths in the range of 6-20 mm, and its initial testing in clinical electron and proton beams. A functioning calorimeter prototype with a total water equivalent thickness of less than 30 mm was constructed in-house and used to obtain measurements in clinical accelerator-based 6 MeV and 8 MeV electron beams and cyclotron-based 60 MeV monoenergetic and modulated proton beams. Corrections for the conductive heat transfer due to dose gradients and non-water materials was also accounted for using a commercial finite element method software package. Absorbed dose to water was measured with an associated type A standard uncertainty of approximately 0.4% and 0.2% for the electron and proton beam experiments, respectively. In terms of thermal stability, drifts were on the order of a couple of hundred µK min-1, with a short-term variation of 5-10 µK. Heat transfer correction factors ranged between 1.021 and 1.049. The overall combined standard uncertainty on the absorbed dose to water was estimated to be 0.6% for the 6 MeV and 8 MeV electron beams, as well as for the 60 MeV monoenergetic protons, and 0.7% for the modulated 60 MeV proton beam. This study establishes the feasibility of developing an absorbed dose transfer standard for short-range clinical electrons and protons and forms the basis for a transportable dose standard for direct calibration of ionization chambers in the user's beam. The largest contributions to the combined standard uncertainty were the positioning (⩽0.5%) and the correction due to conductive heat transfer (⩽0.4%). This is the first time that water calorimetry has been used in such a low energy proton beam.

LET dependence of the response of a PTW-60019 microDiamond detector in a 62MeV proton beam.

This study was initiated following conclusions from earlier experimental work, performed in a low-energy carbon ion beam, indicating a significant LET dependence of the response of a PTW-60019 microDiamond detector. The purpose of this paper is to present a comparison between the response of the same PTW-60019 microDiamond detector and an IBA Roos-type ionization chamber as a function of depth in a 62MeV proton beam. Even though proton beams are considered as low linear energy transfer (LET) beams, the LET value increases slightly in the Bragg peak region. Contrary to the observations made in the carbon ion beam, in the 62MeV proton beam good agreement is found between both detectors in both the plateau and the distal edge region. No significant LET dependent response of the PTW-60019 microDiamond detector is observed consistent with other findings for proton beams in the literature, despite this particular detector exhibiting a substantial LET dependence in a carbon ion beam.

Ion recombination correction in carbon ion beams.

PURPOSE: In this work, ion recombination is studied as a function of energy and depth in carbon ion beams. METHODS: Measurements were performed in three different passively scattered carbon ion beams with energies of 62 MeV/n, 135 MeV/n, and 290 MeV/n using various types of plane-parallel ionization chambers. Experimental results were compared with two analytical models for initial recombination. One model is generally used for photon beams and the other model, developed by Jaffé, takes into account the ionization density along the ion track. An investigation was carried out to ascertain the effect on the ion recombination correction with varying ionization chamber orientation with respect to the direction of the ion tracks. The variation of the ion recombination correction factors as a function of depth was studied for a Markus ionization chamber in the 62 MeV/n nonmodulated carbon ion beam. This variation can be related to the depth distribution of linear energy transfer. RESULTS: Results show that the theory for photon beams is not applicable to carbon ion beams. On the other hand, by optimizing the value of the ionization density and the initial mean-square radius, good agreement is found between Jaffé's theory and the experimental results. As predicted by Jaffé's theory, the results confirm that ion recombination corrections strongly decrease with an increasing angle between the ion tracks and the electric field lines. For the Markus ionization chamber, the variation of the ion recombination correction factor with depth was modeled adequately by a sigmoid function, which is approximately constant in the plateau and strongly increasing in the Bragg peak region to values of up to 1.06. Except in the distal edge region, all experimental results are accurately described by Jaffé's theory. CONCLUSIONS: Experimental results confirm that ion recombination in the investigated carbon ion beams is dominated by initial recombination. Ion recombination corrections are found to be significant and cannot be neglected for reference dosimetry and for the determination of depth dose curves in carbon ion beams.

Under-response of a PTW-60019 microDiamond detector in the Bragg peak of a 62 MeV/n carbon ion beam.

To investigate the linear energy transfer (LET) dependence of the response of a PTW-60019 Freiburg microDiamond detector, its response was compared to the response of a plane-parallel Markus chamber in a 62 MeV/n mono-energetic carbon ion beam. Results obtained with two different experimental setups are in agreement. As recommended by IAEA TRS-398, the response of the Markus chamber was corrected for temperature, pressure, polarity effects and ion recombination. No correction was applied to the response of the microDiamond detector. The ratio of the response of the Markus chamber to the response of the microDiamond is close to unity in the plateau region. In the Bragg peak region, a significant increase of the ratio is observed, which increases to 1.2 in the distal edge region. Results indicate a correlation between the under-response of the microDiamond detector and high LET values. The combined relative standard uncertainty of the results is estimated to be 2.38% in the plateau region and 12% in the distal edge region. These values are dominated by the uncertainty of alignment in the non-uniform beam and the uncertainty of range determination.

Analytical computation of prompt gamma ray emission and detection for proton range verification.

A prompt gamma (PG) slit camera prototype recently demonstrated that Bragg Peak position in a clinical proton scanned beam could be measured with 1-2 mm accuracy by comparing an expected PG detection profile to a measured one. The computation of the expected PG detection profile in the context of a clinical framework is challenging but must be solved before clinical implementation. Obviously, Monte Carlo methods (MC) can simulate the expected PG profile but at prohibitively long calculation times. We implemented a much faster method that is based on analytical processing of precomputed MC data that would allow practical evaluation of this range monitoring approach in clinical conditions. Reference PG emission profiles were generated with MC simulations (PENH) in targets consisting of either (12)C, (14)N, (16)O, (31)P or (40)Ca, with 10% of (1)H. In a given geometry, the local PG emission can then be derived by adding the contribution of each element, according to the local energy of the proton obtained by continuous slowing down approximation and the local composition. The actual incident spot size is taken into account using an optical model fitted to measurements and by super sampling the spot with several rays (up to 113). PG transport in the patient/camera geometries and the detector response are modelled by convolving the PG production profile with a transfer function. The latter is interpolated from a database of transfer functions fitted to MC data (PENELOPE) generated for a photon source in a cylindrical phantom with various radiuses and a camera placed at various positions. As a benchmark, the analytical model was compared to MC and experiments in homogeneous and heterogeneous phantoms. Comparisons with MC were also performed in a thoracic CT. For all cases, the analytical model reproduced the prediction of the position of the Bragg peak computed with MC within 1 mm for the camera in nominal configuration. When compared to measurements, the shape of the profiles was well reproduced and agreement for the estimation of the position of the Bragg peak was within 2.7 mm on average (1.4 mm standard deviation). On a non-optimized MATLAB code, computation time with the analytical model is between 0.3 to 10 s depending on the number of rays simulated per spot. The analytical model can be further used to determine which spots are the best candidates to evaluate the range in clinical conditions and eventually correct for over- and under-shoots depending on the acquired PG profiles.

Influence of the physical data to calibrate TEPCs.

Tissue-equivalent proportional counters (TEPCs) measure distributions of ionisations, produced in the gas cavity by the radiation field which are afterwards converted into distributions of energy imparted by applying a calibration factor. To calibrate the pulse-height spectra, first, a marker point must be identified in the measured spectrum. Then, an accurate value of lineal energy must be assigned to this marker. A common marker that is often used for calibration is the so-called proton-edge (p-edge). It is a distinctive feature of a proton or neutron spectrum which corresponds to the maximum amount of energy that a proton can deposit in the active volume of the detector. A precise method to identify the marker point was applied to identify the p-edge with an uncertainty below 1 %. To evaluate the final uncertainty of the calibration, the uncertainty of the energy value assigned to the p-edge must also be considered. This value can be evaluated using different energy-range tables. This study investigates how the choice of different input databases for calibration purposes influences the calibration. The effect of three different frequently used sets of input data was analysed for pure propane gas and for propane-TE gas mixture.

Reference dosimetry for light-ion beams based on graphite calorimetry.

Developments in hadron therapy require efforts to improve the accuracy of the dose delivered to a target volume. Here, the determination of the absorbed dose under reference conditions was analysed. Based on the International Atomic Energy Agency TRS-398 code of practice, for hadron beams, the combined standard uncertainty on absorbed dose to water under reference conditions, derived from ionisation chambers, is too large. This uncertainty is dominated by the beam quality correction factors, [Formula: see text], mainly due to the mean energy to produce one ion pair in air, wair. A method to reduce this uncertainty is to carry out primary dosimetry, using calorimetry. A [Formula: see text]-value can be derived from a direct comparison between calorimetry and ionometry. Here, this comparison is performed using a graphite calorimeter in an 80-MeV A(-1) carbon ion beam. Assuming recommended TRS-398 values of water-to-graphite stopping power ratio and the perturbation factor for an ionisation chamber, preliminary results indicate a wair-value of 35.5 ± 0.9 J C(-1).

Monte Carlo tools to supplement experimental microdosimetric spectra.

Tissue-equivalent proportional counters (TEPCs) are widely used in experimental microdosimetry for characterising the radiation quality in radiation protection and radiation therapy environments. Generally, TEPCs are filled with tissue-equivalent gas mixtures, at low gas pressure, to simulate tissue site sizes similar to the cell nucleus (1 or 2 µm). The TEPC response using Monte Carlo (MC) codes can be applied to supplement experimental measurements. Most of general-purpose MC codes currently available recourse to the condensed-history approach to model the electron transport and do not transport low-energy electrons (<1 keV), which can lead to systematic errors, especially in thin layers and in gas-condensed medium interfaces. In this work, a comparison between experimental microdosimetric spectra of (60)Co and (137)Cs radiation at different simulated sizes (from 1.0 to 3.0 µm) in pure propane versus simulated spectra obtained with two general-purpose codes FLUKA and PENELOPE, which include a detailed simulation of electron-photon transport in arbitrary materials, including gases, is presented.

Extension of PENELOPE to protons: simulation of nuclear reactions and benchmark with Geant4.

PURPOSE: Describing the implementation of nuclear reactions in the extension of the Monte Carlo code (MC) PENELOPE to protons (PENH) and benchmarking with Geant4. METHODS: PENH is based on mixed-simulation mechanics for both elastic and inelastic electromagnetic collisions (EM). The adopted differential cross sections for EM elastic collisions are calculated using the eikonal approximation with the Dirac-Hartree-Fock-Slater atomic potential. Cross sections for EM inelastic collisions are computed within the relativistic Born approximation, using the Sternheimer-Liljequist model of the generalized oscillator strength. Nuclear elastic and inelastic collisions were simulated using explicitly the scattering analysis interactive dialin database for (1)H and ICRU 63 data for (12)C, (14)N, (16)O, (31)P, and (40)Ca. Secondary protons, alphas, and deuterons were all simulated as protons, with the energy adapted to ensure consistent range. Prompt gamma emission can also be simulated upon user request. Simulations were performed in a water phantom with nuclear interactions switched off or on and integral depth-dose distributions were compared. Binary-cascade and precompound models were used for Geant4. Initial energies of 100 and 250 MeV were considered. For cases with no nuclear interactions simulated, additional simulations in a water phantom with tight resolution (1 mm in all directions) were performed with FLUKA. Finally, integral depth-dose distributions for a 250 MeV energy were computed with Geant4 and PENH in a homogeneous phantom with, first, ICRU striated muscle and, second, ICRU compact bone. RESULTS: For simulations with EM collisions only, integral depth-dose distributions were within 1%/1 mm for doses higher than 10% of the Bragg-peak dose. For central-axis depth-dose and lateral profiles in a phantom with tight resolution, there are significant deviations between Geant4 and PENH (up to 60%/1 cm for depth-dose distributions). The agreement is much better with FLUKA, with deviations within 3%/3 mm. When nuclear interactions were turned on, agreement (within 6% before the Bragg-peak) between PENH and Geant4 was consistent with uncertainties on nuclear models and cross sections, whatever the material simulated (water, muscle, or bone). CONCLUSIONS: A detailed and flexible description of nuclear reactions has been implemented in the PENH extension of PENELOPE to protons, which utilizes a mixed-simulation scheme for both elastic and inelastic EM collisions, analogous to the well-established algorithm for electrons/positrons. PENH is compatible with all current main programs that use PENELOPE as the MC engine. The nuclear model of PENH is realistic enough to give dose distributions in fair agreement with those computed by Geant4.

Conversion from dose-to-graphite to dose-to-water in an 80 MeV/A carbon ion beam.

Based on experiments and numerical simulations, a study is carried out pertaining to the conversion of dose-to-graphite to dose-to-water in a carbon ion beam. This conversion is needed to establish graphite calorimeters as primary standards of absorbed dose in these beams. It is governed by the water-to-graphite mass collision stopping power ratio and fluence correction factors, which depend on the particle fluence distributions in each of the two media. The paper focuses on the experimental and numerical determination of this fluence correction factor for an 80 MeV/A carbon ion beam. Measurements have been performed in the nuclear physics laboratory INFN-LNS in Catania (Sicily, Italy). The numerical simulations have been made with a Geant4 Monte Carlo code through the GATE simulation platform. The experimental data are in good agreement with the simulated results for the fluence correction factors and are found to be close to unity. The experimental values increase with depth reaching 1.010 before the Bragg peak region. They have been determined with an uncertainty of 0.25%. Different numerical results are obtained depending on the level of approximation made in calculating the fluence correction factors. When considering carbon ions only, the difference between measured and calculated values is maximal just before the Bragg peak, but its value is less than 1.005. The numerical value is close to unity at the surface and increases to 1.005 near the Bragg peak. When the fluence of all charged particles is considered, the fluence correction factors are lower than unity at the surface and increase with depth up to 1.025 before the Bragg peak. Besides carbon ions, secondary particles created due to nuclear interactions have to be included in the analysis: boron ions ((10)B and (11)B), beryllium ions ((7)Be), alpha particles and protons. At the conclusion of this work, we have the conversion of dose-to-graphite to dose-to-water to apply to the response of a graphite calorimeter in an 80 MeV/A carbon ion beam. This conversion consists of the product of two contributions: the water-to-graphite electronic mass collision stopping power ratio, which is equal to 1.115, and the fluence correction factor which varies linearly with depth, as k(fl, all) = 0.9995 + 0.0048(zw-eq). The latter has been determined on the basis of experiments and numerical simulations.

Reproducibility of three-dimensional cephalometric landmarks in cone-beam and low-dose computed tomography.

OBJECTIVES: The purpose of this study is to compare the reproducibility of three-dimensional cephalometric landmarks on three-dimensional computed tomography (3D-CT) surface rendering using clinical protocols based on low-dose (35-mAs) spiral CT and cone-beam CT (I-CAT). The absorbed dose levels for radiosensitive organs in the maxillofacial region during exposure in both 3D-CT protocols were also assessed. MATERIALS AND METHODS: The study population consisted of ten human dry skulls examined with low-dose CT and cone-beam CT. Two independent observers identified 24 cephalometric anatomic landmarks at 13 sites on the 3D-CT surface renderings using both protocols, with each observer repeating the identification 1 month later. A total of 1,920 imaging measurements were performed. Thermoluminescent dosimeters were placed at six sites around the thyroid gland, the submandibular glands, and the eyes in an Alderson phantom to measure the absorbed dose levels. RESULTS: When comparing low-dose CT and cone-beam CT protocols, the cone-beam CT protocol proved to be significantly more reproducible for four of the 13 anatomical sites. There was no significant difference between the protocols for the other nine anatomical sites. Both low-dose and cone-beam CT protocols were equivalent in dose absorption to the eyes and submandibular glands. However, thyroid glands were more irradiated with low-dose CT. CONCLUSIONS: Cone-beam CT was more reproducible and procured less irradiation to the thyroid gland than low-dose CT. CLINICAL RELEVANCE: Cone-beam CT should be preferred over low-dose CT for developing three-dimensional bony cephalometric analyses.

Validation of XiO Electron Monte Carlo-based calculations by measurements in a homogeneous phantom and by EGSnrc calculations in a heterogeneous phantom.

The purpose of the present study is to perform a clinical validation of a new commercial Monte Carlo (MC) based treatment planning system (TPS) for electron beams, i.e. the XiO 4.60 electron MC (XiO eMC). Firstly, MC models for electron beams (4, 8, 12 and 18 MeV) have been simulated using BEAMnrc user code and validated by measurements in a homogeneous water phantom. Secondly, these BEAMnrc models have been set as the reference tool to evaluate the ability of XiO eMC to reproduce dose perturbations in the heterogeneous phantom. In the homogeneous phantom calculations, differences between MC computations (BEAMnrc, XiO eMC) and measurements are less than 2% in the homogeneous dose regions and less than 1 mm shifting in the high dose gradient regions. As for the heterogeneous phantom, the accuracy of XiO eMC has been benchmarked with predicted BEAMnrc models. In the lung tissue, the overall agreement between the two schemes lies under 2.5% for the most tested dose distributions at 8, 12 and 18 MeV and is better than the 4 MeV one. In the non-lung tissue, a good agreement has been found between BEAMnrc simulation and XiO eMC computation for 8, 12 and 18 MeV. Results are worse in the case of 4 MeV calculations (discrepancies ≈ 4%). XiO eMC can predict dose perturbation induced by high-density heterogeneities for 8, 12 and 18 MeV. However, significant deviations found in the case of 4 MeV demonstrate that caution is necessary in using XiO eMC at lower electron energies.

Evaluation of Gafchromic® EBT3 films characteristics in therapy photon, electron and proton beams.

PURPOSE: To evaluate the uncertainties and characteristics of radiochromic film-based dosimetry system using the EBT3 model Gafchromic(®) film in therapy photon, electron and proton beams. MATERIAL AND METHODS: EBT3 films were read using an EPSON Expression 10000XL/PRO scanner. They were irradiated in five beams, an Elekta SL25 6 MV and 18 MV photon beam, an IBA 100 MeV 5 × 5 cm(2) proton beam delivered by pencil-beam scanning, a 60 MeV fixed proton beam and an Elekta SL25 6 MeV electron beam. Reference dosimetry was performed using a FC65-G chamber (Elekta beam), a PPC05 (IBA beam) and both Markus 1916 and PPC40 Roos ion-chambers (60 MeV proton beam). Calibration curves of the radiochromic film dosimetry system were acquired and compared within a dose range of 0.4-10 Gy. An uncertainty budget was estimated on films irradiated by Elekta SL25 by measuring intra-film and inter-film reproducibility and uniformity; scanner uniformity and reproducibility; room light and film reading delay influences. RESULTS: The global uncertainty on acquired optical densities was within 0.55% and could be reduced to 0.1% by placing films consistently at the center of the scanner. For all beam types, the calibration curves are within uncertainties of measured dose and optical densities. The total uncertainties on calibration curve due to film reading and fitting were within 1.5% for photon and proton beams. For electrons, the uncertainty was within 2% for dose superior to 0.8 Gy. CONCLUSIONS: The low combined uncertainty observed and low beam and energy-dependence make EBT3 suitable for dosimetry in various applications.

Impact of the number of discrete angles used during dose computation for TomoTherapy treatments.

PURPOSE: To quantify systematically the effect on accuracy of discretizing gantry rotation during the dose calculation process of TomoTherapy treatments. METHODS: Up to version 4.0.x included, TomoTherapy treatment planning system (TPS) approximates gantry rotation by computing dose from 51 discrete angles corresponding to the center of the projections used to control the binary multileaf collimator. Potential effects on dose computation accuracy for off-axis targets and low modulation factors have been shown previously for a few treatment configurations. In versions 4.1.x and later, TomoTherapy oversamples the projections to better account for gantry rotation, but only during full scatter optimization and final calculation (i.e., not during optimization in "beamlet" mode). The effect on accuracy of changing the number of angles was quantified with the following framework: (1) predict the impact of the discretization of gantry rotation for various modulation factors, target sizes, and off-axis positions using a simplified analytical algorithm; (2) perform regular quality assurance using measurements with EDR2 radiographic films; (3) isolating the effect of changing the number of discretized angles only (51, 153, and 459) using a previously validated Monte Carlo model (TomoPen). The diameters of the targets were 2, 3, and 5 cm; off-axis central positions of target volumes were 5, 10 and 15, and 17 cm (when accepted by the treatment unit); planned modulation factors were 1.3 and 2.0. RESULTS: For extreme configurations (3 cm tumor, 1.3 modulation factor, 15 cm off-axis position), effects on dose distributions were significant with 89.3% and 95.4% of the points passing gamma tests with 2%∕2 mm and 3%∕3 mm criteria, respectively, for TPS software version 4.0.x (51 gantry angles). The passing rate was 100% for both gamma criteria for the 4.1.x version (153 gantry angles). Those differences could be attributed almost completely to gantry motion discretization using TomoPen. Using 51 gantry angles for dose computation, TomoPen reproduced within statistical uncertainties (<1% standard deviation) dose distributions computed with version 4.0.x. Using 153 and 459 gantry angles, TomoPen reproduced within statistical uncertainties measurements and dose distributions computed with version 4.1.x. CONCLUSIONS: When low modulation factors and significant off-axis positions are used, accounting for gantry rotation during dose computation using at least 153 gantry angles is required to ensure optimal accuracy.

Monte Carlo computed machine-specific correction factors for reference dosimetry of TomoTherapy static beam for several ion chambers.

PURPOSE: To determine k(Q(msr),Q(o) ) (f(msr),f(o) ) correction factors for machine-specific reference (msr) conditions by Monte Carlo (MC) simulations for reference dosimetry of TomoTherapy static beams for ion chambers Exradin A1SL, A12; PTW 30006, 31010 Semiflex, 31014 PinPoint, 31018 microLion; NE 2571. METHODS: For the calibration of TomoTherapy units, reference conditions specified in current codes of practice like IAEA∕TRS-398 and AAPM∕TG-51 cannot be realized. To cope with this issue, Alfonso et al. [Med. Phys. 35, 5179-5186 (2008)] described a new formalism introducing msr factors k(Q(msr),Q(o) ) (f(msr),f(o) ) for reference dosimetry, applicable to static TomoTherapy beams. In this study, those factors were computed directly using MC simulations for Q(0) corresponding to a simplified (60)Co beam in TRS-398 reference conditions (at 10 cm depth). The msr conditions were a 10 × 5 cm(2) TomoTherapy beam, source-surface distance of 85 cm and 10 cm depth. The chambers were modeled according to technical drawings using the egs++ package and the MC simulations were run with the egs_chamber user code. Phase-space files used as the source input were produced using PENELOPE after simulation of a simplified (60)Co beam and the TomoTherapy treatment head modeled according to technical drawings. Correlated sampling, intermediate phase-space storage, and photon cross-section enhancement variance reduction techniques were used. The simulations were stopped when the combined standard uncertainty was below 0.2%. RESULTS: Computed k(Q(msr),Q(o) ) (f(msr),f(o) ) values were all close to one, in a range from 0.991 for the PinPoint chamber to 1.000 for the Exradin A12 with a statistical uncertainty below 0.2%. Considering a beam quality Q defined as the TPR(20,10) for a 6 MV Elekta photon beam (0.661), the additional correction k(Q(msr,)Q) (f(msr,)f(ref) ) to k(Q,Q(o) ) defined in Alfonso et al. [Med. Phys. 35, 5179-5186 (2008)] formalism was in a range from 0.997 to 1.004. CONCLUSION: The MC computed factors in this study are in agreement with measured factors for chamber types already studied in literature. This work provides msr correction factors for additional chambers used in reference dosimetry. All of them were close to one (within 1%).

SU-E-T-118: Characterization of EBT3 Films in Photon and Proton Beams.

PURPOSE: To measure the calibration curves of EBT3 dosimetry films in photon and proton beams and to quantify the related uncertainties from one beam type to another. METHODS: EBT3 Gafchromic films have similar properties than EBT2 with a symmetric construction and a matte polyester substrate to prevent Newton's ring artefacts. Films from a same batch were exposed in three different beam qualities, an Elekta SL25 6 MV photon beam, a 100 MeV 5×5cm2 proton beam delivered by pencil-beam scanning dedicated system from IBA and a 60 MeV fixed proton beam (2.5cm in diameter) at Clatterbridge Center for Oncology (CCO), UK. The films were read using an EPSON 10000 XL/PRO scanner. Film calibration curves were acquired for all modalities within a range of 0.05 to 20 Gy. Influence of increasing linear-energy transfer (LET) on film response was investigated by comparing dose measured by EBT3 to a silicon diode detector in depth for a fully-modulated beam using the CCO beam line (homogeneous dose with distal end at 3.1cm in water). A comprehensive uncertainty budget (reproducibility, uniformity'¦) was estimated on films irradiated by Elekta SL25. RESULTS: The main source of uncertainty was the non-uniformity of the scanner response. By placing all the irradiated films at the center of the scanner, the uncertainty could be reduced from 5.8% to 1.9% (1 sigma). For all beams and energies, the calibration curves were matched within uncertainties. Along the fully-modulated depth dose curve, diode and EBT3 measurement were in a 4% agreement point-to-point, indicating films weak dependence with LET. CONCLUSIONS: The weak influence of LET, beam type and energy on film response as well as its small uncertainty make EBT3 suitable for relative dosimetry and a promising candidate for measuring correction factors (quality, recombination,'¦) for reference dosimetry with ion chambers of non-standard beams (e.g pencil-beam scanning proton-therapy). â€Å“This work is supported by the Walloon Region under the project name InVivoIGT, convention number 1017266.â€.

SU-E-T-146: Reference Dosimetry for Protons and Light-Ion Beams Based on Graphite Calorimetry.

PURPOSE: The IAEA TRS-398 code of practice can be applied for the measurement of absorbed dose to water under reference conditions with an ionization chamber. For protons, the combined relative standard uncertainty on those measurements is less than 2% while for light-ion beams, it is considerably larger, i.e. 3.2%, mainly due to the higher uncertainty contributions for the water to air stopping power ration and the W air-value on the beam quality correction factors kQ,Q0 . To decrease this uncertainty, a quantification of kQ,Q0 is proposed using a primary standard level graphite calorimeter. This work includes numerical and experimental determinations of dose conversion factors to derive dose to water from graphite calorimetry. It also reports on the first experimental data obtained with the graphite calorimeter in proton, alpha and carbon ion beams. METHODS: Firstly, the dose conversion has been calculated with by Geant4 Monte-Carlo simulations through the determination of the water to graphite stopping power ratio and the fluence correction factor. The latter factor was also derived by comparison of measured ionization curves in graphite and water. Secondly, kQ,Q0 was obtained by comparison of the dose response of ionization chambers with that of the calorimeter. RESULTS: Stopping power ratios are found to vary by no more than 0.35% up to the Bragg peak, while fluence correction factors are shown to increase slightly above unity close to the Bragg peak. The comparison of the calorimeter with ionization chambers is currently under analysis. For the modulated proton beam, preliminary results on W air confirm the value recommended in TRS-398. Data in both the non-modulated proton and light-ion beams indicate higher values but further investigation of heat loss corrections is needed. CONCLUSIONS: The application of graphite calorimetry to proton, alpha and carbon ion beams has been demonstrated successfully. Other experimental campaigns will be held in 2012. This work is supported by the BioWin program of the Wallon Government.