Monoenergetic photon pencil beam
From UCL HEP PBT Wiki
m |
m |
||
Line 76: | Line 76: | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
- | Layers : x[mm] Edep | + | Layers : x[mm] Edep Edep/Ebeam[%] Dose Dose/MaxDose[%] |
- | layer 1: 10 | + | layer 1: 5 119.065 MeV 0.0992211 9.53818e-11 Gy 15.1483 |
- | layer 2: 20 | + | layer 2: 10 206.505 MeV 0.172088 1.65429e-10 Gy 26.2731 |
- | layer 3: 30 | + | layer 3: 15 259.584 MeV 0.21632 2.07949e-10 Gy 33.026 |
- | layer 4: 40 | + | layer 4: 20 311.898 MeV 0.259915 2.49858e-10 Gy 39.6819 |
- | layer 5: 50 | + | layer 5: 25 395.255 MeV 0.329379 3.16634e-10 Gy 50.2871 |
- | layer | + | layer 6: 30 439.052 MeV 0.365876 3.51719e-10 Gy 55.8592 |
- | layer | + | layer 7: 35 503.311 MeV 0.419426 4.03197e-10 Gy 64.0348 |
- | layer | + | layer 8: 40 573.742 MeV 0.478118 4.59618e-10 Gy 72.9954 |
- | layer | + | layer 9: 45 624.707 MeV 0.520589 5.00446e-10 Gy 79.4796 |
- | layer 10: 100 | + | layer 10: 50 678.383 MeV 0.565319 5.43445e-10 Gy 86.3086 |
- | layer | + | layer 11: 55 694.602 MeV 0.578835 5.56437e-10 Gy 88.3721 |
- | layer | + | layer 12: 60 710.771 MeV 0.592309 5.6939e-10 Gy 90.4292 |
- | layer | + | layer 13: 65 744.826 MeV 0.620689 5.96672e-10 Gy 94.762 |
- | layer | + | layer 14: 70 742.436 MeV 0.618697 5.94757e-10 Gy 94.4579 |
- | layer | + | layer 15: 75 771.713 MeV 0.643094 6.1821e-10 Gy 98.1827 |
- | layer | + | layer 16: 80 767.22 MeV 0.63935 6.14611e-10 Gy 97.611 |
- | layer | + | layer 17: 85 775.608 MeV 0.64634 6.21331e-10 Gy 98.6783 |
- | layer | + | layer 18: 90 762.779 MeV 0.635649 6.11053e-10 Gy 97.046 |
- | layer | + | layer 19: 95 785.997 MeV 0.654997 6.29653e-10 Gy 100 |
- | layer | + | 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 | + | 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= | + | 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 | + | 5 15.1483 |
- | 20 50. | + | 10 26.2731 |
- | 30 | + | 15 33.026 |
- | 40 | + | 20 39.6819 |
- | 50 | + | 25 50.2871 |
- | 60 | + | 30 55.8592 |
- | 70 | + | 35 64.0348 |
- | 80 | + | 40 72.9954 |
- | 90 | + | 45 79.4796 |
- | + | 50 86.3086 | |
- | 110 | + | 55 88.3721 |
- | 120 | + | 60 90.4292 |
- | 130 | + | 65 94.762 |
- | 140 | + | 70 94.4579 |
- | 150 | + | 75 98.1827 |
- | 160 | + | 80 97.611 |
- | 170 | + | 85 98.6783 |
- | 180 | + | 90 97.046 |
- | 190 | + | 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/ | + | 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/ | + | 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 | + | /photonPB/det/setSliceSizeYZ 20 cm |
- | /photonPB/det/sliceNumber | + | /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) |