Monoenergetic photon pencil beam
From UCL HEP PBT Wiki
m |
m |
||
Line 67: | Line 67: | ||
water box along the beam line. It also produces text files: '''DoseFile.txt''' with | water box along the beam line. It also produces text files: '''DoseFile.txt''' with | ||
energy and dose deposited in each slice and '''PlotDose.txt''' with dose deposited in each slice. | energy and dose deposited in each slice and '''PlotDose.txt''' with dose deposited in each slice. | ||
- | This is an example output with | + | |
+ | === <span style="color:#000080"> Text files </span> === | ||
+ | |||
+ | This is an example output for '''DoseFile.txt''' with physics process '''emstandard_opt0''' and incident photon energy of '''20 MeV'''. | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 101: | Line 104: | ||
MaxDose is the highest dose value from all slices. | MaxDose is the highest dose value from all slices. | ||
</pre> | </pre> | ||
+ | |||
+ | This is PlotDose.txt. These values can be directly imported in matlab. | ||
+ | |||
+ | <pre style="color: #800000; background-color: #dcdcdc"> | ||
+ | 10 32.1047 | ||
+ | 20 50.0432 | ||
+ | 30 62.3066 | ||
+ | 40 77.3614 | ||
+ | 50 80.9917 | ||
+ | 60 89.7831 | ||
+ | 70 99.7318 | ||
+ | 80 100 | ||
+ | 90 95.7783 | ||
+ | 100 95.3955 | ||
+ | 110 93.4478 | ||
+ | 120 91.2579 | ||
+ | 130 90.4998 | ||
+ | 140 87.7295 | ||
+ | 150 90.0244 | ||
+ | 160 88.9634 | ||
+ | 170 80.4985 | ||
+ | 180 79.051 | ||
+ | 190 77.9945 | ||
+ | 200 0 | ||
+ | </pre> | ||
+ | |||
+ | === <span style="color:#000080"> Root file </span> === | ||
Open Gamma.root file in the following way: | Open Gamma.root file in the following way: | ||
Line 110: | Line 140: | ||
Select ROOT files and Gamma.root | Select ROOT files and Gamma.root | ||
- | |||
</pre> | </pre> | ||
- | + | The histogram inside Gamma.root shows the energy deposition in water box: | |
http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/PhotonPB/Edep_PhotonB.png | http://www.hep.ucl.ac.uk/pbt/RadiotherapyWorkbook/skins/common/images/PhotonPB/Edep_PhotonB.png | ||
+ | |||
+ | === <span style="color:#000080"> Changing macro gamma.mac </span> === | ||
You can change the physics process, incident photon energy, phantom material, number of slices etc. by | You can change the physics process, incident photon energy, phantom material, number of slices etc. by | ||
Line 124: | Line 155: | ||
</pre> | </pre> | ||
- | This is | + | This is the content of the macro: |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 176: | Line 207: | ||
# | # | ||
# number of events | # number of events | ||
- | /run/beamOn | + | /run/beamOn 6000 |
</pre> | </pre> | ||
- | + | '''Change the physics process''' | |
The default physics process is '''emstandard_opt0'''. This package is used in high energy experiments. In gamma.mac change | The default physics process is '''emstandard_opt0'''. This package is used in high energy experiments. In gamma.mac change | ||
Line 199: | Line 230: | ||
</pre> | </pre> | ||
- | + | '''Change the incident particle energy''' | |
- | The default energy is 20 MeV. | + | The default energy is 20 MeV. It is a typical energy used in radiotherapy. In gamma.mac you can change the value of 20 MeV |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 246: | Line 277: | ||
means that only one particle is generated, incident from (x0,y0,z0). The default number of | means that only one particle is generated, incident from (x0,y0,z0). The default number of | ||
- | events is set to | + | events is set to 6000. Therefore, 6000 particles are incident to the water box. |
- | + | '''Change the material''' | |
- | + | In this example we compute the energy deposition of photons in water box. However, | |
there is an option to change the box material from water to lead. | there is an option to change the box material from water to lead. | ||
Line 271: | Line 302: | ||
</pre> | </pre> | ||
- | + | '''Change the number of slices''' | |
- | You can change the number of | + | You can change the number of slices. The default number is 20. Keep in mind that |
- | if you want to have | + | if you want to have bigger number of slices you need to modify the file DetectorConstruction.hh in |
/PhotonPB_source/include/. | /PhotonPB_source/include/. | ||
- | In DetectorConstruction.hh | + | In DetectorConstruction.hh set MaxLayer to a value which is bigger then the number of your slices. The default number is MaxLayer=30. For example, if you want your box to be divided to 35 slices you need to set the value of MaxLayer to ,for example, 40 then in gamma.mac you need to change the number of slices: |
- | number of your slices. | + | |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
- | /photonPB/det/sliceNumber | + | /photonPB/det/sliceNumber 35 |
</pre> | </pre> | ||
- | + | Every time you modify files in PhotonPB_source you need to compile your code again. In directory PhotonPB_build do | |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 302: | Line 326: | ||
</pre> | </pre> | ||
+ | === <span style="color:#000080"> Visualisation </span> === | ||
== <span style="color:#000080"> Files </span> == | == <span style="color:#000080"> Files </span> == | ||
[[List of files with brief description]] | [[List of files with brief description]] |