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 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) |
---|
25 | |
---|
26 | ########################## |
---|
27 | # Initialisation procedure |
---|
28 | # |
---|
29 | /run/initialize |
---|
30 | |
---|
31 | ########################## |
---|
32 | # Visualisation |
---|
33 | # |
---|
34 | /vis/scene/create |
---|
35 | #/vis/open OGLIQt # only if QT library are installed |
---|
36 | /vis/open OGLIX |
---|
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 | /tracking/storeTrajectory 1 |
---|
42 | #/vis/scene/endOfEventAction accumulate |
---|
43 | /vis/scene/endOfEventAction accumulate -1 |
---|
44 | /vis/viewer/update |
---|
45 | |
---|
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 |
---|
52 | |
---|
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 |
---|
61 | |
---|
62 | ######################### |
---|
63 | # Display the event number |
---|
64 | # during the run |
---|
65 | # |
---|
66 | /event/printEventNumber 10 |
---|
67 | |
---|
68 | ######################### |
---|
69 | # Start of the run |
---|
70 | # |
---|
71 | /run/beamOn 10 |
---|