SND@LHC Software
Loading...
Searching...
No Matches
sndLHC_HXgeom_config.py
Go to the documentation of this file.
1import ROOT as r
2import shipunit as u
3from ShipGeoConfig import AttrDict, ConfigRegistry
4
5if "tb_2023_mc" in globals():
6 tb_2023_mc = True
7else: tb_2023_mc = False
8
9with ConfigRegistry.register_config("basic") as c:
10# cave parameters
11 c.cave = AttrDict(z=0*u.cm)
12
13 # Antonia, 482000mm (FASER+2, P3) + 1017mm (DZ) + 245mm (centre emulsion),z=483262./10.*u.cm
14 # centre emulsion now 326.2cm downstream from origin.
15 c.EmulsionDet = AttrDict(z=326.2*u.cm)
16 c.EmulsionDet.PassiveOption = 1 #0 makes emulsion volumes active, 1 makes all emulsion volumes passive
17 c.EmulsionDet.row = 2
18 c.EmulsionDet.col = 2
19 c.EmulsionDet.wall= 5
20 c.EmulsionDet.target = 1 #number of neutrino target volumes
21 c.EmulsionDet.n_plates = 59
22 c.EmulsionDet.EmTh = 0.0070 * u.cm
23 c.EmulsionDet.EmX = 19.2 * u.cm
24 c.EmulsionDet.EmY = 19.2 * u.cm
25 c.EmulsionDet.PBTh = 0.0175 * u.cm
26 c.EmulsionDet.PassiveTh = 0.1 * u.cm
27 c.EmulsionDet.EPlW = 2* c.EmulsionDet.EmTh + c.EmulsionDet.PBTh
28 c.EmulsionDet.AllPW = c.EmulsionDet.PassiveTh + c.EmulsionDet.EPlW
29
30 c.EmulsionDet.BrPackZ = 0.*u.cm
31 c.EmulsionDet.BrPackX = 2*0.05*u.cm
32 c.EmulsionDet.BrPackY = 2*0.05*u.cm
33
34 c.EmulsionDet.BrX = c.EmulsionDet.EmX + c.EmulsionDet.BrPackX
35 c.EmulsionDet.BrY = c.EmulsionDet.EmY + c.EmulsionDet.BrPackY
36 c.EmulsionDet.BrZ = c.EmulsionDet.n_plates * c.EmulsionDet.AllPW + c.EmulsionDet.EPlW + c.EmulsionDet.BrPackZ
37
38 c.EmulsionDet.xdim = 42.2 *u.cm #external wall dimensions
39 c.EmulsionDet.ydim = 42.2 *u.cm
40 c.EmulsionDet.TotalWallZDim = 10.0 *u.cm #wall dimension along z, including border
41 c.EmulsionDet.WallXDim = 38.6 *u.cm #internal wall dimensions
42 c.EmulsionDet.WallYDim = 38.6 *u.cm
43 c.EmulsionDet.WallZDim = 8.15 *u.cm
44 c.EmulsionDet.WallZBorder_offset = 4.75 * u.mm
45 c.EmulsionDet.TTz = 3.0*u.cm
46 c.EmulsionDet.zdim = c.EmulsionDet.wall* c.EmulsionDet.TotalWallZDim + c.EmulsionDet.wall*c.EmulsionDet.TTz
47 c.EmulsionDet.ShiftX = -8.0*u.cm - c.EmulsionDet.xdim/2.
48 c.EmulsionDet.ShiftY = 15.5*u.cm + c.EmulsionDet.ydim/2.
49
50 c.EmulsionDet.startpos = -25.4750 * u.cm + c.EmulsionDet.z
51 c.EmulsionDet.zC = c.EmulsionDet.startpos + c.EmulsionDet.zdim/2.
52
53 # survey points in survey coordinate system!
54 c.EmulsionDet.Xpos0,c.EmulsionDet.Ypos0,c.EmulsionDet.Zpos0 = 53.5*u.mm,2889.2*u.mm,172.0*u.mm
55 c.EmulsionDet.Xpos1,c.EmulsionDet.Ypos1,c.EmulsionDet.Zpos1 = 53.4*u.mm,3019.3*u.mm,172.1*u.mm
56 c.EmulsionDet.Xpos2,c.EmulsionDet.Ypos2,c.EmulsionDet.Zpos2 = 53.3*u.mm,3149.2*u.mm,172.3*u.mm
57 c.EmulsionDet.Xpos3,c.EmulsionDet.Ypos3,c.EmulsionDet.Zpos3 = 53.3*u.mm,3279.2*u.mm,172.4*u.mm
58 c.EmulsionDet.Xpos4,c.EmulsionDet.Ypos4,c.EmulsionDet.Zpos4 = 53.2*u.mm,3409.2*u.mm,172.5*u.mm
59
60 #SciFi parameters
61 c.Scifi = AttrDict(z=0*u.cm)
62 c.Scifi.xdim = 39.0 * u.cm #sensitive only
63 c.Scifi.ydim = 39.0 * u.cm
64 c.Scifi.zdim = c.EmulsionDet.TTz
65 c.Scifi.DZ = c.EmulsionDet.BrZ
66 c.Scifi.nmats = 3
67 c.Scifi.nscifi = 5
68
69 #SiPM parameters
70 c.Scifi.channel_width = 0.25 *u.mm
71 c.Scifi.sipm_edge = 0.17*u.mm # on both ends
72 c.Scifi.charr_gap = 0.2 *u.mm
73 c.Scifi.charr_width = 64 * c.Scifi.channel_width
74 c.Scifi.sipm_diegap = 0.06*u.mm
75 c.Scifi.SiPMarray_width = c.Scifi.sipm_edge+c.Scifi.charr_width+c.Scifi.charr_gap+c.Scifi.charr_width+c.Scifi.sipm_edge
76 c.Scifi.nsipm_channels = 128
77 c.Scifi.nsipm_mat = 4
78 c.Scifi.nsipms = c.Scifi.nsipm_mat * c.Scifi.nmats # 12 per SciFi plane
79 c.Scifi.sipmarr_width = c.Scifi.charr_gap + 2.*c.Scifi.charr_width
80 c.Scifi.firstChannelX = -19.528*u.cm # to agree with SiPM positions of engineering design
81
82 c.Scifi.nfibers_shortrow = 471
83 c.Scifi.nfibers_longrow = 472
84 c.Scifi.nfibers_z = 6
85
86# Guido: Fibre mat total width 500um less wide than 4 SiPM arrays mounted on a single bar, therefore 32.6mm x 4 -0.5=129.9mm
87 c.Scifi.scifimat_width = (c.Scifi.SiPMarray_width+c.Scifi.sipm_diegap)*c.Scifi.nsipm_mat -0.5*u.mm
88 c.Scifi.scifimat_length = c.Scifi.ydim
89 c.Scifi.scifimat_z = 0.135 *u.cm
90 c.Scifi.epoxymat_z = 0.162 *u.cm
91 c.Scifi.scifimat_gap = 0.05 *u.cm
92
93 c.Scifi.fiber_length = c.Scifi.scifimat_length
94 c.Scifi.scintcore_rmax = 0.011 *u.cm
95 c.Scifi.clad1_rmin = c.Scifi.scintcore_rmax
96 c.Scifi.clad1_rmax = 0.01175 *u.cm
97 c.Scifi.clad2_rmin = c.Scifi.clad1_rmax
98 c.Scifi.clad2_rmax = 0.0125 *u.cm
99
100 c.Scifi.horizontal_pitch = 0.0275 *u.cm
101 c.Scifi.vertical_pitch = 0.022 *u.cm
102 c.Scifi.rowlong_offset = 0.035 *u.cm
103 c.Scifi.rowshort_offset = 0.0215 *u.cm
104
105 c.Scifi.carbonfiber_z = 0.02 *u.cm
106 c.Scifi.honeycomb_z = 0.5 *u.cm
107
108 c.Scifi.plastbar_x = 1.5 *u.cm
109 c.Scifi.plastbar_y = c.Scifi.ydim
110 c.Scifi.plastbar_z = 0.195 *u.cm
111
112 c.Scifi.scifi_separation = c.Scifi.zdim + c.EmulsionDet.BrZ
113 c.Scifi.offset_z = - c.EmulsionDet.zdim/2 + c.EmulsionDet.BrZ #SciFi starts at the end of the first ECC
114
115 c.Scifi.timeResol = 150.*u.picosecond
116
117 # absolute edge point positions in survey coordinate system
118 c.Scifi.Xpos0,c.Scifi.Ypos0,c.Scifi.Zpos0 = 39.3*u.mm,2989.7*u.mm,158.2*u.mm
119 c.Scifi.Xpos1,c.Scifi.Ypos1,c.Scifi.Zpos1 = 39.2*u.mm,3119.7*u.mm,158.4*u.mm
120 c.Scifi.Xpos2,c.Scifi.Ypos2,c.Scifi.Zpos2 = 39.1*u.mm,3249.7*u.mm,158.5*u.mm
121 c.Scifi.Xpos3,c.Scifi.Ypos3,c.Scifi.Zpos3 = 39.1*u.mm,3379.7*u.mm,158.6*u.mm
122 c.Scifi.Xpos4,c.Scifi.Ypos4,c.Scifi.Zpos4 = 39.0*u.mm,3509.7*u.mm,158.8*u.mm
123# Scifi technical drawing, distance from first channel on vertical / horizontal plane to edge point 225,225,0, xy plane z perpendicularc.Scifi.ydim
124 c.Scifi.EdgeAX, c.Scifi.EdgeAY, c.Scifi.EdgeAZ = 225*u.mm, 225*u.mm, 0*u.mm
125
126 c.Scifi.FirstChannelVX,c.Scifi.FirstChannelVY,c.Scifi.FirstChannelVZ = -195.28*u.mm, -200.0*u.mm, -12.92*u.mm
127
128 c.Scifi.FirstChannelHX,c.Scifi.FirstChannelHY,c.Scifi.FirstChannelHZ = -200.0*u.mm, -195.28*u.mm, -7.07*u.mm
129
130# Local Scifi position of first channel on a vertical scifi plane in software
131 c.Scifi.LfirstChannelVX,c.Scifi.LfirstChannelVY,c.Scifi.LfirstChannelVZ = -195.135*u.mm, 195.0*u.mm, 11.85*u.mm # sipm on top
132
133# Local Scifi position of last channel (=first channel on technical drawing) on a horizontal scifi plane in software
134 c.Scifi.LfirstChannelHX,c.Scifi.LfirstChannelHY,c.Scifi.LfirstChannelHZ = -195.0*u.mm, 195.178*u.mm, 6.25*u.mm # sipm on the right side direction away from IP1
135
136# Time alignment Scifi, T0 = station 0, mat 0
137 c.Scifi.signalSpeed = 15 * u.cm/u.ns
138
139 c.MuFilter = AttrDict(z=0*u.cm)
140 #coordinates in local gravity based system
141 c.MuFilter.Veto1Dx,c.MuFilter.Veto1Dy,c.MuFilter.Veto1Dz = 40.8*u.mm, 2798.3*u.mm, 192.1*u.mm
142 c.MuFilter.Veto2Dx,c.MuFilter.Veto2Dy,c.MuFilter.Veto2Dz = 40.6*u.mm, 2839.3*u.mm, 172.1*u.mm
143 c.MuFilter.Iron1Dx, c.MuFilter.Iron1Dy, c.MuFilter.Iron1Dz = -22.1*u.mm, 3579.6*u.mm, 146.6*u.mm
144 # US1
145 c.MuFilter.Muon1Dx,c.MuFilter.Muon1Dy,c.MuFilter.Muon1Dz = -46.6*u.mm, 3760.2*u.mm, 128.6 *u.mm
146 c.MuFilter.Iron2Dx, c.MuFilter.Iron2Dy, c.MuFilter.Iron2Dz = -22.1*u.mm, 3804.6*u.mm, 136.6*u.mm
147 # US2
148 c.MuFilter.Muon2Dx,c.MuFilter.Muon2Dy,c.MuFilter.Muon2Dz = -45.7*u.mm, 3984.1*u.mm, 127.6 *u.mm
149 c.MuFilter.Iron3Dx, c.MuFilter.Iron3Dy, c.MuFilter.Iron3Dz = -22.1*u.mm, 4029.6*u.mm, 126.6*u.mm
150 # US3
151 c.MuFilter.Muon3Dx,c.MuFilter.Muon3Dy,c.MuFilter.Muon3Dz = -44.6*u.mm, 4209.5*u.mm, 128.0 *u.mm
152 c.MuFilter.Iron4Dx, c.MuFilter.Iron4Dy, c.MuFilter.Iron4Dz = -22.1*u.mm, 4254.6*u.mm, 116.7*u.mm
153 # US4
154 c.MuFilter.Muon4Dx,c.MuFilter.Muon4Dy,c.MuFilter.Muon4Dz = -45.1*u.mm, 4435.6*u.mm, 128.6 *u.mm
155 c.MuFilter.Iron5Dx, c.MuFilter.Iron5Dy, c.MuFilter.Iron5Dz = -22.1*u.mm, 4479.6*u.mm, 127.7*u.mm
156 # US5
157 c.MuFilter.Muon5Dx,c.MuFilter.Muon5Dy,c.MuFilter.Muon5Dz = -46.8*u.mm, 4663.0*u.mm, 129.9 *u.mm
158 c.MuFilter.Iron6Dx, c.MuFilter.Iron6Dy, c.MuFilter.Iron6Dz = -22.1*u.mm, 4704.6*u.mm, 127.7*u.mm
159 # DS1
160 c.MuFilter.Muon6Dx,c.MuFilter.Muon6Dy,c.MuFilter.Muon6Dz = -45.1*u.mm, 4889.6*u.mm - 1*u.mm, 129.8 *u.mm
161 c.MuFilter.Iron7Dx, c.MuFilter.Iron7Dy, c.MuFilter.Iron7Dz = -22.1*u.mm, 4943.6*u.mm, 127.7*u.mm
162 # DS2
163 c.MuFilter.Muon7Dx,c.MuFilter.Muon7Dy,c.MuFilter.Muon7Dz = -45.2*u.mm, 5125.9*u.mm, 132.8 *u.mm
164 c.MuFilter.Iron8Dx, c.MuFilter.Iron8Dy, c.MuFilter.Iron8Dz = -22.1*u.mm, 5183.6*u.mm, 127.7*u.mm
165 # DS3
166 c.MuFilter.Muon8Dx,c.MuFilter.Muon8Dy,c.MuFilter.Muon8Dz = -7.9*u.mm, 5396.7*u.mm, 132.5 *u.mm
167
168 c.MuFilter.DS4ZGap = 8.82*u.cm
169 # DS4V
170 c.MuFilter.Muon9Dx,c.MuFilter.Muon9Dy,c.MuFilter.Muon9Dz = c.MuFilter.Muon8Dx, c.MuFilter.Muon8Dy + c.MuFilter.DS4ZGap, c.MuFilter.Muon8Dz
171 c.MuFilter.Iron9Dx, c.MuFilter.Iron9Dy, c.MuFilter.Iron9Dz = 177.9*u.mm, 5529.7*u.mm + 1*u.cm, 127.7*u.mm # move downstream by 1cm to avoid overlap
172
173 # relation between edge and bottom bar for VETO
174 c.MuFilter.VETOLocX,c.MuFilter.VETOLocY,c.MuFilter.VETOLocZ = 20.0*u.mm,20.0*u.mm,46.7*u.mm
175
176 # relation between edge and bottom bar for US and DS
177 c.MuFilter.DSHLocX,c.MuFilter.DSHLocY,c.MuFilter.DSHLocZ = 10.5*u.mm, 32.0*u.mm, 11.1*u.mm
178 # relation between edge and right bar
179 c.MuFilter.DSVLocX,c.MuFilter.DSVLocY,c.MuFilter.DSVLocZ = 623.0*u.mm, 47.0*u.mm, 641.3*u.mm
180
181 # offsets in Z of first US bar
182 c.MuFilter.USOffZ1 = 4.35*u.mm
183 c.MuFilter.USOffZ2 = 5.0*u.mm
184 c.MuFilter.USOffZ3 = 4.2*u.mm
185 c.MuFilter.USOffZ4 = 5.21*u.mm
186 c.MuFilter.USOffZ5 = 4.74*u.mm
187
188 #Veto station parameters
189 c.MuFilter.VetonSiPMs = 8
190 c.MuFilter.VetonSides = 2
191 c.MuFilter.NVetoPlanes = 2
192 c.MuFilter.NVetoBars = 7
193
194 c.MuFilter.VetoBarX,c.MuFilter.VetoBarY,c.MuFilter.VetoBarZ = 42 *u.cm, 6 * u.cm, 1 * u.cm
195 c.MuFilter.VetoBarGap = 2*30*u.um # wrapping material
196
197 c.MuFilter.FeX,c.MuFilter.FeY,c.MuFilter.FeZ = 80*u.cm, 60*u.cm, 20*u.cm
198 c.MuFilter.FeEndX,c.MuFilter.FeEndY,c.MuFilter.FeEndZ = 40*u.cm, 40*u.cm, 20*u.cm
199 c.MuFilter.FeBotX,c.MuFilter.FeBotY,c.MuFilter.FeBotZ = 80*u.cm, 9*u.cm, 40*u.cm
200
201 c.MuFilter.UpstreamDetZ = 2.6*u.cm
202 c.MuFilter.UpstreamnSiPMs = 8
203 c.MuFilter.UpstreamnSides = 2
204 c.MuFilter.NUpstreamPlanes = 5
205 c.MuFilter.DownstreamDetZ = 3.9*u.cm
206 c.MuFilter.DownstreamnSiPMs = 1
207 c.MuFilter.DownstreamnSides = 2 # only for horizontal, vertical only one side
208 c.MuFilter.NDownstreamPlanes = 4
209 #upstream bars configuration
210 c.MuFilter.NUpstreamBars = 10
211 c.MuFilter.UpstreamBarX,c.MuFilter.UpstreamBarY,c.MuFilter.UpstreamBarZ = 82.525*u.cm, 6.0*u.cm, 1.0*u.cm
212 c.MuFilter.UpstreamBarGap = 0.1*u.mm
213
214 #downstream bars configuration
215 c.MuFilter.NDownstreamBars = 60 #n.d.r. both for x and y in this case
216 c.MuFilter.DownstreamBarX,c.MuFilter.DownstreamBarY,c.MuFilter.DownstreamBarZ = 82.525*u.cm, 1*u.cm, 1*u.cm
217 c.MuFilter.DownstreamBarGap = 0.1*u.mm
218 c.MuFilter.DownstreamBarX_ver,c.MuFilter.DownstreamBarY_ver,c.MuFilter.DownstreamBarZ_ver = 1*u.cm, 63.525*u.cm, 1*u.cm
219
220 # DS and US support box, inner, Z pointing upward
221 c.MuFilter.SupportBoxD = 0.5*u.mm # empty space between bars and box
222 c.MuFilter.SupportBoxW = 2*u.mm
223 c.MuFilter.DSBoxX1 = c.MuFilter.DSHLocX - c.MuFilter.SupportBoxD
224 c.MuFilter.DSBoxX2 = c.MuFilter.DSHLocX + c.MuFilter.DownstreamBarX + c.MuFilter.SupportBoxD
225 c.MuFilter.DSBoxZ1 = c.MuFilter.DSHLocZ - c.MuFilter.DownstreamBarY/2 - c.MuFilter.SupportBoxD
226 c.MuFilter.DSBoxZ2 = c.MuFilter.DSVLocZ + c.MuFilter.SupportBoxD
227 c.MuFilter.DSBoxY1 = c.MuFilter.DSHLocY - c.MuFilter.DownstreamBarZ/2 - c.MuFilter.SupportBoxD
228 c.MuFilter.DSBoxY2 = c.MuFilter.DSVLocY + c.MuFilter.DownstreamBarZ/2 + c.MuFilter.SupportBoxD
229
230 c.MuFilter.USBoxY1 = c.MuFilter.DSHLocY - c.MuFilter.DownstreamBarZ/2 - c.MuFilter.SupportBoxD
231 c.MuFilter.USBoxY2 = c.MuFilter.DSHLocY + c.MuFilter.DownstreamBarZ/2 + c.MuFilter.SupportBoxD
232
233 # VETO support box
234 c.MuFilter.SupportBoxVW = 4*u.mm
235 c.MuFilter.VETOBoxX1 = c.MuFilter.VETOLocX - c.MuFilter.SupportBoxD
236 c.MuFilter.VETOBoxX2 = c.MuFilter.VETOLocX + c.MuFilter.VetoBarX + c.MuFilter.SupportBoxD
237 c.MuFilter.VETOBoxZ1 = c.MuFilter.VETOLocZ - c.MuFilter.VetoBarY/2 - c.MuFilter.SupportBoxD
238 c.MuFilter.VETOBoxZ2 = c.MuFilter.VETOLocZ + (c.MuFilter.NVetoBars-1)*(c.MuFilter.VetoBarY+c.MuFilter.VetoBarGap) + c.MuFilter.VetoBarY/2 + c.MuFilter.SupportBoxD
239 c.MuFilter.VETOBoxY1 = c.MuFilter.VETOLocY - c.MuFilter.VetoBarZ/2 - c.MuFilter.SupportBoxD
240 c.MuFilter.VETOBoxY2 = c.MuFilter.VETOLocY + c.MuFilter.VetoBarZ/2 + c.MuFilter.SupportBoxD
241
242 # VETO/US/DS plane alignment
243 c.MuFilter.Veto1ShiftY = 0.11 * u.cm
244 c.MuFilter.Veto2ShiftY = -0.04 * u.cm
245 c.MuFilter.US1ShiftY = 0.10 * u.cm
246 c.MuFilter.US2ShiftY = 0.26 * u.cm
247 c.MuFilter.US3ShiftY = 0.24 * u.cm
248 c.MuFilter.US4ShiftY = 0.31 * u.cm
249 c.MuFilter.US5ShiftY = 0.34 * u.cm
250 c.MuFilter.DS1ShiftY = 0.43 * u.cm
251 c.MuFilter.DS1ShiftX = 1.13 * u.cm
252 c.MuFilter.DS2ShiftY = 0.53 * u.cm
253 c.MuFilter.DS2ShiftX = 1.31 * u.cm
254 c.MuFilter.DS3ShiftY = 0.61 * u.cm
255 c.MuFilter.DS3ShiftX = 1.35 * u.cm
256 c.MuFilter.DS4ShiftX = 1.39 * u.cm
257
258 #digitization parameters
259 c.MuFilter.DsAttenuationLength = 350 * u.cm # values between 300 cm and 400cm observed for H6 testbeam
260 c.MuFilter.DsTAttenuationLength = 700 * u.cm # top readout with mirror on bottom
261 c.MuFilter.VandUpAttenuationLength = 999 * u.cm # no significante attenuation observed for H6 testbeam
262 c.MuFilter.VandUpSiPMcalibrationL = 25.*1000. # 1.65 MeV = 41 qcd
263 c.MuFilter.VandUpSiPMcalibrationS = 25.*1000.
264 c.MuFilter.DsSiPMcalibration = 25.*1000.
265 c.MuFilter.timeResol = 150.*u.picosecond
266 c.MuFilter.VandUpPropSpeed = 12.5*u.cm/u.nanosecond
267 c.MuFilter.DsPropSpeed = 14.3*u.cm/u.nanosecond
268
269 c.Floor = AttrDict(z=48000.*u.cm) # to place tunnel in SND_@LHC coordinate system
270 c.Floor.DX = 1.0*u.cm
271 c.Floor.DY = -4.5*u.cm # subtract 4.5cm to avoid overlaps
272 c.Floor.DZ = 0.
273
274 #COLDBOX configuration
275 c.Floor.Acrylic_width = 5.0*u.cm
276 c.Floor.BPoly_width = 4.0*u.cm
277 c.Floor.CBFrontWall_xdim = 219.*u.cm
278 c.Floor.CBFrontWall_ydim = 170.72*u.cm-c.Floor.Acrylic_width
279 c.Floor.CBLatWall_zdim = 176.0*u.cm
280 c.Floor.CBTiny_zdim = 17.0*u.cm
281 c.Floor.CBExtra_zdim = 41.0*u.cm
282 c.Floor.CBExtra_xdim = 67.5 *u.cm
283 c.Floor.SlopedWall_zproj = 110.0*u.cm
284 c.Floor.MFeBlockX = c.MuFilter.FeX
285 c.Floor.MFeBlockY = c.MuFilter.FeY
286 c.Floor.MFeBlockZ = c.MuFilter.FeZ
287
288# for H6 / H8 / HX testbeam and commissioning
289 HX = True
290 if HX:
291 c.Floor.z = 0 # no tunnel, no slope
292
293 # no Veto
294 c.MuFilter.NVetoPlanes = 0
295
296 # only 1 DS station
297 c.MuFilter.NDownstreamPlanes = 1
298 # remove most downstream iron block and its bottom part
299 c.MuFilter.FeEndX,c.MuFilter.FeEndY,c.MuFilter.FeEndZ = 0*u.cm, 0*u.cm, 0*u.cm
300 c.MuFilter.FeBotX,c.MuFilter.FeBotY,c.MuFilter.FeBotZ = 0*u.cm, 0*u.cm, 0*u.cm
301
302 # removing the emulsion+W target
303 c.EmulsionDet.wall = 0
304
305 # set baby SciFi modules
306 # they have much simpler internal structure compared to the TI18 SciFi
307 # most notably they have 7 layer of fibers per mat
308 c.Scifi.xdim = 13.0*u.cm #sensitive only
309 c.Scifi.ydim = 13.0*u.cm
310 c.Scifi.zdim = 3*u.cm # maybe not needed
311 c.Scifi.nmats = 1
312 c.Scifi.nscifi = 4
313 c.Scifi.nfibers_z = 7
314 c.Scifi.scifimat_z = 0.16*u.cm # thickness of a mat of 7 fiber layers
315 c.Scifi.scifimat_length = c.Scifi.ydim
316 c.Scifi.fiber_length = c.Scifi.scifimat_length
317 c.Scifi.plane_gap = 6*u.mm # an air gap btw X an Y planes, in TI18 case this is controlled via c.Scifi.plastbar_z
318 c.Scifi.tedlar_to_plane = 5.42*u.mm # an air gap btw protective tedlar sheet and a sensitive plane
319 #c.Scifi.tedlar_z = 50*u.um # not included in the sw detector model
320 # offset between the edge of a baby module frame and the upstream tedlar inside it
321 c.Scifi.frame_offset = 0.8*u.cm
322 # offsets between the upstream tedlar sheet and the upstream passive block
323 c.Scifi.station_offset1 = c.Scifi.frame_offset
324 c.Scifi.station_offset2 = c.Scifi.frame_offset+1.2*u.cm
325 c.Scifi.station_offset3 = c.Scifi.frame_offset+0.6*u.cm
326
327 # add 3 variable size iron blocks upstream of SciFi stations
328 # complete removal of a wall goes with commenting the respective line below
329 c.Scifi.FeTargetX2, c.Scifi.FeTargetY2, c.Scifi.FeTargetZ2 = 30.*u.cm, 10.*u.cm, 30.*u.cm
330 c.Scifi.FeTargetX3, c.Scifi.FeTargetY3, c.Scifi.FeTargetZ3 = 30.*u.cm, 10.*u.cm, 30.*u.cm
331 c.Scifi.FeTargetX4, c.Scifi.FeTargetY4, c.Scifi.FeTargetZ4 = 30.*u.cm, 10.*u.cm, 30.*u.cm
332
333# absolute edge point positions in survey coordinate system (survey is 'by eye' for now)
334 c.Scifi.Xpos0,c.Scifi.Ypos0,c.Scifi.Zpos0 = 156.6*u.mm,3136.48*u.mm,218.2*u.mm
335 c.Scifi.Xpos1,c.Scifi.Ypos1,c.Scifi.Zpos1 = 156.6*u.mm,3264.52*u.mm,218.4*u.mm
336 c.Scifi.Xpos2,c.Scifi.Ypos2,c.Scifi.Zpos2 = 156.6*u.mm,3404.56*u.mm,218.5*u.mm
337 c.Scifi.Xpos3,c.Scifi.Ypos3,c.Scifi.Zpos3 = 156.6*u.mm,3538.6*u.mm,218.6*u.mm
338# and MuFilter - survey measurements
339 c.MuFilter.Iron1Dx, c.MuFilter.Iron1Dy, c.MuFilter.Iron1Dz = -22.1*u.mm, 3559.6*u.mm, 128.6*u.mm
340 # US1
341 c.MuFilter.Muon1Dx,c.MuFilter.Muon1Dy,c.MuFilter.Muon1Dz = -46.6*u.mm, 3740.2*u.mm, 128.6 *u.mm
342 c.MuFilter.Iron2Dx, c.MuFilter.Iron2Dy, c.MuFilter.Iron2Dz = -22.1*u.mm, 3784.6*u.mm, 128.6*u.mm
343 # US2
344 c.MuFilter.Muon2Dx,c.MuFilter.Muon2Dy,c.MuFilter.Muon2Dz = -46.6*u.mm, 3964.1*u.mm, 128.6 *u.mm
345 c.MuFilter.Iron3Dx, c.MuFilter.Iron3Dy, c.MuFilter.Iron3Dz = -22.1*u.mm, 4009.6*u.mm, 128.6*u.mm
346 # US3
347 c.MuFilter.Muon3Dx,c.MuFilter.Muon3Dy,c.MuFilter.Muon3Dz = -46.6*u.mm, 4189.5*u.mm, 128.6 *u.mm
348 c.MuFilter.Iron4Dx, c.MuFilter.Iron4Dy, c.MuFilter.Iron4Dz = -22.1*u.mm, 4234.6*u.mm, 128.6*u.mm
349 # US4
350 c.MuFilter.Muon4Dx,c.MuFilter.Muon4Dy,c.MuFilter.Muon4Dz = -46.1*u.mm, 4415.6*u.mm, 128.6 *u.mm
351 c.MuFilter.Iron5Dx, c.MuFilter.Iron5Dy, c.MuFilter.Iron5Dz = -22.1*u.mm, 4459.6*u.mm, 128.6*u.mm
352 # US5
353 c.MuFilter.Muon5Dx,c.MuFilter.Muon5Dy,c.MuFilter.Muon5Dz = -46.6*u.mm, 4643.0*u.mm, 128.6*u.mm
354 c.MuFilter.Iron6Dx, c.MuFilter.Iron6Dy, c.MuFilter.Iron6Dz = -22.1*u.mm, 4684.6*u.mm, 128.6*u.mm
355 # DS1
356 c.MuFilter.Muon6Dx,c.MuFilter.Muon6Dy,c.MuFilter.Muon6Dz = -46.6*u.mm, 4869.6*u.mm - 1*u.mm, 128.6 *u.mm
357 c.MuFilter.Iron7Dx, c.MuFilter.Iron7Dy, c.MuFilter.Iron7Dz = -22.1*u.mm, 4923.6*u.mm, 128.6*u.mm
358
359# fine adjustements
360 #c.MuFilter.Muon8Dx = -4.7*u.mm - 4*u.cm
361 #c.MuFilter.Muon7Dx = -47.1*u.mm + 0.74*u.cm
362 #c.MuFilter.Muon1Dz = 146.6*u.mm - 2.3*u.cm
363 #c.MuFilter.Muon2Dz = 136.6*u.mm - 1.5*u.cm
364 #c.MuFilter.Muon5Dz = 127.7*u.mm + 0.25*u.cm
365 #c.MuFilter.Muon7Dz = 127.7*u.mm + 0.25*u.cm
366 #c.MuFilter.Muon8Dz = 127.7*u.mm + 0.56*u.cm
367
368# from Scifi track alignment
369 if tb_2023_mc :
370# spatial
371 c.Scifi.LocM100, c.Scifi.LocM110 = 0.000*u.um, 0.000*u.um
372 c.Scifi.LocM200, c.Scifi.LocM210 = 0.000*u.um, 0.000*u.um
373 c.Scifi.LocM300, c.Scifi.LocM310 = 0.000*u.um, 0.000*u.um
374 c.Scifi.LocM400, c.Scifi.LocM410 = 0.000*u.um, 0.000*u.um
375 c.Scifi.RotPhiS10,c.Scifi.RotPsiS10,c.Scifi.RotThetaS10 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
376 c.Scifi.RotPhiS11,c.Scifi.RotPsiS11,c.Scifi.RotThetaS11 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
377 c.Scifi.RotPhiS20,c.Scifi.RotPsiS20,c.Scifi.RotThetaS20 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
378 c.Scifi.RotPhiS21,c.Scifi.RotPsiS21,c.Scifi.RotThetaS21 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
379 c.Scifi.RotPhiS30,c.Scifi.RotPsiS30,c.Scifi.RotThetaS30 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
380 c.Scifi.RotPhiS31,c.Scifi.RotPsiS31,c.Scifi.RotThetaS31 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
381 c.Scifi.RotPhiS40,c.Scifi.RotPsiS40,c.Scifi.RotThetaS40 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
382 c.Scifi.RotPhiS41,c.Scifi.RotPsiS41,c.Scifi.RotThetaS41 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
383# Time alignment Scifi, T0 = station 0, mat 0
384 c.Scifi.station1t,c.Scifi.station1H0t,c.Scifi.station1V0t = 0.000*u.ns, 0.000*u.ns, 0.000*u.ns
385 c.Scifi.station2t,c.Scifi.station2H0t,c.Scifi.station2V0t = 0.000*u.ns, 0.000*u.ns, 0.000*u.ns
386 c.Scifi.station3t,c.Scifi.station3H0t,c.Scifi.station3V0t = 0.000*u.ns, 0.000*u.ns, 0.000*u.ns
387 c.Scifi.station4t,c.Scifi.station4H0t,c.Scifi.station4V0t = 0.000*u.ns, 0.000*u.ns, 0.000*u.ns
388 else :
389# spatial
390 c.Scifi.LocM100, c.Scifi.LocM110 = 667.33*u.um, 429.73*u.um
391 c.Scifi.LocM200, c.Scifi.LocM210 = 0.000*u.um, 51.80*u.um
392 c.Scifi.LocM300, c.Scifi.LocM310 = 32.36*u.um, 803.10*u.um
393 c.Scifi.LocM400, c.Scifi.LocM410 = 105.10*u.um, 371.07*u.um
394 c.Scifi.RotPhiS10,c.Scifi.RotPsiS10,c.Scifi.RotThetaS10 = 0.00*u.mrad, 3.29*u.mrad, 0.00*u.mrad
395 c.Scifi.RotPhiS11,c.Scifi.RotPsiS11,c.Scifi.RotThetaS11 = 0.00*u.mrad, -0.89*u.mrad, 0.00*u.mrad
396 c.Scifi.RotPhiS20,c.Scifi.RotPsiS20,c.Scifi.RotThetaS20 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
397 c.Scifi.RotPhiS21,c.Scifi.RotPsiS21,c.Scifi.RotThetaS21 = 0.00*u.mrad, 0.25*u.mrad, 0.00*u.mrad
398 c.Scifi.RotPhiS30,c.Scifi.RotPsiS30,c.Scifi.RotThetaS30 = 0.00*u.mrad, 0.93*u.mrad, 0.00*u.mrad
399 c.Scifi.RotPhiS31,c.Scifi.RotPsiS31,c.Scifi.RotThetaS31 = 0.00*u.mrad, 2.86*u.mrad, 0.00*u.mrad
400 c.Scifi.RotPhiS40,c.Scifi.RotPsiS40,c.Scifi.RotThetaS40 = 0.00*u.mrad, 0.15*u.mrad, 0.00*u.mrad
401 c.Scifi.RotPhiS41,c.Scifi.RotPsiS41,c.Scifi.RotThetaS41 = 0.00*u.mrad, 3.13*u.mrad, 0.00*u.mrad
402# Time alignment Scifi, T0 = station 0, mat 0
403 c.Scifi.station1t,c.Scifi.station1H0t,c.Scifi.station1V0t = 0.000*u.ns, 0.000*u.ns, 1.535*u.ns
404 c.Scifi.station2t,c.Scifi.station2H0t,c.Scifi.station2V0t = 1.246*u.ns, 0.000*u.ns, -0.281*u.ns
405 c.Scifi.station3t,c.Scifi.station3H0t,c.Scifi.station3V0t = 3.240*u.ns, 0.000*u.ns, -0.913*u.ns
406 c.Scifi.station4t,c.Scifi.station4H0t,c.Scifi.station4V0t = 1.750*u.ns, 0.000*u.ns, -0.145*u.ns