Handcrafted firearm suicide: About the unusual use of a mole gun.

The body of a man was found in a workshop at the back of his home, a metal pipe embedded in his chest. In the immediate vicinity of the body, among various tools, was a machine containing a holster hit of lead shot evoking a mole-trap-type "taupe gun" (mole gun). The findings at the scene and the presence of a metal tube, sinkers, and a flock in the chest are in favour of the use of a homemade firearm. The autopsy and CT scan revealed atypical ballistic thoracic trauma responsible for a dilaceration of both the heart and thoracic aorta that caused the death. The autopsy also found a contused wound at the back of the skull without intracranial lesion, which may be due to the fall of the victim after the thoracic trauma. There is no other traumatic lesion and in particular no lesion of seizure, constraint, or defense. Toxicological analyses revealed the presence of alcoholic impregnation at the time of death. These findings, confronted with the expertise of a balistician, made it possible to understand how the victim used a mole gun to cause these lesions. A metal tube was used as a cannon to aim more surely at the heart but also as a means of triggering the propulsion mechanism of the trap without having to stretch out the arm. The tube, thinner than the barrel, had disengaged from its tip to be embedded with the lead shot in the chest of the victim.

[Duration of seclusion and exit of the seclusion room, study of the impact of anamnestic and institutional factors]. / Durée d'isolement et levée de la mesure ; étude de l'influence des facteurs anamnestiques et institutionnels.

INTRODUCTION: In France, the use of seclusion in psychiatric hospitals is regulated by the Act of January 26th 2016 which enforces a reduction of the use of coercive measures and limiting their duration. Criteria that are unrelated to the patient's symptoms might affect this decision and extend the duration of seclusion. The goal of the current observational study is, firstly, to determine which factors - unrelated to the patient's symptoms - influence the length of stay in seclusion. Secondly, it is to compare the composition of the medical and nursing teams at two times, the beginning and the end of the seclusion time period. METHOD: We conducted this study in a La Rochelle regional hospital from October 2017 to July 2018. There were one hundred and twenty-four episodes of seclusion that occurred in the different psychiatric inpatient wards during this time. The episodes were divided into two groups: short-term and long-term seclusion, defined by the median duration of seclusion. Data were collected using a survey, completed by the nursing teams, based on the recommendations of good clinical practice published in 1998 and revised in 2017 by the French health authority. For each episode of seclusion, we collected the following data: socio-demographic information, history of psychiatric care of the patient and of violent acting-out, admission status, medical indication for seclusion, day of the week and time at the beginning and the end of seclusion, prior administration of a sedative before seclusion, exceptional events that might justify the end of seclusion, transfer to a protected room after seclusion, as well as the composition of the medical and nursing team on duty at the beginning and the end of the seclusion period. We compared the anamnesis between the short-term and the long-term seclusion groups, and we analyzed the composition of the medical and nursing teams at the two time points of seclusion. Statistical analyses were performed using R software (v. 3.5.1). RESULTS: The mean duration of seclusion was 4.7 days and the median was 1.9 days. The average age was 37-years-old, with a ratio of 0.6 for females to males. Variables associated with a long-duration of seclusion were: the male gender (P=0.005), Compulsory Admission at the Request of a State Representative (P=0.008), a prevention measure of any hetero-aggressive action (P=0.007), the lack of psychiatric care (P=0.004), previous medico-legal issues (P=0.006), violent behavior during a previous hospitalization (P=0.022) and the use of seclusion on the weekend (P=0.01). The composition of the medical and nursing teams related to the period of the end of seclusion were: the presence of the physician in charge of the patient (P<0.001), a male caregiver in the team (P<0.001), a specialized psychiatric nurse (P<0.001) and the training of caregivers in the management of violence (P<0.001). The presence of nurses who do not usually work in the psychiatric ward was associated with the period of the beginning of seclusion (P<0.001). CONCLUSION: Our findings showed a strong relationship between several anamnestic factors and the duration of seclusion. Caregivers lacking information about patients, potential violent behaviors and the beginnings of seclusion on weekends are associated with a long-duration of seclusion. Our study also highlights the roles played by the caregivers according to their composition and level of training to determine the ending of patient's seclusion.

