G4EMMA/BeamSetup3.mac at master · g4emmaStudent/G4EMMA · GitHub

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
# Macro for EMMA GEANT4 simulation

# Prepare beam
#---------------------------------------------------------------
/mydet/beamZ 36
/mydet/beamA 78
/mydet/beamCharge 24
/mydet/energy 153 MeV
/mydet/sigmaEnergy 0.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.5 mm
#/mydet/transEmittance 2.0 mm
#---------------------------------------------------------------


# Add nuclear-reaction process:
# (Z1,A1)+(Z2,A2)->(Z3,A3)+(Z4,A4)
#-----------------------------------------
/twoBodyReaction/Z1 36
/twoBodyReaction/A1 78
/twoBodyReaction/Z2 1
/twoBodyReaction/A2 1
/twoBodyReaction/Z3 37
/twoBodyReaction/A3 79
/twoBodyReaction/Z4 0
/twoBodyReaction/A4 0
/twoBodyReaction/qmin 0.0   # minimum scattering angle in C.M. (deg)
/twoBodyReaction/qmax 1.0   # maximum ...
/twoBodyReaction/Charge3 25 # charge state of fragment 3
#-----------------------------------------

# Specify central trajectory
#----------------------------
/mydet/centralZ 36
/mydet/centralA 79
/mydet/centralQ 25
/mydet/centralE 143
/mydet/updategeo           # necessary for changes to take effect
#----------------------------

# Amount of output info
#-----------------------
/run/verbose      0
/event/verbose    0
/tracking/verbose 0
#-----------------------

# Run recoil simulation
#-----------------------
/mydet/nEvents 100
/mydet/doPrepare
/mydet/doBeam
/mydet/doReaction
#-----------------------