Warning, /geant4/examples/advanced/hadrontherapy/macro/laserDrivenBeamline.mac is written in an unsupported language. File is not indexed.
0001 # A. Tramontana
0002 # tramontana@lns.infn.it
0003 #
0004 # Macro file recommended for the use with laser-driven beam line and built-in physic list.
0005 # This macro uses the physics list already contained in the Geant4 distribution
0006
0007 #######################
0008 # Set of the verboses #
0009 #######################
0010 #
0011 /control/verbose 1
0012 /tracking/verbose 0
0013 /run/verbose 1
0014 /event/verbose 0
0015
0016 ###########################################################################################################
0017 # ACTIVATED PHYSICS MODELS
0018
0019
0020 /Physics/addPhysics standard_opt4
0021
0022 /geometrySetup/selectGeometry LaserDriven
0023
0024 # Initialize geometry and physic
0025 /run/initialize
0026
0027 # Set a very high time threshold to allow all decays to happen
0028 /process/had/rdm/thresholdForVeryLongDecayTime 1.0e+60 year
0029
0030 ##########################################################################################################
0031
0032 /gps/particle proton
0033
0034
0035 /gps/pos/type Beam
0036 /gps/pos/shape Circle
0037
0038 /gps/pos/centre -171.3 0.0 0.0 cm
0039 /gps/pos/radius 0.0 mm
0040 /gps/pos/sigma_r 0.1 mm
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 /gps/ang/rot1 0 0 1
0050 /gps/ang/rot2 0 1 0
0051 /gps/ang/type beam1d
0052 /gps/ang/sigma_r 0.028 deg
0053
0054 # ###########the beam energy is in gaussian profile
0055 #/gps/ene/type Gauss
0056 #/gps/ene/mono 7 MeV
0057 #/gps/ene/sigma 0.0 MeV
0058 #/gps/ang/type iso
0059 #/gps/ang/mintheta 0 deg
0060 #/gps/ang/maxtheta .1 deg
0061
0062 ################### the beam energy is uniform profile from 4.0 MeV to 12 MeV
0063 #/gps/ene/type User
0064 #/gps/ene/type User
0065 #/gps/hist/type energy
0066 #/gps/hist/point 4.000001 0
0067 #/gps/hist/point 12. 1
0068 #/gps/direction 1 0 0
0069
0070 #####################################################################
0071 # TARANIS energy spectrum without angular dipendence: exp(-.44E) for energy between 1. to 12 MeV
0072 /gps/ene/type Exp
0073 /gps/ene/min 1. MeV
0074 /gps/ene/max 12. MeV
0075 /gps/ene/ezero 2.272727273
0076
0077 ######################################################################
0078 /vis/open OGLI 600x600-0+0
0079
0080 # Disable auto refresh and quieten vis messages whilst scene and
0081 # trajectories are established:
0082 /vis/viewer/set/autoRefresh false
0083 /vis/verbose errors
0084
0085 /vis/drawVolume
0086 /vis/viewer/set/background white
0087 /vis/viewer/set/viewpointThetaPhi 60 130
0088 /vis/viewer/zoom 1.0
0089 /vis/viewer/panTo -.27 -.0
0090 /vis/scene/add/axes 0 0 0 1 m
0091 /vis/scene/add/trajectories
0092 /vis/scene/endOfRunAction accumulate
0093 /vis/modeling/trajectories/create/drawByCharge
0094 /vis/scene/endOfEventAction accumulate -1
0095 /vis/viewer/update
0096
0097
0098
0099 ##########################################################################################################
0100 ############### Define the entrance beam: TARANIS angular distribution
0101 ## beam->1-2
0102 #/gps/verbose 2
0103 #/gps/source/intensity 3.6
0104 #/gps/particle proton
0105 #/gps/pos/type Beam
0106 #/gps/pos/shape Circle
0107 #/gps/pos/centre -171.3 0.0 0.0 cm
0108 #/gps/pos/radius 0.007 mm
0109 #/gps/pos/sigma_r 0. mm
0110 #/gps/ang/rot1 0 0 1
0111 #/gps/ang/rot2 0 1 0
0112 #/gps/pos/rot1 0 0 1
0113 #/gps/pos/rot2 0 1 0
0114 #/gps/ang/type iso
0115 #/gps/ang/mintheta 0 deg
0116 #/gps/ang/maxtheta 26.69 deg
0117 #/gps/ene/type Arb
0118 #/gps/hist/type arb
0119 #/gps/hist/point 1.1 3.5
0120 #/gps/hist/point 3.1399999999999999 1.2
0121 #/gps/hist/inter Exp
0122
0123
0124 ################# change first collimator #################
0125 #/LaserDriven/EnergySelector/FirstCollimator/Radius 10.0 mm
0126 ##/LaserDriven/EnergySelector/FirstCollimator/Thickness 20.0 mm
0127 ##/LaserDriven/EnergySelector/FirstCollimator/zPosizion 42. mm
0128 ################# change second collimator #################
0129 #/LaserDriven/EnergySelector/SecondCollimator/Radius 20.0 mm
0130 #/LaserDriven/EnergySelector/SecondCollimator/Thickness 20.0 mm
0131 #/LaserDriven/EnergySelector/SecondCollimator/zPosizion 42. mm
0132 ######################## change slit ########################
0133 ##/LaserDriven/EnergySelector/Slit/thickness 0.1 mm
0134 /LaserDriven/EnergySelector/Slit/HoleDimensionY 8. mm
0135 /LaserDriven/EnergySelector/Slit/HoleDimensionZ 1. mm
0136 /LaserDriven/EnergySelector/Slit/HolePositionZ 35. mm
0137 #
0138 #/LaserDriven/EnergySelector/Disable
0139 #/LaserDriven/Quadrupoles/DisableQuads
0140 #############################################################################################
0141 # CUT AND STEP MAX
0142 # Suggested values of cut and step:
0143 #
0144 # Set cuts OUTSIDE the detector
0145 #/run/setCut 1 mm
0146 # Set cuts ONLY inside the detector
0147 #
0148 #/run/setCutForRegion DetectorLog 0.5 mm
0149
0150 /Step/waterPhantomStepMax 1 mm
0151
0152
0153 ##################################################
0154 # DETECTOR DESCRIPTION
0155 ###################################################
0156 # USE THE FOLLOWING FOR A VOXELLIZED DETECTOR
0157 ###################################################
0158 #
0159 /changePhantom/position 100 100 0 cm
0160 /changePhantom/size 40 40 40 cm
0161 /changeDetector/size 4 4 4 cm
0162 /changeDetector/voxelSize 1 40 40 mm
0163
0164 # Put the detector in the lower left corner of the phantom
0165 /changeDetector/displacement 0 18 18 cm
0166 /changePhantom/update
0167 /event/printEventNumber 100
0168
0169 /run/beamOn 500
