Monoenergetic proton pencil beam
From UCL HEP PBT Wiki
m |
m |
||
Line 78: | Line 78: | ||
pico DoseFile.txt | pico DoseFile.txt | ||
</pre> | </pre> | ||
- | + | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
- | Layers : | + | Layers : x[mm] Edep Edep/Ebeam[%] Dose Dose/MaxDose[%] |
- | layer 1: | + | layer 1: 0.8 5.17382 GeV 1.39081 6.47608e-07 Gy 17.3315 |
- | layer | + | layer 2: 1.6 5.24175 GeV 1.40907 6.5611e-07 Gy 17.5591 |
- | layer | + | layer 3: 2.4 5.29862 GeV 1.42436 6.63229e-07 Gy 17.7496 |
- | layer 4: 5. | + | layer 4: 3.2 5.41481 GeV 1.45559 6.77772e-07 Gy 18.1388 |
- | layer 5: | + | layer 5: 4 5.43787 GeV 1.46179 6.80658e-07 Gy 18.216 |
- | layer 6: 8 | + | layer 6: 4.8 5.52101 GeV 1.48414 6.91065e-07 Gy 18.4946 |
- | layer | + | layer 7: 5.6 5.64633 GeV 1.51783 7.06751e-07 Gy 18.9144 |
- | layer 8: | + | layer 8: 6.4 5.63969 GeV 1.51605 7.0592e-07 Gy 18.8921 |
- | layer | + | layer 9: 7.2 5.71744 GeV 1.53695 7.15652e-07 Gy 19.1526 |
- | layer | + | layer 10: 8 5.78086 GeV 1.55399 7.2359e-07 Gy 19.365 |
- | layer | + | layer 11: 8.8 5.94371 GeV 1.59777 7.43975e-07 Gy 19.9106 |
- | layer | + | layer 12: 9.6 6.02518 GeV 1.61967 7.54172e-07 Gy 20.1835 |
- | layer | + | layer 13: 10.4 6.10292 GeV 1.64057 7.63903e-07 Gy 20.4439 |
- | layer | + | layer 14: 11.2 6.18071 GeV 1.66148 7.7364e-07 Gy 20.7045 |
- | layer | + | layer 15: 12 6.2621 GeV 1.68336 7.83827e-07 Gy 20.9771 |
- | layer | + | layer 16: 12.8 6.37762 GeV 1.71441 7.98286e-07 Gy 21.3641 |
- | layer | + | layer 17: 13.6 6.52458 GeV 1.75392 8.16682e-07 Gy 21.8564 |
- | layer | + | layer 18: 14.4 6.67805 GeV 1.79517 8.35892e-07 Gy 22.3705 |
- | layer | + | layer 19: 15.2 6.86252 GeV 1.84476 8.58982e-07 Gy 22.9884 |
- | layer | + | layer 20: 16 6.95226 GeV 1.86889 8.70215e-07 Gy 23.289 |
- | layer | + | layer 21: 16.8 7.11679 GeV 1.91312 8.90809e-07 Gy 23.8402 |
- | layer | + | layer 22: 17.6 7.15125 GeV 1.92238 8.95122e-07 Gy 23.9556 |
- | layer | + | layer 23: 18.4 7.43274 GeV 1.99805 9.30356e-07 Gy 24.8986 |
- | layer | + | layer 24: 19.2 7.58811 GeV 2.03981 9.49804e-07 Gy 25.419 |
- | layer | + | layer 25: 20 7.8156 GeV 2.10097 9.78279e-07 Gy 26.1811 |
- | layer | + | layer 26: 20.8 7.94754 GeV 2.13643 9.94794e-07 Gy 26.6231 |
- | layer | + | layer 27: 21.6 8.35363 GeV 2.2456 1.04562e-06 Gy 27.9834 |
- | layer | + | layer 28: 22.4 8.44564 GeV 2.27033 1.05714e-06 Gy 28.2917 |
- | layer | + | layer 29: 23.2 8.74817 GeV 2.35166 1.09501e-06 Gy 29.3051 |
- | layer | + | layer 30: 24 9.08194 GeV 2.44138 1.13679e-06 Gy 30.4232 |
+ | layer 31: 24.8 9.50886 GeV 2.55615 1.19022e-06 Gy 31.8533 | ||
+ | layer 32: 25.6 9.93302 GeV 2.67017 1.24332e-06 Gy 33.2741 | ||
+ | layer 33: 26.4 10.5627 GeV 2.83943 1.32213e-06 Gy 35.3834 | ||
+ | layer 34: 27.2 11.1563 GeV 2.99902 1.39644e-06 Gy 37.3721 | ||
+ | layer 35: 28 12.0025 GeV 3.22647 1.50235e-06 Gy 40.2065 | ||
+ | layer 36: 28.8 13.1124 GeV 3.52485 1.64128e-06 Gy 43.9247 | ||
+ | layer 37: 29.6 14.5158 GeV 3.90209 1.81694e-06 Gy 48.6257 | ||
+ | layer 38: 30.4 16.8932 GeV 4.54119 2.11452e-06 Gy 56.5898 | ||
+ | layer 39: 31.2 20.9827 GeV 5.64051 2.62641e-06 Gy 70.289 | ||
+ | layer 40: 32 29.8521 GeV 8.02475 3.73658e-06 Gy 100 | ||
+ | layer 41: 32.8 14.922 GeV 4.0113 1.86779e-06 Gy 49.9866 | ||
+ | layer 42: 33.6 330.005 MeV 0.088711 4.13067e-08 Gy 1.10547 | ||
+ | layer 43: 34.4 2.92265 MeV 0.000785658 3.65828e-10 Gy 0.00979044 | ||
+ | layer 44: 35.2 864.413 keV 0.000232369 1.08199e-10 Gy 0.00289566 | ||
+ | layer 45: 36 673.958 keV 0.000181172 8.43594e-11 Gy 0.00225766 | ||
+ | layer 46: 36.8 1.18851 MeV 0.000319491 1.48765e-10 Gy 0.00398132 | ||
+ | layer 47: 37.6 3.81319 MeV 0.00102505 4.77297e-10 Gy 0.0127736 | ||
+ | layer 48: 38.4 4.58219 MeV 0.00123177 5.73552e-10 Gy 0.0153496 | ||
+ | layer 49: 39.2 7.29449 MeV 0.00196088 9.13052e-10 Gy 0.0244355 | ||
+ | layer 50: 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 | + | The run consists of 6000 proton of 62 MeV through 4 cm of Water (density: 1 g/cm3 ) divided into 50 slices. |
Edep is the deposited energy in every slice. | Edep is the deposited energy in every slice. | ||
Total incident energy(Ebeam)= 372 GeV | Total incident energy(Ebeam)= 372 GeV | ||
- | Total energy deposit= 367. | + | Total energy deposit= 367.368 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 125: | Line 145: | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
- | 1. | + | 0.8 17.3315 |
- | 2. | + | 1.6 17.5591 |
- | 4 | + | 2.4 17.7496 |
- | 5. | + | 3.2 18.1388 |
- | 6. | + | 4 18.216 |
- | 8 | + | 4.8 18.4946 |
- | 9. | + | 5.6 18.9144 |
- | 10. | + | 6.4 18.8921 |
- | 12 | + | 7.2 19.1526 |
- | 13. | + | 8 19.365 |
- | 14. | + | 8.8 19.9106 |
- | 16 | + | 9.6 20.1835 |
- | 17. | + | 10.4 20.4439 |
- | 18. | + | 11.2 20.7045 |
- | 20 | + | 12 20.9771 |
- | 21. | + | 12.8 21.3641 |
- | 22. | + | 13.6 21.8564 |
- | 24 | + | 14.4 22.3705 |
- | 25. | + | 15.2 22.9884 |
- | 26. | + | 16 23.289 |
- | 28 | + | 16.8 23.8402 |
- | 29. | + | 17.6 23.9556 |
- | 30. | + | 18.4 24.8986 |
+ | 19.2 25.419 | ||
+ | 20 26.1811 | ||
+ | 20.8 26.6231 | ||
+ | 21.6 27.9834 | ||
+ | 22.4 28.2917 | ||
+ | 23.2 29.3051 | ||
+ | 24 30.4232 | ||
+ | 24.8 31.8533 | ||
+ | 25.6 33.2741 | ||
+ | 26.4 35.3834 | ||
+ | 27.2 37.3721 | ||
+ | 28 40.2065 | ||
+ | 28.8 43.9247 | ||
+ | 29.6 48.6257 | ||
+ | 30.4 56.5898 | ||
+ | 31.2 70.289 | ||
32 100 | 32 100 | ||
- | 33. | + | 32.8 49.9866 |
- | 34. | + | 33.6 1.10547 |
- | 36 0. | + | 34.4 0.00979044 |
- | 37. | + | 35.2 0.00289566 |
- | 38. | + | 36 0.00225766 |
- | 40 0 | + | 36.8 0.00398132 |
+ | 37.6 0.0127736 | ||
+ | 38.4 0.0153496 | ||
+ | 39.2 0.0244355 | ||
+ | 40 0 | ||
</pre> | </pre> | ||
Line 401: | Line 441: | ||
The scripts '''PlotDataAndSim.C''' and '''PlotDataAndSimMesh.C''' in folder '''RootScripts''' compare data with simulation. | The scripts '''PlotDataAndSim.C''' and '''PlotDataAndSimMesh.C''' in folder '''RootScripts''' compare data with simulation. | ||
- | You can use script '''PlotDataAndSim.C''' to compare data ( | + | You can use script '''PlotDataAndSim.C''' to compare data (ClatterbridgeData.txt) and simulation ('''PlotDose.txt'''). Both text files must be in the folder where you run the script. |
- | + | <pre style="color: #800000; background-color: #dcdcdc"> | |
+ | cp /home/username/ProtonPBFolder/ProtonPB_source/RootScripts/PlotDataAndSim.C . | ||
- | You can also compare data ( | + | root -l |
+ | |||
+ | .x PlotDataAndSim.C | ||
+ | </pre> | ||
+ | |||
+ | This will create '''BraggPeakComparison.root''' file with the following plot: | ||
+ | |||
+ | http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/DataSimulation.png | ||
+ | |||
+ | You can also compare data (ClatterbridgeData.txt) with simulation done with scoring mesh ('''DoseLongitudinalMesh.txt'''). | ||
This can be done with script '''PlotDataAndSimMesh.C'''. This script works only if before running it you substitude the commas | This can be done with script '''PlotDataAndSimMesh.C'''. This script works only if before running it you substitude the commas | ||
- | in '''DoseLongitudinalMesh.txt''' with spaces. Remove also the header in the text file. | + | in '''DoseLongitudinalMesh.txt''' with spaces. Remove also the header in the text file. |
+ | |||
+ | <pre style="color: #800000; background-color: #dcdcdc"> | ||
+ | cp /home/username/ProtonPBFolder/ProtonPB_source/RootScripts/PlotDataAndSimMesh.C . | ||
+ | |||
+ | root -l | ||
+ | |||
+ | .x PlotDataAndSimMesh.C | ||
+ | </pre> | ||
- | + | This creates '''BraggPeakComparison_Mesh.root''' file with the following plot: | |
- | + | http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/DataSimulationMesh.png | |
== <span style="color:#000080"> Files </span> == | == <span style="color:#000080"> Files </span> == | ||
[[List of monoenergetic proton pencil beam files with brief description]] | [[List of monoenergetic proton pencil beam files with brief description]] |