Monoenergetic photon pencil beam

From UCL HEP PBT Wiki

Jump to: navigation, search
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 (pico, vi, emacs etc.). For example open the macro with editor pico:
+
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. Change
+
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 for low energy physics experiments. Now run the code:
+
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 energy used in radiotherapy. You can change the value of 20 MeV  
+
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 your desired value.   
+
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 your code therefore you need to comppile it. In directory PhotonPB_build do
+
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">

Revision as of 07:46, 17 July 2014

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