[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 |
---|