D Sawkey's research while affiliated with Siemens Healthineers and other places

Publications (60)

Article
Full-text available
For the new Geant4 series 11.X, the electromagnetic (EM) physics sub-libraries were revised and reorganized in view of requirements for simulation of Phase-2 LHC experiments. EM physics simulation takes a significant fraction of the available CPU during massive production of Monte Carlo events for LHC experiments. We present the recent evolution of...
Preprint
Purpose: To develop a Monte Carlo (MC) therapeutic dose calculation toolkit of a recently released ring gantry linac in Geant4 (Version 10.7) for secondary dose validation of radiotherapy plan. Methods: For the Halcyon (Varian Medical Systems), the DSMLC was modeled and radiation transport in DSMLC and patient phantom was simulated using Geant4. Ra...
Article
Full-text available
The Geant4 electromagnetic (EM) physics sub-packages are a component of LHC experiment simulations. During long shutdown 2 for LHC, these packages are under intensive development and we report progress of EM physics in Geant4 versions 10.5 and 10.6, which includes faster computation, more accurate EM models, and extensions to the validation suite....
Article
Purpose: DEAR is a novel dynamic technique that achieves highly conformal dose through simultaneous couch and gantry motion during delivery. The purpose of this study is to develop a framework integrating a Monte Carlo dose engine (VirtuaLinac) to a treatment planning system (TPS, Eclipse) for DEAR. A quality assurance (QA) procedure is also devel...
Article
Full-text available
We report on developments of the Geant4 electromagnetic physics sub-libraries of Geant4 release 10.4 and beyond. Modifications are introduced to the models of photoelectric effect, bremsstrahlung, gamma conversion, single and multiple scattering. The theory-based Goudsmit-Saunderson model of electron/positron multiple scattering has been recently r...
Article
Purpose: A good Monte Carlo model with an accurate head shielding model is important in estimating the long term risks of unwanted radiation exposure during radiation therapy. The aim of this paper is to validate the Monte Carlo simulation of a TrueBeam linear accelerator (linac) head shielding model. We approach this by evaluating the accuracy of...
Article
Full-text available
We report on the recent progress within the Geant4 electromagnetic physics subpackages. Several new interfaces and models recently introduced are already used in LHC applications and may be useful for any type of simulation. Significant developments were carried out to improve the user interface, develop models of single and multiple scattering, an...
Article
Full-text available
Geant4 is a software toolkit for the simulation of the passage of particles through matter. It is used by a large number of experiments and projects in a variety of application domains, including high energy physics, astrophysics and space science, medical physics and radiation protection. Over the past several years, major changes have been made t...
Article
Purpose:Patient dose far from the treatment field is comprised of scatter from within the patient, and treatment head leakage. We quantify the treatment head leakage for TrueBeam linear accelerator for 6X and 6X-FFF beams by comparing measurements to Monte Carlo simulations for a variety of jaw sizes and collimator rotations. This work is conceptua...
Article
Purpose:DEAR is a radiation therapy technique utilizing synchronized motion of gantry and couch during delivery to optimize dose distribution homogeneity and penumbra for treatment of superficial disease. Dose calculation for DEAR is not yet supported by commercial TPSs. The purpose of this study is to demonstrate the feasibility of using a web-bas...
Article
Full-text available
In this work we report on recent improvements in the electromagnetic (EM) physics models of Geant4 and new validations of EM physics. Improvements have been made in models of the photoelectric effect, Compton scattering, gamma conversion to electron and muon pairs, fluctuations of energy loss, multiple scattering, synchrotron radiation, and high en...
Article
Full-text available
Purpose: Electronic portal imagers (EPIDs) with high detective quantum efficiencies (DQEs) are sought to facilitate the use of the megavoltage (MV) radiotherapy treatment beam for image guidance. Potential advantages include high quality (treatment) beam’s eye view imaging, and improved cone-beam computed tomography (CBCT) generating images with m...
Article
Calculate neutron production and transport in the TrueBeam treatment head, as input for vault design and phantom dose calculations. A detailed model of the treatment head, including shielding components off the beam axis, was created from manufacturer's engineering drawings. Simulations were done with Geant4 for the 18X, 15X, 10X and 10FFF beams, t...
Article
Full-text available
To develop a framework for accurate electron Monte Carlo dose calculation. In this study, comprehensive validations of vendor provided electron beam phase space files for Varian TrueBeam Linacs against measurement data are presented. In this framework, the Monte Carlo generated phase space files were provided by the vendor and used as input to the...
Conference Paper
Purpose: To validate the physics parameters of a Monte Carlo model for patient plane leakage calculations on the 6MV Unique linac by comparing the simulations against IEC patient plane leakage measurements. The benchmarked model can further be used for shielding design optimization, to predict leakage in the proximity of intended treatment fields,...
Article
Full-text available
In this work we present recent progress in Geant4 electromagnetic physics modelling, with an emphasis on the new refinements for the processes of multiple and single scattering, ionisation, high energy muon interactions, and gamma induced processes. The future LHC upgrade to 13 TeV will bring new requirements regarding the quality of electromagneti...
Conference Paper
Full-text available
The most recent upgrades of the electromagnetic (EM) physics “standard” and “low energy” sub-libraries of the general purpose Geant4 Monte Carlo simulation toolkit are described. These upgrades are relevant to different application domains including high energy physics, medical physics and space science. Validation results are presented and discuss...
Article
Purpose: Dynamic radiotherapy involving electron beams such as Dynamic Electron Arc Radiotherapy (DEAR) requires accurate dose modelling of small field sizes, similar to the requirement of IMRT field on the small photon field. The current commercial electron Monte Carlo algorithms such as eMC v11 in Eclipse were developed for standard field sizes a...
Article
Purpose: For TrueBeam Monte Carlo simulations, Varian does not distribute linac head geometry and material compositions, instead providing a phase space file (PSF) for the users. The PSF has a finite number of particle histories and can have very large file size, yet still contains inherent statistical noises. The purpose of this study is to charac...
Article
Purpose To simulate and measure radiation backscattered into the monitor chamber of a TrueBeam linac; establish a rigorous framework for absolute dose calculations for TrueBeam Monte Carlo (MC) simulations through a novel approach, taking into account the backscattered radiation and the actual machine output during beam delivery; improve agreement...
Article
Purpose: The design of the linac head is different for TrueBeam than Clinac, and there are differences in measured dose distributions in water phantoms between TrueBeam and Clinac for electron beams. Therefore, MC models for Clinac may not be applied directly to the Truebeam linac. The purpose of this study is to validate a Monte Carlo (MC) dose ca...
Article
Full-text available
The focus of this work was the demonstration and validation of VirtuaLinac with clinical photon beams and to investigate the implementation of low-Z targets in a TrueBeam linear accelerator (Linac) using Monte Carlo modeling. VirtuaLinac, a cloud based web application utilizing Geant4 Monte Carlo code, was used to model the Linac treatment head com...
Article
Purpose: Lateral dimensions of the incident electron beam (spot sizes) are essential for accurate modelling of x-ray fields, especially for small field output factors and penumbrae. We present spot sizes measured on the Varian TrueBeam linac and propose an explanation for the observed shapes. Methods: Spot were measured on 3 TrueBeam linacs usin...
Article
Purpose: Field‐independent phase‐space files are an effective means to reduce the time to simulate particle transport through medical linac treatment head components. The purpose of this work is to: (a) introduce the methodology used to newly create the full TrueBeam phase‐space library; (b) generate ∼50GB of data for each of seven photon beams; (c...
Article
Purpose: Monte Carlo simulations of linear accelerator treatment heads have traditionally involved building the simulated treatment head from specifications provided by the manufacturer. Here, we make pre-compiled applications available through a web interface, with processors running in the cloud. Methods: Simulations are run on computers on Am...
Presentation
Full-text available
Purpose: The focus of this work was to investigate the implementation of low-Z targets in a TrueBeam linear accelerator (linac) using Monte Carlo simulations. Methods: Simulation of a TrueBeam linac was accomplished by using Varian Virtualinac. Virtualinac is a cloud based web application utilizing GEANT4 Monte Carlo code to model TrueBeam compo...
Article
Purpose: The detective quantum efficiency (DQE(f)) is a critical metric to evaluate detector performance. For diagnostic imagers operating at keV energies, there are clear standards from the International Electrotechnical Commission (IEC) using patient equivalent filters to appropriately harden the beam. Such standards do not exist for MeV imaging...
Article
Full-text available
An overview of the current status of electromagnetic physics (EM) of the Geant4 toolkit is presented. Recent improvements are focused on the performance of large scale production for LHC and on the precision of simulation results over a wide energy range. Significant efforts have been made to improve the accuracy without compromising of CPU speed f...
Article
Full-text available
Strategies for delivering radiation to a moving lesion each require a margin to compensate for uncertainties in treatment. These motion margins have been determined here by separating the total uncertainty into components. Probability density functions for the individual sources of uncertainty were calculated for ten motion traces obtained from the...
Article
Purpose: To validate the generic phase space files for Varian TrueBeam linac head simulations. Methods: The generic phase space files include the simulation results of 6MV, 10MV, 6MV FFF, and 10MV FFF (flattening-filter free) operating modes of TrueBeam for patient-independent linac head components. Using the generic phase space files as the radiat...
Article
Full-text available
Electron scattering algorithms in Geant4 versions 9.4 and 9.5 were benchmarked by comparing scattered distributions against previously measured values at 13 and 20 MeV, for low, intermediate, and high atomic number materials. Several scattering models were used: Versions 93 and 95 of the Urban model, with different step size limits near boundaries;...
Article
Full-text available
The assumption of cylindrical symmetry in radiotherapy accelerator models can pose a challenge for precise Monte Carlo modeling. This assumption makes it difficult to account for measured asymmetries in clinical dose distributions. We have performed a sensitivity study examining the effect of varying symmetric and asymmetric beam and geometric para...
Article
Purpose: To develop a patient‐independent IAEA‐compliant phase space database using accurate geometry representation of the treatment head components imported from CAD (Computer‐Aided‐Design) drawings and a physics list exploration for the specific materials of the linac treatment head. Methods: Geant4.9.4 was employed to simulate the TrueBeam 6X,...
Article
Purpose: Massive computing resources are needed to generate an extensive linac phase‐space (PHSP) database for different radiotherapy beams that could serve as a starting point for several user‐applications, such as patient‐dependent dose calculations or Monte‐Carlo‐based quality assurance procedures. The goal of this project is to disseminate our...
Article
Full-text available
Monte Carlo simulation can accurately calculate electron fluence at the patient surface and the resultant dose deposition if the initial source electron beam and linear accelerator treatment head geometry parameters are well characterized. A recent approach used large electron fields to extract these simulation parameters. This method took advantag...
Article
Purpose: To account for a fringe magnetic field from the bending magnet in Monte Carlo treatment head simulation. Method and Materials: A Tesla meter was used to map the magnetic field between the exit window and secondary scattering foil of a clinical linac for 6–21 MeV electron beams. The foil and monitor chamber were removed to access this regio...
Article
Purpose: Write general?purpose, modular code for phase space input and output in Geant4 including latch and npass. Method and Materials: Geant4.9.3 does not include the ability to read or write phase space files. These files are useful for dividing the simulation into parts, for example, into patient?independent and patient?specific parts. Stand?al...
Article
Purpose: Verify that Monte Carlo simulations of clinically?relevant x?ray fields based on open?field simulations agree with measured dose distributions. Incident beam energy and spot?sized were determined independently, previously; the only simulation parameters varied here were the jaw/MLC positions. Methods & Materials:Monte Carlo simulations usi...
Article
To use an imaging beam line (IBL) to obtain the first megavoltage cone-beam computed tomography (MV CBCT) images of patients with a low atomic number (Z) target, and to compare these images to those taken of the same patients with the 6 MV flattened beam from the treatment beam line (TBL). The IBL, which produces a 4.2 MV unflattened beam from a ca...
Article
To determine the properties of a megavoltage cone-beam CT system using the unflattened beam from a sintered diamond target at 4 and 6 MV. A sintered diamond target was used in place of a graphite target as part of an imaging beam line (an unflattened beam from a graphite target) installed on a linear accelerator. The diamond target, with a greater...
Chapter
This paper reports on improved simulation of small electron fields from a Siemens Primus linear accelerator. Accelerator simulation was performed using the Monte Carlo user code BEAMnrc. Source and geometry parameters used were obtained from previous measurements and simulation of the linear accelerator with maximum electron field size (40 × 40 cm2...
Article
Obtain an accurate simulation of the dose from the 6 and 18 MV x-ray beams from a Siemens Oncor linear accelerator by comparing simulation to measurement. Constrain the simulation by independently determining parameters of the treatment head and incident beam, in particular, the energy and spot size. Measurements were done with the treatment head i...
Article
The purposes of this study are to improve the accuracy of source and geometry parameters used in the simulation of large electron fields from a clinical linear accelerator and to evaluate improvement in the accuracy of the calculated dose distributions. The monitor chamber and scattering foils of a clinical machine not in clinical service were remo...
Article
Purpose: To validate the use of improved source and geometry parameters, derived from Monte Carlo simulation of a Siemens Primus accelerator at maximum field size (40 × 40 cm2), for clinically relevant small electron fields. Method and Materials: Measurements were performed on a Siemens Primus linear accelerator for electron energies 6–21 MeV with...
Article
Purpose : Measure and simulate dose distributions from a 6 and 18 MV x‐ray linac with open jaws, including asymmetry, with independent measurements of the source and geometry details. Materials and Methods : Energy of the incident electron beam was measured by removing the target and flattening filter and measuring the percent‐depth ionization curv...
Article
Purpose/Objective(s): Digital Tomosynthesis (DT) is an imaging technique for reconstructing multiple tomograms from 2D projections acquired over a limited gantry arc (20°–40°). This study evaluates the use of megavoltage cone‐beam (MVCB) DT images acquired with a high contrast imaging beam line for patient setup. The objective is to determine the o...
Article
Purpose. Compare images and spectra using a dual energy diamond x?ray target for megavoltage cone?beam CT (MVCBCT) with those of carbon and tungsten targets. Materials and Methods. A Siemens Oncor linac with megavoltage cone?beam CT was modified to use either a tungsten,carbon, or diamond target for imaging. The carbon target was used with an accel...
Article
Purpose: To validate the use of improved source and geometry parameters, derived from Monte Carlo simulation of a Siemens Primus accelerator at maximum field size (40 × 40 cm2), for clinically relevant small electron fields. Method and Materials: Measurements were performed on a Siemens Primus linear accelerator for electron energies 6–21 MeV with...
Article
Purpose: High atomic number (Z) materials such as dental fillings complicate the use of conventional CTimages and application of the superposition dose calculation algorithm. This work investigates the accuracy of current planning methods for patients with dental work and proposes to combine the strengths of two imaging modalities in order to obtai...
Article
Monte Carlo simulations of x-ray beams typically take parameters of the electron beam in the accelerating waveguide to be free parameters. In this paper, a methodology is proposed and implemented to determine the energy, spectral width, and beam divergence of the electron source. All treatment head components were removed from the beam path, leavin...
Article
Purpose: To facilitate the development and commissioning of radiotherapy electron beam models by studying the dosimetric response of a Varian Trilogy electron accelerator to variations of the geometric and electron beam parameters of the treatment head, including lateral offsets of linac components and lateral and angular offsets of the incident el...
Article
Purpose: Determine fluence, differential in energy and angle, of Siemens Oncor 6/18 MV x?ray beams, a necessary step towards an accurate, easily?commissionable beam model for Monte Carlo treatment planning.Method and Materials: Dose to water is measured under various conditions, including: various field sizes including the largest; flattening filte...
Article
Purpose: Design a retrofit to Siemens electron applicators meeting IEC leakage standards with minimal effect on treatment field. Method and Materials: Dose to air measured on a Siemens Oncor linac with electron applicators, both at the patient plane and outside the applicator enclosure, by scanning an ion chamber and with film. Material added to ap...
Article
Purpose: To facilitate the development and commissioning of radiotherapy accelerator beam models by studying the impact of the variation of geometric and electron beam parameters on photon dose and fluence output for a Siemens Primus accelerator. In particular, the impact of lateral offsets of linac components and lateral and angular offsets of the...
Article
Purpose: Determine the characteristics, especially spot size and position, of a linac's primary electron beam, required as input to high‐precision Monte Carlo calculations of the X‐ray beam. Method and Materials: A parallel slit beam spot camera is placed on the block tray and a diode (oriented with normal to flat surface perpendicular to the beam)...
Conference Paper
Methodology is presented with the goal to determine fluence to 1-2% for electron and x-ray beams used in external beam radiotherapy. Fluence was determined by combining accurate and detailed treatment head simulation of a modern clinical accelerator with accurate measurements done over a wide range of operating parameters. The treatment head was di...

Citations

... MaximumEpsilonStep Maximum relative tolerance of the position and momentum inaccuracy for the integration of a step standard_opt4 physics list in TOPAS, corresponding to G4EMStandardPhysics_option4 7 in Geant4, is used to define the electromagnetic processes, models and its parameters for all particles with settings which are appropriate for radiotherapy with photons and charged particle beams (Ivanchenko et al 2020). For photons, the pair production is handled by the BetheHeitler5D model, the Compton scattering by the G4LowEPComptonModel (<20 MeV) and both the photo-electric effect and the Rayleigh scattering by the Livermore models. ...
... Muon interactions with materials in Geant4 are simulated using the Wentzel-VI Combined Scattering Model for electromagnetic interactions. 28,29 The drift radius is the distance of the closest approach between muon tracks and the central anode wire as the muon travels through a drift tube. While muon tracks can be reconstructed from experimental drift radius measurements, secondary electrons (delta rays) produced by muons' interactions with electrons in the detector gas also deposit energy in the drift tube gas. ...
... Feng et al. used BEAMnrc and DOSXYZnrc codes to simulate the beam characteristics from small fields, based on the TrueBeam linear particle accelerator (LINAC) (Varian Medical Systems, Palo Alto, CA, USA) flattening-filter-free (FFF) mode [3]. The Geant4based MC model is used to evaluate the accuracy of out-of-field dose predictions at extended distances, which are comprised of scattering from within the patient and treatment head leakage, thus reflecting the accuracy of the head shielding model proposed by Wijsooriya et al. [12]. The MCNPX code was used by Chibani et al. to simulate three different photon beams from two machines: the Siemens 18 MV, Varian 15 MV (Varian Medical Systems) and Varian 18 MV (Varian Medical Systems) [7]. ...
... For our simulation, we used a GEANT4-based [17] framework where different configurations have been investigated as shown in figure 3. The standard emstandard with opt4 physics constructor was employed as it is suggested to be used for applications that require higher accuracy of electrons, hadrons and ion tracking. For more details on the implementation with a particular focus on the MS see [18,19] and references therein. The first simulation regarded the estimation of the MS effect in the case of a full Fe blast ( = 7.87 g/cm 3 ) furnace as the worst-case scenario (left side of figure 3). ...
... To understand what happens during these collisions, sophisticated simulations of the detectors inside the LHC are generated. The majority of these simulations are founded on Monte Carlo methods [2], which yield results of high quality, albeit with a significant computational expense. In response, researchers investigate strategies to optimize resource use, including the application of generative machine learning methods to lessen CPU consumption [3], [4]. ...
... Multiple general-purpose MC codes are currently available, such as EGS4 [4], MCNP [5], PENELOPE [6] and Geant4 [7]. Nowadays, Geant4, a software toolkit originally created for high energy physics applications, is considered by many the main reference in the field of medical physics [8]. It is the most common choice for nuclear medicine applications due to its accurate physics models [9]. ...
... The tortuous path of the energetic electron emitted by the gamma was modeled with a realistic inelastic mean free path in BGO at 511 keV in the order of a fraction of microns [19]. The optical photons produced by scintillation and Cerenkov emission from energetic electrons were simultaneously generated using the Livermore Model [28]. Optical photons transport was conducted using the LUT Davis model [15,18]. ...
... [23][24][25][26] Prior research indicates that a DQE(0) of at least 20% is needed for sufficient image quality without delivering excessive dose to the patient, with non-negligible dose concerns emerging due to daily imaging of patients for modern IGRT. 2,43 However, conventional linear accelerators typically acquire MV images using a cone-beam CT geometry and flat-panel electronic portal imaging device (EPID) with a much lower detective efficiency of 1% to 2%. 44 Thus, it is of interest to consider the applicability of our work in the context of conventional cone-beam CT. ...
... Rodrigues et al. previously validated the vendor-provided electron phase-space files for conventional electron therapy, with the MLCs fully retracted. 23 The number of particles in the vendor-provided phase space is not large enough for these simulations, so we assume rotational symmetry for the upper part of the treatment head (and hence the phase-space file), and each particle is randomly rotated about the central axis (CAX) of the beam prior to calculation. ...
... This approach allows for a precise and simultaneous evaluation of absorbed dose and average LETs, for a voxel size of arbitrary dimensions, with negligible dependence with the cut variations [18]. Electromagnetic and hadronic interactions were simulated with the G4EmStandardPhysics-option4 and the QGSP-BIC-HP physics list, respectively [19,20]. More details on the Monte Carlo simulations of absorbed dose and LET are given in [17]. ...