Clatterbridge: Difference between revisions
| Line 32: | Line 32: | ||
== Changing parameters == | == Changing parameters == | ||
=== Scoring mesh === | |||
''' | '''Longitudinal scoring mesh''' | ||
''Beam radius'' | A longitudinal scoring mesh extends along the length of the beamline from the source to the water volume. The mesh utilises a filter to detect the flux of protons per cm<sup>2</sup> and writes the data to the text file <code> FluxLongitudinal.txt </code>. The location of the mesh centre can be changed in <code> proton.mac </code>, in addition to the dimensions of the mesh and the number of bins. | ||
'''Lateral scoring mesh''' | |||
A lateral scoring mesh | |||
=== Initial beam parameters === | |||
Initial parameters of the proton beam can be modified in <code> proton.mac </code> | |||
'''Beam radius''' | |||
<pre> | <pre> | ||
| Line 42: | Line 52: | ||
</pre> | </pre> | ||
'''Beam energy''' | |||
This simulation models the proton beam source with a Gaussian distribution. | |||
<pre> | <pre> | ||
| Line 51: | Line 62: | ||
</pre> | </pre> | ||
''Source position'' | '''Source position''' | ||
The proton source is positioned at <code> z = -420 cm </code> relative to the centre of the inner room (the mother volume), which translates as the wall surface of the inner room. | |||
<pre> | <pre> | ||
Revision as of 15:12, 25 July 2016
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.
The protons are generated using the G4ParticleGun class, and the physics list used is QGSP_BIC_HP, standard for simulating clinical proton beams.
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].
Running the simulation
Run macro proton.mac
[username@plus1 ProtonPB_build]$ ./protonPB proton.mac
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
[username@plus1 ProtonPB_build]$ root -l root [1] .x simulation_analysis.C
Changing parameters
Scoring mesh
Longitudinal scoring mesh
A longitudinal scoring mesh extends along the length of the beamline from the source to the water volume. The mesh utilises a filter to detect the flux of protons per cm2 and writes the data to the text file FluxLongitudinal.txt . The location of the mesh centre can be changed in proton.mac , in addition to the dimensions of the mesh and the number of bins.
Lateral scoring mesh
A lateral scoring mesh
Initial beam parameters
Initial parameters of the proton beam can be modified in proton.mac
Beam radius
/gps/pos/radius 3 mm
Beam energy
This simulation models the proton beam source with a Gaussian distribution.
/gps/ene/type Gauss /gps/ene/mono 62.5 MeV /gps/ene/sigma 0.082 MeV
Source position
The proton source is positioned at z = -420 cm relative to the centre of the inner room (the mother volume), which translates as the wall surface of the inner room.
/gps/pos/type Plane /gps/pos/shape Circle /gps/pos/centre 0.0 0.0 -420 cm