Off-axis x-ray spectra: a comparison of Monte Carlo simulated and computed x-ray spectra with measured spectra

Simulations of X-ray spectra, half value layer, and mean energy from mammography using EGSnrc Monte Carlo and SpekPy

INTRODUCTION

This study presents the simulation results of X-ray spectra, half value layers (HVLs), and mean energies (E_mean) of two mammography units using EGSnrc Monte Carlo (MC) and SpekPy computer codes.

METHODS

The spectra caused by different combinations of targets/filters at various tube voltages (kVps) of two mammography units were simulated using two different computer codes. The EGSnrc MC simulated data of spectra and E_mean were compared with those obtained from SpekPy. The simulated values of two units' HVLs obtained from two computer codes were compared with those from physical measurements from mammography machines used in clinical practice.

RESULTS

The maximum discrepancies in E_mean simulated from two codes were less than 4.1% and 1.5% for the target/filter combination of W/Rh and Mo/Mo, respectively. The HVLs of the SpekPy were well matched to the physical measurements. The percentage differences were within 3.3% and 6.8% for two units, respectively. The EGSnrc MC simulated values of HVLs show the percentage differences within 8.9% and 7.0% with those from physical measurements for two units, respectively. All methods of HVLs determination comply with the requirements of IAEA Human Health Series No.17.

CONCLUSIONS

The HVLs, E_mean, spectra varied depending on the target/filter combinations and composition of the mammography tubes. The simulation results verify that the HVLs evaluation using the EGSnrc MC and SpekPy can be validated for mammography standard beam qualities and provide prediction almost immediately compared with physical experiments.

IMPLICATIONS FOR PRACTICE

EGSnrc MC and SpekPy have been considered powerful toolkits to simulate typical X-ray tubes used in mammography due to the good agreement between the calculation of E_mean, physical measurements and simulated HVLs.

A model for the energy and angular distribution of x rays emitted from an x-ray tube. Part II. Validation of x-ray spectra from 20 to 300 kV

PURPOSE

To present and validate a complete x-ray emission model (bremsstrahlung and characteristic x-ray emission) for the energy range 20-300 kV.

METHODS

An analytical x-ray spectrum model that combines the bremsstrahlung emission model developed in Part I with a previously developed characteristic x-ray emission model is validated by comparison with Monte Carlo calculations, published measured spectra, and models developed by other authors. Furthermore, the assumptions and limitations of previous spectrum models are summarized, and their predictions are compared with results obtained by Monte Carlo simulations of x rays emitted from tungsten and molybdenum targets.

RESULTS

The model is able to reproduce narrow-beam Monte Carlo calculations to within 0.5% in terms of the first and second aluminum half-value layer thickness (HVL). Compared with measured spectra, the difference in HVL is < 2% for typical diagnostic and therapeutic beam qualities available at primary standard laboratories. Compared with previous spectrum models, the present model performs especially well for low kilovoltage x-ray beams (below 50 kV), and is reliable for a wider range of take-off angles, that is, the angle between the target surface and the direction of emission. The difference in model and Monte Carlo predictions of the energy-fluence weighted air kerma (i.e., the photon energy absorption in air) is < 0.5% using the present model, while previous spectrum models can differ by more than 10%.

CONCLUSIONS

The x-ray emission model developed in this work has been validated against Monte Carlo calculations and measured results. The model provides an efficient alternative to comprehensive Monte Carlo simulations and is an improvement over previous models. The model can be used to predict both central- and off-axis spectra, as well as off-axis effects such as the (anode) heel effect.

Validation of mammographic x-ray spectra generated using Particle and Heavy Ion Transport code System

