Monoenergetic proton pencil beam

From UCL HEP PBT Wiki

Jump to: navigation, search
m
m
Line 70: Line 70:
<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.7788 GeV 2.35989 1.05489e-05 Gy 20.404
+
  layer 1: 1.33333 8.84429 GeV 2.3775 1.06276e-05 Gy 20.6613
-
  layer 2: 2.66667 8.95366 GeV 2.4069 1.0759e-05 Gy 20.8104
+
  layer 2: 2.66667 8.91694 GeV 2.39703 1.07149e-05 Gy 20.831
-
  layer 3: 4 9.20046 GeV 2.47324 1.10556e-05 Gy 21.384
+
  layer 3: 4 9.17299 GeV 2.46586 1.10226e-05 Gy 21.4291
-
  layer 4: 5.33333 9.33632 GeV 2.50976 1.12188e-05 Gy 21.6998
+
  layer 4: 5.33333 9.38711 GeV 2.52342 1.12799e-05 Gy 21.9294
-
  layer 5: 6.66667 9.58072 GeV 2.57546 1.15125e-05 Gy 22.2678
+
  layer 5: 6.66667 9.62373 GeV 2.58702 1.15642e-05 Gy 22.4821
-
  layer 6: 8 9.76978 GeV 2.62628 1.17397e-05 Gy 22.7073
+
  layer 6: 8 9.75588 GeV 2.62255 1.1723e-05 Gy 22.7908
-
  layer 7: 9.33333 10.0149 GeV 2.69217 1.20342e-05 Gy 23.277
+
  layer 7: 9.33333 9.98464 GeV 2.68404 1.19979e-05 Gy 23.3252
-
  layer 8: 10.6667 10.248 GeV 2.75484 1.23143e-05 Gy 23.8187
+
  layer 8: 10.6667 10.3128 GeV 2.77225 1.23922e-05 Gy 24.0918
-
  layer 9: 12 10.5163 GeV 2.82695 1.26367e-05 Gy 24.4422
+
  layer 9: 12 10.5176 GeV 2.82732 1.26383e-05 Gy 24.5704
-
  layer 10: 13.3333 10.7679 GeV 2.89459 1.2939e-05 Gy 25.0271
+
  layer 10: 13.3333 10.8275 GeV 2.91061 1.30107e-05 Gy 25.2942
-
  layer 11: 14.6667 11.2918 GeV 3.03543 1.35686e-05 Gy 26.2448
+
  layer 11: 14.6667 11.164 GeV 3.00107 1.3415e-05 Gy 26.0803
-
  layer 12: 16 11.5998 GeV 3.11822 1.39387e-05 Gy 26.9606
+
  layer 12: 16 11.6363 GeV 3.12804 1.39825e-05 Gy 27.1837
-
  layer 13: 17.3333 11.9594 GeV 3.21489 1.43708e-05 Gy 27.7964
+
  layer 13: 17.3333 12.0232 GeV 3.23203 1.44474e-05 Gy 28.0874
-
  layer 14: 18.6667 12.6008 GeV 3.38732 1.51416e-05 Gy 29.2873
+
  layer 14: 18.6667 12.6138 GeV 3.3908 1.51571e-05 Gy 29.4672
-
  layer 15: 20 13.0799 GeV 3.51611 1.57173e-05 Gy 30.4008
+
  layer 15: 20 12.9498 GeV 3.48112 1.55609e-05 Gy 30.2521
-
  layer 16: 21.3333 13.853 GeV 3.72393 1.66463e-05 Gy 32.1977
+
  layer 16: 21.3333 13.6288 GeV 3.66365 1.63768e-05 Gy 31.8383
-
  layer 17: 22.6667 14.4719 GeV 3.89029 1.73899e-05 Gy 33.636
+
  layer 17: 22.6667 14.5047 GeV 3.89912 1.74294e-05 Gy 33.8847
-
  layer 18: 24 15.3391 GeV 4.12341 1.8432e-05 Gy 35.6517
+
  layer 18: 24 15.4695 GeV 4.15847 1.85887e-05 Gy 36.1385
-
  layer 19: 25.3333 16.7448 GeV 4.50128 2.01211e-05 Gy 38.9188
+
  layer 19: 25.3333 16.5838 GeV 4.45801 1.99276e-05 Gy 38.7416
-
  layer 20: 26.6667 18.2636 GeV 4.90957 2.19461e-05 Gy 42.4489
+
  layer 20: 26.6667 18.3271 GeV 4.92663 2.20224e-05 Gy 42.8141
-
  layer 21: 28 20.5916 GeV 5.53539 2.47436e-05 Gy 47.8598
+
  layer 21: 28 20.6828 GeV 5.5599 2.48532e-05 Gy 48.3174
-
  layer 22: 29.3333 24.3455 GeV 6.54449 2.92543e-05 Gy 56.5847
+
  layer 22: 29.3333 24.4737 GeV 6.57896 2.94084e-05 Gy 57.1734
-
  layer 23: 30.6667 32.4192 GeV 8.71484 3.8956e-05 Gy 75.3499
+
  layer 23: 30.6667 32.4582 GeV 8.72532 3.90028e-05 Gy 75.826
-
  layer 24: 32 43.0249 GeV 11.5658 5.17001e-05 Gy 100
+
  layer 24: 32 42.8062 GeV 11.507 5.14373e-05 Gy 100
