Proton beam with realistic geometry

From UCL HEP PBT Wiki

Jump to: navigation, search
m
Line 81: Line 81:
<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.80235 GeV 2.366 1.058e-05 Gy 21.8
+
  layer 1: 0.8 5.20624 GeV 1.399 6.517e-07 Gy 18.49
-
  layer 2: 2.667 8.898 GeV 2.392 1.069e-05 Gy 22.03
+
layer 2: 1.6 5.255 GeV 1.413 6.578e-07 Gy 18.67
-
  layer 3: 4 9.213 GeV 2.476 1.107e-05 Gy 22.81
+
  layer 3: 2.4 5.349 GeV 1.438 6.696e-07 Gy 19
-
  layer 4: 5.333 9.374 GeV 2.52 1.126e-05 Gy 23.21
+
layer 4: 3.2 5.464 GeV 1.469 6.839e-07 Gy 19.41
-
  layer 5: 6.667 9.594 GeV 2.579 1.153e-05 Gy 23.76
+
  layer 5: 4 5.533 GeV 1.487 6.926e-07 Gy 19.66
-
  layer 6: 8 9.91 GeV 2.664 1.191e-05 Gy 24.54
+
  layer 6: 4.8 5.546 GeV 1.491 6.942e-07 Gy 19.7
-
  layer 7: 9.333 10.07 GeV 2.708 1.211e-05 Gy 24.95
+
layer 7: 5.6 5.592 GeV 1.503 7e-07 Gy 19.87
-
  layer 8: 10.67 10.26 GeV 2.758 1.233e-05 Gy 25.41
+
  layer 8: 6.4 5.586 GeV 1.502 6.993e-07 Gy 19.84
-
  layer 9: 12 10.55 GeV 2.835 1.267e-05 Gy 26.11
+
layer 9: 7.2 5.748 GeV 1.545 7.194e-07 Gy 20.42
-
  layer 10: 13.33 10.88 GeV 2.925 1.308e-05 Gy 26.95
+
  layer 10: 8 5.827 GeV 1.566 7.293e-07 Gy 20.7
-
  layer 11: 14.67 11.21 GeV 3.012 1.347e-05 Gy 27.75
+
  layer 11: 8.8 5.922 GeV 1.592 7.412e-07 Gy 21.04
-
  layer 12: 16 11.64 GeV 3.13 1.399e-05 Gy 28.83
+
layer 12: 9.6 5.949 GeV 1.599 7.446e-07 Gy 21.13
-
  layer 13: 17.33 11.96 GeV 3.215 1.437e-05 Gy 29.62
+
  layer 13: 10.4 6.05 GeV 1.626 7.573e-07 Gy 21.49
-
  layer 14: 18.67 12.5 GeV 3.361 1.502e-05 Gy 30.96
+
layer 14: 11.2 6.195 GeV 1.665 7.755e-07 Gy 22.01
-
  layer 15: 20 12.99 GeV 3.492 1.561e-05 Gy 32.17
+
  layer 15: 12 6.306 GeV 1.695 7.894e-07 Gy 22.4
-
  layer 16: 21.33 13.71 GeV 3.686 1.648e-05 Gy 33.95
+
layer 16: 12.8 6.365 GeV 1.711 7.967e-07 Gy 22.61
-
  layer 17: 22.67 14.5 GeV 3.899 1.743e-05 Gy 35.91
+
  layer 17: 13.6 6.499 GeV 1.747 8.135e-07 Gy 23.09
-
  layer 18: 24 15.37 GeV 4.131 1.847e-05 Gy 38.06
+
  layer 18: 14.4 6.634 GeV 1.783 8.303e-07 Gy 23.56
-
  layer 19: 25.33 16.58 GeV 4.457 1.993e-05 Gy 41.06
+
layer 19: 15.2 6.731 GeV 1.809 8.425e-07 Gy 23.91
-
  layer 20: 26.67 18.18 GeV 4.887 2.185e-05 Gy 45.02