A Monte Carlo (MC) code is a robust method to generate a mammographic x-ray spectrum because the geometry of a mammography system can be flexible and directly modeled in MC simulation. However, simulations from MC code need to be validated before it can be reliably used for specific applications. This study aimed to generate and validate the x-ray spectra of relevant anodes used in mammography and breast tomosynthesis using Particle and Heavy Ion Transport code System (PHITS). PHITS version 3.08 was used to generate the x-ray spectra of molybdenum (Mo), rhodium (Rh), and tungsten (W) anodes. The Mo anode spectrum derived using PHITS was compared with those obtained using other MC codes. The generated spectra of all anodes were compared with the literature. Parameters including spectral shape, K characteristic x-ray yield, heel effect, and half-value layer (HVL) were used for a comparative assessment. The differences in these assessment parameters conducted by PHITS and PHITSEGS5 simulations were studied. Regarding the comparative parameters, PHITSEGS5 simulation improved the accuracy of mammographic x-ray generation compared to PHITS simulation; K x-ray and bremsstrahlung yields of the Mo anode spectrum generated by PHITSEGS5 simulation were a better agreement with those generated by other MC code simulations. The PHITSEGS5 spectra overestimated K x-ray and low-energy bremsstrahlung photons in comparison with measured spectra. Subsequently, HVLs calculated from PHITSEGS5 spectra were 1.0% (Mo/Mo) and 7.0% (W/Al) lower than those derived from measured spectra. For Mo and Rh anodes, relative difference of HVLs calculated from PHITSEGS5 spectra and those obtained from literature and measurement were within the TRS 457 acceptance criteria (±0.02 mm Al). The observed difference exceeded the acceptance criteria for W anode. Regarding existed consistency in HVL between simulation and measurement, PHITSEGS5 simulation can be reliably used to generate x-ray spectra of Mo and Rh anodes. However, its accuracy should be improved for generating W anode spectrum.

VALIDATION OF A BEAMNRC MONTE CARLO SIMULATION OF A BROAD BEAM DIAGNOSTIC X-RAY UNIT

A BEAMnrc Monte Carlo simulation of a diagnostic X-ray unit in a broad beam geometry intended for the derivation of effective linear attenuation coefficients was validated against measurements made on the X-ray unit on which the simulation was based. The validation process assessed the size of the beam, the first and second half value layer, the variation in kerma with anode-heel effect and the accuracy of values of effective linear attenuation coefficient for water and solid water attenuators. The agreement between the simulated and measured results for all tests was consistently within the experimental uncertainty for the measurements. The average deviation between values of effective linear attenuation coefficient for simulated and measured results is 3.8%, with the highest individual deviation 7.1%. The simulation can be used to produce values of effective linear attenuation coefficient for a range of kVp, field size and attenuator thickness that are close to those measured.

Experimental measurements and Monte Carlo modelling of the XSTRAHL 150 superficial X-ray therapy unit

Background

Superficial X-ray therapy units are used for the treatment of certain types of skin cancer and some severe dermatological conditions. The performance assessment and beam characteristics of the superficial unit are very important to ensure accurate dose delivery during patient treatment. Both experimental measurements and Monte Carlo calculations can be used for this purpose.

Purpose

This study aims to investigate whether it is possible to reproduce experimentally measured data for the XSTRAHL 150 superficial X-ray unit with simulations using the BEAMnrc Monte Carlo code.

Materials and Methods

The experimental procedure applied in this study included the following: experimental measurements of different X-ray spectra, half-value layers, percentage depth dose and beam profiles. Monte Carlo modelling of the XSTRAHL 150 unit was performed with the BEAMnrc code. The validity of the model was checked by comparing the theoretical calculations with experimental measurements.

Results

There was good agreement (∼1%) between experimentally measured and simulated X-ray spectra. Results of half-value layers obtained from simulated and measured spectra showed that there was a maximum of 3·6% difference between BEAMnrc and measurements and a minimum of 2·3%. In addition, simulated percentage depth dose and profile curves have been compared against experimental measurements and show good agreement (within 2% for the depth dose curves and 3–5% for beam profile curves, depending on the applicator size).

Conclusion

The results of this study provide information about particles’ interaction in different kilovoltage and filter combinations. This information is useful for X-ray tube design and development of new target/filter combinations to improve beam quality in superficial X-ray radiotherapy. The data presented here may provide a base for comparison and a reference for other or potential new users of the XSTRAHL 150 X-ray unit.

Monte Carlo simulator of realistic x-ray beam for diagnostic applications

PURPOSE

Monte Carlo simulation is a very useful tool for radiotherapy and diagnostic radiology. Yet even with the latest PCs, simulation of photon spectra emitted by an x-ray tube is a time-consuming task, potentially reducing the possibility to obtain relevant data such as dose evaluations, simulation of geometric settings, or monitor detector efficiency. This study developed and validated a method to generate random numbers for realistic beams in terms of photon spectrum and intensity to simulate x-ray tubes via Monte Carlo algorithms.

METHODS

Starting from literature data, the most common semiempirical models of bremsstrahlung are analyzed and implemented, adjusting their formulation to describe a large irradiation area (i.e., large field of view) and to take account of the heel effect as in common practice during patient examinations.

