Proton beam with realistic geometry

From UCL HEP PBT Wiki

Jump to: navigation, search
Line 63: Line 63:
== <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  
+
The macro produces a root file '''Proton.root''' with two histograms. The first histogram shows the energy deposition in water box along the beam line, the second histogram shows zoomed Bragg peak. The macro also produces two text files '''DoseFile.txt''' and '''PlotDose.txt'''. The file '''DoseFile.txt''' contains energy and dose deposition at each water box layer. The file '''PlotDose.txt''' contains depth vs dose values. These values can be imported in matlab. This is an example output for '''DoseFile.txt''' with physics process '''QGSP_BIC_EMY''' and incident proton energy of '''62 MeV'''.  
-
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">
<pre style="color: #800000; background-color: #dcdcdc">
-
  Cumulated Doses : X[mm]   Edep  Edep/Ebeam Dose
+
  Layers x[mm] Edep  Edep/Ebeam[%] Dose Dose/MaxDose[%]
-
  layer 1: 3.33333 22.8684 GeV 6.14768 % 2.74795e-06 Gy
+
  layer 1: 1.33333 8.83908 GeV 2.376 1.062e-05 Gy 21.88
-
  layer 2: 6.66667 24.4078 GeV 6.56151 % 2.93293e-06 Gy
+
layer 2: 2.667 8.928 GeV 2.4 1.073e-05 Gy 22.11
-
  layer 3: 10 25.8277 GeV 6.94321 % 3.10354e-06 Gy
+
  layer 3: 4 9.214 GeV 2.477 1.107e-05 Gy 22.81
-
  layer 4: 13.3333 27.87 GeV 7.49224 % 3.34895e-06 Gy
+
layer 4: 5.333 9.468 GeV 2.545 1.138e-05 Gy 23.44
-
  layer 5: 16.6667 30.3243 GeV 8.15203 % 3.64387e-06 Gy
+
layer 5: 6.667 9.577 GeV 2.575 1.151e-05 Gy 23.71
-
  layer 6: 20 34.0826 GeV 9.16235 % 4.09547e-06 Gy
+
layer 6: 8 9.767 GeV 2.626 1.174e-05 Gy 24.18
-
  layer 7: 23.3333 39.1645 GeV 10.5285 % 4.70613e-06 Gy
+
layer 7: 9.333 10 GeV         2.689 1.202e-05 Gy 24.77
-
  layer 8: 26.6667 50.4216 GeV 13.5547 % 6.05883e-06 Gy
+
  layer 8: 10.67 10.27 GeV 2.761 1.234e-05 Gy 25.43
-
  layer 9: 30 85.9984 GeV 23.1188 % 1.03338e-05 Gy
+
layer 9: 12 10.48 GeV 2.817 1.259e-05 Gy 25.94
-
  layer 10: 33.3333 3.84312 GeV 1.03314 % 4.61801e-07 Gy
+
  layer 10: 13.33 10.79 GeV 2.901 1.297e-05 Gy 26.72
-
  layer 11: 36.6667 1.15485 MeV 0.000310456 % 1.38771e-10 Gy
+
  layer 11: 14.67 11.11 GeV 2.986 1.335e-05 Gy 27.5
-
  layer 12: 40 0 eV 0 % 0 Gy
+
layer 12: 16 11.58 GeV 3.114 1.392e-05 Gy 28.68
 +
layer 13: 17.33 12.04 GeV 3.236 1.447e-05 Gy 29.81
 +
layer 14: 18.67 12.4 GeV 3.334 1.49e-05 Gy 30.7
 +
  layer 15: 20 12.98 GeV 3.489 1.56e-05 Gy 32.14
 +
layer 16: 21.33 13.67 GeV 3.675 1.643e-05 Gy 33.84
 +
  layer 17: 22.67 14.39 GeV 3.868 1.729e-05 Gy 35.63
 +
layer 18: 24 15.4 GeV 4.139 1.85e-05 Gy 38.12
 +
layer 19: 25.33 16.57 GeV 4.454 1.991e-05 Gy 41.02
 +
  layer 20: 26.67 18.19 GeV 4.89 2.186e-05 Gy 45.04
 +
layer 21: 28 20.69 GeV 5.562 2.486e-05 Gy 51.22
 +