+
  layer 20: 16 6.915 GeV 1.859 8.655e-07 Gy 24.56
-
  layer 21: 28 20.45 GeV 5.496 2.457e-05 Gy 50.63
+
layer 21: 16.8 7.099 GeV 1.908 8.886e-07 Gy 25.22
-
  layer 22: 29.33 24.39 GeV 6.555 2.93e-05 Gy 60.38
+
  layer 22: 17.6 7.168 GeV 1.927 8.972e-07 Gy 25.46
-
  layer 23: 30.67 32.75 GeV 8.803 3.935e-05 Gy 81.09
+
  layer 23: 18.4 7.301 GeV 1.963 9.138e-07 Gy 25.93
-
  layer 24: 32 40.39 GeV 10.86 4.853e-05 Gy 100
+
  layer 24: 19.2 7.593 GeV 2.041 9.505e-07 Gy 26.97
-
  layer 25: 33.33 3.931 GeV 1.057 4.724e-06 Gy 9.734
+
layer 25: 20 7.718 GeV 2.075 9.66e-07 Gy 27.42
-
  layer 26: 34.67 1.427 MeV 0.0003836 1.715e-09 Gy 0.003534
+
  layer 26: 20.8 7.952 GeV 2.138 9.954e-07 Gy 28.25
-
  layer 27: 36 903.1 keV 0.0002427 1.085e-09 Gy 0.002236
+
layer 27: 21.6 8.149 GeV 2.191 1.02e-06 Gy 28.95
-
  layer 28: 37.33 7.17 MeV 0.001927 8.616e-09 Gy 0.01775
+
  layer 28: 22.4 8.446 GeV 2.27 1.057e-06 Gy 30
-
  layer 29: 38.67 972.7 keV 0.0002615 1.169e-09 Gy 0.002408
+
layer 29: 23.2 8.784 GeV 2.361 1.1e-06 Gy 31.2
-
  layer 30: 40 0 eV 0 0 Gy 0
+
  layer 30: 24 9.082 GeV 2.441 1.137e-06 Gy 32.26
 +
layer 31: 24.8 9.47 GeV 2.546 1.185e-06 Gy 33.64
 +
  layer 32: 25.6 9.929 GeV 2.669 1.243e-06 Gy 35.27
 +
  layer 33: 26.4 10.54 GeV 2.833 1.319e-06 Gy 37.44
 +
layer 34: 27.2 11.25 GeV 3.023 1.408e-06 Gy 39.95
 +
  layer 35: 28 12.11 GeV 3.256 1.516e-06 Gy 43.03
 +
  layer 36: 28.8 13.21 GeV 3.55 1.653e-06 Gy 46.91
 +
layer 37: 29.6 14.65 GeV 3.938 1.834e-06 Gy 52.05
 +
layer 38: 30.4 16.98 GeV 4.564 2.125e-06 Gy 60.31
 +
  layer 39: 31.2 21.42 GeV 5.759 2.682e-06 Gy 76.1
 +
  layer 40: 32 28.15 GeV 7.567 3.524e-06 Gy 100
 +
  layer 41: 32.8 15.14 GeV 4.07 1.895e-06 Gy 53.78
 +
  layer 42: 33.6 1.256 GeV 0.3376 1.572e-07 Gy 4.462
 +
layer 43: 34.4 18.93 MeV 0.005087  2.369e-09 Gy 0.06723
 +
  layer 44: 35.2 597.5 keV 0.0001606  7.479e-11 Gy 0.002122
 +
layer 45: 36 665.8 keV 0.000179  8.334e-11 Gy 0.002365
 +
layer 46: 36.8 130.1 keV 3.498e-05  1.629e-11 Gy 0.0004623
 +
  layer 47: 37.6 457.2 keV 0.0001229  5.722e-11 Gy 0.001624
 +
  layer 48: 38.4 247.8 keV 6.661e-05  3.102e-11 Gy 0.0008802
 +
