-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathmyQweakCerenkovOnly.mac
141 lines (111 loc) · 4.33 KB
/
myQweakCerenkovOnly.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
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
#=================================================
# Buddhini W. 01-27-2015
# Macro file myQweakCerenkovOnly.mac
# For simulating PMT DD effect
#
# Defaults to production LH2 simulation
#=================================================
#==========================
# Store Tracks
#==========================
/tracking/storeTrajectory 1
#==========================
# Cerenkov Detector related
#==========================
# Air :without Cerenkov light production
# Quartz :with Cerenkov light production
#
# Angle is definedagaing vertical (y) axis
# - a tilt angle of 0.0*deg means vertical
# - negative angle will tilt detector towards target
/Cerenkov/SetCerenkovMaterial Quartz
/Cerenkov/SetTiltingAngle 0.0 degree
/Cerenkov/SetPreradiatorMaterial PBA
# Define the number of Cereknov Bars
/Cerenkov/SetNumberOfDetectors 8
# Move Cerenkov Detectors to actual position
/Cerenkov/Cerenkov1/SetCenterPositionInY 334.724 cm
/Cerenkov/Cerenkov1/SetCenterPositionInZ 576.665 cm
/Cerenkov/Cerenkov2/SetCenterPositionInY 334.742 cm
/Cerenkov/Cerenkov2/SetCenterPositionInZ 576.705 cm
/Cerenkov/Cerenkov3/SetCenterPositionInY 334.738 cm
/Cerenkov/Cerenkov3/SetCenterPositionInZ 577.020 cm
/Cerenkov/Cerenkov4/SetCenterPositionInY 335.036 cm
/Cerenkov/Cerenkov4/SetCenterPositionInZ 577.425 cm
/Cerenkov/Cerenkov5/SetCenterPositionInY 335.074 cm
/Cerenkov/Cerenkov5/SetCenterPositionInZ 577.515 cm
/Cerenkov/Cerenkov6/SetCenterPositionInY 335.252 cm
/Cerenkov/Cerenkov6/SetCenterPositionInZ 577.955 cm
/Cerenkov/Cerenkov7/SetCenterPositionInY 335.224 cm
/Cerenkov/Cerenkov7/SetCenterPositionInZ 577.885 cm
/Cerenkov/Cerenkov8/SetCenterPositionInY 334.888 cm
/Cerenkov/Cerenkov8/SetCenterPositionInZ 577.060 cm
# ELOG Detector/74: Positions of background detectors
#
# For practical reasons (MD9 is in the same octant as MD1) the background
# detector is here in the opposite octant than it actually was in (behind MD5).
#
# PRELIMINARY: More accurate numbers will require both the final survey
# and the local survey of the monument with respect to the detector itself.
#
# "Det. Origin" is (assumed to be) the center of the back of the bar, outer edge.
#
# Z position of Det. Origin = 737.86 cm
# - 1 cm (half thickness) = 736.86 cm (Z in MC)
# X position of Det. Origin = 377.84 cm
# - 9 cm (half width bar) = 368.84 cm (Y in MC)
# Y position of Det. Origin = 72.86 cm (X in MC)
#/Cerenkov/Cerenkov9/SetCenterPositionInX -72.86 cm
#/Cerenkov/Cerenkov9/SetCenterPositionInY 368.84 cm
#/Cerenkov/Cerenkov9/SetCenterPositionInZ 736.86 cm
#====================================
# Event generator related
#====================================
# Usually this will be changed by a user macro that calls this
# configuration macro.
# List of Generators (Included in src/QweakSimEPEvent.cc):
# 0 - Combination of all reactions
# 1 - LH2 elastic (default)
# 2 - Al elastic
# 3 - Al quasi-elastic (proton)
# 4 - Al quasi-elastic (neutron)
# 5 - LH2 inelastic (delta resonance)
# 6 - Moller scattering
# 7 - LH2 radiative lookup table (3.35 GeV)
# 8 - Al quasi-elastic, Bosted
# 88 - pion
# 9 - Al nuclear inelastics SPS
# 10 - Al GDR
# 11 - Al inelastic (Bosted)
# 12-17 - Al alloy (Zn, Mg, Cu, Cr, Fe, Si, respectively)
/EventGen/SelectReactionType 1
#==========================
# Octant Selection
#==========================
# Valid values range from 0 to 12
# default 0: events thrown in all octants
# 1-8: events thrown in corresponding octant
# 9-12: events thrown in two octants (9 for 1+5, 12 for 4+8)
#/EventGen/SelectOctant 11
#==========================
# Primary generator
#==========================
# If below set points are not given, the beam starts
# at (0,0,0) (inside QTOR)
#
# place beam after pre-radiator of MD3
/PrimaryEvent/SetBeamPositionX 0.0 cm
/PrimaryEvent/SetBeamPositionY 335.0 cm
/PrimaryEvent/SetBeamPositionZ 575.0 cm
# Set energy of primary particle
/EventGen/SetBeamEnergy 50 MeV
#====================================
# Random Seed Handling
#====================================
# Inform the RunManager to systematically save the seed
/random/setSavingFlag 0
# Copy the following command line to your user macro and uncomment it
# to use the seeds of the last event in last run
#/random/resetEngineFrom currentEvent.rndm
# Or use specific seeds from file 'myseed.rndm', for example
#/random/resetEngineFrom myseed.rndm