How I can calculate and plot the spread out Bragg Peak, or SOBP in Geant4?

Hi.
I am working on simulating an ion beam in the water phantom. I need to plot spread out Bragg Peak (SOBP) for my project. But unfortunately, I do not know what should I do.
I would be grateful for any help or suggestion on this issue.
Thank you

Dear,

You can calculate and plot the spread out Bragg Peak using the advanced example “hadrontherapy” in Geant4.

best wishes,
Batmunkh.