layer 49: 39.2 236.6 keV 6.361e-05  2.962e-11 Gy 0.0008406
 +
  layer 50: 40 0 eV 0 0 Gy 0
-
  The run consists of 6000  protons of 62 MeV through 4 cm  of Water (density: 1 g/cm3 ) divided into 30 slices.
+
  The run consists of 6000  protons 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= 366.7 GeV
+
  Total energy deposit= 367.3 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 126: Line 146:
<pre style="color: #800000; background-color: #dcdcdc">
<pre style="color: #800000; background-color: #dcdcdc">
-
1.33333 21.7953
+
0.8 18.4941
-
2.66667 22.033
+
1.6 18.668
-
4 22.8122
+
2.4 19.0022
-
5.33333 23.2101
+
3.2 19.4096
-
6.66667 23.7561
+
4 19.6559
-
8 24.5372
+
4.8 19.7012
-
9.33333 24.9451
+
5.6 19.8661
-
10.6667 25.4071
+
6.4 19.8446
-
12 26.114
+
7.2 20.4168
-
13.3333 26.9454
+
8 20.6987
-
14.6667 27.7468
+
8.8 21.0351
-
16 28.8325
+
9.6 21.1327
-
17.3333 29.6189
+
10.4 21.493
-
18.6667 30.9594
+
11.2 22.0079
-
20 32.1666
+
12 22.4025
-
21.3333 33.9549
+
12.8 22.6104
-
22.6667 35.9146
+
13.6 23.0872
-
24 38.0563
+
14.4 23.5642
-
25.3333 41.0607
+
15.2 23.9102
-
26.6667 45.0231
+
16 24.564
-
28 50.6268
+
16.8 25.2168
-
29.3333 60.3833
+
17.6 25.4611
-
30.6667 81.0937
+
18.4 25.9348
 +
19.2 26.9738
 +
20 27.4152
 +
20.8 28.2496
 +
21.6 28.9492
 +
22.4 30.0033
 +
23.2 31.2049
 +
24 32.2636
 +
24.8 33.6414
 +
25.6 35.2707
 +
26.4 37.4412
 +
27.2 39.9539
 +
28 43.0264
 +
28.8 46.9139
 +
29.6 52.0461
 +
30.4 60.3067
 +
31.2 76.1028
32 100
32 100
-
33.3333 9.73398
+
32.8 53.7805
-
34.6667 0.00353384
+
33.6 4.46152
-
36 0.00223614
+
34.4 0.0672278
-
37.3333 0.0177538
+
35.2 0.00212242
-
38.6667 0.00240845
+
36 0.00236509
 +
36.8 0.000462313
 +
37.6 0.00162397
 +
38.4 0.00088024
 +
