Monoenergetic photon pencil beam
From UCL HEP PBT Wiki
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 | + | The macro produces root file '''Gamma.root''' with a histogram showing the energy deposition in |
| - | water box along the beam line. It also produces | + | water box along the beam line. It also produces text files: '''DoseFile.txt''' with |
| - | energy and dose deposited | + | energy and dose deposited in each slice and '''PlotDose.txt''' with dose deposited in each slice. |
| - | with the default settings: physics process emstandard_opt0 and incident gamma energy of 20 MeV. | + | 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"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
| - | + | Layers : x[mm] Edep Edep/Ebeam[%] Dose Dose/MaxDose[%] | |
| - | layer 1: | + | layer 1: 10 513.656 MeV 0.428046 5.14354e-09 Gy 32.1047 |
| - | layer 2: | + | layer 2: 20 800.659 MeV 0.667216 8.01749e-09 Gy 50.0432 |
| - | layer 3: | + | layer 3: 30 996.867 MeV 0.830722 9.98223e-09 Gy 62.3066 |
| - | layer 4: | + | layer 4: 40 1.23773 GeV 1.03144 1.23942e-08 Gy 77.3614 |
| - | layer 5: | + | layer 5: 50 1.29582 GeV 1.07985 1.29758e-08 Gy 80.9917 |
| - | layer 6: | + | layer 6: 60 1.43647 GeV 1.19706 1.43843e-08 Gy 89.7831 |
| - | layer 7: | + | layer 7: 70 1.59565 GeV 1.32971 1.59782e-08 Gy 99.7318 |
| - | layer 8: | + | layer 8: 80 1.59994 GeV 1.33328 1.60211e-08 Gy 100 |
| - | layer 9: | + | layer 9: 90 1.53239 GeV 1.27699 1.53448e-08 Gy 95.7783 |
| - | layer 10: | + | layer 10: 100 1.52627 GeV 1.27189 1.52834e-08 Gy 95.3955 |
| - | layer 11: | + | layer 11: 110 1.49511 GeV 1.24592 1.49714e-08 Gy 93.4478 |
| - | layer 12: | + | layer 12: 120 1.46007 GeV 1.21672 1.46206e-08 Gy 91.2579 |
| - | layer 13: | + | layer 13: 130 1.44794 GeV 1.20662 1.44991e-08 Gy 90.4998 |
| - | layer 14: | + | layer 14: 140 1.40362 GeV 1.16968 1.40553e-08 Gy 87.7295 |
| - | layer 15: | + | layer 15: 150 1.44033 GeV 1.20028 1.44229e-08 Gy 90.0244 |
| - | layer 16: | + | layer 16: 160 1.42336 GeV 1.18613 1.42529e-08 Gy 88.9634 |
| - | layer 17: | + | layer 17: 170 1.28793 GeV 1.07327 1.28968e-08 Gy 80.4985 |
| - | layer 18: | + | layer 18: 180 1.26477 GeV 1.05397 1.26649e-08 Gy 79.051 |
| - | layer 19: | + | layer 19: 190 1.24786 GeV 1.03989 1.24956e-08 Gy 77.9945 |
| - | layer 20: | + | layer 20: 200 5.30499e-307 eV 4.42082e-316 0 Gy 0 |
| - | The run consists of | + | The run consists of 6000 gamma of 20 MeV through 20 cm of Water (density: 1 g/cm3 ) divided into 20 slices. |
| + | |||
| + | Edep is the deposited energy in every slice. | ||
| + | Total incident energy(Ebeam)= 120 GeV | ||
| + | Total energy deposit= 25.169 GeV | ||
| + | Dose is the deposited dose in every slice. | ||
| + | MaxDose is the highest dose value from all slices. | ||
</pre> | </pre> | ||
| - | + | Open Gamma.root file in the following way: | |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
| Line 109: | Line 115: | ||
Inside is 1D histogram showing the energy deposition in water box: | Inside is 1D histogram showing the energy deposition in water box: | ||
| - | http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/PhotonPB/ | + | http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/PhotonPB/Edep_PhotonB.png |
You can change the physics process, incident photon energy, phantom material, number of slices etc. by | You can change the physics process, incident photon energy, phantom material, number of slices etc. by | ||
