A probabilistic approach for determining Monte Carlo beam source parameters: I. Modeling of a CyberKnife M6 unit

Objective. During Monte Carlo modeling of external radiotherapy beams, models must be adjusted to reproduce the experimental measurements of the linear accelerator being considered. The aim of this work is to propose a new method for the determination of the energy and spot size of the electron beam incident on the target of a linear accelerator using a maximum likelihood estimation. Approach. For that purpose, the method introduced by Francescon et al (2008 Med. Phys. 35 504-13) is expanded upon in this work. Simulated tissue-phantom ratios and uncorrected output factors using a set of different detector models are compared to experimental measurements. A probabilistic formalism is developed and a complete uncertainty budget, which includes a detailed simulation of positioning errors, is evaluated. The method is applied to a CyberKnife M6 unit using four detectors (PTW 60012, PTW 60019, Exradin A1SL and IBA CC04), with simulations being performed using the EGSnrc suite. Main results. The likelihood distributions of the electron beam energy and spot size are evaluated, leading to (E) over cap = 7.42 +/- 0.17 MeV and (F) over cap = 2.15 + 0.06 mm. Using these results and a 95% confidence region, simulations reproduce measurements in 13 out of the 14 considered setups. Significance. The proposed method allows an accurate beam parameter optimization and uncertainty evaluation during the Monte Carlo modeling of a radiotherapy unit.