RESULTS

Simulation results show that Birch and Marshall's model is the fastest and most accurate for the aims of this work. Correction of the geometric size of the beam and validation of the intensity variation (heel effect) yielded excellent results with differences between experimental and simulated data of less than 6%.

CONCLUSIONS

The results of validation and execution time showed that the tube simulator calculates the x-ray photons quickly and efficiently and is perfectly capable of considering all the phenomena occurring in a real beam (total filtration, focal spot size, and heel effect), so it can be used in a wide range of applications such as industry, medical physics, or quality assurance.

A comparison of mammographic x-ray spectra: simulation with EGSnrc and experiment with CdTe detector

Mammographic x-ray spectra simulated by BEAMnrc/EGSnrc Monte Carlo code were qualitatively compared with the results obtained from the direct measurement using a cadmium telluride x-ray spectroscopy system and from the generation of IPEM report number 78. Generally, there is good agreement between the simulated and measured spectra, though there are slight differences at low energy in which the K-characteristic x-ray intensity is relatively higher for IPEM spectra. In addition, transmission curves were measured and simulated using a breast tissue-equivalent phantom (BR-12) as filtration. Comparison of the transmission curves shows good agreement. Moreover, the first half value layer (HVL) from direct measurement using ion chamber was consistent with the first HVL calculated by simulated spectra. Therefore, Monte Carlo may be used as an alternative tool for obtaining x-ray spectra when direct measurement is not available.

DXRaySMCS: a user-friendly interface developed for prediction of diagnostic radiology X-ray spectra produced by Monte Carlo (MCNP-4C) simulation

In this work, the general purpose Monte Carlo N-particle radiation transport computer code (MCNP-4C) was used for the simulation of X-ray spectra in diagnostic radiology. The electron's path in the target was followed until its energy was reduced to 10 keV. A user-friendly interface named 'diagnostic X-ray spectra by Monte Carlo simulation (DXRaySMCS)' was developed to facilitate the application of MCNP-4C code for diagnostic radiology spectrum prediction. The program provides a user-friendly interface for: (i) modifying the MCNP input file, (ii) launching the MCNP program to simulate electron and photon transport and (iii) processing the MCNP output file to yield a summary of the results (relative photon number per energy bin). In this article, the development and characteristics of DXRaySMCS are outlined. As part of the validation process, output spectra for 46 diagnostic radiology system settings produced by DXRaySMCS were compared with the corresponding IPEM78. Generally, there is a good agreement between the two sets of spectra. No statistically significant differences have been observed between IPEM78 reported spectra and the simulated spectra generated in this study.

Dosimetric characterisation of bismuth shields in CT: measurements and Monte Carlo simulations

Although bismuth shields are frequently used in radiology to reduce radiation dose, its mechanism has not been fully investigated. Dosimetric characteristics of bismuth shields in computed tomography (CT) were assessed with ion chamber and Monte Carlo (MC) simulations. Primary attenuation and backscatter effects of paediatric (2-ply) and adult (4-ply) bismuth shields were measured. Simulated CT beams were used for ion chamber measurements. Radiation doses were measured free-in-air and in the tissue-equivalent slabs. MC simulations for the same settings were also performed. Mean dose reductions from primary attenuation were 23% (2-ply) and 40% (4-ply). The dose increase from backscatter was 2% for both shields. MC simulations for primary beam dose reduction were 20% (2-ply) and 38% (4-ply); the backscatter dose increase was around 6% for both shields. In summary, primary attenuation is the major factor that introduces the dose reduction in bismuth and the dose increase from backscatter is negligible.

Calculation of x-ray spectra emerging from an x-ray tube. Part II. X-ray production and filtration in x-ray targets

