Monoenergetic photon pencil beam

From UCL HEP PBT Wiki

Jump to: navigation, search
m
Line 1: Line 1:
== <span style="color:#000080"> Introduction </span> ==
== <span style="color:#000080"> Introduction </span> ==
-
This example shows the dose distribution in water along the incident photon beam. The beam hits the water cube surface and deposits a dose under the surface of the water. The volume of the water cube is divided into voxels. At each voxel the dose contribution from the pencil beams is summed up and yields the total dose along the beam axis.  
+
This example shows the dose distribution in water along the incident photon beam. The beam hits the water cube surface and deposits a dose under the surface of the water. The volume of the water cube is divided into slices perpendicular to the incident beam. The dose and energy deposited from the pencil beams at each slice is computed.  
 +
 
== <span style="color:#000080"> Setting up the environment </span> ==  
== <span style="color:#000080"> Setting up the environment </span> ==  
Line 21: Line 22:
<pre style="color: #800000; background-color: #dcdcdc">
<pre style="color: #800000; background-color: #dcdcdc">
-
source   
+
source /unix/pbt/software/dev/bin/pbt-dev.sh  
</pre>
</pre>
== <span style="color:#000080"> How to get the code </span> ==
== <span style="color:#000080"> How to get the code </span> ==
-
; Copy the code to your area
+
; Copy the code to your working directory
<pre style="color: #800000; background-color: #dcdcdc">
<pre style="color: #800000; background-color: #dcdcdc">
Line 34: Line 35:
== <span style="color:#000080"> How to run the code </span> ==  
== <span style="color:#000080"> How to run the code </span> ==  
-
; Inside /home/username/.../ create directory
+
; Inside /home/username/.../ create a directory
<pre style="color: #800000; background-color: #dcdcdc">
<pre style="color: #800000; background-color: #dcdcdc">
Line 40: Line 41:
</pre>
</pre>
-
; Change to this directory and run cmake and make to build the PhotonPB tutorial
+
; To compile the code enter this directory and run cmake and make  
<pre style="color: #800000; background-color: #dcdcdc">
<pre style="color: #800000; background-color: #dcdcdc">
Line 50: Line 51:
</pre>
</pre>
 +
; Run the macro gamma.mac. The macro generates 10000 events.
 +
<pre style="color: #800000; background-color: #dcdcdc">
 +
./photonPB gamma.mac
 +
</pre>
 +
== <span style="color:#000080"> How to analyze data </span> ==
== <span style="color:#000080"> How to analyze data </span> ==
 +
 +
The macro produces a root file with 1D histogram showing the energy deposition in
 +
the water box along the beam line. It also produces a text file with energy and dose
 +
deposited at each slice of the water box. This is an example output with the default settings:
 +
physics process-emstandard_opt0 and incident gamma energy of 20 MeV.
 +
 +
<pre style="color: #800000; background-color: #dcdcdc">
 +
Cumulated Doses : X[mm] Edep      Edep/Ebeam Dose
 +
layer 1: 10 695.299 MeV 0.347649 % 6.96245e-09 Gy
 +
layer 2: 20 1.17329 GeV 0.586646 % 1.17489e-08 Gy
 +
layer 3: 30 1.60782 GeV 0.80391 % 1.61001e-08 Gy
 +
layer 4: 40 1.89698 GeV 0.94849 % 1.89956e-08 Gy
 +
layer 5: 50 2.23625 GeV 1.11812 % 2.23929e-08 Gy
 +
layer 6: 60 2.41072 GeV 1.20536 % 2.414e-08 Gy
 +
layer 7: 70 2.5396 GeV 1.2698 % 2.54306e-08 Gy
 +
layer 8: 80 2.66291 GeV 1.33145 % 2.66653e-08 Gy
 +
layer 9: 90 2.5538 GeV 1.2769 % 2.55727e-08 Gy
 +
layer 10: 100 2.47482 GeV 1.23741 % 2.47819e-08 Gy
 +
layer 11: 110 2.45013 GeV 1.22507 % 2.45346e-08 Gy
 +
layer 12: 120 2.34969 GeV 1.17484 % 2.35288e-08 Gy
 +
layer 13: 130 2.3329 GeV 1.16645 % 2.33608e-08 Gy
 +
layer 14: 140 2.37302 GeV 1.18651 % 2.37625e-08 Gy
 +
layer 15: 150 2.4226 GeV 1.2113 % 2.4259e-08 Gy
 +
layer 16: 160 2.4147 GeV 1.20735 % 2.41799e-08 Gy
 +
layer 17: 170 2.38873 GeV 1.19437 % 2.39198e-08 Gy
 +
layer 18: 180 2.36667 GeV 1.18333 % 2.36989e-08 Gy
 +
layer 19: 190 2.28832 GeV 1.14416 % 2.29143e-08 Gy
 +
layer 20: 200 8.08046e+167 J  2.52171e+177 % 5.05029e+169 Gy
 +
 +
 +
The run consists of 10000 gamma of 20 MeV through 20 cm  of Water (density: 1 g/cm3 )
 +
</pre>

Revision as of 19:59, 15 July 2014

Personal tools