Proton beam with realistic geometry
From UCL HEP PBT Wiki
for
Proton beam with realistic geometry
Jump to:
navigation
,
search
== <span style="color:#000080"> Introduction </span> == This example shows the dose distribution in water along the incident proton beam. The difference with the monoenergetic proton pencil beam example is that the beam is defined with realistic geometry. Instead of particle gun, we use general particle source for the generation of the proton beam. For more details on general particle source look [https://geant4.web.cern.ch/geant4/UserDocumentation/UsersGuides/ForApplicationDeveloper/html/ch02s07.html here]. == <span style="color:#000080"> Setting up the environment </span> == ; Connect to HEP cluster and create folder ProtonGBFolder in your area <pre style="color: #800000; background-color: #dcdcdc"> ssh -X username@plus1.hep.ucl.ac.uk password: type your password here cd /home/username/ mkdir ProtonGBFolder cd ProtonGBFolder </pre> ; Setup your environment <pre style="color: #800000; background-color: #dcdcdc"> source /unix/pbt/software/dev/bin/pbt-dev.sh </pre> == <span style="color:#000080"> How to get the code </span> == ; Copy the code to your working directory and rename it <pre style="color: #800000; background-color: #dcdcdc"> cp -r /unix/pbt/tutorials/basic/ProtonGB . mv ProtonGB ProtonGB_source </pre> == <span style="color:#000080"> How to run the code </span> == ; Inside /home/username/ProtonGBFolder/ create a directory <pre style="color: #800000; background-color: #dcdcdc"> mkdir ProtonGB_build </pre> ; To compile the code enter this directory and run cmake and make <pre style="color: #800000; background-color: #dcdcdc"> cd ProtonGB_build cmake -DGeant4_DIR=/unix/pbt/software/dev /home/username/ProtonGBFolder/ProtonGB_source make </pre> ; Run macro proton.mac. <pre style="color: #800000; background-color: #dcdcdc"> ./protonGB proton.mac </pre> == <span style="color:#000080"> How to analyze data </span> == The macro produces a root file '''Proton.root''' with a 1D histogram showing the energy deposition in water box along the beam line and another histogram showing zoomed Bragg Peak. It also produces a text file '''DoseFile.txt''' with energy and dose deposited at each slice of the water box. This is an example output with the default settings: physics process '''QGSP_BIC_EMY''' and incident proton energy of 62 MeV. <pre style="color: #800000; background-color: #dcdcdc"> Cumulated Doses : X[mm] Edep Edep/Ebeam Dose layer 1: 3.33333 22.8684 GeV 6.14768 % 2.74795e-06 Gy layer 2: 6.66667 24.4078 GeV 6.56151 % 2.93293e-06 Gy layer 3: 10 25.8277 GeV 6.94321 % 3.10354e-06 Gy layer 4: 13.3333 27.87 GeV 7.49224 % 3.34895e-06 Gy layer 5: 16.6667 30.3243 GeV 8.15203 % 3.64387e-06 Gy layer 6: 20 34.0826 GeV 9.16235 % 4.09547e-06 Gy layer 7: 23.3333 39.1645 GeV 10.5285 % 4.70613e-06 Gy layer 8: 26.6667 50.4216 GeV 13.5547 % 6.05883e-06 Gy layer 9: 30 85.9984 GeV 23.1188 % 1.03338e-05 Gy layer 10: 33.3333 3.84312 GeV 1.03314 % 4.61801e-07 Gy layer 11: 36.6667 1.15485 MeV 0.000310456 % 1.38771e-10 Gy layer 12: 40 0 eV 0 % 0 Gy The run consists of 6000 protons of 62 MeV through 4 cm of Water (density: 1 g/cm3 ) Total incident energy(Ebeam)= 371.985 GeV Total energy deposit= 366.469 GeV </pre> You can open Proton.root file in the following way: <pre style="color: #800000; background-color: #dcdcdc"> root -l Proton.root new TBrowser Select ROOT Files and Proton.root </pre> You can change the physics process, incident proton energy and number of slices etc. by modifying the macro proton.mac. In addition, you can configure the particle source by using '''/gps/''' commands. In this macro the proton beam has Guassian energy distribution. The beam characteristics are the same as the ones used at the Laboratori Nazionali del Sud (INFN) in Catania, Italy. Have a look [https://geant4.web.cern.ch/geant4/UserDocumentation/UsersGuides/ForApplicationDeveloper/html/ch02s07.html here] to find out the meaning of all /gps/ commands. Now, use your favorite editor '''pico''', '''vi''', '''emacs''' etc. to open the macro proton.mac: <pre style="color: #800000; background-color: #dcdcdc"> pico proton.mac </pre> You will see: <pre style="color: #800000; background-color: #dcdcdc"> # proton.mac # /control/verbose 2 /run/verbose 2 /tracking/verbose 0 /run/particle/verbose 1 /run/particle/dumpList # # set geometry /protonGB/det/setSizeX 4 cm /protonGB/det/setSizeYZ 4 cm /protonGB/det/setSliceSizeYZ 2 cm /protonGB/det/sliceNumber 12 # # set physics process /protonGB/phys/addPhysics QGSP_BIC_EMY #/protonGB/phys/addPhysics emlivermore #/protonGB/phys/addPhysics empenelope # # production tresholds (recommended range #cut off not bigger than 10% of slice thickness) /protonGB/phys/setCuts 0.5 mm #/protonGB/phys/setGCut 1 um #/protonGB/phys/setECut 1 um #/protonGB/phys/setPCut 1 um # # initialize /run/initialize # # General particle source # proton circle source /gps/pos/shape Circle /gps/pos/centre -4. 0. 0. cm /gps/pos/radius 0. mm /gps/pos/sigma_r 2. mm /gps/particle proton /gps/pos/type Beam # # the incident surface is in the y-z plane /gps/pos/rot1 0 1 0 /gps/pos/rot2 0 0 1 # # the beam is travelling along the x-axis without any angular #dispersion (angular despersion set to 0.0) /gps/ang/rot1 0 0 1 /gps/ang/rot2 0 1 0 /gps/ang/type beam1d /gps/ang/sigma_r 0. deg # # the beam energy is in gaussian profile /gps/ene/type Gauss /gps/ene/mono 62 MeV /gps/ene/sigma 0.3 MeV # # step limit (recommended not bigger than 5% of # slice thickness) /protonGB/stepMax 0.2 mm # /protonGB/event/printModulo 50 # # output file /analysis/setFileName Proton # /analysis/h1/set 2 50 25 35 mm # number of events /run/beamOn 6000 # </pre> == <span style="color:#000080"> Files </span> == [[List of files with brief description]]
Return to
Proton beam with realistic geometry
.
Views
Page
Discussion
View source
History
Personal tools
Log in
Navigation
Main page
Community portal
Current events
Recent changes
Random page
Help
Search
Toolbox
What links here
Related changes
Special pages