layer 22: 29.33 24.52 GeV 6.592 2.946e-05 Gy 60.71
 +
  layer 23: 30.67 33.19 GeV 8.922 3.988e-05 Gy 82.17
 +
  layer 24: 32 40.39 GeV 10.86 4.853e-05 Gy 100
 +
layer 25: 33.33 3.536 GeV 0.9507 4.249e-06 Gy 8.756
 +
  layer 26: 34.67 1.466 MeV 0.0003942 1.762e-09 Gy 0.00363
 +
layer 27: 36 10.16 MeV 0.002731  1.221e-08 Gy 0.02515
 +
layer 28: 37.33 126 keV         3.387e-05 1.514e-10 Gy 0.0003119
 +
  layer 29: 38.67 18.63 keV 5.009e-06 2.239e-11 Gy 4.613e-05
 +
layer 30: 40 0 eV         0           0 Gy 0
-
  The run consists of 6000  protons of 62 MeV through 4 cm  of Water (density: 1 g/cm3 )
+
  The run consists of 6000  protons of 62 MeV through 4 cm  of Water (density: 1 g/cm3 ) divided into 30 slices.
-
  Total incident energy(Ebeam)= 371.985 GeV
+
 
-
  Total energy deposit= 366.469 GeV
+
Edep is the deposited energy in every slice.
 +
  Total incident energy(Ebeam)= 372 GeV
 +
  Total energy deposit= 366.6 GeV
 +
Dose is the deposited dose in every slice.
 +
MaxDose is the highest dose value from all slices.
</pre>
</pre>
 +
 +
The corresponding '''PlotDose.txt''' is
 +
 +
<pre style="color: #800000; background-color: #dcdcdc">
 +
1.33333 21.8846
 +
2.66667 22.106
 +
4 22.8118
 +
5.33333 23.4411
 +
6.66667 23.7111
 +
8 24.1827
 +
9.33333 24.7673
 +
10.6667 25.4289
 +
12 25.9402
 +
13.3333 26.7199
 +
14.6667 27.4969
 +
16 28.681
 +
17.3333 29.8068
 +
18.6667 30.704
 +
20 32.136
 +
21.3333 33.8436
 +
22.6667 35.6258
 +
24 38.1195
 +
25.3333 41.0227
 +
26.6667 45.0388
 +
28 51.2228
 +
29.3333 60.709
 +
30.6667 82.172
 +
32 100
 +
33.3333 8.75568
 +
34.6667 0.00363025
 +
36 0.0251539
 +
37.3333 0.000311934
 +
38.6667 4.61342e-05
 +
40 0
 +
</pre>
You can open Proton.root file in the following way:
You can open Proton.root file in the following way:
Line 98: Line 153:
</pre>
</pre>
   
   
-
 
You can change the physics process, incident proton energy and number of slices etc. by
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.
modifying the macro proton.mac. In addition, you can configure the particle source by using '''/gps/''' commands.
Line 123: Line 177:
/protonGB/det/setSizeX  4 cm
/protonGB/det/setSizeX  4 cm
/protonGB/det/setSizeYZ 4 cm
/protonGB/det/setSizeYZ 4 cm
-
/protonGB/det/setSliceSizeYZ 2 cm
+
/protonGB/det/setSliceSizeYZ 1 cm
-
/protonGB/det/sliceNumber 12
+
/protonGB/det/sliceNumber 30
#
#
# set physics process
# set physics process
Line 133: Line 187:
# production tresholds (recommended range  
# production tresholds (recommended range  
#cut off not bigger than 10% of slice thickness)
#cut off not bigger than 10% of slice thickness)
-
/protonGB/phys/setCuts 0.5 mm
+
/protonGB/phys/setCuts 0.2 mm
#/protonGB/phys/setGCut 1 um
#/protonGB/phys/setGCut 1 um
#/protonGB/phys/setECut 1 um
#/protonGB/phys/setECut 1 um
Line 168: Line 222:
# step limit (recommended not bigger than 5% of  
# step limit (recommended not bigger than 5% of  
# slice thickness)
# slice thickness)
-
/protonGB/stepMax 0.2 mm
+
/protonGB/stepMax 0.1 mm
#
#
/protonGB/event/printModulo 50
/protonGB/event/printModulo 50

Revision as of 19:56, 21 July 2014

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