Biomedical subjects

F Verhaegen

Publications and source records attributed to F Verhaegen.

At least 19 recordsLinked to original sources

PubMed2006

Perturbation correction factors for the NACP-02 plane-parallel ionization chamber in water in high-energy electron beams.

Recent dosimetry protocols for clinical high-energy electron beams recommend measurements of absorbed dose-to-water with a plane-parallel or cylindrical ionization chamber. For well-guarded plane-parallel ionization chambers, the ionization chamber perturbation factor in water, p(Q), has a recommended value of unity in all protocols. This assumption was investigated in detail in this study for one of the recommended ionization chambers in the protocols: the Scanditronix NACP-02 plane-parallel ionization chamber. Monte Carlo (MC) simulations of the NACP-02 ionization chamber with the EGSnrc code were validated against backscatter experiments. MC simulations were then used to calculate p(wall), p(cav) and p(Q) perturbation factors and water-to-air Spencer-Attix stopping powers in 4-19 MeV electron beams of a calibration laboratory (NPL), and in 6-22 MeV clinical electron beams from a Varian CL2300 accelerator. Differences between calculated and the currently recommended (Burns et al 1996 Med. Phys. 23 383-8) stopping powers, water-to-air, were found to be limited to 0.9% at depths between the reference depth z(ref) and the depth where the dose has decreased to 50% of the maximum dose, R50. p(wall) was found to exceed unity by 2.3% in the 4 MeV NPL calibration beam at z(ref). For higher energy electron beams p(wall) decreased to a value of about 1%. Combined with a p(cav) about 1% below unity for all energies at z(ref), this was found to cause p(Q) to exceed unity significantly for all energies. In clinical electron beams all three perturbation factors were found to increase with depth. Our findings indicate that the perturbation factors have to be taken into account in calibration procedures and for clinical depth dose measurements with the NACP-02 ionization chamber.

Aluminum↗

PubMed2005

Electron beam treatment verification using measured and Monte Carlo predicted portal images.

Electron beam treatments may benefit from techniques to verify patient positioning and dose delivery. This is particularly so for complex techniques such as mixed photon and electron beam radiotherapy and electron beam modulated therapy. This study demonstrates that it is possible to use the bremsstrahlung photons in an electron beam from a dual scattering foil linear accelerator to obtain portal images of electron beam treatments. The possibility of using Monte Carlo (MC) simulations to predict the electron beam treatment portal images was explored. The MC code EGSnrc was used to model a Varian CL21EX linear accelerator (linac) and to characterize the bremsstrahlung photon production in the linac head. It was found that the main sources of photons in the electron beam are the scattering foils, the applicator and the beam-shaping cut-out. Images were acquired using the Varian CL21EX linac and the Varian aS500 electronic portal imager (EPI); four electron energies (6, 9, 12, 16 MeV), and different applicator and cut-out sizes were used. It was possible to acquire images with as little as 10.7 MU per image. The contrast, the contrast-to-noise ratio (CNR), the signal-to-noise ratio (SNR), the resolution and an estimate of the modulated transfer function (MTF) of the electron beam portal images were computed using a quality assurance (QA) phantom and were found to be comparable to those of a 6 MV photon beam. Images were also acquired using a Rando anthropomorphic phantom. MC simulations were used to model the aS500 EPID and to obtain predicted portal images of the QA and Rando phantom. The contrast in simulated and measured portal images agrees within +/-5% for both the QA and the Rando phantom. The measured and simulated images allow for a verification of the phantom positioning by making sure that the structure edges are well aligned. This study suggests that the Varian aS500 portal imager can be used to obtain patient portal images of electron beams in the scattering foil linacs.

Computer Simulation↗

PubMed2005

Head-and-neck IMRT treatments assessed with a Monte Carlo dose calculation engine.

