Warning, /geant4/examples/advanced/hadrontherapy/macro/defaultMacro.mac is written in an unsupported language. File is not indexed.
0001 # Macro file for the initialization of example B1
0002 # in interactive session
0003 #
0004 # Set some default verbose
0005 /control/verbose 2
0006 /control/saveHistory
0007 /run/verbose 2
0008 #
0009 # Change the default number of threads (in multi-threaded mode)
0010 /run/numberOfThreads 4
0011
0012 ##########################
0013 # Set of the physic models
0014 #
0015 #/Physics/addPhysics HADRONTHERAPY_1
0016 /Physics/addPhysics QGSP_BIC_HP
0017
0018 #/run/numberOfThreads 1
0019
0020 # Initialize kernel
0021 /run/initialize
0022
0023 # Set a very high time threshold to allow all decays to happen
0024 #/process/had/rdm/thresholdForVeryLongDecayTime 1.0e+60 year
0025
0026 ####################################################
0027 # Set here the CUT and the STEP MAX for the tracking.
0028 # Suggested values of cut and step
0029 #
0030 # Set cuts OUTSIDE the detector
0031 #
0032 /run/setCut 1 m
0033
0034 # Set cuts ONLY inside the detector
0035 #
0036 /run/setCutForRegion DetectorLog 0.1 mm
0037
0038
0039 /run/setCutForAGivenParticle e- 0.1 mm
0040 /run/setCutForAGivenParticle e+ 0.1 mm
0041 /run/setCutForAGivenParticle proton 0.1 mm
0042 #/Step/waterPhantomStepMax 0.1 mm
0043
0044 #########################
0045 # Set the primary particle type,
0046 # energy and position along the X direction
0047 #
0048
0049 #---------------------------gps-----------------
0050 /gps/pos/shape Circle
0051 /gps/pos/centre -310. 0. 0. cm
0052 /gps/pos/radius 0. mm
0053 /gps/pos/sigma_r 2. mm
0054 /gps/particle proton
0055 #/gps/particle ion
0056 #/gps/ion 6 12 6
0057 /gps/pos/type Beam
0058
0059 # the incident surface is in the y-z plane
0060 #
0061 /gps/pos/rot1 0 1 0
0062 /gps/pos/rot2 0 0 1
0063
0064 # the beam is travelling along the x-axis without any angular dispersion (angular despersion set to 0.0)
0065 #
0066 /gps/ang/rot1 0 0 1
0067 /gps/ang/rot2 0 1 0
0068 /gps/ang/type beam1d
0069 /gps/ang/sigma_r 0.028 deg
0070
0071 # the beam energy is in gaussian profile
0072 #
0073 /gps/ene/type Gauss
0074 /gps/ene/mono 62.3 MeV
0075 /gps/ene/sigma 0.25 MeV
0076 #/gps/ene/mono 744 MeV
0077 #/gps/ene/sigma 0.740 MeV
0078
0079 ###################################################
0080 # DETECTOR DESCRIPTION
0081 ########################################################
0082 # USE THE FOLLOWING FOUR LINES FOR DEFINITION OF PHANTOM
0083 ########################################################
0084 #
0085 #/changePhantom/position 20. 0. 0. cm
0086 #/changePhantom/update
0087
0088 /changePhantom/size 40 40 40 cm
0089 /changePhantom/position 20 0 0 cm
0090
0091 ###################################################
0092 # USE THE FOLLOWING FOR A VOXELLIZED DETECTOR
0093 ###################################################
0094 #
0095 /changeDetector/size 4 4 4 cm
0096 /changeDetector/voxelSize .1 40 40 mm
0097
0098 # Put the detector in the lower left corner of the phantom
0099 #
0100 /changeDetector/displacement 0 18 18 cm
0101 /changePhantom/update
0102
0103 #########################
0104 # Start of the run
0105 #
0106 # If secondary particles dose & fluence are needed
0107 /analysis/secondary true
0108 /analysis/computeLet
0109
0110
0111 ### RBE
0112 /rbe/verbose 2
0113 /rbe/loadLemTable data/rbe/lem1.csv
0114 /rbe/calculation 1
0115 /rbe/accumulate 1
0116
0117 ## U87
0118 #/rbe/dcut 8.0 gray
0119 #/rbe/alphaX 0.11 gray(-1)
0120 #/rbe/betaX 0.06 gray(-2)
0121 #/rbe/lemTable data/rbe/U87_LEM.dat
0122 /rbe/cellLine U87
0123
0124 ## AG01522
0125 #/rbe/dcut 9.5 gray
0126 #/rbe/alphaX 0.54 gray(-1)
0127 #/rbe/betaX 0.062 gray(-2)
0128 #/rbe/lemTable data/rbe/AG01522_LEM.dat
0129
0130 /rbe/doseScale 7777770
0131 # Dose scale should be: 2e7 / nEvents * voxelSizeInCm2
0132
0133
0134 # Default material is water liquid
0135 /changePhantom/material G4_WATER
0136
0137 /run/printProgress 1000
0138 /run/beamOn 100000
0139
0140 # Visualization setting (in case you wish it)
0141 #
0142 # /control/execute vis.mac