Comparison of IMRT and Rapidarc treatment plans using AAPM task group test suites.

The purpose of this study is to examine the plan quality and monitor unit with sliding window IMRT and RapidArc (RA) treatment plans using American Association Physicists in Medicine TG119 test suite DICOM-RT images and structure sets. The structure set includes multi-target (superior, central, inferior), prostate, head and neck and C-shape. Plans were performed with Eclipse planning system using AAA algorithm with the plan goals specified in TG119. The plan results for multitarget shows that the D99 is greater than the plan goal for all the targets. The D10 is less than the plan goal for superior and inferior targets in both IMRT and RA plans. The D10 is 5% more with IMRT plan and 7% more with RA plan for central target in comparison with plan goal. The plan results for prostate shows that D95 is greater than the plan goal for both IMRT and RA plans. The D5 is less than the plan goal for IMRT plan and almost equal to plan goal for RA plan. The D30 is less than the plan goal for bladder and rectum in both the plans. The D10 is higher than the plan goal by 1.9% and 2.5% in IMRT and RA plan for rectum. The plan results for head and neck shows that the D99 and D90 were greater than the plan goal for PTV. The spinal cord and parotid doses were less than the plan goal in both the plans. The plan results for C-shape shows that the D95 was greater than the plan goal and D10 was less than the plan goal for PTV. The dose to central core was less than the plan goal in both IMRT and RA plans. Both the IMRT and RapidArc plans have met the plan goal for all the target and normal structures. RapidArc optimization and treatment planning requires more time than the IMRT plan. The monitor unit calculated by the RapidArc plan is less compared to IMRT plan, which reduces the treatment error caused by patient motion during treatment and integral dose.

A quality assurance phantom for IMRT dose verification.

This paper investigates a quality assurance (QA) phantom specially designed to verify the accuracy of dose distributions and monitor units (MU) calculated by clinical treatment planning optimization systems and by the Monte Carlo method for intensity-modulated radiotherapy (IMRT). The QA phantom is a PMMA cylinder of 30 cm diameter and 40 cm length with various bone and lung inserts. A procedure (and formalism) has been developed to measure the absolute dose to water in the PMMA phantom. Another cylindrical phantom of the same dimensions, but made of water, was used to confirm the results obtained with the PMMA phantom. The PMMA phantom was irradiated by 4, 6 and 15 MV photon beams and the dose was measured using an ionization chamber and compared to the results calculated by a commercial inverse planning system (CORVUS, NOMOS, Sewickley, PA) and by the Monte Carlo method. The results show that the dose distributions calculated by both CORVUS and Monte Carlo agreed to within 2% of dose maximum with measured results in the uniform PMMA phantom for both open and intensity-modulated fields. Similar agreement was obtained between Monte Carlo calculations and measured results with the bone and lung heterogeneity inside the PMMA phantom while the CORVUS results were 4% different. The QA phantom has been integrated as a routine QA procedure for the patient's IMRT dose verification at Stanford since 1999.

Role of beam orientation optimization in intensity-modulated radiation therapy.