IMRT is frequently used in the head-and-neck region, which contains materials of widely differing densities (soft tissue, bone, air-cavities). Conventional methods of dose computation for these complex, inhomogeneous IMRT cases involve significant approximations. In the present work, a methodology for the development, commissioning and implementation of a Monte Carlo (MC) dose calculation engine for intensity modulated radiotherapy (MC-IMRT) is proposed which can be used by radiotherapy centres interested in developing MC-IMRT capabilities for research or clinical evaluations. The method proposes three levels for developing, commissioning and maintaining a MC-IMRT dose calculation engine: (a) development of a MC model of the linear accelerator, (b) validation of MC model for IMRT and (c) periodic quality assurance (QA) of the MC-IMRT system. The first step, level (a), in developing an MC-IMRT system is to build a model of the linac that correctly predicts standard open field measurements for percentage depth-dose and off-axis ratios. Validation of MC-IMRT, level (b), can be performed in a rando phantom and in a homogeneous water equivalent phantom. Ultimately, periodic quality assurance of the MC-IMRT system is needed to verify the MC-IMRT dose calculation system, level (c). Once the MC-IMRT dose calculation system is commissioned it can be applied to more complex clinical IMRT treatments. The MC-IMRT system implemented at the Royal Marsden Hospital was used for IMRT calculations for a patient undergoing treatment for primary disease with nodal involvement in the head-and-neck region (primary treated to 65 Gy and nodes to 54 Gy), while sparing the spinal cord, brain stem and parotid glands. Preliminary MC results predict a decrease of approximately 1-2 Gy in the median dose of both the primary tumour and nodal volumes (compared with both pencil beam and collapsed cone). This is possibly due to the large air-cavity (the larynx of the patient) situated in the centre of the primary PTV and the approximations present in the dose calculation.

Algorithms↗

PubMed2004

The dose distribution of low dose rate Cs-137 in intracavitary brachytherapy: comparison of Monte Carlo simulation, treatment planning calculation and polymer gel measurement.

In this study, the dose distribution delivered by low dose rate Cs-137 brachytherapy sources was investigated using Monte Carlo (MC) techniques and polymer gel dosimetry. The results obtained were compared with a commercial treatment planning system (TPS). The 20 mm and the 30 mm diameter Selectron vaginal applicator set (Nucletron) were used for this study. A homogeneous and a heterogeneous-with an air cavity-polymer gel phantom was used to measure the dose distribution from these sources. The same geometrical set-up was used for the MC calculations. Beyond the applicator tip, differences in dose as large as 20% were found between the MC and TPS. This is attributed to the presence of stainless steel in the applicator and source set, which are not considered by the TPS calculations. Beyond the air cavity, differences in dose of around 5% were noted, due to the TPS assuming a homogeneous water medium. The polymer gel results were in good agreement with the MC calculations for all the cases investigated.

Body Burden↗

PubMed2004

Theoretical analysis of microdosimetric spectra and cluster formation for 103Pd and 125I photon emitters.

In this work we have compared 125I or 103Pd from a microdosimetric point of view. The photon spectra at different positions around the seeds have first been calculated using EGSnrc Monte Carlo (MC) code. These photon spectra are used as input for the event-by-event MC code TRION to calculate the microdosimetric lineal energy (y) distribution for each isotope. The microdosimetric dose average lineal energy, yD, calculated in a sphere of 1 microm is 3.5 keV microm(-1) for 125I and 4.0 keV microm(-1) for 103Pd, agreeing well with values reported in the literature. yD in a 1 microm sphere diminishes slightly with the distance from the seed for 103Pd. This is due to the spectral hardening caused by the presence of a gamma-ray of 357.5 keV in the initial spectrum of 103Pd. In parallel with the calculation of the microdosimetric spectra, we have analysed the distribution of the size of the energy deposition clusters generated by these low energy photons in structures of 2 and 10 nm of radius. Due to Compton interactions, the fraction of very low energy electrons (<5 keV) generated by 125I photons is 51%, whereas it is only 27% for 103Pd. As these electrons deposit their energy very locally, the pattern of energy depositions contains more clusters of a few nm of radius for 125I than for 103Pd; the mean cluster orders are respectively 3.3 and 3.0 for 10 nm clusters. This is in opposition with the prediction based on the microdosimetric spectrum and the parameter yD and could be of importance for the damage to the cells.

Algorithms↗

PubMed2004

The radial dose function of low-energy brachytherapy seeds in different solid phantoms: comparison between calculations with the EGSnrc and MCNP4C Monte Carlo codes and measurements.