A new approach to the calculation of the x-ray spectrum emerging from an x-ray tube is proposed. Theoretical results for the bremsstrahlung cross section appearing in the literature are summarized. Four different treatments of electron penetration, based on the work presented in Part I, are then used to generate bremsstrahlung spectra. These spectra are compared to experimental data at 50, 80 and 100 kVp tube potentials. The most sophisticated treatment of electron penetration was required to obtain good agreement. With this treatment both the National Institute of Standards and Technology bremsstrahlung cross sections, based on accurate partial wave calculations, and the Bethe-Heitler cross section [H. A. Bethe and W. Heitler, Proc R. Soc. London, Ser. A. 146, 83-112 (1934)] corrected by a modified Elwert factor [G. Elwert, Ann. Phys. (Leipzig) 426, 178-208 (1939)], provided good agreement to measured data. An approximate treatment of the characteristic spectrum is suggested. The dependencies of the bremsstrahlung and characteristic outputs of an x-ray tube on tube potential are compared to experimentally derived data for 70-140 kVp potentials. Agreement is to within a few percent of the total output over the entire range. The spectral predictions of the semiempirical models of Birch and Marshall [R. Birch and M. Marshall, Phys. Med. Biol. 24, 505-513 (1979)] (IPEM Report 78) and of Tucker et al. [D. M. Tucker, G. T. Barnes, and D. P. Chakraborty, Med. Phys. 18, 211-218 (1991).] are also assessed. The predictions of Tucker et al. are very close to the model developed here. The predictions of IPEM Report 78 are similar, but consistently harder for the range of tube potentials examined (50-100 kV). Unlike the semiempirical models, the model proposed here requires the introduction of no empirical and unphysical parameters in the differential bremsstrahlung cross section, bar an overall normalization factor which is close to unity.

Efficiency improvements of x-ray simulations in EGSnrc user-codes using bremsstrahlung cross-section enhancement (BCSE)

This paper presents the implementation of the bremsstrahlung cross-section enhancement (BCSE) variance-reduction technique into the EGSnrc/BEAMnrc system. BCSE makes the simulation of x-ray production from bremsstrahlung targets more efficient; it does so by artificially making the rare event of bremsstrahlung emission more abundant, which increases the number of statistically-independent photons that contribute to reducing the variance of the quantity of interest without increasing the CPU time appreciably. BCSE does not perturb the charged-particle transport in EGSnrc and it is made compatible with all other variance-reduction techniques already used in EGSnrc and BEAMnrc, including range rejection, uniform bremsstrahlung splitting, and directional bremsstrahlung splitting. When optimally combining BCSE with splitting to simulate typical situations of interest in medical physics research and in clinical practice, efficiencies can be up to five orders of magnitude larger than those obtained with analog simulations, and up to a full order of magnitude larger than those obtained with optimized splitting alone (which is the state-of-the-art of the EGSnrc/BEAMnrc system before this study was carried out). This study recommends that BCSE be combined with the existing splitting techniques for all EGSnrc/BEAMnrc simulations that involve bremsstrahlung targets, both in the kilovoltage and megavoltage range. Optimum crosssection enhancement factors for typical situations in diagnostic x-ray imaging and in radiotherapy are recommended, along with an easy algorithm for simulation optimization.

Dose delivered from Varian's CBCT to patients receiving IMRT for prostate cancer

With the increased use of cone beam CT (CBCT) for daily patient setup, the accumulated dose from CBCT may be significantly higher than that from simulation CT or portal imaging. The objective of this work is to measure the dose from daily pelvic scans with fixed technical settings and collimations. CBCT scans were acquired in half-fan mode using a half bowtie and x-rays were delivered in pulsed-fluoro mode. The skin doses for seven prostate patients were measured on an IRB-approved protocol. TLD capsules were placed on the patient's skin at the central axis of three beams: AP, left lateral (Lt Lat) and right lateral (Rt Lat). To avoid the ring artefacts centred in the prostate, the treatment couch was dropped 3 cm from the patient's tattoo (central axis). The measured AP skin doses ranged 3-6 cGy for 20-33 cm separation. The larger the patient size the less the AP skin dose. Lateral doses did not change much with patient size. The Lt Lat dose was approximately 4.0 cGy, which was approximately 40% higher than the Rt Lat dose of approximately 2.6 cGy. To verify this dose asymmetry, surface doses on an IMRT QA phantom (oval shaped, 30 cm x 20 cm) were measured at the same three sites using TLD capsules with 3 cm table-drop. The dose asymmetry was due to: (1) kV source rotation which always starts from the patient's Lt Lat and ends at Lt Lat. Gantry rotation gets much slower near the end of rotation but dose rate stays constant and (2) 370 degrees scan rotation (10 degrees scan overlap on the Lt Lat side). In vivo doses were measured inside a Rando pelvic heterogeneous phantom using TLDs. The left hip (femoral head and neck) received the highest doses of approximately 10-11 cGy while the right hip received approximately 6-7 cGy. The surface and in vivo doses were also measured for phantoms at the central-axis setup. The difference was less than approximately 12% to the table-drop setup.

