Bragg Peak Analysis

The function performs numerical analysis of Bragg Peak curve using three methods

The function performs analysis of Bragg Peak curve giving the parameters of the depth dose distribution. Three methods of analysis are possible using spline fit, 3rd polynomial fit to the peak [1] and Bortfeld fit [2]. All parameters of Bragg curve are calculated using normalisation to maximum of spline, poly3 fit or bortfeld fit for spline, poly3 and Bortfeld analysis method respectively. For Borthelfd analysis the user must assure that depth is in [mm] unit. Each analysis type gives following parameters:
D100 - dose at maximum
R100 - position of maximum
R90P - position of 90% proximal build-up
R90D - position of 90% distal fall-off
R80P - position of 80% proximal build-up
R80D - position of 80% distal fall-off
R20 - position of 20% distal fall-off
R10 - position of 10% distal fall-off
R50P - position of 50% proximal build-up
R50D - position of 50% distal fall-off
FWHM - full width at half maximum
DFO2080 - distal fall-off 20%-80%
DFO1090 - distal fall-off 10%-90%
[1] Parodi, K., et al. (2012). "Monte Carlo simulations to support start-up and treatment planning of scanned proton and carbon ion therapy at a synchrotron-based facility", Physics in Medicine and Biology, 57(12), 3759–84
[2] Bortfeld, T. (1997). "An analytical approximation of the Bragg curve for therapeutic proton beams", Medical Physics, 24(12), 2024.