Monoenergetic proton pencil beam

From UCL HEP PBT Wiki

Jump to: navigation, search
m
m
Line 64: Line 64:
== <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 energy deposition around the peak. The macro also produces two text files '''DoseFile.txt''' and '''PlotDose.txt'''. The file '''DoseFile.txt''' contains energy and dose deposition for every layer. The file '''PlotDose.txt''' contains only depth vs dose values for each layer. These values can be imported in matlab.
-
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.  
+
=== <span style="color:#000080"> Text files </span> ===
 +
 
 +
This is an example output for '''DoseFile.txt''' with 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 23.1611 GeV 6.22611 % 2.78311e-06 Gy
+
  layer 1: 1.33333 8.7788 GeV 2.35989 1.05489e-05 Gy 20.404
-
  layer 2: 6.66667 24.3512 GeV 6.54602 % 2.92612e-06 Gy
+
layer 2: 2.66667 8.95366 GeV 2.4069 1.0759e-05 Gy 20.8104
-
  layer 3: 10 25.9655 GeV 6.97998 % 3.1201e-06 Gy
+
  layer 3: 4 9.20046 GeV 2.47324 1.10556e-05 Gy 21.384
-
  layer 4: 13.3333 27.6521 GeV 7.43336 % 3.32277e-06 Gy
+
layer 4: 5.33333 9.33632 GeV 2.50976 1.12188e-05 Gy 21.6998
-
  layer 5: 16.6667 30.1527 GeV 8.10556 % 3.62325e-06 Gy
+
layer 5: 6.66667 9.58072 GeV 2.57546 1.15125e-05 Gy 22.2678
-
  layer 6: 20 33.6152 GeV 9.03635 % 4.03931e-06 Gy
+
layer 6: 8 9.76978 GeV 2.62628 1.17397e-05 Gy 22.7073
-
  layer 7: 23.3333 39.396 GeV 10.5903 % 4.73395e-06 Gy
+
layer 7: 9.33333 10.0149 GeV 2.69217 1.20342e-05 Gy 23.277
-
  layer 8: 26.6667 50.3177 GeV 13.5263 % 6.04634e-06 Gy
+
  layer 8: 10.6667 10.248 GeV 2.75484 1.23143e-05 Gy 23.8187
-
  layer 9: 30 88.6343 GeV 23.8264 % 1.06506e-05 Gy
+
layer 9: 12 10.5163 GeV 2.82695 1.26367e-05 Gy 24.4422
-
  layer 10: 33.3333 2.09587 GeV 0.563405 % 2.51846e-07 Gy
+
  layer 10: 13.3333 10.7679 GeV 2.89459 1.2939e-05 Gy 25.0271
-
  layer 11: 36.6667 4.59086 MeV 0.0012341 % 5.51653e-10 Gy
+
layer 11: 14.6667 11.2918 GeV 3.03543 1.35686e-05 Gy 26.2448
-
  layer 12: 40 0 eV 0 % 0 Gy
+
  layer 12: 16 11.5998 GeV 3.11822 1.39387e-05 Gy 26.9606
 +
layer 13: 17.3333 11.9594 GeV 3.21489 1.43708e-05 Gy 27.7964
 +
layer 14: 18.6667 12.6008 GeV 3.38732 1.51416e-05 Gy 29.2873
 +
  layer 15: 20 13.0799 GeV 3.51611 1.57173e-05 Gy 30.4008
 +
  layer 16: 21.3333 13.853 GeV 3.72393 1.66463e-05 Gy 32.1977
 +
layer 17: 22.6667 14.4719 GeV 3.89029 1.73899e-05 Gy 33.636
 +
layer 18: 24 15.3391 GeV 4.12341 1.8432e-05 Gy 35.6517
 +
  layer 19: 25.3333 16.7448 GeV 4.50128 2.01211e-05 Gy 38.9188
 +
layer 20: 26.6667 18.2636 GeV 4.90957 2.19461e-05 Gy 42.4489
 +
layer 21: 28 20.5916 GeV 5.53539 2.47436e-05 Gy 47.8598
 +
layer 22: 29.3333 24.3455 GeV 6.54449 2.92543e-05 Gy 56.5847
 +
  layer 23: 30.6667 32.4192 GeV 8.71484 3.8956e-05 Gy 75.3499
 +
layer 24: 32 43.0249 GeV 11.5658 5.17001e-05 Gy 100
 +
  layer 25: 33.3333 1.97182 GeV 0.530059 2.3694e-06 Gy 4.58297
 +
  layer 26: 34.6667 0 eV         0 0 Gy         0
 +
layer 27: 36 0 eV         0 0 Gy         0
 +
layer 28: 37.3333 442.484 keV 0.000118947 5.31703e-10 Gy 0.00102844
 +
  layer 29: 38.6667 0 eV         0     0 Gy 0
 +
layer 30: 40 0 eV         0     0 Gy 0
 +
 
 +
 
 +
The run consists of 6000 proton of 62 MeV through 4 cm  of Water (density: 1 g/cm3 ) divided into 30 slices.
 +
 
 +
Edep is the deposited energy in every slice.
 +
Total incident energy(Ebeam)= 372 GeV
 +
Total energy deposit= 367.309 GeV
 +
Dose is the deposited dose in every slice.
 +
MaxDose is the highest dose value from all slices.
</pre>
</pre>
-
You can open Proton.root file in the following way:
+
 
 +
Open Proton.root file in the following way:
<pre style="color: #800000; background-color: #dcdcdc">
<pre style="color: #800000; background-color: #dcdcdc">
Line 95: Line 125:
Inside is 1D histogram showing the energy deposition in water box:
Inside is 1D histogram showing the energy deposition in water box:
-
+
 
-
[[File:/pbt/RadiotherapyWorkbook/skins/common/images/PhotonPB/Edep.png|Energy deposit per event along the beam]]
+
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

Revision as of 10:31, 22 July 2014

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