Monoenergetic photon pencil beam

From UCL HEP PBT Wiki

Jump to: navigation, search
m
m
Line 76: Line 76:
<pre style="color: #800000; background-color: #dcdcdc">
<pre style="color: #800000; background-color: #dcdcdc">
-
  Layers : x[mm] Edep         Edep/Ebeam[%]     Dose Dose/MaxDose[%]
+
  Layers : x[mm] Edep       Edep/Ebeam[%] Dose       Dose/MaxDose[%]
-
  layer 1: 10 513.656 MeV 0.428046 5.14354e-09 Gy 32.1047
+
  layer 1: 5 119.065 MeV 0.0992211 9.53818e-11 Gy 15.1483
-
  layer 2: 20 800.659 MeV 0.667216 8.01749e-09 Gy 50.0432
+
layer 2: 10 206.505 MeV 0.172088 1.65429e-10 Gy 26.2731
-
  layer 3: 30 996.867 MeV 0.830722 9.98223e-09 Gy 62.3066
+
  layer 3: 15 259.584 MeV 0.21632         2.07949e-10 Gy 33.026
-
  layer 4: 40 1.23773 GeV 1.03144         1.23942e-08 Gy 77.3614
+
layer 4: 20 311.898 MeV 0.259915 2.49858e-10 Gy 39.6819
-
  layer 5: 50 1.29582 GeV 1.07985         1.29758e-08 Gy 80.9917
+
  layer 5: 25 395.255 MeV 0.329379 3.16634e-10 Gy 50.2871
-
  layer 6: 60 1.43647 GeV 1.19706         1.43843e-08 Gy 89.7831
+
layer 6: 30 439.052 MeV 0.365876 3.51719e-10 Gy 55.8592
-
  layer 7: 70 1.59565 GeV 1.32971         1.59782e-08 Gy 99.7318
+
  layer 7: 35 503.311 MeV 0.419426 4.03197e-10 Gy 64.0348
-
  layer 8: 80 1.59994 GeV 1.33328         1.60211e-08 Gy 100
+
layer 8: 40 573.742 MeV 0.478118 4.59618e-10 Gy 72.9954
-
  layer 9: 90 1.53239 GeV 1.27699         1.53448e-08 Gy 95.7783
+
  layer 9: 45 624.707 MeV 0.520589 5.00446e-10 Gy 79.4796
-
  layer 10: 100 1.52627 GeV 1.27189         1.52834e-08 Gy 95.3955
+
layer 10: 50 678.383 MeV 0.565319 5.43445e-10 Gy 86.3086
-
  layer 11: 110 1.49511 GeV 1.24592         1.49714e-08 Gy 93.4478
+
  layer 11: 55 694.602 MeV 0.578835 5.56437e-10 Gy 88.3721
-
  layer 12: 120 1.46007 GeV 1.21672         1.46206e-08 Gy 91.2579
+
layer 12: 60 710.771 MeV 0.592309 5.6939e-10 Gy 90.4292
-
  layer 13: 130 1.44794 GeV 1.20662         1.44991e-08 Gy 90.4998
+
  layer 13: 65 744.826 MeV 0.620689 5.96672e-10 Gy 94.762
-
  layer 14: 140 1.40362 GeV 1.16968         1.40553e-08 Gy 87.7295
+
layer 14: 70 742.436 MeV 0.618697 5.94757e-10 Gy 94.4579
-
  layer 15: 150 1.44033 GeV 1.20028         1.44229e-08 Gy 90.0244
+
  layer 15: 75 771.713 MeV 0.643094 6.1821e-10 Gy 98.1827
-
  layer 16: 160 1.42336 GeV 1.18613         1.42529e-08 Gy 88.9634
+
layer 16: 80 767.22 MeV 0.63935         6.14611e-10 Gy 97.611
-
  layer 17: 170 1.28793 GeV 1.07327         1.28968e-08 Gy 80.4985
+
  layer 17: 85 775.608 MeV 0.64634         6.21331e-10 Gy 98.6783
-
  layer 18: 180 1.26477 GeV 1.05397         1.26649e-08 Gy 79.051
+
  layer 18: 90 762.779 MeV 0.635649 6.11053e-10 Gy 97.046
-
  layer 19: 190 1.24786 GeV 1.03989         1.24956e-08 Gy 77.9945
+
layer 19: 95 785.997 MeV 0.654997 6.29653e-10 Gy 100
-
  layer 20: 200 5.30499e-307 eV 4.42082e-316 0 Gy         0
+
layer 20: 100 735.186 MeV 0.612655 5.88949e-10 Gy 93.5355
 +
  layer 21: 105 761.414 MeV 0.634512 6.0996e-10 Gy 96.8724
 +
layer 22: 110 721.836 MeV 0.60153         5.78254e-10 Gy 91.837
 +
  layer 23: 115 730.726 MeV 0.608939 5.85376e-10 Gy 92.9681
 +
layer 24: 120 728.394 MeV 0.606995 5.83508e-10 Gy 92.6714
 +
  layer 25: 125 744.904 MeV 0.620753 5.96734e-10 Gy 94.7719
 +
layer 26: 130 731.53 MeV 0.609608 5.8602e-10 Gy 93.0703
 +
  layer 27: 135 702.85 MeV 0.585709 5.63045e-10 Gy 89.4215
 +
