forked from EMMA-Collaboration/G4EMMA
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathBeamSetup.mac
More file actions
executable file
·120 lines (102 loc) · 3.53 KB
/
Copy pathBeamSetup.mac
File metadata and controls
executable file
·120 lines (102 loc) · 3.53 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
# Macro for EMMA GEANT4 simulation
# Prepare beam
#---------------------------------------------------------------
/mydet/beamZ 8
/mydet/beamA 16
/mydet/beamCharge 8
/mydet/energy 90.1 MeV
/mydet/sigmaEnergy 1 # beam-energy resolution (sigma) in percent of beam energy
#/mydet/Angle 0.036 deg # beam angle sampled uniformly between 0 and theta=sigmaAngle
/mydet/beamSpotDiameter 1.0 mm
#/mydet/transEmittance 2.0 mm
#---------------------------------------------------------------
# Add nuclear-reaction process:
# (Z1,A1)+(Z2,A2)->(Z3,A3)+(Z4,A4)
#-----------------------------------------
/twoBodyReaction/Z1 8
/twoBodyReaction/A1 17
/twoBodyReaction/Z2 3
/twoBodyReaction/A2 6
/twoBodyReaction/Z3 10
/twoBodyReaction/A3 21
/twoBodyReaction/Z4 1
/twoBodyReaction/A4 2
/twoBodyReaction/qmin 0.0 # minimum scattering angle in C.M. (deg)
/twoBodyReaction/qmax 1.0 # maximum ...
/twoBodyReaction/Charge3 10 # charge state of fragment 3
#-----------------------------------------
# Specify central trajectory
#----------------------------
/mydet/centralZ 10
/mydet/centralA 21
/mydet/centralQ 10
/mydet/centralE 37.5
/mydet/updategeo # necessary for changes to take effect
#----------------------------
# Amount of output info
#-----------------------
/run/verbose 0
/event/verbose 0
/tracking/verbose 0
#-----------------------
# Run simulation
#-----------------------
/mydet/nEvents 10
/mydet/doPrepare
/mydet/doBeam
/mydet/doReaction
#-----------------------
#-----------------------------------------
/twoBodyReaction/Z1 8
/twoBodyReaction/A1 17
/twoBodyReaction/Z2 3
/twoBodyReaction/A2 6
/twoBodyReaction/Z3 10
/twoBodyReaction/A3 20
/twoBodyReaction/Z4 1
/twoBodyReaction/A4 3
/twoBodyReaction/qmin 0.0 # minimum scattering angle in C.M. (deg)
/twoBodyReaction/qmax 1.0 # maximum ...
/twoBodyReaction/Charge3 10 # charge state of fragment 3
#-----------------------------------------
#-----------------------
#/mydet/doPrepare
#/mydet/doBeam
#/mydet/doReaction
#-----------------------
#-----------------------------------------
/twoBodyReaction/Z1 8
/twoBodyReaction/A1 17
/twoBodyReaction/Z2 3
/twoBodyReaction/A2 6
/twoBodyReaction/Z3 9
/twoBodyReaction/A3 20
/twoBodyReaction/Z4 2
/twoBodyReaction/A4 3
/twoBodyReaction/qmin 0.0 # minimum scattering angle in C.M. (deg)
/twoBodyReaction/qmax 1.0 # maximum ...
/twoBodyReaction/Charge3 9 # charge state of fragment 3
#-----------------------------------------
#-----------------------
#/mydet/doPrepare
#/mydet/doBeam
#/mydet/doReaction
#-----------------------
#-----------------------------------------
/twoBodyReaction/Z1 8
/twoBodyReaction/A1 17
/twoBodyReaction/Z2 3
/twoBodyReaction/A2 6
/twoBodyReaction/Z3 10
/twoBodyReaction/A3 20
/twoBodyReaction/Z4 1
/twoBodyReaction/A4 3
/twoBodyReaction/qmin 0.0 # minimum scattering angle in C.M. (deg)
/twoBodyReaction/qmax 1.0 # maximum ...
/twoBodyReaction/Charge3 9 # charge state of fragment 3
#-----------------------------------------
#-----------------------
#/mydet/doPrepare
#/mydet/doBeam
#/mydet/doReaction
#-----------------------