39.2 0.000840567
40 0
40 0
</pre>  
</pre>  
-
[http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/DoseLongitudinalMesh.txt '''DoseLongitudinalMesh.txt'''] and [http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/EnergyLongitudinalMesh.txt '''EnergyLongitudinalMesh.txt'''] contain information about the dose and energy deposition in 30 voxels along the beam.   
+
[http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/DoseLongitudinalMesh.txt '''DoseLongitudinalMesh.txt'''] and [http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/EnergyLongitudinalMesh.txt '''EnergyLongitudinalMesh.txt'''] contain information about the dose and energy deposition in 50 voxels along the beam.   
-
[http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/DoseLateralMesh.txt '''DoseLateralMesh.txt'''] and [http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/EnergyLateralMesh.txt '''EnergyLateralMesh.txt'''] contain information about the dose and energy deposition in 30 voxels in direction perpendicular to the beam at its peak location along the beam.  
+
[http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/DoseLateralMesh.txt '''DoseLateralMesh.txt'''] and [http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/EnergyLateralMesh.txt '''EnergyLateralMesh.txt'''] contain information about the dose and energy deposition in 50 voxels in direction perpendicular to the beam at its peak location along the beam.  
=== <span style="color:#000080"> Root file </span> ===
=== <span style="color:#000080"> Root file </span> ===
Line 192: Line 232:
This will create root file with the following plot:
This will create root file with the following plot:
-
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/Simulation.png   
+
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/ProtonGB/Simulation1.png   
=== <span style="color:#000080"> Changes in proton.mac </span> ===
=== <span style="color:#000080"> Changes in proton.mac </span> ===
Line 219: Line 259:
/protonGB/det/setSizeX  4 cm
/protonGB/det/setSizeX  4 cm
/protonGB/det/setSizeYZ 4 cm
/protonGB/det/setSizeYZ 4 cm
-
/protonGB/det/setSliceSizeYZ 1 cm
+
/protonGB/det/setSliceSizeYZ 4 cm
-
/protonGB/det/sliceNumber 30
+
/protonGB/det/sliceNumber 50
#
#
# define longitudinal scoring mesh
# define longitudinal scoring mesh
Line 226: Line 266:
/score/create/boxMesh waterMeshlongitudinal
/score/create/boxMesh waterMeshlongitudinal
/score/mesh/boxSize 2. 2. 2. cm
/score/mesh/boxSize 2. 2. 2. cm
-
/score/mesh/nBin 30 1 1
+
/score/mesh/nBin 50 1 1
/score/mesh/translate/xyz 0. 0. 0. cm
/score/mesh/translate/xyz 0. 0. 0. cm
/score/quantity/energyDeposit energyDeposit  
/score/quantity/energyDeposit energyDeposit  
Line 235: Line 275:
# centered at the Bragg peak
# centered at the Bragg peak
/score/create/boxMesh waterMeshlateral
/score/create/boxMesh waterMeshlateral
-
/score/mesh/boxSize 1. 2. 2. cm
+
/score/mesh/boxSize 0.1 2. 2. cm
-
/score/mesh/nBin 1 30 1
+
/score/mesh/nBin 1 1 50
/score/mesh/translate/xyz 1.2 0. 0. cm
/score/mesh/translate/xyz 1.2 0. 0. cm
/score/quantity/energyDeposit energyDeposit  
/score/quantity/energyDeposit energyDeposit  
Line 256: Line 296:
# initialize
# initialize
/run/initialize
/run/initialize
 +
#
 +
# visualisation
 +
#/control/execute visualisation.mac
 +
#
# General particle source
# General particle source
# proton circle source   
# proton circle source   
/gps/pos/shape Circle
/gps/pos/shape Circle
-
/gps/pos/centre -4. 0. 0. cm
+
/gps/pos/centre -2. 0. 0. cm
/gps/pos/radius 0. mm
/gps/pos/radius 0. mm
/gps/pos/sigma_r 2. mm
/gps/pos/sigma_r 2. mm
Line 294: Line 338:
/run/beamOn 6000
/run/beamOn 6000
#
#
 +
# drawing projections
 +
#/score/drawProjection waterMeshlongitudinal doseDeposit
 +
#/score/drawProjection waterMeshlateral doseDeposit
 +
# dump scores to a file
# dump scores to a file
/score/dumpQuantityToFile waterMeshlongitudinal doseDeposit DoseLongitudinalMesh.txt
/score/dumpQuantityToFile waterMeshlongitudinal doseDeposit DoseLongitudinalMesh.txt
Line 302: Line 350:
You can modify '''proton.mac''' as it is done in the tutorial '''Monoenergetic proton pencil beam'''.   
You can modify '''proton.mac''' as it is done in the tutorial '''Monoenergetic proton pencil beam'''.   
 +
In addition to that you can modify 
=== <span style="color:#000080"> Visualisation </span> ===
=== <span style="color:#000080"> Visualisation </span> ===
-
Under construction (problems with OpenGL installed on plus1)
+
 
== <span style="color:#000080"> Comparison with data from The Clatterbridge Cancer Centre </span> ==
== <span style="color:#000080"> Comparison with data from The Clatterbridge Cancer Centre </span> ==

Revision as of 11:33, 24 August 2014

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