Monoenergetic proton pencil beam
From UCL HEP PBT Wiki
m |
m |
||
Line 70: | Line 70: | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
- | Layers : x[mm] | + | Layers : x[mm] Edep Edep/Ebeam[%] Dose Dose/MaxDose[%] |
- | layer 1: 1.33333 8. | + | layer 1: 1.33333 8.84429 GeV 2.3775 1.06276e-05 Gy 20.6613 |
- | layer 2: 2.66667 8. | + | layer 2: 2.66667 8.91694 GeV 2.39703 1.07149e-05 Gy 20.831 |
- | layer 3: 4 9. | + | layer 3: 4 9.17299 GeV 2.46586 1.10226e-05 Gy 21.4291 |
- | layer 4: 5.33333 9. | + | layer 4: 5.33333 9.38711 GeV 2.52342 1.12799e-05 Gy 21.9294 |
- | layer 5: 6.66667 9. | + | layer 5: 6.66667 9.62373 GeV 2.58702 1.15642e-05 Gy 22.4821 |
- | layer 6: 8 9. | + | layer 6: 8 9.75588 GeV 2.62255 1.1723e-05 Gy 22.7908 |
- | layer 7: 9.33333 | + | layer 7: 9.33333 9.98464 GeV 2.68404 1.19979e-05 Gy 23.3252 |
- | layer 8: 10.6667 10. | + | layer 8: 10.6667 10.3128 GeV 2.77225 1.23922e-05 Gy 24.0918 |
- | layer 9: 12 10. | + | layer 9: 12 10.5176 GeV 2.82732 1.26383e-05 Gy 24.5704 |
- | layer 10: 13.3333 10. | + | layer 10: 13.3333 10.8275 GeV 2.91061 1.30107e-05 Gy 25.2942 |
- | layer 11: 14.6667 11. | + | layer 11: 14.6667 11.164 GeV 3.00107 1.3415e-05 Gy 26.0803 |
- | layer 12: 16 11. | + | layer 12: 16 11.6363 GeV 3.12804 1.39825e-05 Gy 27.1837 |
- | layer 13: 17.3333 | + | layer 13: 17.3333 12.0232 GeV 3.23203 1.44474e-05 Gy 28.0874 |
- | layer 14: 18.6667 12. | + | layer 14: 18.6667 12.6138 GeV 3.3908 1.51571e-05 Gy 29.4672 |
- | layer 15: 20 | + | layer 15: 20 12.9498 GeV 3.48112 1.55609e-05 Gy 30.2521 |
- | layer 16: 21.3333 13. | + | layer 16: 21.3333 13.6288 GeV 3.66365 1.63768e-05 Gy 31.8383 |
- | layer 17: 22.6667 14. | + | layer 17: 22.6667 14.5047 GeV 3.89912 1.74294e-05 Gy 33.8847 |
- | layer 18: 24 15. | + | layer 18: 24 15.4695 GeV 4.15847 1.85887e-05 Gy 36.1385 |
- | layer 19: 25.3333 16. | + | layer 19: 25.3333 16.5838 GeV 4.45801 1.99276e-05 Gy 38.7416 |
- | layer 20: 26.6667 18. | + | layer 20: 26.6667 18.3271 GeV 4.92663 2.20224e-05 Gy 42.8141 |
- | layer 21: 28 20. | + | layer 21: 28 20.6828 GeV 5.5599 2.48532e-05 Gy 48.3174 |
- | layer 22: 29.3333 24. | + | layer 22: 29.3333 24.4737 GeV 6.57896 2.94084e-05 Gy 57.1734 |
- | layer 23: 30.6667 32. | + | layer 23: 30.6667 32.4582 GeV 8.72532 3.90028e-05 Gy 75.826 |
- | layer 24: 32 | + | layer 24: 32 42.8062 GeV 11.507 5.14373e-05 Gy 100 |
- | layer 25: 33.3333 1. | + | layer 25: 33.3333 1.99409 GeV 0.536046 2.39617e-06 Gy 4.65842 |
- | layer 26: 34.6667 0 | + | layer 26: 34.6667 1.01076 MeV 0.000271708 1.21456e-09 Gy 0.00236124 |
- | layer 27: 36 0 | + | layer 27: 36 503.954 keV 0.000135472 6.05568e-10 Gy 0.00117729 |
- | layer 28: 37.3333 | + | layer 28: 37.3333 33.8047 keV 9.08729e-06 4.06209e-11 Gy 7.89716e-05 |
- | layer 29: 38.6667 0 | + | layer 29: 38.6667 2.03563 MeV 0.000547213 2.44608e-09 Gy 0.00475547 |
- | layer 30: 40 0 eV 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. | + | 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. | + | x[mm] Dose/MaxDose[%] |
- | 2.66667 20. | + | 1.33333 20.6613 |
- | 4 21. | + | 2.66667 20.831 |
- | 5.33333 21. | + | 4 21.4291 |
- | 6.66667 22. | + | 5.33333 21.9294 |
- | 8 22. | + | 6.66667 22.4821 |
- | 9.33333 23. | + | 8 22.7908 |
- | 10.6667 | + | 9.33333 23.3252 |
- | 12 24. | + | 10.6667 24.0918 |
- | 13.3333 25. | + | 12 24.5704 |
- | 14.6667 26. | + | 13.3333 25.2942 |
- | 16 | + | 14.6667 26.0803 |
- | 17.3333 | + | 16 27.1837 |
- | 18.6667 29. | + | 17.3333 28.0874 |
- | 20 30. | + | 18.6667 29.4672 |
- | 21.3333 | + | 20 30.2521 |
- | 22.6667 33. | + | 21.3333 31.8383 |
- | 24 | + | 22.6667 33.8847 |
- | 25.3333 38. | + | 24 36.1385 |
- | 26.6667 42. | + | 25.3333 38.7416 |
- | 28 | + | 26.6667 42.8141 |
- | 29.3333 | + | 28 48.3174 |
- | 30.6667 75. | + | 29.3333 57.1734 |
+ | 30.6667 75.826 | ||
32 100 | 32 100 | ||
- | 33.3333 4. | + | 33.3333 4.65842 |
- | 34.6667 0 | + | 34.6667 0.00236124 |
- | 36 0 | + | 36 0.00117729 |
- | 37.3333 | + | 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 | ||
- | #/ | + | #/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> | ||