Range verification based on prompt gamma detection is an important step to improve dose control for proton therapy. To deduce the proton range from the detected prompt gamma emission, a prediction of the measured profile is required. This study introduces the Boltzmann solver as a faster alternative to the Monte Carlo simulations to produce dose distributions and prompt gamma source terms from proton therapy treatment plans.