layer 28: 140 671.716 MeV 0.559763 5.38104e-10 Gy 85.4604
 +
layer 29: 145 676.297 MeV 0.56358         5.41773e-10 Gy 86.0432
 +
  layer 30: 150 653.849 MeV 0.544875 5.23791e-10 Gy 83.1873
 +
  layer 31: 155 674.152 MeV 0.561793 5.40055e-10 Gy 85.7703
 +
layer 32: 160 687.008 MeV 0.572507 5.50354e-10 Gy 87.406
 +
layer 33: 165 706.2 MeV 0.5885         5.65729e-10 Gy 89.8478
 +
  layer 34: 170 700.841 MeV 0.584034 5.61435e-10 Gy 89.1659
 +
  layer 35: 175 686.009 MeV 0.571674 5.49553e-10 Gy 87.2788
 +
layer 36: 180 658.891 MeV 0.549076 5.2783e-10 Gy 83.8287
 +
  layer 37: 185 644.999 MeV 0.537499 5.16701e-10 Gy 82.0613
 +
  layer 38: 190 655.81 MeV 0.546508 5.25362e-10 Gy 83.4367
 +
layer 39: 195 645.691 MeV 0.538075 5.17255e-10 Gy 82.1493
 +
layer 40: 200 0 eV         0         0 Gy         0
-
  The run consists of 6000 gamma of 20 MeV through 20 cm  of Water (density: 1 g/cm3 ) divided into 20 slices.
+
  The run consists of 6000 gamma of 20 MeV through 20 cm  of Water (density: 1 g/cm3 )  
 +
divided into 40 slices.
  Edep is the deposited energy in every slice.
  Edep is the deposited energy in every slice.
  Total incident energy(Ebeam)= 120 GeV
  Total incident energy(Ebeam)= 120 GeV
-
  Total energy deposit= 25.169 GeV
+
  Total energy deposit= 24.8344 GeV
  Dose is the deposited dose in every slice.
  Dose is the deposited dose in every slice.
  MaxDose is the highest dose value from all slices.
  MaxDose is the highest dose value from all slices.
Line 111: Line 132:
<pre style="color: #800000; background-color: #dcdcdc">   
<pre style="color: #800000; background-color: #dcdcdc">   
-
10 32.1047
+
5 15.1483
-
20 50.0432
+
10 26.2731
-
30 62.3066
+
15 33.026
-
40 77.3614
+
20 39.6819
-
50 80.9917
+
25 50.2871
-
60 89.7831
+
30 55.8592
-
70 99.7318
+
35 64.0348
-
80 100
+
40 72.9954
-
90 95.7783
+
45 79.4796
-
100 95.3955
+
50 86.3086
-
110 93.4478
+
55 88.3721
-
120 91.2579
+
60 90.4292
-
130 90.4998
+
65 94.762
-
140 87.7295
+
70 94.4579
-
150 90.0244
+
75 98.1827
-
160 88.9634
+
80 97.611
-
170 80.4985
+
85 98.6783
-
180 79.051
+
90 97.046
-
190 77.9945
+
95 100
 +
100 93.5355
 +
105 96.8724
 +
110 91.837
 +
115 92.9681
 +
120 92.6714
 +
125 94.7719
 +
130 93.0703
 +
135 89.4215
 +
140 85.4604
 +
145 86.0432
 +
150 83.1873
 +
155 85.7703
 +
160 87.406
 +
165 89.8478
 +
170 89.1659
 +
175 87.2788
 +
180 83.8287
 +
185 82.0613
 +
190 83.4367
 +
195 82.1493
200 0
200 0
</pre>
</pre>
Line 147: Line 188:
The histogram inside Gamma.root shows the energy deposition in water box:
The histogram inside Gamma.root shows the energy deposition in water box:
   
   
-
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/PhotonPB/Edep_PhotonB.png  
+
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/PhotonPB/Edep_PhotonB1.png  
To exit the ROOT terminal type  
To exit the ROOT terminal type  
Line 182: Line 223:
This is the result:
This is the result:
-
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/PhotonPB/DoseDeposition.png   
+
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/PhotonPB/DoseDeposition1.png   
=== <span style="color:#000080"> Changes in macro gamma.mac </span> ===
=== <span style="color:#000080"> Changes in macro gamma.mac </span> ===
Line 209: Line 250:
/photonPB/det/setSizeX  20 cm
/photonPB/det/setSizeX  20 cm
/photonPB/det/setSizeYZ 20 cm
/photonPB/det/setSizeYZ 20 cm
-
/photonPB/det/setSliceSizeYZ 4 cm
+
/photonPB/det/setSliceSizeYZ 20 cm
-
/photonPB/det/sliceNumber 20
+
/photonPB/det/sliceNumber 40
#
#
# set physics process
# set physics process
/photonPB/phys/addPhysics emstandard_opt0
/photonPB/phys/addPhysics emstandard_opt0
-
#photonPB/phys/addPhysics emlivermore
+
#/photonPB/phys/addPhysics emlivermore
#/photonPB/phys/addPhysics empenelope
#/photonPB/phys/addPhysics empenelope
#
#
Line 227: Line 268:
/run/initialize
/run/initialize
#
#
 +
# visualisation
 +
/control/execute visualisation.mac
 +
# particle gun properties (type of  
# particle gun properties (type of  
#particle and energy)
#particle and energy)

Revision as of 16:05, 24 August 2014

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