Generation of calibrated tungsten target x-ray spectra: modified TBC model

In spite of the recent advances in the experimental detection of x-ray spectra, theoretical or semi-empirical approaches for determining realistic x-ray spectra in the range of diagnostic energies are important tools for planning experiments, estimating radiation doses in patients, and formulating radiation shielding models. The TBC model is one of the most useful approaches since it allows for straightforward computer implementation, and it is able to accurately reproduce the spectra generated by tungsten target x-ray tubes. However, as originally presented, the TBC model fails in situations where the determination of x-ray spectra produced by an arbitrary waveform or the calculation of realistic values of air kerma for a specific x-ray system is desired. In the present work, the authors revisited the assumptions used in the original paper published by . They proposed a complementary formulation for taking into account the waveform and the representation of the calculated spectra in a dosimetric quantity. The performance of the proposed model was evaluated by comparing values of air kerma and first and second half value layers from calculated and measured spectra by using different voltages and filtrations. For the output, the difference between experimental and calculated data was better then 5.2%. First and second half value layers presented differences of 23.8% and 25.5% in the worst case. The performance of the model in accurately calculating these data was better for lower voltage values. Comparisons were also performed with spectral data measured using a CZT detector. Another test was performed by the evaluation of the model when considering a waveform distinct of a constant potential. In all cases the model results can be considered as a good representation of the measured data. The results from the modifications to the TBC model introduced in the present work reinforce the value of the TBC model for application of quantitative evaluations in radiation physics.

Efficient x-ray tube simulations

This article describes an efficiency study of directional bremsstrahlung splitting (DBS) for x-ray tube modeling. DBS is shown to be up to five or six orders of magnitude more efficient at 50 or 135 kV tube potential than a simulation without splitting, and 60 times more efficient compared to uniform bremsstrahlung splitting. A methodology is presented to determine the optimum splitting number for a given situation using a second degree polynomial expression derived from theoretical considerations. Very large optimum splitting numbers are found for small fields (5 mm radius) at 1 m from the x-ray source, which are relevant for half-value layer (HVL) calculations and for simulations related to primary air kerma standards. Two approaches for the calculation of kerma at a plane and inside a volume using track-length estimation are implemented in BEAMnrc, a user-code from the EGSnrc Monte Carlo simulation system for photon and electron transport. A practical application of DBS to HVL calculations for a Comet MXR-320 x-ray tube is reported. The agreement with measured HVLs at different constant tube potentials is found to be better than 2.3%.

Assessment of different computational models for generation of x-ray spectra in diagnostic radiology and mammography

Different computational methods based on empirical or semi-empirical models and sophisticated Monte Carlo calculations have been proposed for prediction of x-ray spectra both in diagnostic radiology and mammography. In this work, the x-ray spectra predicted by various computational models used in the diagnostic radiology and mammography energy range have been assessed by comparison with measured spectra and their effect on the calculation of absorbed dose and effective dose (ED) imparted to the adult ORNL hermaphroditic phantom quantified. This includes empirical models (TASMIP and MASMIP), semi-empirical models (X-rayb&m, X-raytbc, XCOMP, IPEM, Tucker et al., and Blough et al.), and Monte Carlo modeling (EGS4, ITS3.0, and MCNP4C). As part of the comparative assessment, the K x-ray yield, transmission curves, and half value layers (HVLs) have been calculated for the spectra generated with all computational models at different tube voltages. The measured x-ray spectra agreed well with the generated spectra when using X-raytbc and IPEM in diagnostic radiology and mammography energy ranges, respectively. Despite the systematic differences between the simulated and reference spectra for some models, the student's t-test statistical analysis showed there is no statistically significant difference between measured and generated spectra for all computational models investigated in this study. The MCNP4C-based Monte Carlo calculations showed there is no discernable discrepancy in the calculation of absorbed dose and ED in the adult ORNL hermaphroditic phantom when using different computational models for generating the x-ray spectra. Nevertheless, given the limited flexibility of the empirical and semi-empirical models, the spectra obtained through Monte Carlo modeling offer several advantages by providing detailed information about the interactions in the target and filters, which is relevant for the design of new target and filter combinations and optimization of radiological imaging protocols.

Monte Carlo simulation of x-ray spectra in diagnostic radiology and mammography using MCNP4C