The use of low-energy photon emitters for brachytherapy applications, as in the treatment of the prostate or of eye tumours, has drastically increased in the last few years. New seed models for 103Pd and 125I have recently been introduced. The American Association of Physicists in Medicine recommends that measurements are made to obtain the dose rate constant, the radial dose function and the anisotropy function. These results must then be compared with Monte Carlo calculations to finally obtain the dosimetric parameters in liquid water. We have used the results obtained during the characterization of the new InterSource (furnished by IBt, Seneffe, Belgium) palladium and iodine sources to compare two Monte Carlo codes against experiment for these low energies. The measurements have been performed in three different media: two solid water plastics, WT1 and RW1, and polymethylmetacrylate. The Monte Carlo calculations were made using two different codes: MCNP4C and EGSnrc. These codes use photon cross-section data of a different origin. Differences were observed between both sets of input data below 100 keV, especially for the photoelectric effect. We obtained differences in the radial dose functions calculated with each code, which can be explained by the difference between the input data. New cross-section data were then tested for both codes. The agreement between the calculations using these new libraries is excellent. The differences are within the statistical uncertainties of the calculations. These results were compared with the experimental data. A good agreement is reached for both isotopes and in the three phantoms when the measured values are corrected for the presence of the TLDs in the phantom.

Algorithms↗

PubMed2004

[Hepatic abscesses in childhood: retrospective study about 33 cases observed in New-Caledonia between 1985 and 2003].

UNLABELLED: Hepatic abscesses in childhood are rarely observed in Europe. The aim of this word was to study how to diagnose and how to treat an hepatic abscess. METHODS: Between 1985 and 2003, we recensed retrospectively 33 cases of hepatic abscesses hospitalised in the paediatric unit of Noumea. RESULTS: Children were mainly melanesians (79%), 7 years old on average, having abdominal pains, a clinical and biological infectious syndrome, and abscesses images on ultrasonography or computed tomography. The identified micro-organisms included Entamoeba histolytica in 30% (10 cases); Staphylococcus aureus in 15% (five cases), Staphylococcus coagulase negative in 6% (two cases), Streptococcus D in 3% (one case); Bartonella henselae in 9% (three cases); ascaris in 6% (two cases); Mycobacterium tuberculosis in 6% (two cases). In eight cases no bacteria was identified (24%) but the good evolution after antibiotics and the negative amoebic serology looked like pyogenic abscesses. Two abscesses were aspirated, two were drained, one child had a surgical intervention. There was no death. Following a mean duration of 1 month for antibiotics treatment, outcome was always favourable. CONCLUSION: Diagnosis of hepatic abscess can be difficult. Ultrasonography shows the abscess but not the causal agent. The amoebic serology is sensible, consequently, its negativity leads to evoke a pyogenic agent. Early antibiotic treatment against pyogenic, anaerobic bacteria, and Entamoeba histolytica is required. Hepatic abscesses in ascaridiosis, tuberculosis and cat-scratch disease are less frequently encountered. If diagnosis remains doubtful or clinical evolution worsens, or if abscess volume increases, a percutaneous aspiration or drainage is needed.

Child↗

PubMed2004

Microdosimetric analysis of various mammography spectra: lineal energy distributions and ionization cluster analysis.

In view of recent recommendations on the frequency and the starting age of mammography screening in healthy women, it is desirable to quantify the enhanced relative biological effectiveness (RBE) of mammography X rays compared to hard X rays. While there is little doubt that the former are more potent in inducing biological damage than the latter, the magnitude of the effect is still hotly debated in the literature. We used Monte Carlo simulations and track structure analysis in micrometer and nanometer volumes to investigate differences in distributions of lineal energy and ionization clusters for a range of mammography X-ray qualities. Dose-averaged lineal energies, (yD), in breast tissue for various mammography qualities were found to result in quality factors about 40% higher than unity. Among the various mammography qualities studied, the popular molybdenum/molybdenum target/filter combination was found to have the highest (yD) in 1-microm spheres (about 5.0 keV/microm near the entrance surface of breast tissue). In 10-nm radius spheres, the mean ionization cluster order was found to be about 35% higher in mammography X rays compared to 300 keV electrons (roughly representing 60Co or 192Ir photon radiation). In even smaller spheres (2 nm radius), no significant differences were observed for the mean ionization cluster order between mammography X rays and 300 keV electrons. We conclude that the potential of mammography X rays to induce biological damage is probably not much higher than a factor of two compared to hard X rays.

Breast↗

PubMed2003

An investigation of entrance surface dose calculations for diagnostic radiology using Monte Carlo simulations and radiotherapy dosimetry formalisms.

