Warning, /geant4/examples/advanced/hadrontherapy/macro/carbon_beamline.mac is written in an unsupported language. File is not indexed.
0001 # G.A.P.Cirrone
0002 #
0003 # Default macro file. It is called if no argument is provided at run
0004 #
0005 # i.e. simply typing $G4WORKDIR/bin/Linux-++/Hadrontherapy <no argument here!>
0006 #
0007 # This macro can be used for a proton beam in water. Both electrmagnetic and
0008 # hadronic models are swiched on
0009
0010 #########################
0011 # Set of the verboses
0012 #
0013 /control/verbose 1
0014 /tracking/verbose 0
0015 /run/verbose 1
0016 /event/verbose 0
0017
0018 ##########################
0019 # Set of the physic models
0020 #
0021 /Physics/addPhysics HADRONTHERAPY_1
0022
0023 /run/numberOfThreads 4
0024
0025 ##########################
0026 # INFN-LNS transport beam line,
0027 # normally used for interdisciplinary
0028 # researches with carbon and
0029 # and other ion beams is activated
0030 # with the following command
0031 #
0032 /geometrySetup/selectGeometry Carbon
0033
0034 ##########################
0035 # Initialisation procedure
0036 #
0037 /run/initialize
0038
0039 # Set a very high time threshold to allow all decays to happen
0040 /process/had/rdm/thresholdForVeryLongDecayTime 1.0e+60 year
0041
0042 ##########################
0043 # Set here the cut and the step max for the tracking.
0044 # Suggested values of cut and step:
0045 # Set cuts OUTSIDE the detector
0046 /run/setCut 10 mm
0047 # Set cuts ONLY inside the detector
0048 #
0049 /run/setCutForRegion DetectorLog 0.1 mm
0050
0051 /Step/waterPhantomStepMax 0.1 mm
0052
0053 ##########################################################################################################
0054 # PRIMARY PARTICLES
0055
0056 # the beam spot is centered at the origin and is
0057 # of 1d gaussian shape with a 3mm central plateau
0058 #
0059
0060 #---------------------------gps-----------------
0061 /gps/pos/shape Circle
0062 /gps/pos/centre -1710. 0. 0. mm
0063 /gps/pos/radius 0. mm
0064 /gps/pos/sigma_r 2. mm
0065 /gps/particle ion
0066 /gps/ion 6 12 6
0067 /gps/pos/type Beam
0068
0069
0070 # the incident surface is in the y-z plane
0071
0072 /gps/pos/rot1 0 1 0
0073 /gps/pos/rot2 0 0 1
0074
0075 # the beam is travelling along the x-axis without any angular dispersion (angular despersion set to 0.0)
0076 #
0077
0078 /gps/ang/rot1 0 0 1
0079 /gps/ang/rot2 0 1 0
0080 /gps/ang/type beam1d
0081 /gps/ang/sigma_r 0. deg
0082
0083
0084 # the beam energy is in gaussian profile
0085 #
0086 /gps/ene/type Gauss
0087 /gps/ene/mono 744 MeV
0088 /gps/ene/sigma 0.740 MeV
0089
0090
0091 # Voxelised detector
0092 /changePhantom/size 40 40 40 cm
0093 /changePhantom/position 20 0 0 cm
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 /changeDetector/displacement 0 18 18 cm
0100 /changePhantom/update
0101
0102
0103 ##########################
0104 # Visualisation
0105 #
0106 #/vis/open OGLI 600x600-0+0
0107 ##/vis/open ## open the default viewer
0108
0109 # Disable auto refresh and quieten vis messages whilst scene and
0110 # trajectories are established:
0111 /vis/viewer/set/autoRefresh false
0112 #/vis/verbose errors
0113
0114 /vis/drawVolume
0115 /vis/viewer/set/viewpointThetaPhi 30 140 deg
0116 #/vis/viewer/zoom 1
0117 #/vis/viewer/pan -10 0 cm
0118 /vis/scene/add/trajectories smooth
0119
0120 /vis/scene/endOfEventAction accumulate
0121
0122 /vis/viewer/set/autoRefresh true
0123 /vis/verbose warnings
0124
0125 /vis/viewer/flush
0126 #########################
0127 # Display the event number
0128 # during the run
0129 #
0130 /run/printProgress 1000
0131
0132 # Compute Dose and Fluence for secondary
0133 #/analysis/secondary true
0134
0135
0136 #########################
0137 # Start of the run
0138 #
0139 #########################
0140
0141 /run/beamOn 10