PURPOSE: To investigate the role of beam orientation optimization in intensity-modulated radiation therapy (IMRT) and to examine the potential benefits of noncoplanar intensity-modulated beams. METHODS AND MATERIALS: A beam orientation optimization algorithm was implemented. For this purpose, system variables were divided into two groups: beam position (gantry and table angles) and beam profile (beamlet weights). Simulated annealing was used for beam orientation optimization and the simultaneous iterative inverse treatment planning algorithm (SIITP) for beam intensity profile optimization. Three clinical cases were studied: a localized prostate cancer, a nasopharyngeal cancer, and a paraspinal tumor. Nine fields were used for all treatments. For each case, 3 types of treatment plan optimization were performed: (1) beam intensity profiles were optimized for 9 equiangular spaced coplanar beams; (2) orientations and intensity profiles were optimized for 9 coplanar beams; (3) orientations and intensity profiles were optimized for 9 noncoplanar beams. RESULTS: For the localized prostate case, all 3 types of optimization described above resulted in dose distributions of a similar quality. For the nasopharynx case, optimized noncoplanar beams provided a significant gain in the gross tumor volume coverage. For the paraspinal case, orientation optimization using noncoplanar beams resulted in better kidney sparing and improved gross tumor volume coverage. CONCLUSION: The sensitivity of an IMRT treatment plan with respect to the selection of beam orientations varies from site to site. For some cases, the choice of beam orientations is important even when the number of beams is as large as 9. Noncoplanar beams provide an additional degree of freedom for IMRT treatment optimization and may allow for notable improvement in the quality of some complicated plans.

Modeling the extrafocal radiation and monitor chamber backscatter for photon beam dose calculation.

A simple analytical approach has been developed to model extrafocal radiation and monitor chamber backscatter for clinical photon beams. Model parameters for both the extrafocal source and monitor chamber backscatter are determined simultaneously using conventional measured data, i.e., in-air output factors for square and rectangular fields defined by the photon jaws. The model has been applied to 6 MV and 15 MV photon beams produced by a Varian Clinac 2300C/D accelerator. Contributions to the in-air output factor from the extrafocal radiation and monitor chamber backscatter, as predicted by the model, are in good agreement with the measurements. The model can be used to calculate the in-air output factors analytically, with an accuracy of 0.2% for symmetric or asymmetric rectangular fields defined by jaws when the calculation point is at the isocenter and 0.5% when the calculation point is at an extended SSD. For MLC-defined fields, with the jaws at the recommended positions, calculated in-air output factors agree with the measured data to within 0.3% at the isocenter and 0.7% at off-axis positions. The model has been incorporated into a Monte Carlo dose algorithm to calculate the absolute dose distributions in patients or phantoms. For three MLC-defined irregular fields (triangle shape, C-shape, and L-shape), the calculations agree with the measurements to about 1% even for points at off-axis positions. The model will be particularly useful for IMRT dose calculations because it accurately predicts beam output and penumbra dose.

Verification of IMRT dose distributions using a water beam imaging system.

A water beam imaging system (WBIS) has been developed and used to verify dose distributions for intensity modulated radiotherapy using dynamic multileaf collimator. This system consisted of a water container, a scintillator screen, a charge-coupled device camera, and a portable personal computer. The scintillation image was captured by the camera. The pixel value in this image indicated the dose value in the scintillation screen. Images of radiation fields of known spatial distributions were used to calibrate the device. The verification was performed by comparing the image acquired from the measurement with a dose distribution from the IMRT plan. Because of light scattering in the scintillator screen, the image was blurred. A correction for this was developed by recognizing that the blur function could be fitted to a multiple Gaussian. The blur function was computed using the measured image of a 10 cm x 10 cm x-ray beam and the result of the dose distribution calculated using the Monte Carlo method. Based on the blur function derived using this method, an iterative reconstruction algorithm was applied to recover the dose distribution for an IMRT plan from the measured WBIS image. The reconstructed dose distribution was compared with Monte Carlo simulation result. Reasonable agreement was obtained from the comparison. The proposed approach makes it possible to carry out a real-time comparison of the dose distribution in a transverse plane between the measurement and the reference when we do an IMRT dose verification.

A dose-volume histogram based optimization algorithm for ultrasound guided prostate implants.

