|
|
| (376 intermediate revisions by 3 users not shown) |
| Line 1: |
Line 1: |
| == Simulation of the Clatterbridge beamline ==
| |
|
| |
|
| This simulation models the monoenergetic 62.5 MeV proton beam at Clatterbridge Cancer Centre, as it traverses the components of the beamline and is then deposited into a volume of water. The beamline components are contained within a geometry modelling the Clatterbridge treatment room. | | This simulation is a model of the monoenergetic 62.5 MeV proton beam at the [http://www.clatterbridgecc.nhs.uk/patients/treatment-and-support/proton-therapy Clatterbridge Cancer Centre] as it traverses the components of the beamline and finally hits a volume of water. |
| | | The simulation was built using: |
| The protons are generated using the G4ParticleGun class, and the physics list used is QGSP_BIC_HP, standard for simulating clinical proton beams. | | * [[/Geant4 | Geant4]] |
| | | ** The simulation was built on the example in <code>examples/extended/electromagnetic/TestEm7</code> supplied with the Geant4 package and detailed [http://www.hep.ucl.ac.uk/pbt/wiki/Software/Geant4/Tutorials/Basic/Monoenergetic_Proton_Pencil_Beam here]. |
| The energy of the beam after travelling through the beamline components is measured by tracking the energy deposition of individual protons within the water volume, using the G4Track class. This simulation produces a post-beamline energy of 60.08 MeV and a Bragg peak at depth [insert depth].
| | * [[/TOPAS | TOPAS]] |
| | |
| == Running the simulation ==
| |
| | |
| '''Run macro proton.mac'''
| |
| | |
| <pre>
| |
| [username@plus1 ProtonPB_build]$ ./protonPB proton.mac
| |
| </pre>
| |
| | |
| == Data Analysis ==
| |
| | |
| The macro produces several output files. | |
| | |
| '''hits.out'''
| |
| | |
| This file contains information about the energy deposition in the water volume.
| |
| | |
| '''Open ROOT and run the simulation analysis file '''
| |
| | |
| <pre>
| |
| [username@plus1 ProtonPB_build]$ root -l
| |
| root [1] .x simulation_analysis.C
| |
| </pre>
| |
| | |
| == Changing parameters ==
| |
| | |
| '''Scoring mesh'''
| |
| | |
| '''Initial beam parameters'''
| |
| | |
| ''Beam radius''
| |
| | |
| <pre> | |
| /gps/pos/radius 3 mm | |
| </pre> | |
| | |
| | |
| ''Beam energy''
| |
| | |
| <pre>
| |
| /gps/ene/type Gauss | |
| /gps/ene/mono 62.5 MeV
| |
| /gps/ene/sigma 0.082 MeV
| |
| </pre>
| |
| | |
| ''Source position''
| |
| | |
| <pre>
| |
| /gps/pos/type Plane | |
| /gps/pos/shape Circle | |
| /gps/pos/centre 0.0 0.0 -420 cm
| |
| </pre>
| |