Monoenergetic photon pencil beam
From UCL HEP PBT Wiki
m |
|||
Line 106: | Line 106: | ||
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 | ||
- | modifying the macro gamma.mac. Use your favorite editor | + | modifying the macro gamma.mac. Use your favorite editor '''pico''', '''vi''', '''emacs''' etc. For example open the macro with editor pico: |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
pico gamma.mac | pico gamma.mac | ||
</pre> | </pre> | ||
+ | |||
+ | <pre style="color: #800000; background-color: #dcdcdc"> | ||
+ | This is what you will see | ||
+ | |||
+ | # | ||
+ | /control/verbose 2 | ||
+ | /run/verbose 2 | ||
+ | /tracking/verbose 0 | ||
+ | /run/particle/verbose 1 | ||
+ | /run/particle/dumpList | ||
+ | # | ||
+ | # set geometry and material | ||
+ | /photonPB/det/setMat Water | ||
+ | #/photonPB/det/setMat Lead | ||
+ | /photonPB/det/setSizeX 20 cm | ||
+ | /photonPB/det/setSizeYZ 20 cm | ||
+ | /photonPB/det/setSliceSizeYZ 4 cm | ||
+ | /photonPB/det/sliceNumber 20 | ||
+ | # | ||
+ | # set physics process | ||
+ | /photonPB/phys/addPhysics emstandard_opt0 | ||
+ | #photonPB/phys/addPhysics emlivermore | ||
+ | #/photonPB/phys/addPhysics empenelope | ||
+ | # | ||
+ | # production tresholds (range cut off- | ||
+ | # not bigger than 10% of slice thickness) | ||
+ | /photonPB/phys/setCuts 1 mm | ||
+ | #/photonPB/phys/setGCut 1 um | ||
+ | #/photonPB/phys/setECut 1 um | ||
+ | #/photonPB/phys/setPCut 1 um | ||
+ | # | ||
+ | # initialize | ||
+ | /run/initialize | ||
+ | # | ||
+ | # particle gun properties (type of | ||
+ | #particle and energy:10keV-25MeV) | ||
+ | /gun/particle gamma | ||
+ | /gun/energy 20 MeV | ||
+ | # | ||
+ | # beam size | ||
+ | #/photonPB/gun/rndm 3 mm | ||
+ | # | ||
+ | # step limit (not bigger than 5% of | ||
+ | # slice thickness) | ||
+ | /photonPB/stepMax 0.5 mm | ||
+ | # | ||
+ | /photonPB/event/printModulo 50 | ||
+ | # | ||
+ | # output root file | ||
+ | /analysis/setFileName Gamma | ||
+ | # | ||
+ | # number of events | ||
+ | /run/beamOn 10000 | ||
+ | </pre> | ||
=== Change the physics process === | === Change the physics process === | ||
- | The default physics process is '''emstandard_opt0''' used in high energy experiments. | + | The default physics process is '''emstandard_opt0'''. This package is used in high energy experiments. In gamma.mac change |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 126: | Line 180: | ||
</pre> | </pre> | ||
- | The process '''emlivermore''' is used | + | The process '''emlivermore''' is used in low energy physics experiments. Now run the code: |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 134: | Line 188: | ||
=== Change the incident particle energy === | === Change the incident particle energy === | ||
- | The default energy is 20 MeV. This is the | + | The default energy is 20 MeV. This is one of the energies 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 140: | Line 194: | ||
</pre> | </pre> | ||
- | to for example 10 MeV | + | to, for example, 10 MeV |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 152: | Line 206: | ||
</pre> | </pre> | ||
- | Keep in mind that the primary particle generation is done at /PhotonPB_source/src/PrimaryGeneratorAction.cc | + | Keep in mind that the primary particle generation is done at /PhotonPB_source/src/PrimaryGeneratorAction.cc. |
+ | This is part of the PrimaryGeneratorAction.cc: | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 209: | Line 264: | ||
/PhotonPB_source/include/. | /PhotonPB_source/include/. | ||
- | In DetectorConstruction.hh change the MaxLayer to | + | In DetectorConstruction.hh change the MaxLayer to a value which is bigger or equal to the |
+ | number of your slices. | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 221: | Line 277: | ||
</pre> | </pre> | ||
- | You modified | + | You modified the code therefore you need to comppile it again. In directory PhotonPB_build do |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> |