The task of treatment planning for prostate implants is to find an optimal seed configuration, comprising the target coverage and dosimetric consideration of critical structures such as the rectum and urethra. An efficient method to accomplish this is to use an inverse planning technique that derives the optimized solution from a prescribed treatment goal. The goal can be specified in the voxel domain as the desired doses to the voxels of the target and critical structures, or in the dose volume representation as the desired dose volume histograms (DVHs) of the target and critical structures. The DVH based optimization has been successfully used in plan optimization for intensity-modulated radiation therapy (IMRT) but little attention has been paid to its application in prostate implants. Clinically, it has long been known that some normal structure tolerances are more accurately assessed by volumetric information. Dose-volume histograms are also widely used for plan evaluation. When working in the DVH domain for optimization one has more control over the final DVHs. We have constructed an objective function sensitive to the DVHs of the target and critical structures. The objective function is minimized using an iterative algorithm, starting from a randomly selected initial seed configuration. At each iteration step, a trial position is given to a randomly selected source and the trial position is accepted if the objective function is decreased. To avoid being trapped in a less optimal local minimum, the optimization process is repeated. The final plan is selected from a pool of optimized plans obtained from a series of randomized initial seed configurations.

Monte Carlo verification of IMRT dose distributions from a commercial treatment planning optimization system.

The purpose of this work was to use Monte Carlo simulations to verify the accuracy of the dose distributions from a commercial treatment planning optimization system (Corvus, Nomos Corp., Sewickley, PA) for intensity-modulated radiotherapy (IMRT). A Monte Carlo treatment planning system has been implemented clinically to improve and verify the accuracy of radiotherapy dose calculations. Further modifications to the system were made to compute the dose in a patient for multiple fixed-gantry IMRT fields. The dose distributions in the experimental phantoms and in the patients were calculated and used to verify the optimized treatment plans generated by the Corvus system. The Monte Carlo calculated IMRT dose distributions agreed with the measurements to within 2% of the maximum dose for all the beam energies and field sizes for both the homogeneous and heterogeneous phantoms. The dose distributions predicted by the Corvus system, which employs a finite-size pencil beam (FSPB) algorithm, agreed with the Monte Carlo simulations and measurements to within 4% in a cylindrical water phantom with various hypothetical target shapes. Discrepancies of more than 5% (relative to the prescribed target dose) in the target region and over 20% in the critical structures were found in some IMRT patient calculations. The FSPB algorithm as implemented in the Corvus system is adequate for homogeneous phantoms (such as prostate) but may result in significant under or over-estimation of the dose in some cases involving heterogeneities such as the air-tissue, lung-tissue and tissue-bone interfaces.

Calculation of x-ray transmission through a multileaf collimator.

A ray tracing based method has been developed to calculate the x-ray transmission through a multileaf collimator (MLC) for beam delivery verification and dose calculation in intensity modulated radiotherapy (IMRT). The path length of a ray line in the MLC is accurately calculated using the exact geometry of the MLC leaves. The fluence distribution of an IMRT field is calculated first using a point source. The fluence distribution for a realistic beam model is obtained, as an approximation, by convolving the point source fluence distribution with the distribution of source strength. Full ray tracing calculations are performed using analytic and Monte Carlo simulated beam models to verify the accuracy of the convolution method. The calculation is in better agreement with measurements using either film or a beam imaging system (BIS) than previous calculations for MLC transmission using a simplified model. This ray tracing calculation can be applied to the problem of verifying dynamic MLC leaf sequences as part of a patient-specific quality assurance process for IMRT.

Energy- and intensity-modulated electron beams for radiotherapy.

This work investigates the feasibility of optimizing energy- and intensity-modulated electron beams for radiation therapy. A multileaf collimator (MLC) specially designed for modulated electron radiotherapy (MERT) was investigated both experimentally and by Monte Carlo simulations. An inverse-planning system based on Monte Carlo dose calculations was developed to optimize electron beam energy and intensity to achieve dose conformity for target volumes near the surface. The results showed that an MLC with 5 mm leaf widths could produce complex field shapes for MERT. Electron intra- and inter-leaf leakage had negligible effects on the dose distributions delivered with the MLC, even at shallow depths. Focused leaf ends reduced the electron scattering contributions to the dose compared with straight leaf ends. As anticipated, moving the MLC position toward the patient surface reduced the penumbra significantly. There were significant differences in the beamlet distributions calculated by an analytic 3-D pencil beam algorithm and the Monte Carlo method. The Monte Carlo calculated beamlet distributions were essential to the accuracy of the MERT dose distribution in cases involving large air gaps, oblique incidence and heterogeneous treatment targets (at the tissue-bone and bone-lung interfaces). To demonstrate the potential of MERT for target dose coverage and normal tissue sparing for treatment of superficial targets, treatment plans for a hypothetical treatment were compared using photon beams and MERT.

