Warning, /geant4/examples/advanced/hadrontherapy/macro/BESTBeamLine.mac is written in an unsupported language. File is not indexed.
0001 ##############################
0002 ####### BEST BEAM LINE #######
0003 ##############################
0004 # Set of the verboses
0005 #
0006 /control/verbose 1
0007 /tracking/verbose 0
0008 /run/verbose 1
0009 /event/verbose 0
0010
0011 ##########################
0012 # Set of the physic models
0013 #
0014 /Physics/addPhysics HADRONTHERAPY_1
0015 /run/numberOfThreads 16
0016
0017 ######## Geometry selection #######
0018 /geometrySetup/selectGeometry BESTBeamLine
0019
0020 ##########################
0021 # Initialisation procedure
0022 #
0023 /run/initialize
0024
0025 # Set a very high time threshold to allow all decays to happen
0026 /process/had/rdm/thresholdForVeryLongDecayTime 1.0e+60 year
0027
0028 #########################
0029 # Set the primary particle type,
0030 # energy and position along the X direction
0031 #
0032 #---------------------------gps-----------------
0033 ####### BEST beam parameters #######
0034 /gps/pos/shape Circle
0035 /gps/pos/centre -310. 0. 0. cm
0036 /gps/pos/radius 0. mm
0037 /gps/pos/sigma_r 3. mm
0038 /gps/particle proton
0039 /gps/pos/type Beam
0040
0041
0042 # the incident surface is in the y-z plane
0043
0044 /gps/pos/rot1 0 1 0
0045 /gps/pos/rot2 0 0 1
0046
0047 # the beam is travelling along the x-axis without any angular dispersion (angular despersion set to 0.0)
0048 #
0049
0050 /gps/ang/rot1 0 0 1
0051 /gps/ang/rot2 0 1 0
0052 /gps/ang/type beam1d
0053 /gps/ang/sigma_r 0.16 deg
0054
0055
0056 # the beam energy is in gaussian profile
0057 #
0058 /gps/ene/type Gauss
0059 /gps/ene/mono 70 MeV
0060 /gps/ene/sigma 0.3 MeV
0061
0062
0063 ######## Range shifter settings #########
0064 #/beamLine/BESTRangeShifter/RSMat G4_PLEXIGLASS
0065 #/beamLine/BESTRangeShifter/thickness 13 mm
0066 /beamLine/BESTFinalCollimator/halfInnerRad 12.5 mm
0067
0068 ##########################
0069 # Set here the cut and the step max for the tracking.
0070 # Suggested values of cut and step:
0071 # Set cuts OUTSIDE the detector
0072 #
0073 /run/setCut 1 m
0074
0075
0076
0077 # Set cuts ONLY inside the detector
0078 #
0079 /run/setCutForRegion DetectorLog 0.05 mm
0080
0081 /Step/waterPhantomStepMax 1 mm
0082
0083 #########################
0084 ##############################################################################################################################
0085 # DETECTOR DESCRIPTION
0086 ###################################################
0087 # USE THE FOLLOWING FOUR LINES FOR DEFINITION OF PHANTOM
0088 ###################################################
0089 #
0090 #/changePhantom/position 20. 0. 0. cm
0091 #/changePhantom/update
0092
0093 /changePhantom/size 40 40 40 cm
0094 /changePhantom/position 20 0 0 cm
0095
0096 ###################################################
0097 # USE THE FOLLOWING FOR A VOXELLIZED DETECTOR
0098 ###################################################
0099 #
0100
0101
0102 /changeDetector/size 4 1 1 cm
0103 /changeDetector/voxelSize .1 10 10 mm
0104
0105 # Put the detector in the lower left corner of the phantom
0106 #
0107 /changeDetector/displacement 0 19.5 19.5 cm
0108 /changePhantom/update
0109
0110
0111 # Visualisation
0112 #
0113 #/vis/scene/create
0114 #/vis/open OGL
0115 ##/vis/open ## open the default viewer
0116 #/vis/viewer/flush
0117 #/vis/viewer/set/viewpointThetaPhi 30 140 deg
0118 #/vis/viewer/zoomTo 1
0119 #/vis/viewer/panTo 100 20 cm
0120 #/vis/viewer/panTo -10 0 cm
0121 #/vis/scene/add/trajectories
0122 #/tracking/storeTrajectory 1
0123 #/vis/scene/endOfEventAction accumulate -1
0124 #/vis/viewer/update
0125
0126 #########################
0127 # Display the event number
0128 # during the run
0129 #
0130 /event/printEventNumber 1
0131 /run/printProgress 1000
0132
0133 # Default material is liquid water (G4_WATER)
0134 #/changePhantom/material G4_Pb
0135
0136 /run/setCutForRegion DetectorLog 0.1 mm
0137
0138
0139 /run/setCutForAGivenParticle e- 0.1 mm
0140 /run/setCutForAGivenParticle e+ 0.1 mm
0141 /run/setCutForAGivenParticle proton 0.1 mm
0142
0143 # For dose/fleunce computation of primary and secondary particles
0144 #/analysis/secondary true
0145
0146 # For dose average LET computations of primary and secondary particles#
0147
0148 /analysis/computeLet
0149
0150 ####### Dose scorer at the isocentre #######
0151 /score/create/boxMesh detector
0152 /score/mesh/boxSize 1 20 20 mm
0153 /score/mesh/nBin 1 80 80
0154 /score/mesh/translate/xyz 0.5 0 0 mm
0155 /score/quantity/doseDeposit dose_detector
0156 /score/close
0157
0158 /run/beamOn 100
0159 #########################
0160
0161 /score/dumpQuantityToFile detector dose_detector dose_70Mev.txt
