| [1230] | 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 | #
|
|---|
| 21 | /physic/addPhysics emstandard_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)
|
|---|
| [807] | 25 |
|
|---|
| [1230] | 26 | ##########################
|
|---|
| 27 | # Initialisation procedure
|
|---|
| 28 | #
|
|---|
| 29 | /run/initialize
|
|---|
| 30 |
|
|---|
| 31 | ##########################
|
|---|
| 32 | # Visualisation
|
|---|
| 33 | #
|
|---|
| [807] | 34 | /vis/scene/create
|
|---|
| [1230] | 35 | #/vis/open OGLIQt # only if QT library are installed
|
|---|
| [807] | 36 | /vis/open OGLIX
|
|---|
| 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
|
|---|
| [807] | 41 | /tracking/storeTrajectory 1
|
|---|
| [1230] | 42 | #/vis/scene/endOfEventAction accumulate
|
|---|
| 43 | /vis/scene/endOfEventAction accumulate -1
|
|---|
| [807] | 44 | /vis/viewer/update
|
|---|
| 45 |
|
|---|
| [1230] | 46 | ##########################
|
|---|
| 47 | # Set here the cut and the step max for the tracking.
|
|---|
| 48 | # Suggested values of cut and step:
|
|---|
| 49 | #
|
|---|
| 50 | /physic/setCuts 0.01 mm
|
|---|
| 51 | /Step/waterPhantomStepMax 0.01 mm
|
|---|
| [807] | 52 |
|
|---|
| [1230] | 53 | #########################
|
|---|
| 54 | # Set the primary particle type,
|
|---|
| 55 | # energy and position along the X direction
|
|---|
| 56 | #
|
|---|
| 57 | /gun/particle proton
|
|---|
| 58 | /beam/energy/meanEnergy 62 MeV
|
|---|
| 59 | /beam/energy/sigmaEnergy 400 keV
|
|---|
| 60 | /beam/position/Xposition -2700 mm
|
|---|
| [807] | 61 |
|
|---|
| [1230] | 62 | #########################
|
|---|
| 63 | # Display the event number
|
|---|
| 64 | # during the run
|
|---|
| 65 | #
|
|---|
| 66 | /event/printEventNumber 10
|
|---|
| [807] | 67 |
|
|---|
| [1230] | 68 | #########################
|
|---|
| 69 | # Start of the run
|
|---|
| 70 | #
|
|---|
| 71 | /run/beamOn 10
|
|---|
