Monoenergetic proton pencil beam
From UCL HEP PBT Wiki
m |
m |
||
Line 64: | Line 64: | ||
== <span style="color:#000080"> How to analyze data </span> == | == <span style="color:#000080"> How to analyze data </span> == | ||
- | The macro produces a root file '''Proton.root''' with | + | The macro produces a root file '''Proton.root''' with two histograms. The first histogram shows the energy deposition in water box along the beam line, the second histogram shows zoomed energy deposition around the peak. The macro also produces two text files '''DoseFile.txt''' and '''PlotDose.txt'''. The file '''DoseFile.txt''' contains energy and dose deposition for every layer. The file '''PlotDose.txt''' contains only depth vs dose values for each layer. These values can be imported in matlab. |
- | water box along the beam line | + | |
- | with | + | === <span style="color:#000080"> Text files </span> === |
+ | |||
+ | This is an example output for '''DoseFile.txt''' with physics process '''QGSP_BIC_EMY''' and incident proton energy of '''62 MeV'''. | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
- | + | Layers : x[mm] Edep Edep/Ebeam[%] Dose Dose/MaxDose[%] | |
- | layer 1: | + | layer 1: 1.33333 8.7788 GeV 2.35989 1.05489e-05 Gy 20.404 |
- | layer 2: | + | layer 2: 2.66667 8.95366 GeV 2.4069 1.0759e-05 Gy 20.8104 |
- | layer | + | layer 3: 4 9.20046 GeV 2.47324 1.10556e-05 Gy 21.384 |
- | layer | + | layer 4: 5.33333 9.33632 GeV 2.50976 1.12188e-05 Gy 21.6998 |
- | layer | + | layer 5: 6.66667 9.58072 GeV 2.57546 1.15125e-05 Gy 22.2678 |
- | layer | + | layer 6: 8 9.76978 GeV 2.62628 1.17397e-05 Gy 22.7073 |
- | layer | + | layer 7: 9.33333 10.0149 GeV 2.69217 1.20342e-05 Gy 23.277 |
- | layer | + | layer 8: 10.6667 10.248 GeV 2.75484 1.23143e-05 Gy 23.8187 |
- | layer | + | layer 9: 12 10.5163 GeV 2.82695 1.26367e-05 Gy 24.4422 |
- | layer | + | layer 10: 13.3333 10.7679 GeV 2.89459 1.2939e-05 Gy 25.0271 |
- | layer | + | layer 11: 14.6667 11.2918 GeV 3.03543 1.35686e-05 Gy 26.2448 |
- | layer | + | layer 12: 16 11.5998 GeV 3.11822 1.39387e-05 Gy 26.9606 |
+ | layer 13: 17.3333 11.9594 GeV 3.21489 1.43708e-05 Gy 27.7964 | ||
+ | layer 14: 18.6667 12.6008 GeV 3.38732 1.51416e-05 Gy 29.2873 | ||
+ | layer 15: 20 13.0799 GeV 3.51611 1.57173e-05 Gy 30.4008 | ||
+ | layer 16: 21.3333 13.853 GeV 3.72393 1.66463e-05 Gy 32.1977 | ||
+ | layer 17: 22.6667 14.4719 GeV 3.89029 1.73899e-05 Gy 33.636 | ||
+ | layer 18: 24 15.3391 GeV 4.12341 1.8432e-05 Gy 35.6517 | ||
+ | layer 19: 25.3333 16.7448 GeV 4.50128 2.01211e-05 Gy 38.9188 | ||
+ | layer 20: 26.6667 18.2636 GeV 4.90957 2.19461e-05 Gy 42.4489 | ||
+ | layer 21: 28 20.5916 GeV 5.53539 2.47436e-05 Gy 47.8598 | ||
+ | layer 22: 29.3333 24.3455 GeV 6.54449 2.92543e-05 Gy 56.5847 | ||
+ | layer 23: 30.6667 32.4192 GeV 8.71484 3.8956e-05 Gy 75.3499 | ||
+ | layer 24: 32 43.0249 GeV 11.5658 5.17001e-05 Gy 100 | ||
+ | layer 25: 33.3333 1.97182 GeV 0.530059 2.3694e-06 Gy 4.58297 | ||
+ | layer 26: 34.6667 0 eV 0 0 Gy 0 | ||
+ | layer 27: 36 0 eV 0 0 Gy 0 | ||
+ | layer 28: 37.3333 442.484 keV 0.000118947 5.31703e-10 Gy 0.00102844 | ||
+ | layer 29: 38.6667 0 eV 0 0 Gy 0 | ||
+ | layer 30: 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. | ||
+ | |||
+ | Edep is the deposited energy in every slice. | ||
+ | Total incident energy(Ebeam)= 372 GeV | ||
+ | Total energy deposit= 367.309 GeV | ||
+ | Dose is the deposited dose in every slice. | ||
+ | MaxDose is the highest dose value from all slices. | ||
</pre> | </pre> | ||
- | + | ||
+ | Open Proton.root file in the following way: | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 95: | Line 125: | ||
Inside is 1D histogram showing the energy deposition in water box: | Inside is 1D histogram showing the energy deposition in water box: | ||
- | + | ||
- | + | ||
You can change the physics process, incident proton energy and number of slices etc. by | You can change the physics process, incident proton energy and number of slices etc. by |