Beam orientation optimization in intensity-modulated radiation treatment planning.

Beam direction optimization is an important problem in radiation therapy. In intensity modulated radiation therapy (IMRT), the difficulty for computer optimization of the beam directions arises from the fact that they are coupled with the intensity profiles of the incident beams. In order to obtain the optimal incident beam directions using iterative or stochastic methods, the beam profiles ought to be optimized after every change of beam configuration. In this paper we report an effective algorithm to optimize gantry angles for IMRT. In our calculation the gantry angles and the beam profiles (beamlet weights) were treated as two separate groups of variables. The gantry angles were sampled according to a simulated annealing algorithm. For each sampled beam configuration, beam profile calculation was done using a fast filtered backprojection (FBP) method. Simulated annealing was also used for beam profile optimization to examine the performance of the FBP for beam orientation optimization. Relative importance factors were incorporated into the objective function to control the relative importance of the target and the sensitive structures. Minimization of the objective function resulted in the best possible beam orientations and beam profiles judged by the given objective function. The algorithm was applied to several model problems and the results showed that the approach has potential for IMRT applications.

Breast-conserving radiation therapy using combined electron and intensity-modulated radiotherapy technique.

BACKGROUND AND PURPOSE: To explore the feasibility of a multi-modality breast-conserving radiation therapy treatment technique to reduce high dose to the ipsilateral lung and the heart when compared with the conventional treatment technique using two tangential fields. MATERIALS AND METHODS: An electron beam with appropriate energy was combined with four intensity modulated photon beams. The direction of the electron beam was chosen to be tilted 10-20 degrees laterally from the anteroposterior direction. Two of the intensity-modulated photon beams had the same gantry angles as the conventional tangential fields, whereas the other two beams were rotated 15-25 degrees toward the anteroposterior directions from the first two photon beams. An iterative algorithm was developed which optimizes the weight of the electron beam as well as the fluence profiles of the photon beams for a given patient. Two breast cancer patients with early-stage breast tumors were planned with the new technique and the results were compared with those from 3D planning using tangential fields as well as 9-field intensity-modulated radiotherapy (IMRT) techniques. RESULTS: The combined electron and IMRT plans showed better dose conformity to the target with significantly reduced dose to the ipsilateral lung and, in the case of the left-breast patient, reduced dose to the heart, than the tangential field plans. In both the right-sided and left-sided breast plans, the dose to other normal structures was similar to that from conventional plans and was much smaller than that from the 9-field IMRT plans. The optimized electron beam provided between 70 to 80% of the prescribed dose at the depth of maximum dose of the electron beam. CONCLUSIONS: The combined electron and IMRT technique showed improvement over the conventional treatment technique using tangential fields with reduced dose to the ipsilateral lung and the heart. The customized beam directions of the four IMRT fields also kept the dose to other critical structures to a minimum.

Dosimetric effects of patient displacement and collimator and gantry angle misalignment on intensity modulated radiation therapy.

