Monoenergetic proton pencil beam

From UCL HEP PBT Wiki

Jump to: navigation, search
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 : x[mm]   Edep     Edep/Ebeam[%]  Dose Dose/MaxDose[%]
+
  Layers : x[mm] Edep     Edep/Ebeam[%]  Dose       Dose/MaxDose[%]
-
  layer 1: 1.33333 8.84429 GeV 2.3775 1.06276e-05 Gy 20.6613
+
  layer 1: 0.8 5.17382 GeV 1.39081 6.47608e-07 Gy 17.3315
-
  layer 2: 2.66667 8.91694 GeV 2.39703 1.07149e-05 Gy 20.831
+
layer 2: 1.6 5.24175 GeV 1.40907 6.5611e-07 Gy 17.5591
-
  layer 3: 4 9.17299 GeV 2.46586 1.10226e-05 Gy 21.4291
+
  layer 3: 2.4 5.29862 GeV 1.42436 6.63229e-07 Gy 17.7496
-
  layer 4: 5.33333 9.38711 GeV 2.52342 1.12799e-05 Gy 21.9294
+
layer 4: 3.2 5.41481 GeV 1.45559 6.77772e-07 Gy 18.1388
-
  layer 5: 6.66667 9.62373 GeV 2.58702 1.15642e-05 Gy 22.4821
+
  layer 5: 4 5.43787 GeV 1.46179 6.80658e-07 Gy 18.216
-
  layer 6: 8 9.75588 GeV 2.62255 1.1723e-05 Gy 22.7908
+
  layer 6: 4.8 5.52101 GeV 1.48414 6.91065e-07 Gy 18.4946
-
  layer 7: 9.33333 9.98464 GeV 2.68404 1.19979e-05 Gy 23.3252
+
layer 7: 5.6 5.64633 GeV 1.51783 7.06751e-07 Gy 18.9144
-
  layer 8: 10.6667 10.3128 GeV 2.77225 1.23922e-05 Gy 24.0918
+
layer 8: 6.4 5.63969 GeV 1.51605 7.0592e-07 Gy 18.8921
-
  layer 9: 12 10.5176 GeV 2.82732 1.26383e-05 Gy 24.5704
+
layer 9: 7.2 5.71744 GeV 1.53695 7.15652e-07 Gy 19.1526
-
  layer 10: 13.3333 10.8275 GeV 2.91061 1.30107e-05 Gy 25.2942
+
  layer 10: 8 5.78086 GeV 1.55399 7.2359e-07 Gy 19.365
-
  layer 11: 14.6667 11.164 GeV 3.00107 1.3415e-05 Gy 26.0803
+
layer 11: 8.8 5.94371 GeV 1.59777 7.43975e-07 Gy 19.9106
-
  layer 12: 16 11.6363 GeV 3.12804 1.39825e-05 Gy 27.1837
+
layer 12: 9.6 6.02518 GeV 1.61967 7.54172e-07 Gy 20.1835
-
  layer 13: 17.3333 12.0232 GeV 3.23203 1.44474e-05 Gy 28.0874
+
layer 13: 10.4 6.10292 GeV 1.64057 7.63903e-07 Gy 20.4439
-
  layer 14: 18.6667 12.6138 GeV 3.3908 1.51571e-05 Gy 29.4672
+
layer 14: 11.2 6.18071 GeV 1.66148 7.7364e-07 Gy 20.7045
-
  layer 15: 20 12.9498 GeV 3.48112 1.55609e-05 Gy 30.2521
+
  layer 15: 12 6.2621 GeV 1.68336 7.83827e-07 Gy 20.9771
-
  layer 16: 21.3333 13.6288 GeV 3.66365 1.63768e-05 Gy 31.8383
+
layer 16: 12.8 6.37762 GeV 1.71441 7.98286e-07 Gy 21.3641
-
  layer 17: 22.6667 14.5047 GeV 3.89912 1.74294e-05 Gy 33.8847
+
layer 17: 13.6 6.52458 GeV 1.75392 8.16682e-07 Gy 21.8564
-
  layer 18: 24 15.4695 GeV 4.15847 1.85887e-05 Gy 36.1385
+
  layer 18: 14.4 6.67805 GeV 1.79517 8.35892e-07 Gy 22.3705
-
  layer 19: 25.3333 16.5838 GeV 4.45801 1.99276e-05 Gy 38.7416
+
layer 19: 15.2 6.86252 GeV 1.84476 8.58982e-07 Gy 22.9884
-
  layer 20: 26.6667 18.3271 GeV 4.92663 2.20224e-05 Gy 42.8141
+
  layer 20: 16 6.95226 GeV 1.86889 8.70215e-07 Gy 23.289
