Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/radiobiology/macro/defaultMacro.mac

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ] 

  1 ################################################################################
  2 # Default macro for radiobiology example.
  3 # Creates a beam of protons and calculates dose, LET, and RBE.
  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