PURPOSE AND OBJECTIVE: The primary goal of this study was to examine systematically the dosimetric effect of small patient movements and linear accelerator angular setting misalignments in the delivery of intensity modulated radiation therapy. We will also provide a method for estimating dosimetric errors for an arbitrary combination of these uncertainties. MATERIALS AND METHODS: Sites in two patients (lumbar-vertebra and nasopharynx) were studied. Optimized intensity modulated radiation therapy treatment plans were computed for each patient using a commercially available inverse planning system (CORVUS, NOMOS Corporation, Sewickley, PA). The plans used nine coplanar beams. For each patient the dose distributions and relevant dosimetric quantities were calculated, including the maximum, minimum, and average doses in targets and sensitive structures. The corresponding dose volumetric information was recalculated by purposely varying the collimator angle or gantry angle of an incident beam while keeping other beams unchanged. Similar calculations were carried out by varying the couch indices in either horizontal or vertical directions. The intensity maps of all the beams were kept the same as those in the optimized plan. The change of a dosimetric quantity, Q, for a combination of collimator and gantry angle misalignments and patient displacements was estimated using Delta=Sigma(DeltaQ/Deltax(i))Deltax(i). Here DeltaQ is the variation of Q due to Deltax(i), which is the change of the i-th variable (collimator angle, gantry angle, or couch indices), and DeltaQ/Deltax(i) is a quantity equivalent to the partial derivative of the dosimetric quantity Q with respect to x(i). RESULTS: While the change in dosimetric quantities was case dependent, it was found that the results were much more sensitive to small changes in the couch indices than to changes in the accelerator angular setting. For instance, in the first example in the paper, a 3-mm movement of the couch in the anterior-posterior direction can cause a 38% decrease in the minimum target dose or a 41% increase in the maximum cord dose, whereas a 5 degrees change in the θ(1)=20 degrees beam only gave rise to a 1.5% decrease in the target minimum or 5.1% in the cord maximum. The effect of systematic positioning uncertainties of the machine settings was more serious than random uncertainties, which tended to smear out the errors in dose distributions. CONCLUSIONS: The dose distribution of an intensity modulated radiation therapy (IMRT) plan changes with patient displacement and angular misalignment in a complex way. A method was proposed to estimate dosimetric errors for an arbitrary combination of uncertainties in these quantities. While it is important to eliminate the angular misalignment, it was found that the couch indices (or patient positioning) played a much more important role. Accurate patient set-up and patient immobilization is crucial in order to take advantage fully of the technological advances of IMRT. In practice, a sensitivity check should be useful to foresee potential IMRT treatment complications and a warning should be given if the sensitivity exceeds an empirical value. Quality assurance action levels for a given plan can be established out of the sensitivity calculation.

Monitor unit calculation for an intensity modulated photon field by a simple scatter-summation algorithm.

An important issue in intensity modulated radiation therapy (IMRT) is the verification of the monitor unit (MU) calculation of the planning system using an independent procedure. Because of the intensity modulation and the dynamic nature of the delivery process, the problem becomes much more involved than that in conventional radiation therapy. In this work, a closed formula for MU calculation is derived. The approach is independent of the specific form of leaf sequence algorithms. It is straightforward to implement the procedure using a simple computer program. The approach is illustrated by a simplified example and is demonstrated by a few CORVUS (NOMOS Corporation, Sewickley, PA) treatment plans. The results indicate that it is robust and suitable for IMRT MU verification.

A dosimetric comparison of fan-beam intensity modulated radiotherapy with Gamma Knife stereotactic radiosurgery for treating intermediate intracranial lesions.

PURPOSE: To compare and evaluate treatment plans for the fan-beam intensity modulated radiotherapy and the Gamma Knife radiosurgery for treating medium-size intracranial lesions (range 4-25 cm3). METHODS AND MATERIALS: Treatment plans were developed for the Leksell Gamma Knife and a fan-beam inverse treatment planning system for intensity modulated radiotherapy. Treatment plan comparisons were carried out using dose-volume histogram (DVH), tissue-volume ratio (TVR), and maximum dose to the prescription dose (MDPD) ratio. The study was carried out for both simulated targets and clinical targets with irregular shapes and at different locations. RESULTS: The MDPD ratio was significantly greater for the Gamma Knife plans than for the fan-beam IMRT plans. The Gamma Knife plans produced equivalent TVR values to the fan-beam IMRT plans. Based on the DVH comparison, the fan-beam IMRT delivered significantly more dose to the normal brain tissue than the Gamma Knife. The results of the comparison were found to be insensitive to the target locations. CONCLUSION: The Gamma Knife is better than the fan-beam IMRT in sparing normal brain tissue while producing equivalent tumor dose conformity for treating medium-size intracranial lesions. However, the target dose homogeneity is significantly better for the fan-beam IMRT than for the Gamma Knife.

