Monte Carlo simulations for angular and spatial distributions in therapeutic-energy proton beams