The general purpose Monte Carlo N-particle radiation transport computer code (MCNP4C) was used for the simulation of x-ray spectra in diagnostic radiology and mammography. The electrons were transported until they slow down and stop in the target. Both bremsstrahlung and characteristic x-ray production were considered in this work. We focus on the simulation of various target/filter combinations to investigate the effect of tube voltage, target material and filter thickness on x-ray spectra in the diagnostic radiology and mammography energy ranges. The simulated x-ray spectra were compared with experimental measurements and spectra calculated by IPEM report number 78. In addition, the anode heel effect and off-axis x-ray spectra were assessed for different anode angles and target materials and the results were compared with EGS4-based Monte Carlo simulations and measured data. Quantitative evaluation of the differences between our Monte Carlo simulated and comparison spectra was performed using student's t-test statistical analysis. Generally, there is a good agreement between the simulated x-ray and comparison spectra, although there are systematic differences between the simulated and reference spectra especially in the K-characteristic x-rays intensity. Nevertheless, no statistically significant differences have been observed between IPEM spectra and the simulated spectra. It has been shown that the difference between MCNP simulated spectra and IPEM spectra in the low energy range is the result of the overestimation of characteristic photons following the normalization procedure. The transmission curves produced by MCNP4C have good agreement with the IPEM report especially for tube voltages of 50 kV and 80 kV. The systematic discrepancy for higher tube voltages is the result of systematic differences between the corresponding spectra.

An alignment method for mammographic X-ray spectroscopy under clinical conditions

This paper describes an alignment method for mammographic X-ray spectroscopy under clinical conditions. A pinhole, a fluorescent screen, a laser device and the case for a detector are used for alignment of the focal spot, a collimator and a detector. The method determines the line between the focal spot and the point of interest in an X-ray field radiographically. The method allows alignment for both central axis and off-axis directions.

Field nonuniformity correction for quantitative analysis of digitized mammograms

Several factors, including the heel effect, variation in distance from the x-ray source to points in the image and path obliquity contribute to the signal nonuniformity of mammograms. To best use digitized mammograms for quantitative image analysis, these field non-uniformities must be corrected. An empirically based correction method, which uses a bowl-shaped calibration phantom, has been developed. Due to the annular spherical shape of the phantom, its attenuation is constant over the entire image. Remaining nonuniformities are due only to the heel and inverse square effects as well as the variable path through the beam filter, compression plate and image receptor. In logarithmic space, a normalized image of the phantom can be added to mammograms to correct for these effects. Then, an analytical correction for path obliquity in the breast can be applied to the images. It was found that the correction causes the errors associated with field nonuniformity to be reduced from 14% to 2% for a 4 cm block of material corresponding to a combination of 50% fibroglandular and 50% fatty breast tissue. A repeatability study has been conducted to show that in regions as far as 20 cm away from the chest wall, variations due to imaging conditions and phantom alignment contribute to <2% of overall corrected signal.

Enhancing the ratio of fluorescence to bremsstrahlung radiation in X-ray tube spectra

This paper describes techniques that can be used to improve the ratio of fluorescence to bremsstrahlung radiation (F/B) in X-ray tube spectra. Firstly, an extension of the EGS4 code system is used to evaluate the impact of the substrate in thin target applications, in terms of the yield of bremsstrahlung photons produced. The choice of materials to filter X-ray tube spectra, and its effect in the F/B and the tube efficiency is discussed. The characteristics of spectra produced in transmission tubes with different target thicknesses, substrates and tube voltages are also presented.

Monte Carlo simulation of x-ray spectra in mammography

A model for generating x-ray spectra in mammography is presented. This model used the ITS version 3 Monte Carlo code for simulating the radiation transport. Various target/filter combinations such as tungsten/aluminium, molybdenum/molybdenum, molybdenum/rhodium and rhodium/rhodium were used in the simulation. Both bremsstrahlung and characteristic x-ray production were included in the model. The simulated x-ray emission spectra were compared with two sets of spectra, those of Boone et al (1997 Med. Phys. 24 1863-74) and IPEM report 78. The chi2 test was used for the overall goodness of fit of the spectral data. There is good agreement between the simulated x-ray spectra and the comparison spectra as the test yielded a probability value of nearly 1. When the transmitted x-ray spectra for specific target/filter combinations were generated and compared with a measured molybdenum/rhodium spectrum and spectra generated in IPEM report 78, close agreement is also observed. This was demonstrated by the probability value for the chi2 test being almost 1 for all the cases. However, minor differences between the simulated spectra and the 'standard' ones are observed.