Proton beam with realistic geometry
From UCL HEP PBT Wiki
m |
|||
Line 62: | Line 62: | ||
== <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 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> == | == <span style="color:#000080"> Files </span> == | ||
[[List of files with brief description]] | [[List of files with brief description]] |