Proton beam with realistic geometry

From UCL HEP PBT Wiki

Jump to: navigation, search
m
Line 62: Line 62:
== <span style="color:#000080"> How to analyze data </span> ==
== <span style="color:#000080"> How to analyze data </span> ==
 +
 +
The macro produces a root file '''Proton.root''' with a 1D histogram showing the energy deposition in
 +
water box along the beam line and another histogram showing zoomed Bragg Peak. It also produces a text file '''DoseFile.txt''' with energy and dose deposited at each slice of the water box. This is an example output
 +
with the default settings: physics process '''QGSP_BIC_EMY''' and incident proton energy of 62 MeV.
 +
 +
<pre style="color: #800000; background-color: #dcdcdc">
 +
Cumulated Doses : X[mm]    Edep  Edep/Ebeam Dose
 +
layer 1: 3.33333 22.8684 GeV 6.14768 % 2.74795e-06 Gy
 +
layer 2: 6.66667 24.4078 GeV 6.56151 % 2.93293e-06 Gy
 +
layer 3: 10 25.8277 GeV 6.94321 % 3.10354e-06 Gy
 +
layer 4: 13.3333 27.87 GeV 7.49224 % 3.34895e-06 Gy
 +
layer 5: 16.6667 30.3243 GeV 8.15203 % 3.64387e-06 Gy
 +
layer 6: 20 34.0826 GeV 9.16235 % 4.09547e-06 Gy
 +
layer 7: 23.3333 39.1645 GeV 10.5285 % 4.70613e-06 Gy
 +
layer 8: 26.6667 50.4216 GeV 13.5547 % 6.05883e-06 Gy
 +
layer 9: 30 85.9984 GeV 23.1188 % 1.03338e-05 Gy
 +
layer 10: 33.3333 3.84312 GeV 1.03314 % 4.61801e-07 Gy
 +
layer 11: 36.6667 1.15485 MeV 0.000310456 % 1.38771e-10 Gy
 +
layer 12: 40 0 eV 0 % 0 Gy
 +
 +
 +
The run consists of 6000  protons of 62 MeV through 4 cm  of Water (density: 1 g/cm3 )
 +
Total incident energy(Ebeam)= 371.985 GeV
 +
Total energy deposit= 366.469 GeV
 +
</pre>
 +
 +
You can open Proton.root file in the following way:
 +
 +
<pre style="color: #800000; background-color: #dcdcdc">
 +
root -l Proton.root
 +
 +
new TBrowser
 +
 +
Select ROOT Files and Proton.root
 +
</pre>
 +
 +
 +
You can change the physics process, incident proton energy and number of slices etc. by
 +
modifying the macro proton.mac. In addition, you can configure the particle source by using '''/gps/''' commands.
 +
In this macro the proton beam has Guassian energy distribution. The beam characteristics are the same as the ones used at the Laboratori Nazionali del Sud (INFN) in Catania, Italy. Have a look [https://geant4.web.cern.ch/geant4/UserDocumentation/UsersGuides/ForApplicationDeveloper/html/ch02s07.html here] to find out the meaning of all /gps/ commands.
 +
 +
Now, use your favorite editor '''pico''', '''vi''', '''emacs''' etc. to open the macro proton.mac:
 +
 +
<pre style="color: #800000; background-color: #dcdcdc">
 +
pico proton.mac
 +
</pre> 
 +
 +
You will see:
 +
 +
<pre style="color: #800000; background-color: #dcdcdc">
 +
# proton.mac
 +
#
 +
/control/verbose 2
 +
/run/verbose 2
 +
/tracking/verbose 0
 +
/run/particle/verbose 1
 +
/run/particle/dumpList
 +
#
 +
# set geometry
 +
/protonGB/det/setSizeX  4 cm
 +
/protonGB/det/setSizeYZ 4 cm
 +
/protonGB/det/setSliceSizeYZ 2 cm
 +
/protonGB/det/sliceNumber 12
 +
#
 +
# set physics process
 +
/protonGB/phys/addPhysics QGSP_BIC_EMY
 +
#/protonGB/phys/addPhysics emlivermore
 +
#/protonGB/phys/addPhysics empenelope
 +
#
 +
# production tresholds (recommended range
 +
#cut off not bigger than 10% of slice thickness)
 +
/protonGB/phys/setCuts 0.5 mm
 +
#/protonGB/phys/setGCut 1 um
 +
#/protonGB/phys/setECut 1 um
 +
#/protonGB/phys/setPCut 1 um
 +
#
 +
# initialize
 +
/run/initialize
 +
#
 +
# General particle source
 +
# proton circle source 
 +
/gps/pos/shape Circle
 +
/gps/pos/centre -4. 0. 0. cm
 +
/gps/pos/radius 0. mm
 +
/gps/pos/sigma_r 2. mm
 +
/gps/particle proton
 +
/gps/pos/type Beam
 +
#
 +
# the incident surface is in the y-z plane
 +
/gps/pos/rot1 0 1 0
 +
/gps/pos/rot2 0 0 1
 +
#
 +
# the beam is travelling along the x-axis without any angular
 +
#dispersion (angular despersion set to 0.0)
 +
/gps/ang/rot1 0 0 1
 +
/gps/ang/rot2 0 1 0
 +
/gps/ang/type beam1d
 +
/gps/ang/sigma_r 0. deg
 +
#
 +
# the beam energy is in gaussian profile
 +
/gps/ene/type Gauss
 +
/gps/ene/mono 62 MeV
 +
/gps/ene/sigma 0.3 MeV
 +
#
 +
# step limit (recommended not bigger than 5% of
 +
# slice thickness)
 +
/protonGB/stepMax 0.2 mm
 +
#
 +
/protonGB/event/printModulo 50
 +
#
 +
# output file
 +
/analysis/setFileName Proton
 +
#
 +
/analysis/h1/set 2 50 25 35 mm
 +
# number of events
 +
/run/beamOn 6000
 +
#
 +
</pre>
 +
== <span style="color:#000080"> Files </span> ==
== <span style="color:#000080"> Files </span> ==
[[List of files with brief description]]
[[List of files with brief description]]

Revision as of 14:35, 21 July 2014

Retrieved from "https://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/index.php/Proton_beam_with_realistic_geometry"
Personal tools