Warning, /geant4/examples/advanced/hadrontherapy/macro/Trento_parameters.mac is written in an unsupported language. File is not indexed.
0001 #
0002 # Macro file recommended for the use with proton beams and built-in physic list.
0003 # This macro uses the physics list already contained in the Geant4 distribution
0004 #
0005
0006 #########################
0007 # Set of the verboses
0008 #
0009 /control/verbose 1
0010 /tracking/verbose 0
0011 /run/verbose 1
0012 /event/verbose 0
0013
0014 #######
0015 # Select geometry
0016 /geometrySetup/selectGeometry TrentoLine
0017
0018
0019
0020 ##########################
0021 # Set of the physic models
0022 #
0023 /Physics/addPhysics HADRONTHERAPY_1
0024
0025 /random/setSeeds 891 61
0026
0027 ##########################
0028 # Initialisation procedure
0029 #
0030 /run/numberOfThreads 1
0031
0032 /run/initialize
0033
0034 # Set a very high time threshold to allow all decays to happen
0035 /process/had/rdm/thresholdForVeryLongDecayTime 1.0e+60 year
0036
0037 #########################
0038 # Set the primary particle type,
0039 # energy and position along the X direction
0040 #
0041 #---------------------------gps-----------------
0042
0043 /gps/pos/type Beam
0044 /gps/pos/shape Circle
0045 /gps/pos/centre -95. 0. 0. cm
0046 #/gps/pos/centre -82. 0. 0. cm
0047 #/gps/pos/radius 0. mm
0048 #/gps/pos/sigma_r 4.515 mm
0049 /gps/pos/sigma_x 1.805 mm
0050 /gps/pos/sigma_y 2.835 mm
0051 /gps/particle proton
0052 /gps/direction 1 0 0
0053
0054 # the incident surface is in the y-z plane
0055
0056 /gps/pos/rot1 0 1 0
0057 /gps/pos/rot2 0 0 1
0058
0059 # the beam is travelling along the x-axis without any angular dispersion (angular despersion set to 0.0)
0060 #
0061
0062 /gps/ang/rot1 0 0 1
0063 /gps/ang/rot2 0 1 0
0064 /gps/ang/type beam2d
0065 #/gps/ang/focuspoint 0 0 0 cm
0066 /gps/ang/sigma_x 0.00321 rad
0067 /gps/ang/sigma_y 0.00364 rad
0068 #/gps/ang/sigma_r 0 deg
0069
0070
0071 # the beam energy is in gaussian profile
0072 #
0073 /gps/ene/type Gauss
0074 /gps/ene/mono 147.1 MeV
0075 /gps/ene/sigma 1.1 MeV
0076
0077
0078 ##########################
0079 # Set here the cut and the step max for the tracking.
0080 # Suggested values of cut and step:
0081 # Set cuts OUTSIDE the detector
0082 #
0083 #/run/setCut 0.1 mm
0084
0085 # Set cuts ONLY inside the detector
0086 #
0087 /run/setCutForRegion DetectorLog 0.01 mm
0088
0089 /Step/waterPhantomStepMax 0.1 mm
0090
0091 ######################### GEOMETRY OF BEAM LINE DESCRIPTION
0092 ####################################################################################################
0093 #/beamLine/ScatteringFoil/thickness 0.25 cm
0094 #/beamLine/ScatteringFoil/material G4_Pb
0095 #/beamLine/ScatteringFoil/material G4_Ta
0096 #/beamLine/ScatteringFoil/material G4_Cu
0097 #/beamLine/preCollimator/ZSize mm
0098 #(cambia la dimensione lungo x dato che รจ ruotato)
0099 #/beamLine/preCollimator/XPosition 150. cm
0100 #/beamLine/preCollimator/innerRadius 2.5 mm
0101 #/beamLine/FirstCollimator/innerRadius 2.5 mm
0102
0103
0104
0105 ######################### DETECTOR DESCRIPTION
0106 ###################################################
0107 # USE THE FOLLOWING FOUR LINES FOR DEFINITION OF PHANTOM
0108 ###################################################
0109 #
0110 #/changePhantom/position 220. 0. 0. cm
0111 #/changePhantom/update
0112
0113 #/changePhantom/size 20 20 40 cm
0114 #/changePhantom/position 20 0 0 cm
0115
0116 ###################################################
0117 # USE THE FOLLOWING FOR A VOXELLIZED DETECTOR
0118 ###################################################
0119 #
0120
0121
0122 /changeDetector/size 20 10 10 cm #size for Giraffe
0123 #/changeDetector/size 20 20 20 cm #large detector for dose calculation
0124 /changeDetector/voxelSize 0.02 1 1 cm
0125
0126 # Put the detector in the lower left corner of the phantom
0127 #
0128 /changeDetector/displacement 0 15 15 cm #displacement when detector is small
0129 #/changeDetector/displacement 0 10 10 cm #displacement when detector is larger
0130 /changePhantom/update
0131
0132
0133 # Visualisation
0134 #
0135 /vis/scene/create
0136 /vis/open OGL
0137 #/vis/open ## open the default viewer
0138 /vis/viewer/flush
0139 /vis/viewer/set/viewpointThetaPhi 30 140 deg
0140 /vis/viewer/zoomTo 1
0141 /vis/viewer/panTo 100 20 cm
0142 /vis/viewer/panTo -10 0 cm
0143 /vis/scene/add/trajectories
0144 /tracking/storeTrajectory 1
0145 /vis/scene/endOfEventAction accumulate 1
0146 /vis/viewer/update
0147
0148 #########################
0149 # Display the event number
0150 # during the run
0151 #
0152 /event/printEventNumber 500
0153
0154
0155 # Default material is liquid water (G4_WATER)
0156
0157 # For dose average LET computations of primary and secondary particles
0158 #/analysis/computeLet
0159
0160 /run/beamOn 10
0161
0162