Our aim in this work was to investigate the methodology used in the determination of the entrance surface dose (ESD) in diagnostic radiology. In kV x-rays for low-energy photons (tube potential up to 160 kV, HVL: 1-8 mm Al), the ESD is based on the use of the ratio of mass-energy absorption coefficients and backscatter factors. A full simulation of the photon and electron transport in a kilovoltage x-ray unit, using the Monte Carlo code BEAM/EGS4, was performed to obtain an accurate beam phase space for use in dose calculation. The modelled phase space was experimentally validated for the beam qualities (measured HVL: 3.3 mm Al-2.2 mm Cu) and showed good agreement between calculated and measured HVLs, air kerma and relative dose distributions. We have computed the conversion factors from air kerma to water or soft tissue absorbed dose at the surface of a phantom for beam qualities (HVL: 3.3-8.35 mm Al). The same model was also used to calculate the ESD in water and in soft tissue for the low-energy photon range considered. The results show that the numerical differences between the air kerma and the water kerma based backscatter factors are insignificant. The same conclusion was reached for the (mu(en)/rho) ratios, for soft tissue to air, evaluated using either the primary photon spectra or the spectra at the surface of a phantom. Furthermore, the good agreement obtained for the computation of the conversion factors with a full BEAM/EGS4 model confirms the previous studies which are based on different sources for the spectral distribution and different beam geometries (pencil beam or point source assumptions). On the other hand, the ESD in water or soft tissue is well described either with the B(air) or the B(w) formalism. Conversion factors from air kerma to ESD in these media are proposed in this work for several beam qualities in diagnostic radiology.

Biophysical Phenomena↗

PubMed2003

Commissioning and quality assurance of the Pinnacle(3) radiotherapy treatment planning system for external beam photons.

The commissioning of a Pinnacle(3) treatment planning system is described. Four Elekta linear accelerators were commissioned for external beam photons. Measured data were used to derive parameter values for the Pinnacle(3) beam model by (1). fitting a Monte Carlo model of the accelerator head to measured data and then extracting the parameters for the Pinnacle(3) beam model, and by (2). using the auto-modelling facility within Pinnacle(3). Both of these methods yielded dose distributions in accord with published recommendations. A separate small-field beam model, customized for an in-house compact blocking system, was also created, which satisfied appropriate acceptance criteria for stereotactically guided conformal brain treatments. Inhomogeneous, oblique, asymmetrical and irregular fields were also assessed, with calculated and measured doses agreeing to within +/-3%. Dose-volume histogram calculation was found to be accurate to within +/-5% dose or volume for a grid size of 4 mm x 4 mm x 4 mm, with better accuracy being achieved for finer grids. Isocentric doses were compared between Pinnacle(3)'s collapsed cone convolution algorithm and the Bentley-Milan algorithm within the Target-2 treatment planning system. Dose differences were generally less than 3% in the dose prescribed, with larger values for breast plans, where the Pinnacle(3) algorithm calculated scatter more accurately. Pelvic and thoracic plans were also verified using an anthropomorphic phantom, with local dose differences between calculated and delivered dose of up to 8%, but mainly less than 3%, and with no systematic difference. Ionization chamber verifications using START and RT-01 trial procedures demonstrated differences between calculated and measured doses of less than 2%. Following satisfactory performance in the commissioning process, Pinnacle(3) has now been introduced into routine clinical use.

Algorithms↗

PubMed2002

Optimization of accelerator target and detector for portal imaging using Monte Carlo simulation and experiment.

Megavoltage portal images suffer from poor quality compared to those produced with kilovoltage x-rays. Several authors have shown that the image quality can be improved by modifying the linear accelerator to generate more low-energy photons. This work addresses the problem of using Monte Carlo simulation and experiment to optimize the beam and detector combination to maximize image quality for a given patient thickness. A simple model of the whole imaging chain was developed for investigation of the effect of the target parameters on the quality of the image. The optimum targets (6 mm thick aluminium and 1.6 mm copper) were installed in an Elekta SL25 accelerator. The first beam will be referred to as A16 and the second as Cu1.6. A tissue-equivalent contrast phantom was imaged with the 6 MV standard photon beam and the experimental beams with standard radiotherapy and mammography film/screen systems. The arrangement with a thin Al target/mammography system improved the contrast from 1.4 cm bone in 5 cm water to 19% compared with 2% for the standard arrangement of a thick, high-Z target/radiotherapy verification system. The linac/phantom/detector system was simulated with the BEAM/EGS4 Monte Carlo code. Contrast calculated from the predicted images was in good agreement with the experiment (to within 2.5%). The use of MC techniques to predict images accurately, taking into account the whole imaging system, is a powerful new method for portal imaging system design optimization.

