A system, method, and computer program product of clinical treatment planning implements complex Monte Carlo (MC) based brachytherapy dose distributions using conventional brachytherapy treatment planning systems (TPS). Dose distributions from complex brachytherapy source configurations determined with MC methods are used as inputs. Radial dose functions and 2D anisotropy functions are obtained by positioning the coordinate system origin along the dose distribution cylindrical axis of symmetry. Origin to tissue distance and active length are chosen to minimize TPS interpolation errors. A 2D anisotropy function is determined, and a brachytherapy dose rate constant is selected. A virtual brachytherapy source dose distribution is calculated based upon the complex treatment configuration. Additional dosimetry parameters may be considered as well, and dose distributions may be calculated and compared to the original MC-derived dose distributions. The present techniques may calculate dose to a specific tissue type instead of dose to water as used in the TG-43 formalism
展开▼