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