| [1313] | 1 | # G.A.P.Cirrone
|
|---|
| [807] | 2 | #
|
|---|
| [1230] | 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!>
|
|---|
| [807] | 6 | #
|
|---|
| [1230] | 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
|
|---|
| [807] | 12 | #
|
|---|
| 13 | /control/verbose 1
|
|---|
| 14 | /tracking/verbose 0
|
|---|
| [1230] | 15 | /run/verbose 1
|
|---|
| [807] | 16 | /event/verbose 0
|
|---|
| 17 |
|
|---|
| [1230] | 18 | ##########################
|
|---|
| 19 | # Set of the physic models
|
|---|
| 20 | #
|
|---|
| [1313] | 21 | /physic/addPhysics standard_opt3 # Electromagnetic model
|
|---|
| [1230] | 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)
|
|---|
| [807] | 25 |
|
|---|
| [1313] | 26 |
|
|---|
| 27 |
|
|---|
| [1230] | 28 | ##########################
|
|---|
| 29 | # Initialisation procedure
|
|---|
| 30 | #
|
|---|
| 31 | /run/initialize
|
|---|
| 32 |
|
|---|
| 33 | ##########################
|
|---|
| 34 | # Visualisation
|
|---|
| 35 | #
|
|---|
| [1313] | 36 | /vis/open OGL
|
|---|
| [807] | 37 | /vis/viewer/flush
|
|---|
| [1230] | 38 | /vis/viewer/set/viewpointThetaPhi 30 140 deg
|
|---|
| 39 | /vis/viewer/zoom 1
|
|---|
| 40 | /vis/viewer/pan -10 0 cm
|
|---|
| [1313] | 41 | /vis/scene/add/trajectories
|
|---|
| 42 | /vis/scene/endOfEventAction accumulate
|
|---|
| [807] | 43 |
|
|---|
| [1230] | 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
|
|---|
| [807] | 50 |
|
|---|
| [1230] | 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
|
|---|
| [807] | 59 |
|
|---|
| [1230] | 60 | #########################
|
|---|
| 61 | # Display the event number
|
|---|
| 62 | # during the run
|
|---|
| 63 | #
|
|---|
| 64 | /event/printEventNumber 10
|
|---|
| [807] | 65 |
|
|---|
| [1230] | 66 | #########################
|
|---|
| 67 | # Start of the run
|
|---|
| 68 | #
|
|---|
| [1313] | 69 | # If secondaries dose/fluence are needed
|
|---|
| 70 | #/analysis/secondaries true
|
|---|
| 71 |
|
|---|
| 72 | /run/beamOn 100
|
|---|
| 73 |
|
|---|
| 74 |
|
|---|