Bone and Bones↗

PubMed2002

Microdosimetric characterisation of 28 kVp Mo/Mo, Rh/Rh, Rh/Al, W/Rh and Mo/Rh mammography X ray spectra.

Microdosimetric characteristics of 28 kVp mammography X ray spectra were studied for several target/added filtration combinations (Mo/Mo, Rh/Rh, Rh/Al, W/Rh, Mo/Rh). Monte Carlo techniques were used to model X ray production from mammography units and to calculate distributions of absorbed dose and energy imparted in breast tissue. The results show that the dose averaged lineal energy is about 5.0 keV.micron-1, about 25% higher than for general diagnostic imaging X ray spectra. Significant differences in lineal energy between the five X ray qualities were noted, with the highest value for the commonly used Mo/Mo combination. Spectral hardening with depth in the tissues causes a 5% decrease in lineal energy over 5 cm. No significant differences were found for the different tissue compositions.

Breast↗

PubMed2002

An investigation of energy spectrum and lineal energy variations in mega-voltage photon beams used for radiotherapy.

Megavoltage photon beams are routinely used for external-beam radiotherapy. Recently, new treatment modalities based on dynamic and intensity modulated (IM), beam delivery systems are increasingly used in clinical practice. The purpose of this work is to investigate the energy spectrum and microdosimetric features of these photon beams. A Monte Carlo technique was first used to simulate beam lines of medical accelerators and to compute photon fluence and spectrum per unit dose-to-water inside the irradiated medium. Subsequently, a track structure code was used to compute the lineal energy and its distribution in a 1 micron sphere based on the individual photon spectrum. Results showed that the low energy photon component varied significantly with field size and location within the field due to the presence of the scattered photons. The calculated dose-mean lineal energy ranged from 2.3 keV.micron-1 at a depth of 1.5 cm along the central axis of a 4 cm x 4 cm field, to 3.7 keV.micron-1 at a depth of 20 cm of the field edge in a 10 cm x 10 cm field. In the tested IM fields at a depth of 10 cm, ranged from 2.6 to 3.5 keV.micron-1, which was inversely related to the dose intensity in the field. The values for the clinical photon beams were also significantly greater than that of a reference 60Co beam (1.8 keV.micron-1 from the calculation). The beam quality factor was estimated to vary within 20% due to the change of the energy spectrum for the radiotherapy photon fields.

Cobalt Radioisotopes↗

PubMed2001

Incorporating dynamic collimator motion in Monte Carlo simulations: an application in modelling a dynamic wedge.

In radiation therapy, new treatment modalities employing dynamic collimation and intensity modulation increase the complexity of dose calculation because a new dimension, time, has to be incorporated into the traditional three-dimensional problem. In this work, we investigated two classes of sampling technique to incorporate dynamic collimator motion in Monte Carlo simulation. The methods were initially evaluated for modelling enhanced dynamic wedges (EDWs) from Varian accelerators (Varian Medical Systems, Palo Alto, USA). In the position-probability-sampling or PPS method, a cumulative probability distribution function (CPDF) was computed for the collimator position, which could then be sampled during simulations. In the static-component-simulation or SCS method, a dynamic field is approximated by multiple static fields in a step-shoot fashion. The weights of the particles or the number of particles simulated for each component field are computed from the probability distribution function (PDF) of the collimator position. The CPDF and PDF were computed from the segmented treatment tables (STTs) for the EDWs. An output correction factor had to be applied in this calculation to account for the backscattered radiation affecting monitor chamber readings. Comparison of the phase-space data from the PPS method (with the step-shoot motion) with those from the SCS method showed excellent agreement. The accuracy of the PPS method was further verified from the agreement between the measured and calculated dose distributions. Compared to the SCS method, the PPS method is more automated and efficient from an operational point of view. The principle of the PPS method can be extended to simulate other dynamic motions, and in particular, intensity-modulated beams using multileaf collimators.

Biophysical Phenomena↗

PubMed2001

Quantifying effects of lead shielding in electron beams: a Monte Carlo study.