-
  layer 21: 28 20.6828 GeV 5.5599 2.48532e-05 Gy 48.3174
+
layer 21: 16.8 7.11679 GeV 1.91312 8.90809e-07 Gy 23.8402
-
  layer 22: 29.3333 24.4737 GeV 6.57896 2.94084e-05 Gy 57.1734
+
  layer 22: 17.6 7.15125 GeV 1.92238 8.95122e-07 Gy 23.9556
-
  layer 23: 30.6667 32.4582 GeV 8.72532 3.90028e-05 Gy 75.826
+
layer 23: 18.4 7.43274 GeV 1.99805 9.30356e-07 Gy 24.8986
-
  layer 24: 32 42.8062 GeV 11.507 5.14373e-05 Gy 100
+
  layer 24: 19.2 7.58811 GeV 2.03981 9.49804e-07 Gy 25.419
-
  layer 25: 33.3333 1.99409 GeV 0.536046     2.39617e-06 Gy 4.65842
+
  layer 25: 20 7.8156 GeV 2.10097 9.78279e-07 Gy 26.1811
-
  layer 26: 34.6667 1.01076 MeV 0.000271708 1.21456e-09 Gy 0.00236124
+
layer 26: 20.8 7.94754 GeV 2.13643 9.94794e-07 Gy 26.6231
-
  layer 27: 36 503.954 keV 0.000135472 6.05568e-10 Gy 0.00117729
+
  layer 27: 21.6 8.35363 GeV 2.2456 1.04562e-06 Gy 27.9834
-
  layer 28: 37.3333 33.8047 keV 9.08729e-06 4.06209e-11 Gy 7.89716e-05
+
  layer 28: 22.4 8.44564 GeV 2.27033 1.05714e-06 Gy 28.2917
-
  layer 29: 38.6667 2.03563 MeV 0.000547213 2.44608e-09 Gy 0.00475547
+
  layer 29: 23.2 8.74817 GeV 2.35166 1.09501e-06 Gy 29.3051
-
  layer 30: 40 0 eV         0     0 Gy         0
+
  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 30 slices.
+
  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.244 GeV
+
  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.33333 20.6613
+
0.8 17.3315
-
2.66667 20.831
+
1.6 17.5591
-
4 21.4291
+
2.4 17.7496
-
5.33333 21.9294
+
3.2 18.1388
-
6.66667 22.4821
+
4 18.216
-
8 22.7908
+
4.8 18.4946
-
9.33333 23.3252
+
5.6 18.9144
-
10.6667 24.0918
+
6.4 18.8921
-
12 24.5704
+
7.2 19.1526
-
13.3333 25.2942
+
8 19.365
-
14.6667 26.0803
+
8.8 19.9106
-
16 27.1837
+
9.6 20.1835
-
17.3333 28.0874
+
10.4 20.4439
-
18.6667 29.4672
+
11.2 20.7045
-
20 30.2521
+
12 20.9771
-
21.3333 31.8383
+
12.8 21.3641
-
22.6667 33.8847
+
13.6 21.8564
-
24 36.1385
+
14.4 22.3705
-
25.3333 38.7416
+
15.2 22.9884
-
26.6667 42.8141
+
16 23.289
-
28 48.3174
+
16.8 23.8402
-
29.3333 57.1734
+
17.6 23.9556
-
30.6667 75.826
+
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.3333 4.65842
+
32.8 49.9866
-
34.6667 0.00236124
+
33.6 1.10547
-
36 0.00117729
+
34.4 0.00979044
-
37.3333 7.89716e-05
+
35.2 0.00289566
-
38.6667 0.00475547
+
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 (from the Clatterbridge Cancer Center) and simulation ('''PlotDose.txt''').  
+
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.  
   
   
-
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/ClatterbridgeGeant4.png
+
<pre style="color: #800000; background-color: #dcdcdc">
 +
cp /home/username/ProtonPBFolder/ProtonPB_source/RootScripts/PlotDataAndSim.C .
-
You can also compare data (from the Clatterbridge Cancer Center) with simulation done with scoring mesh ('''DoseLongitudinalMesh.txt''').
+
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. This is the result:
+
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> 
-
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonPB/ClatterbridgeGeant4Mesh.png
+
This creates '''BraggPeakComparison_Mesh.root''' file with the following plot:
-
There is discrepancy between the data in '''PlotDose.txt''' and '''DoseLongitudinalMesh.txt''' (two different ways of retrieving the dose deposition). Probably wrong position of the scoring mesh. This is will be fixed after fixing the problem with the visualisation.
+
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]]

Revision as of 15:22, 23 August 2014

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