| Geant4 Cross Reference |
1 ############################################## 1 ################################################################################
2 # Default macro for radiobiology example. 2 # Default macro for radiobiology example.
3 # Creates a beam of protons and calculates dos 3 # Creates a beam of protons and calculates dose, LET, and RBE.
4 # May 2023 4 # May 2023
5 5
6 6
7 ############################################## 7 ################################################################################
8 # Set of the verboses 8 # Set of the verboses
9 # 9 #
10 /control/verbose 0 10 /control/verbose 0
11 /tracking/verbose 0 11 /tracking/verbose 0
12 /run/verbose 0 12 /run/verbose 0
13 /event/verbose 0 13 /event/verbose 0
14 14
15 ########################## 15 ##########################
16 # Set of the physic models 16 # Set of the physic models
17 # 17 #
18 /Physics/addPhysics HADRONTHERAPY_1 18 /Physics/addPhysics HADRONTHERAPY_1
19 /run/numberOfThreads 4 19 /run/numberOfThreads 4
20 20
21 ############################################## 21 ################################################################################
22 # Geometry 22 # Geometry
23 # 23 #
24 /detectorGeom/setMat G4_WATER 24 /detectorGeom/setMat G4_WATER
25 /detectorGeom/setBoxSizes 4 4 4 cm 25 /detectorGeom/setBoxSizes 4 4 4 cm
26 26
27 ############################################## 27 ################################################################################
28 # Voxelization procedure 28 # Voxelization procedure
29 # 29 #
30 /voxels/setVoxelSizes 1 40 40 mm 30 /voxels/setVoxelSizes 1 40 40 mm
31 31
32 ############################################## 32 ################################################################################
33 # Initialisation procedure 33 # Initialisation procedure
34 # 34 #
35 /run/initialize 35 /run/initialize
36 36
37 /run/setCutForAGivenParticle e- 0.1 mm 37 /run/setCutForAGivenParticle e- 0.1 mm
38 /run/setCutForAGivenParticle e+ 0.1 mm 38 /run/setCutForAGivenParticle e+ 0.1 mm
39 /run/setCutForAGivenParticle proton 0.1 mm 39 /run/setCutForAGivenParticle proton 0.1 mm
40 40
41 ############################################## 41 ################################################################################
42 # Primary definitions 42 # Primary definitions
43 # 43 #
44 /gps/particle proton 44 /gps/particle proton
45 45
46 /gps/pos/centre -2. 0. 0. cm 46 /gps/pos/centre -2. 0. 0. cm
47 /gps/pos/type Point 47 /gps/pos/type Point
48 /gps/pos/shape Circle 48 /gps/pos/shape Circle
49 /gps/pos/radius 0.2 cm 49 /gps/pos/radius 0.2 cm
50 50
51 /gps/direction 1 0 0 51 /gps/direction 1 0 0
52 #/gps/ang/type planar 52 #/gps/ang/type planar
53 #/gps/ang/sigma_r 0.005 rad 53 #/gps/ang/sigma_r 0.005 rad
54 54
55 /gps/ene/type Gauss 55 /gps/ene/type Gauss
56 /gps/ene/mono 58. MeV 56 /gps/ene/mono 58. MeV
57 /gps/ene/sigma 0.65 MeV 57 /gps/ene/sigma 0.65 MeV
58 58
59 ############################################## 59 ################################################################################
60 # Dose 60 # Dose
61 # 61 #
62 /dose/reset 62 /dose/reset
63 /dose/verbose 1 63 /dose/verbose 1
64 /dose/calculate true 64 /dose/calculate true
65 /dose/fileName dose.out 65 /dose/fileName dose.out
66 /dose/print 66 /dose/print
67 67
68 ############################################## 68 ################################################################################
69 # LET 69 # LET
70 # 70 #
71 /LET/reset 71 /LET/reset
72 /LET/verbose 1 72 /LET/verbose 1
73 /LET/calculate true 73 /LET/calculate true
74 /LET/fileName LET.out 74 /LET/fileName LET.out
75 /LET/print 75 /LET/print
76 76
77 ############################################## 77 ################################################################################
78 # RBE 78 # RBE
79 # 79 #
80 /rbe/reset 80 /rbe/reset
81 /rbe/verbose 1 81 /rbe/verbose 1
82 /rbe/loadLemTable data/rbe/lem1.csv 82 /rbe/loadLemTable data/rbe/lem1.csv
83 /rbe/calculate true 83 /rbe/calculate true
84 84
85 ## U87 85 ## U87
86 /rbe/cellLine U87 86 /rbe/cellLine U87
87 /rbe/print 87 /rbe/print
88 88
89 ############################################## 89 ################################################################################
90 # Start of the run 90 # Start of the run
91 # 91 #
92 /run/printProgress 10000 92 /run/printProgress 10000
93 /run/beamOn 100000 93 /run/beamOn 100000
94 94
95 95