Modeling the primary source intensity distribution: reconstruction and inter-comparison of six Varian TrueBeam sources

The primary source size is one of the most important beam model parameters in small photon fields. In this work we apply a recently suggested reconstruction technique to characterize the primary source of 6 Varian TrueBeam (TB) linacs. A series of photon fluence profile measurements were performed on 6 Varian TB linacs in the crossplane and inplane orientation using radiochromic film in air and a 2 mm Pb foil as a build-up layer. An image reconstruction algorithm was then applied, based on the maximum likelihood expectation-maximization (MLEM) algorithm, to estimate the source distribution. The method iteratively ray-traces photons from the source plane to the measurement plane to extract source profile corrections. The technique was first benchmarked using a Monte Carlo (MC) model of a Varian TrueBeam with known input Gaussian source sizes. The robustness of the suggested technique was also tested by randomly sampling different combinations of source and field size values and repeating the reconstruction. At the MC benchmarking stage the MLEM reconstruction algorithm was capable of reproducing the Gaussian shape with a RMSE less than 4.0%, while the reconstructed source size (FWHM) and field size were determined with an accuracy level of 0.14 mm and 0.10 mm respectively. Experimentally, the reconstructed TB sources presented FWHM values between 1.02-1.5 mm (-0.18 mm) and 1.08-1.42 mm (-0.13 mm) in the crossplane and inplane orientations respectively. All TB sources studied in this work can be considered symmetric within uncertainties with the exception of one. The source distribution presented systematic deviations from a Gaussian distribution mostly in the lower tail region. Multi-parameter functional forms, such as Pearson VII or double Gaussian presented improvements in modeling the source in this region, but increase the model complexity. The reconstructed sources measured in this work can serve as reference values for commissioning beam models in small fields and set upper and lower thresholds values of the expected source size for a TB linac.