forked from JeffersonLab/g4sbs
-
Notifications
You must be signed in to change notification settings - Fork 0
/
run_example.mac
103 lines (82 loc) · 2.4 KB
/
run_example.mac
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
/g4sbs/filename test.root
/g4sbs/beamcur 10.5 muA
# Have LH2, LD2 (gives scattering chamber)
# GMn config
/g4sbs/target LD2
/g4sbs/kine elastic #Quasielastic
/g4sbs/exp gmn
/g4sbs/targlen 10.0 cm
#or
# GEn config
# H2, 3He (gaseous cell)
#/g4sbs/target 3He
#/g4sbs/kine elastic # Quasielastic
#/g4sbs/targpres 10 atmosphere
#/g4sbs/targlen 55.0 cm
#/g4sbs/gen
# GEp config
#/g4sbs/target LH2
#/g4sbs/kine elastic
#/g4sbs/exp gep # <---- This is important for GEp
#/g4sbs/targlen 40.0 cm
# A1n
#/g4sbs/target 3He
#/g4sbs/kine inelastic <-- Christy/Bosted
#/g4sbs/kine dis <-- CTEQ6.6m
#/g4sbs/targlen 40.0 cm
#/g4sbs/exp a1n
# Wiser pions
#/g4sbs/kine wiser
#/g4sbs/hadron pi-
#/g4sbs/kine flat <-- Flat in phase space generation
/g4sbs/gemres 0.070 mm #GEM resolution
# Raster settings
/g4sbs/rasterx 4.0 mm
/g4sbs/rastery 4.0 mm
/g4sbs/cerdist 0.57 m
/g4sbs/cerdepth 0.75 m
/g4sbs/gemsep 1.34 m
/g4sbs/bbcaldist 1.4 m
# gemconfig = 1 for 4 gems, 2 x 5cm apart <- gemsep -> 2 x 5cm apart
# gemconfig = 2 for 5 gems, 1 <- gemsep -> 4 x 10 cm apart
# gemconfig = 3 for 3 gems, 1 <- gemsep -> 1 <- gemsep -> 1
/g4sbs/gemconfig 2
/g4sbs/bbfield 1 // Use standard BB field
/g4sbs/48d48field 1 // Use constant 48D48field
# If you wanted to use a tosca field map
#/g4sbs/tosfield GEN-map1.table
# Shoot geantinos instead of electrons (good for testing no interactions)
#/g4sbs/shootgeantino
# Set magnet materials to total absorbers (otherwise you get tons of showers
# when just looking at things like acceptances) Default is true
#/g4sbs/totalabs false
# Use lead shielding or not 1 = yes
# * This is only optmized for the GEp 12 GeV point *
#Otherwise use 0 for no shielidng
#/g4sbs/uselead 1
# SBS Clamp Option
#
#/g4sbs/sbsclampopt 1
# 1 for Bigbite (neutron experiments)
# 2 for GEp
/g4sbs/runtime 1.0 s
/g4sbs/beamE 4.4 GeV
/g4sbs/bbang 32.5 deg
/g4sbs/bbdist 1.55 m
/g4sbs/hcalang 31.1 deg
/g4sbs/hcaldist 6.2 m
/g4sbs/48D48dist 1.6 m
/g4sbs/thmin 22.0 deg
/g4sbs/thmax 42.0 deg
/g4sbs/phmin 155.0 deg
/g4sbs/phmax 205.0 deg
/g4sbs/run 100000
/g4sbs/filename test2.root
/g4sbs/beamE 4.4 GeV
/g4sbs/bbang 41.9 deg
/g4sbs/hcalang 24.7 deg
/g4sbs/hcaldist 6.2 m
/g4sbs/48D48dist 1.6 m
/g4sbs/thmin 32.0 deg
/g4sbs/thmax 52.0 deg
/g4sbs/run 100000