Back to home page

EIC code displayed by LXR

 
 

    


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