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