| 1 | # G.A.P.Cirrone
|
|---|
| 2 | #
|
|---|
| 3 | # Default macro file. It is called if no argument is provided at run
|
|---|
| 4 | #
|
|---|
| 5 | # i.e. simply typing $G4WORKDIR/bin/Linux-++/Hadrontherapy <no argument here!>
|
|---|
| 6 | #
|
|---|
| 7 | # This macro can be used for a proton beam in water. Both electrmagnetic and
|
|---|
| 8 | # hadronic models are swiched on
|
|---|
| 9 |
|
|---|
| 10 | #########################
|
|---|
| 11 | # Set of the verboses
|
|---|
| 12 | #
|
|---|
| 13 | /control/verbose 1
|
|---|
| 14 | /tracking/verbose 0
|
|---|
| 15 | /run/verbose 1
|
|---|
| 16 | /event/verbose 0
|
|---|
| 17 |
|
|---|
| 18 | ##########################
|
|---|
| 19 | # Set of the physic models
|
|---|
| 20 | #
|
|---|
| 21 | /physic/addPhysics standard_opt3 # Electromagnetic model
|
|---|
| 22 | /physic/addPhysics elastic # Hadronic elastic model
|
|---|
| 23 | /physic/addPhysics binary # Hadronic inelastic model
|
|---|
| 24 | /physic/addPhysics local_ion_ion_inelastic # Hadronic inelastic model for ions (local physic list)
|
|---|
| 25 |
|
|---|
| 26 |
|
|---|
| 27 |
|
|---|
| 28 | ##########################
|
|---|
| 29 | # Initialisation procedure
|
|---|
| 30 | #
|
|---|
| 31 | /run/initialize
|
|---|
| 32 |
|
|---|
| 33 | ##########################
|
|---|
| 34 | # Visualisation
|
|---|
| 35 | #
|
|---|
| 36 | /vis/open OGL
|
|---|
| 37 | /vis/viewer/flush
|
|---|
| 38 | /vis/viewer/set/viewpointThetaPhi 30 140 deg
|
|---|
| 39 | /vis/viewer/zoom 1
|
|---|
| 40 | /vis/viewer/pan -10 0 cm
|
|---|
| 41 | /vis/scene/add/trajectories
|
|---|
| 42 | /vis/scene/endOfEventAction accumulate
|
|---|
| 43 |
|
|---|
| 44 | ##########################
|
|---|
| 45 | # Set here the cut and the step max for the tracking.
|
|---|
| 46 | # Suggested values of cut and step:
|
|---|
| 47 | #
|
|---|
| 48 | /physic/setCuts 0.01 mm
|
|---|
| 49 | /Step/waterPhantomStepMax 0.01 mm
|
|---|
| 50 |
|
|---|
| 51 | #########################
|
|---|
| 52 | # Set the primary particle type,
|
|---|
| 53 | # energy and position along the X direction
|
|---|
| 54 | #
|
|---|
| 55 | /gun/particle proton
|
|---|
| 56 | /beam/energy/meanEnergy 62 MeV
|
|---|
| 57 | /beam/energy/sigmaEnergy 400 keV
|
|---|
| 58 | /beam/position/Xposition -2700 mm
|
|---|
| 59 |
|
|---|
| 60 | #########################
|
|---|
| 61 | # Display the event number
|
|---|
| 62 | # during the run
|
|---|
| 63 | #
|
|---|
| 64 | /event/printEventNumber 10
|
|---|
| 65 |
|
|---|
| 66 | #########################
|
|---|
| 67 | # Start of the run
|
|---|
| 68 | #
|
|---|
| 69 | # If secondaries dose/fluence are needed
|
|---|
| 70 | #/analysis/secondaries true
|
|---|
| 71 |
|
|---|
| 72 | /run/beamOn 100
|
|---|
| 73 |
|
|---|
| 74 |
|
|---|