Estimation theory and model parameter selection for therapeutic treatment plan optimization.

Treatment optimization is usually formulated as an inverse problem, which starts with a prescribed dose distribution and obtains an optimized solution under the guidance of an objective function. The solution is a compromise between the conflicting requirements of the target and sensitive structures. In this paper, the treatment plan optimization is formulated as an estimation problem of a discrete and possibly nonconvex system. The concept of preference function is introduced. Instead of prescribing a dose to a structure (or a set of voxels), the approach prioritizes the doses with different preference levels and reduces the problem into selecting a solution with a suitable estimator. The preference function provides a foundation for statistical analysis of the system and allows us to apply various techniques developed in statistical analysis to plan optimization. It is shown that an optimization based on a quadratic objective function is a special case of the formalism. A general two-step method for using a computer to determine the values of the model parameters is proposed. The approach provides an efficient way to include prior knowledge into the optimization process. The method is illustrated using a simplified two-pixel system as well as two clinical cases. The generality of the approach, coupled with promising demonstrations, indicates that the method has broad implications for radiotherapy treatment plan optimization.

Matching photon and electron fields with dynamic intensity modulation.

A technique was developed to reduce the size and magnitude of the hot and cold spots in the abutting regions of photon and electron fields. The photon and electron fields were set up such that the photon field extended approximately 2 cm into the electron field in the abutting region. The region of the photon beam that overlapped the electron field was modulated using a multileaf collimator, effectively broadening the photon penumbra to make it complimentary to the electron penumbra. The computer calculations were verified using film measurements for abutting a 6 MV photon beam with a 9 MeV electron beam. A uniform dose was achieved at a prespecified depth of 2 cm, and dose uniformity was improved at the specified depth and beyond compared with unmodulated photon beams. A slight increase in dose inhomogeneity was seen at shallower depths. The overall areas of the hot and cold spots were significantly reduced. The technique also reduced the sensitivity of dose homogeneity to setup errors such that the magnitudes of the hot and cold spots were about half of those produced with unmodulated photon beam when an overlap or gap of 4 mm was introduced. The technique was applied to the treatment of a head and neck cancer and a lymphoma involving the right pleura with markedly reduced dose inhomogeneity in the abutting regions.

Clinical implementation of a Monte Carlo treatment planning system.

The purpose of this study was to implement the Monte Carlo method for clinical radiotherapy dose calculations. We used the EGS4/BEAM code to obtain the phase-space data for 6-20 MeV electron beams and 4, 6, and 15 MV photon beams for Varian Clinac 1800, 2100C, and 2300CD accelerators. A multiple-source model was used to reconstruct the phase-space data for both electron and photon beams, which retained the accuracy of the Monte Carlo beam data. The multiple-source model reduced the phase-space data storage requirement by a factor of 1000 and the accelerator simulation time by a factor of 10 or more. Agreement within 2% was achieved between the Monte Carlo calculations and measurements of the dose distributions in homogeneous and heterogeneous phantoms for various field sizes, source-surface distances, and beam modulations. The Monte Carlo calculated electron output factors were within 2% of the measured values for various treatment fields while the heterogeneity correction factors for various lung and bone phantoms were within 1% for photon beams and within 2% for electron beams. The EGS4/DOSXYZ Monte Carlo code was used for phantom and patient dose calculations. The results were compared to the dose distributions produced by a conventional treatment planning system and an intensity-modulated radiotherapy inverse-planning system. Significant differences (>5% in dose and >5 mm shift in isodose lines) were found between Monte Carlo calculations and the analytical calculations implemented in the commercial systems. Treatment sites showing the largest dose differences were for head and neck, lung, and breast cases.