Lead shielding in contact with the patient's skin is often encountered in radiotherapy with electron beams. The influence of the lead shielding on dose distributions in the patient cannot fully be assessed using modern treatment planning systems. In this work the problem of quantifying the effect of lead shielding on dose distributions is addressed. Monte Carlo dose calculations were performed in a half-blocked water phantom shielded by lead, using a realistic model for the fluence of an electron linear accelerator. Electron beam energies of 6-20 MeV and lead thicknesses of 1-7 mm are used for 10 x 10 cm2 and 5 x 5 cm2 fields. The perturbation of the particle fluence and dose distributions in water introduced by the lead shielding is quantified. The effect of oblique electron beams on the dose perturbation is shown. A fictitious clinical example, the shielding of an eye in electron beam treatment, is used to demonstrate the usefulness of Monte Carlo based treatment planning algorithms that can incorporate the effects of lead shielding.

Algorithms↗

PubMed2001

Experimental p(wall) and p(cel) correction factors for ionization chambers in low-energy clinical proton beams.

Current dosimetry protocols for clinical protons using ionization chambers do not take into account ionization chamber-dependent perturbation correction factors. In the present investigation, the relative response of 17 cylindrical ionization chambers was evaluated at three proton beam qualities: at two points in a modulated beam and one point in a non-modulated beam, all with an incident energy of 75 MeV. Thirteen of the ionization chambers had a Farmer-type geometry but consisted of different combinations of wall and central electrode materials. All ionization chambers were calibrated in terms of air kerma as well as in terms of absorbed dose to water in a 60Co beam. The relative response of the ionization chambers was compared with results of Monte Carlo simulations of proton and secondary electron transport in the phantom and the ionization chamber geometry. The results of the measurements for cylindrical ionization chambers show relative perturbation effects that are limited to 0.5-1%, resulting in perturbation correction factors that are larger than unity compared with an NE2571 ionization chamber. The experimental relative wall and total perturbation correction factors agree with Monte Carlo calculated values, indicating that the differences between the responses of different ionization chambers are due to secondary electron effects. This conclusion is supported by the comparison of our results with those from other investigators after re-analysis of their data. The central electrode perturbation correction factor for an aluminium electrode in a Farmer-type geometry was found to be unity within the experimental uncertainties.

Electrons↗

PubMed2001

A method of simulating dynamic multileaf collimators using Monte Carlo techniques for intensity-modulated radiation therapy.

A method of modelling the dynamic motion of multileaf collimators (MLCs) for intensity-modulated radiation therapy (IMRT) was developed and implemented into the Monte Carlo simulation. The simulation of the dynamic MLCs (DMLCs) was based on randomizing leaf positions during a simulation so that the number of particle histories being simulated for each possible leaf position was proportional to the monitor units delivered to that position. This approach was incorporated into an EGS4 Monte Carlo program, and was evaluated in simulating the DMLCs for Varian accelerators (Varian Medical Systems, Palo Alto. CA, USA). The MU index of each segment, which was specified in the DMLC-control data, was used to compute the cumulative probability distribution function (CPDF) for the leaf positions. This CPDF was then used to sample the leaf positions during a real-time simulation, which allowed for either the step-shoot or sweeping-leaf motion in the beam delivery. Dose intensity maps for IMRT fields were computed using the above Monte Carlo method, with its accuracy verified by film measurements. The DMLC simulation improved the operational efficiency by eliminating the need to simulate multiple segments individually. More importantly, the dynamic motion of the leaves could be simulated more faithfully by using the above leaf-position sampling technique in the Monte Carlo simulation.

Biophysical Phenomena↗

PubMed2001

Monte Carlo calculation of output factors for circular, rectangular, and square fields of electron accelerators (6-20 MeV).

Monte Carlo (MC) techniques can be used to build a simulation model of an electron accelerator to calculate output factors for electron fields. This can be useful during commissioning of electron beams from a linac and in clinical practice where irregular fields are also encountered. The Monte Carlo code BEAM/EGS4 was used to model electron beams (6-20 MeV) from a Varian 2100C linear accelerator. After optimization of the Monte Carlo simulation model, agreement within 1% to 2% was obtained between calculated and measured (with a Si diode) lateral and depth dose distributions or within 1 mm in the penumbral regions. Output factors for square, rectangular, and circular fields were measured using two different plane-parallel ion chambers (Markus and NACP) and compared to MC simulations. The agreement was usually within 1% to 2%. This study was not primarily concerned with minimizing the simulation time required to obtain output factors but some considerations with respect to this are presented. It would be particularly useful if the MC model could also be used to calculate output factors for other, similar linacs. To see if this was possible, the primary electron energies in the MC model were retuned to model a recently commissioned similar linac. Good agreement between calculated and measured output factors was obtained for most field sizes for this second accelerator.

Biophysical Phenomena↗