-
  layer 25: 33.3333 1.97182 GeV 0.530059 2.3694e-06 Gy 4.58297
+
  layer 25: 33.3333 1.99409 GeV 0.536046    2.39617e-06 Gy 4.65842
-
  layer 26: 34.6667 0 eV         0 0 Gy         0
+
  layer 26: 34.6667 1.01076 MeV 0.000271708 1.21456e-09 Gy 0.00236124
-
  layer 27: 36 0 eV         0 0 Gy         0
+
  layer 27: 36 503.954 keV 0.000135472 6.05568e-10 Gy 0.00117729
-
  layer 28: 37.3333 442.484 keV 0.00011 5.31703e-10 Gy 0.00102844
+
  layer 28: 37.3333 33.8047 keV 9.08729e-06 4.06209e-11 Gy 7.89716e-05
-
  layer 29: 38.6667 0 eV         0     0 Gy 0
+
  layer 29: 38.6667 2.03563 MeV 0.000547213  2.44608e-09 Gy 0.00475547
-
  layer 30: 40 0 eV         0     0 Gy 0
+
  layer 30: 40 0 eV         0     0 Gy         0
Line 107: Line 107:
  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.309 GeV
+
  Total energy deposit= 367.244 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 115: Line 115:
<pre style="color: #800000; background-color: #dcdcdc">
<pre style="color: #800000; background-color: #dcdcdc">
-
1.33333 20.404
+
x[mm]  Dose/MaxDose[%]
-
2.66667 20.8104
+
1.33333 20.6613
-
4 21.384
+
2.66667 20.831
-
5.33333 21.6998
+
4 21.4291
-
6.66667 22.2678
+
5.33333 21.9294
-
8 22.7073
+
6.66667 22.4821
-
9.33333 23.277
+
8 22.7908
-
10.6667 23.8187
+
9.33333 23.3252
-
12 24.4422
+
10.6667 24.0918
-
13.3333 25.0271
+
12 24.5704
-
14.6667 26.2448
+
13.3333 25.2942
-
16 26.9606
+
14.6667 26.0803
-
17.3333 27.7964
+
16 27.1837
-
18.6667 29.2873
+
17.3333 28.0874
-
20 30.4008
+
18.6667 29.4672
-
21.3333 32.1977
+
20 30.2521
-
22.6667 33.636
+
21.3333 31.8383
-
24 35.6517
+
22.6667 33.8847
-
25.3333 38.9188
+
24 36.1385
-
26.6667 42.4489
+
25.3333 38.7416
-
28 47.8598
+
26.6667 42.8141
-
29.3333 56.5847
+
28 48.3174
-
30.6667 75.3499
+
29.3333 57.1734
 +
30.6667 75.826
32 100
32 100
-
33.3333 4.58297
+
33.3333 4.65842
-
34.6667 0
+
34.6667 0.00236124
-
36 0
+
36 0.00117729
-
37.3333 0.00102844
+
37.3333 7.89716e-05
-
38.6667 0
+
38.6667 0.00475547
40 0
40 0
</pre>  
</pre>  
Line 192: Line 193:
/protonPB/det/setSliceSizeYZ 1 cm
/protonPB/det/setSliceSizeYZ 1 cm
/protonPB/det/sliceNumber 30  
/protonPB/det/sliceNumber 30  
 +
#
 +
# define longitudinal scoring mesh
 +
# along the beam
 +
/score/create/boxMesh waterMeshlongitudinal
 +
/score/mesh/boxSize 2. 2. 2. cm
 +
/score/mesh/nBin 30 1 1
 +
/score/quantity/energyDeposit energyDeposit
 +
/score/quantity/doseDeposit doseDeposit
 +
/score/close
 +
#
 +
# define lateral scoring mesh
 +
# centered at the Bragg peak
 +
/score/create/boxMesh waterMeshlateral
 +
/score/mesh/boxSize 1. 2. 2. cm
 +
/score/mesh/nBin 1 30 1
 +
/score/mesh/translate/xyz 1.2 0. 0. cm
 +
/score/quantity/energyDeposit energyDeposit
 +
/score/quantity/doseDeposit doseDeposit
 +
/score/close
#
#
# set physics process
# set physics process
Line 213: Line 233:
#
#
# beam size
# beam size
-
#/protonPB/gun/rndm 3 mm
+
#/photonPB/gun/rndm 3 mm
#
#
# step limit (recommended not bigger than 5% of  
# step limit (recommended not bigger than 5% of  
Line 224: Line 244:
/analysis/setFileName Proton
/analysis/setFileName Proton
#
#
-
#  
+
# histogram
/analysis/h1/set 2 50 25 35 mm
/analysis/h1/set 2 50 25 35 mm
# number of events
# number of events
/run/beamOn 6000
/run/beamOn 6000
-
#
+
# dump scores to a file
 +
/score/dumpQuantityToFile waterMeshlongitudinal doseDeposit DoseLongitudinalMesh.txt
 +
/score/dumpQuantityToFile waterMeshlongitudinal energyDeposit EnergyLongitudinalMesh.txt
 +
/score/dumpQuantityToFile waterMeshlateral doseDeposit DoseLateralMesh.txt
 +
/score/dumpQuantityToFile waterMeshlateral energyDeposit EnergyLateralMesh.txt
</pre>
</pre>

Revision as of 11:56, 24 July 2014

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