Fast Monte Carlo simulation for total body irradiation using a [Formula: see text] teletherapy unit

Our institution delivers TBI using a modified Theratron 780 [Formula: see text] unit. Due to limitations of our treatment planning system in calculating dose for this treatment, we have developed a fast Monte Carlo code to calculate dose distributions within the patient. The algorithm is written in C and uses voxel density information from CT images to calculate dose in heterogeneous media. To test the algorithm, film‐based dose measurements were made separately in a simple water phantom with a high‐density insert and a RANDO phantom and then compared to doses calculated by the Monte Carlo algorithm. In addition, a separate simulation in GEANT4 was run for the RANDO phantom and compared to both film and the in‐house simulation. All results were analyzed using RIT113 film analysis software. Simulations in the water phantom accurately predict the depth of maximum dose in the phantom at 0.5 cm. The measured PDD along the central axis of the beam closely matches the PDD generated from the Monte Carlo code, deviating on average by only 3% along the depth of the water phantom. Dose measured at planes inside the high‐density insert had a mean difference of 4.9% on cross‐profile measurement. In the RANDO phantom, gamma pass rates vary between 91% and 99% at 3 mm, 3%, and were [Formula: see text] at 5 mm, 5% for the four film planes measured. Profiles taken across the film and both simulations resulted in mean relative differences of [Formula: see text] for all profiles in each slice measured. The Monte Carlo algorithm presented here is potentially a viable method for calculating dose distributions delivered in TBI treatments at our center. While not yet refined enough to be the primary method of treatment planning, the algorithm at its current resolution determines the dose distribution for one patient within a few hours, and provides clinically useful information in planning TBI. With appropriate optimization, the Monte Carlo method presented here could potentially be implemented as a first‐line treatment planning option for [Formula: see text] TBI. PACS number: 87.10.Rt