SND@LHC Software
Loading...
Searching...
No Matches
sndLHC_H4geom_config.py
Go to the documentation of this file.
1import ROOT as r
2import shipunit as u
3from ShipGeoConfig import AttrDict, ConfigRegistry
4
5with_tungsten = False
6if "tb_2024_mc" in globals():
7 tb_2024_mc = True
8 if globals()["target_material"]=="W":
9 with_tungsten = True
10else: tb_2024_mc = False
11
12with ConfigRegistry.register_config("basic") as c:
13# cave parameters
14 c.cave = AttrDict(z=0*u.cm)
15
16 # Antonia, 482000mm (FASER+2, P3) + 1017mm (DZ) + 245mm (centre emulsion),z=483262./10.*u.cm
17 # centre emulsion now 326.2cm downstream from origin.
18 c.EmulsionDet = AttrDict(z=326.2*u.cm)
19 c.EmulsionDet.PassiveOption = 1 #0 makes emulsion volumes active, 1 makes all emulsion volumes passive
20 c.EmulsionDet.row = 2
21 c.EmulsionDet.col = 2
22 c.EmulsionDet.wall= 5
23 c.EmulsionDet.target = 1 #number of neutrino target volumes
24 c.EmulsionDet.n_plates = 59
25 c.EmulsionDet.EmTh = 0.0070 * u.cm
26 c.EmulsionDet.EmX = 19.2 * u.cm
27 c.EmulsionDet.EmY = 19.2 * u.cm
28 c.EmulsionDet.PBTh = 0.0175 * u.cm
29 c.EmulsionDet.PassiveTh = 0.1 * u.cm
30 c.EmulsionDet.EPlW = 2* c.EmulsionDet.EmTh + c.EmulsionDet.PBTh
31 c.EmulsionDet.AllPW = c.EmulsionDet.PassiveTh + c.EmulsionDet.EPlW
32
33 c.EmulsionDet.BrPackZ = 0.*u.cm
34 c.EmulsionDet.BrPackX = 2*0.05*u.cm
35 c.EmulsionDet.BrPackY = 2*0.05*u.cm
36
37 c.EmulsionDet.BrX = c.EmulsionDet.EmX + c.EmulsionDet.BrPackX
38 c.EmulsionDet.BrY = c.EmulsionDet.EmY + c.EmulsionDet.BrPackY
39 c.EmulsionDet.BrZ = c.EmulsionDet.n_plates * c.EmulsionDet.AllPW + c.EmulsionDet.EPlW + c.EmulsionDet.BrPackZ
40
41 c.EmulsionDet.xdim = 42.2 *u.cm #external wall dimensions
42 c.EmulsionDet.ydim = 42.2 *u.cm
43 c.EmulsionDet.TotalWallZDim = 10.0 *u.cm #wall dimension along z, including border
44 c.EmulsionDet.WallXDim = 38.6 *u.cm #internal wall dimensions
45 c.EmulsionDet.WallYDim = 38.6 *u.cm
46 c.EmulsionDet.WallZDim = 8.15 *u.cm
47 c.EmulsionDet.WallZBorder_offset = 4.75 * u.mm
48 c.EmulsionDet.TTz = 3.0*u.cm
49 c.EmulsionDet.zdim = c.EmulsionDet.wall* c.EmulsionDet.TotalWallZDim + c.EmulsionDet.wall*c.EmulsionDet.TTz
50 c.EmulsionDet.ShiftX = -8.0*u.cm - c.EmulsionDet.xdim/2.
51 c.EmulsionDet.ShiftY = 15.5*u.cm + c.EmulsionDet.ydim/2.
52
53 c.EmulsionDet.startpos = -25.4750 * u.cm + c.EmulsionDet.z
54 c.EmulsionDet.zC = c.EmulsionDet.startpos + c.EmulsionDet.zdim/2.
55
56 # survey points in survey coordinate system!
57 c.EmulsionDet.Xpos0,c.EmulsionDet.Ypos0,c.EmulsionDet.Zpos0 = 53.5*u.mm,2889.2*u.mm,172.0*u.mm
58 c.EmulsionDet.Xpos1,c.EmulsionDet.Ypos1,c.EmulsionDet.Zpos1 = 53.4*u.mm,3019.3*u.mm,172.1*u.mm
59 c.EmulsionDet.Xpos2,c.EmulsionDet.Ypos2,c.EmulsionDet.Zpos2 = 53.3*u.mm,3149.2*u.mm,172.3*u.mm
60 c.EmulsionDet.Xpos3,c.EmulsionDet.Ypos3,c.EmulsionDet.Zpos3 = 53.3*u.mm,3279.2*u.mm,172.4*u.mm
61 c.EmulsionDet.Xpos4,c.EmulsionDet.Ypos4,c.EmulsionDet.Zpos4 = 53.2*u.mm,3409.2*u.mm,172.5*u.mm
62
63 #SciFi parameters
64 c.Scifi = AttrDict(z=0*u.cm)
65 c.Scifi.xdim = 39.0 * u.cm #sensitive only
66 c.Scifi.ydim = 39.0 * u.cm
67 c.Scifi.zdim = c.EmulsionDet.TTz
68 c.Scifi.DZ = c.EmulsionDet.BrZ
69 c.Scifi.nmats = 3
70 c.Scifi.nscifi = 5
71
72 #SiPM parameters
73 c.Scifi.channel_width = 0.25 *u.mm
74 c.Scifi.sipm_edge = 0.17*u.mm # on both ends
75 c.Scifi.charr_gap = 0.2 *u.mm
76 c.Scifi.charr_width = 64 * c.Scifi.channel_width
77 c.Scifi.sipm_diegap = 0.06*u.mm
78 c.Scifi.SiPMarray_width = c.Scifi.sipm_edge+c.Scifi.charr_width+c.Scifi.charr_gap+c.Scifi.charr_width+c.Scifi.sipm_edge
79 c.Scifi.nsipm_channels = 128
80 c.Scifi.nsipm_mat = 4
81 c.Scifi.nsipms = c.Scifi.nsipm_mat * c.Scifi.nmats # 12 per SciFi plane
82 c.Scifi.sipmarr_width = c.Scifi.charr_gap + 2.*c.Scifi.charr_width
83 c.Scifi.firstChannelX = -19.528*u.cm # to agree with SiPM positions of engineering design
84
85 c.Scifi.nfibers_shortrow = 471
86 c.Scifi.nfibers_longrow = 472
87 c.Scifi.nfibers_z = 6
88
89# 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
90 c.Scifi.scifimat_width = (c.Scifi.SiPMarray_width+c.Scifi.sipm_diegap)*c.Scifi.nsipm_mat -0.5*u.mm
91 c.Scifi.scifimat_length = c.Scifi.ydim
92 c.Scifi.scifimat_z = 0.135 *u.cm
93 c.Scifi.epoxymat_z = 0.162 *u.cm
94 c.Scifi.scifimat_gap = 0.05 *u.cm
95
96 c.Scifi.fiber_length = c.Scifi.scifimat_length
97 c.Scifi.scintcore_rmax = 0.011 *u.cm
98 c.Scifi.clad1_rmin = c.Scifi.scintcore_rmax
99 c.Scifi.clad1_rmax = 0.01175 *u.cm
100 c.Scifi.clad2_rmin = c.Scifi.clad1_rmax
101 c.Scifi.clad2_rmax = 0.0125 *u.cm
102
103 c.Scifi.horizontal_pitch = 0.0275 *u.cm
104 c.Scifi.vertical_pitch = 0.022 *u.cm
105 c.Scifi.rowlong_offset = 0.035 *u.cm
106 c.Scifi.rowshort_offset = 0.0215 *u.cm
107
108 c.Scifi.carbonfiber_z = 0.02 *u.cm
109 c.Scifi.honeycomb_z = 0.5 *u.cm
110
111 c.Scifi.plastbar_x = 1.5 *u.cm
112 c.Scifi.plastbar_y = c.Scifi.ydim
113 c.Scifi.plastbar_z = 0.195 *u.cm
114
115 c.Scifi.scifi_separation = c.Scifi.zdim + c.EmulsionDet.BrZ
116 c.Scifi.offset_z = - c.EmulsionDet.zdim/2 + c.EmulsionDet.BrZ #SciFi starts at the end of the first ECC
117
118 c.Scifi.timeResol = 150.*u.picosecond
119
120 # absolute edge point positions in survey coordinate system
121 c.Scifi.Xpos0,c.Scifi.Ypos0,c.Scifi.Zpos0 = 39.3*u.mm,2989.7*u.mm,158.2*u.mm
122 c.Scifi.Xpos1,c.Scifi.Ypos1,c.Scifi.Zpos1 = 39.2*u.mm,3119.7*u.mm,158.4*u.mm
123 c.Scifi.Xpos2,c.Scifi.Ypos2,c.Scifi.Zpos2 = 39.1*u.mm,3249.7*u.mm,158.5*u.mm
124 c.Scifi.Xpos3,c.Scifi.Ypos3,c.Scifi.Zpos3 = 39.1*u.mm,3379.7*u.mm,158.6*u.mm
125 c.Scifi.Xpos4,c.Scifi.Ypos4,c.Scifi.Zpos4 = 39.0*u.mm,3509.7*u.mm,158.8*u.mm
126# Scifi technical drawing, distance from first channel on vertical / horizontal plane to edge point 225,225,0, xy plane z perpendicularc.Scifi.ydim
127 c.Scifi.EdgeAX, c.Scifi.EdgeAY, c.Scifi.EdgeAZ = 225*u.mm, 225*u.mm, 0*u.mm
128
129 c.Scifi.FirstChannelVX,c.Scifi.FirstChannelVY,c.Scifi.FirstChannelVZ = -195.28*u.mm, -200.0*u.mm, -12.92*u.mm
130
131 c.Scifi.FirstChannelHX,c.Scifi.FirstChannelHY,c.Scifi.FirstChannelHZ = -200.0*u.mm, -195.28*u.mm, -7.07*u.mm
132
133# Local Scifi position of first channel on a vertical scifi plane in software
134 c.Scifi.LfirstChannelVX,c.Scifi.LfirstChannelVY,c.Scifi.LfirstChannelVZ = -195.135*u.mm, 195.0*u.mm, 11.85*u.mm # sipm on top
135
136# Local Scifi position of last channel (=first channel on technical drawing) on a horizontal scifi plane in software
137 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
138
139# Time alignment Scifi, T0 = station 0, mat 0
140 c.Scifi.signalSpeed = 15 * u.cm/u.ns
141
142 c.MuFilter = AttrDict(z=0*u.cm)
143 #coordinates in local gravity based system
144 c.MuFilter.Veto1Dx,c.MuFilter.Veto1Dy,c.MuFilter.Veto1Dz = 40.8*u.mm, 2798.3*u.mm, 192.1*u.mm
145 c.MuFilter.Veto2Dx,c.MuFilter.Veto2Dy,c.MuFilter.Veto2Dz = 40.6*u.mm, 2839.3*u.mm, 172.1*u.mm
146 c.MuFilter.Iron1Dx, c.MuFilter.Iron1Dy, c.MuFilter.Iron1Dz = -22.1*u.mm, 3579.6*u.mm, 146.6*u.mm
147 # US1
148 c.MuFilter.Muon1Dx,c.MuFilter.Muon1Dy,c.MuFilter.Muon1Dz = -46.6*u.mm, 3760.2*u.mm, 128.6 *u.mm
149 c.MuFilter.Iron2Dx, c.MuFilter.Iron2Dy, c.MuFilter.Iron2Dz = -22.1*u.mm, 3804.6*u.mm, 136.6*u.mm
150 # US2
151 c.MuFilter.Muon2Dx,c.MuFilter.Muon2Dy,c.MuFilter.Muon2Dz = -45.7*u.mm, 3984.1*u.mm, 127.6 *u.mm
152 c.MuFilter.Iron3Dx, c.MuFilter.Iron3Dy, c.MuFilter.Iron3Dz = -22.1*u.mm, 4029.6*u.mm, 1318.6*u.mm
153 # US3
154 c.MuFilter.Muon3Dx,c.MuFilter.Muon3Dy,c.MuFilter.Muon3Dz = -44.6*u.mm, 4209.5*u.mm, 128.0 *u.mm
155 c.MuFilter.Iron4Dx, c.MuFilter.Iron4Dy, c.MuFilter.Iron4Dz = -22.1*u.mm, 4254.6*u.mm, 116.7*u.mm
156 # US4
157 c.MuFilter.Muon4Dx,c.MuFilter.Muon4Dy,c.MuFilter.Muon4Dz = -45.1*u.mm, 4435.6*u.mm, 128.6 *u.mm
158 c.MuFilter.Iron5Dx, c.MuFilter.Iron5Dy, c.MuFilter.Iron5Dz = -22.1*u.mm, 4479.6*u.mm, 127.7*u.mm
159 # US5
160 c.MuFilter.Muon5Dx,c.MuFilter.Muon5Dy,c.MuFilter.Muon5Dz = -46.8*u.mm, 4663.0*u.mm, 129.9 *u.mm
161 c.MuFilter.Iron6Dx, c.MuFilter.Iron6Dy, c.MuFilter.Iron6Dz = -22.1*u.mm, 4704.6*u.mm, 127.7*u.mm
162 # DS1
163 c.MuFilter.Muon6Dx,c.MuFilter.Muon6Dy,c.MuFilter.Muon6Dz = -45.1*u.mm, 4889.6*u.mm - 1*u.mm, 129.8 *u.mm
164 c.MuFilter.Iron7Dx, c.MuFilter.Iron7Dy, c.MuFilter.Iron7Dz = -22.1*u.mm, 4943.6*u.mm, 127.7*u.mm
165 # DS2
166 c.MuFilter.Muon7Dx,c.MuFilter.Muon7Dy,c.MuFilter.Muon7Dz = -45.2*u.mm, 5125.9*u.mm, 132.8 *u.mm
167 c.MuFilter.Iron8Dx, c.MuFilter.Iron8Dy, c.MuFilter.Iron8Dz = -22.1*u.mm, 5183.6*u.mm, 127.7*u.mm
168 # DS3
169 c.MuFilter.Muon8Dx,c.MuFilter.Muon8Dy,c.MuFilter.Muon8Dz = -7.9*u.mm, 5396.7*u.mm, 132.5 *u.mm
170
171 c.MuFilter.DS4ZGap = 8.82*u.cm
172 # DS4V
173 c.MuFilter.Muon9Dx,c.MuFilter.Muon9Dy,c.MuFilter.Muon9Dz = c.MuFilter.Muon8Dx, c.MuFilter.Muon8Dy + c.MuFilter.DS4ZGap, c.MuFilter.Muon8Dz
174 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
175
176 # relation between edge and bottom bar for VETO
177 c.MuFilter.VETOLocX,c.MuFilter.VETOLocY,c.MuFilter.VETOLocZ = 20.0*u.mm,20.0*u.mm,46.7*u.mm
178
179 # relation between edge and bottom bar for US and DS
180 c.MuFilter.DSHLocX,c.MuFilter.DSHLocY,c.MuFilter.DSHLocZ = 10.5*u.mm, 32.0*u.mm, 11.1*u.mm
181 # relation between edge and right bar
182 c.MuFilter.DSVLocX,c.MuFilter.DSVLocY,c.MuFilter.DSVLocZ = 623.0*u.mm, 47.0*u.mm, 641.3*u.mm
183
184 # offsets in Z of first US bar
185 c.MuFilter.USOffZ1 = 4.35*u.mm
186 c.MuFilter.USOffZ2 = 5.0*u.mm
187 c.MuFilter.USOffZ3 = 4.2*u.mm
188 c.MuFilter.USOffZ4 = 5.21*u.mm
189 c.MuFilter.USOffZ5 = 4.74*u.mm
190
191 #Veto station parameters
192 c.MuFilter.VetonSiPMs = 8
193 c.MuFilter.VetonSides = 2
194 c.MuFilter.NVetoPlanes = 2
195 c.MuFilter.NVetoBars = 7
196
197 c.MuFilter.VetoBarX,c.MuFilter.VetoBarY,c.MuFilter.VetoBarZ = 42 *u.cm, 6 * u.cm, 1 * u.cm
198 c.MuFilter.VetoBarGap = 2*30*u.um # wrapping material
199
200 c.MuFilter.FeX,c.MuFilter.FeY,c.MuFilter.FeZ = 80*u.cm, 60*u.cm, 20*u.cm
201 c.MuFilter.FeEndX,c.MuFilter.FeEndY,c.MuFilter.FeEndZ = 40*u.cm, 40*u.cm, 20*u.cm
202 c.MuFilter.FeBotX,c.MuFilter.FeBotY,c.MuFilter.FeBotZ = 80*u.cm, 9*u.cm, 40*u.cm
203
204 c.MuFilter.UpstreamDetZ = 2.6*u.cm
205 c.MuFilter.UpstreamnSiPMs = 8
206 c.MuFilter.UpstreamnSides = 2
207 c.MuFilter.NUpstreamPlanes = 5
208 c.MuFilter.DownstreamDetZ = 3.9*u.cm
209 c.MuFilter.DownstreamnSiPMs = 1
210 c.MuFilter.DownstreamnSides = 2 # only for horizontal, vertical only one side
211 c.MuFilter.NDownstreamPlanes = 4
212 #upstream bars configuration
213 c.MuFilter.NUpstreamBars = 10
214 c.MuFilter.UpstreamBarX,c.MuFilter.UpstreamBarY,c.MuFilter.UpstreamBarZ = 82.525*u.cm, 6.0*u.cm, 1.0*u.cm
215 c.MuFilter.UpstreamBarGap = 0.1*u.mm
216
217 #downstream bars configuration
218 c.MuFilter.NDownstreamBars = 60 #n.d.r. both for x and y in this case
219 c.MuFilter.DownstreamBarX,c.MuFilter.DownstreamBarY,c.MuFilter.DownstreamBarZ = 82.525*u.cm, 1*u.cm, 1*u.cm
220 c.MuFilter.DownstreamBarGap = 0.1*u.mm
221 c.MuFilter.DownstreamBarX_ver,c.MuFilter.DownstreamBarY_ver,c.MuFilter.DownstreamBarZ_ver = 1*u.cm, 63.525*u.cm, 1*u.cm
222
223 # DS and US support box, inner, Z pointing upward
224 c.MuFilter.SupportBoxD = 0.5*u.mm # empty space between bars and box
225 c.MuFilter.SupportBoxW = 2*u.mm
226 c.MuFilter.DSBoxX1 = c.MuFilter.DSHLocX - c.MuFilter.SupportBoxD
227 c.MuFilter.DSBoxX2 = c.MuFilter.DSHLocX + c.MuFilter.DownstreamBarX + c.MuFilter.SupportBoxD
228 c.MuFilter.DSBoxZ1 = c.MuFilter.DSHLocZ - c.MuFilter.DownstreamBarY/2 - c.MuFilter.SupportBoxD
229 c.MuFilter.DSBoxZ2 = c.MuFilter.DSVLocZ + c.MuFilter.SupportBoxD
230 c.MuFilter.DSBoxY1 = c.MuFilter.DSHLocY - c.MuFilter.DownstreamBarZ/2 - c.MuFilter.SupportBoxD
231 c.MuFilter.DSBoxY2 = c.MuFilter.DSVLocY + c.MuFilter.DownstreamBarZ/2 + c.MuFilter.SupportBoxD
232
233 c.MuFilter.USBoxY1 = c.MuFilter.DSHLocY - c.MuFilter.DownstreamBarZ/2 - c.MuFilter.SupportBoxD
234 c.MuFilter.USBoxY2 = c.MuFilter.DSHLocY + c.MuFilter.DownstreamBarZ/2 + c.MuFilter.SupportBoxD
235
236 # VETO support box
237 c.MuFilter.SupportBoxVW = 4*u.mm
238 c.MuFilter.VETOBoxX1 = c.MuFilter.VETOLocX - c.MuFilter.SupportBoxD
239 c.MuFilter.VETOBoxX2 = c.MuFilter.VETOLocX + c.MuFilter.VetoBarX + c.MuFilter.SupportBoxD
240 c.MuFilter.VETOBoxZ1 = c.MuFilter.VETOLocZ - c.MuFilter.VetoBarY/2 - c.MuFilter.SupportBoxD
241 c.MuFilter.VETOBoxZ2 = c.MuFilter.VETOLocZ + (c.MuFilter.NVetoBars-1)*(c.MuFilter.VetoBarY+c.MuFilter.VetoBarGap) + c.MuFilter.VetoBarY/2 + c.MuFilter.SupportBoxD
242 c.MuFilter.VETOBoxY1 = c.MuFilter.VETOLocY - c.MuFilter.VetoBarZ/2 - c.MuFilter.SupportBoxD
243 c.MuFilter.VETOBoxY2 = c.MuFilter.VETOLocY + c.MuFilter.VetoBarZ/2 + c.MuFilter.SupportBoxD
244
245 # VETO/US/DS plane alignment
246 c.MuFilter.Veto1ShiftY = 0.11 * u.cm
247 c.MuFilter.Veto2ShiftY = -0.04 * u.cm
248 c.MuFilter.US1ShiftY = 0.10 * u.cm
249 c.MuFilter.US2ShiftY = 0.26 * u.cm
250 c.MuFilter.US3ShiftY = 0.24 * u.cm
251 c.MuFilter.US4ShiftY = 0.31 * u.cm
252 c.MuFilter.US5ShiftY = 0.34 * u.cm
253 c.MuFilter.DS1ShiftY = 0.43 * u.cm
254 c.MuFilter.DS1ShiftX = 1.13 * u.cm
255 c.MuFilter.DS2ShiftY = 0.53 * u.cm
256 c.MuFilter.DS2ShiftX = 1.31 * u.cm
257 c.MuFilter.DS3ShiftY = 0.61 * u.cm
258 c.MuFilter.DS3ShiftX = 1.35 * u.cm
259 c.MuFilter.DS4ShiftX = 1.39 * u.cm
260
261 #digitization parameters
262 c.MuFilter.DsAttenuationLength = 350 * u.cm # values between 300 cm and 400cm observed for H6 testbeam
263 c.MuFilter.DsTAttenuationLength = 700 * u.cm # top readout with mirror on bottom
264 c.MuFilter.VandUpAttenuationLength = 999 * u.cm # no significante attenuation observed for H6 testbeam
265 c.MuFilter.VandUpSiPMcalibrationL = 25.*1000. # 1.65 MeV = 41 qcd
266 c.MuFilter.VandUpSiPMcalibrationS = 25.*1000.
267 c.MuFilter.DsSiPMcalibration = 25.*1000.
268 c.MuFilter.timeResol = 150.*u.picosecond
269 c.MuFilter.VandUpPropSpeed = 12.5*u.cm/u.nanosecond
270 c.MuFilter.DsPropSpeed = 14.3*u.cm/u.nanosecond
271
272 c.Floor = AttrDict(z=48000.*u.cm) # to place tunnel in SND_@LHC coordinate system
273 c.Floor.DX = 1.0*u.cm
274 c.Floor.DY = -4.5*u.cm # subtract 4.5cm to avoid overlaps
275 c.Floor.DZ = 0.
276
277 #COLDBOX configuration
278 c.Floor.Acrylic_width = 5.0*u.cm
279 c.Floor.BPoly_width = 4.0*u.cm
280 c.Floor.CBFrontWall_xdim = 219.*u.cm
281 c.Floor.CBFrontWall_ydim = 170.72*u.cm-c.Floor.Acrylic_width
282 c.Floor.CBLatWall_zdim = 176.0*u.cm
283 c.Floor.CBTiny_zdim = 17.0*u.cm
284 c.Floor.CBExtra_zdim = 41.0*u.cm
285 c.Floor.CBExtra_xdim = 67.5 *u.cm
286 c.Floor.SlopedWall_zproj = 110.0*u.cm
287 c.Floor.MFeBlockX = c.MuFilter.FeX
288 c.Floor.MFeBlockY = c.MuFilter.FeY
289 c.Floor.MFeBlockZ = c.MuFilter.FeZ
290
291# for H4 / H6 / H8 / HX testbeam and commissioning
292 H4 = True
293 if H4:
294 c.Floor.z = 0 # no tunnel, no slope
295
296 # no Veto
297 c.MuFilter.NVetoPlanes = 0
298 # no US
299 c.MuFilter.NUpstreamPlanes = 0
300
301 # only 1 DS station
302 c.MuFilter.NDownstreamPlanes = 1
303 # remove most downstream iron block and its bottom part
304 c.MuFilter.FeEndX,c.MuFilter.FeEndY,c.MuFilter.FeEndZ = 0*u.cm, 0*u.cm, 0*u.cm
305 c.MuFilter.FeBotX,c.MuFilter.FeBotY,c.MuFilter.FeBotZ = 0*u.cm, 0*u.cm, 0*u.cm
306 # remove the iron block before the DS1
307 c.MuFilter.FeX,c.MuFilter.FeY,c.MuFilter.FeZ = 0*u.cm, 0*u.cm, 0*u.cm
308
309 # removing the emulsion+W target for the second part of the TB,
310 # where we have Fe blocks btw SciFi stations
311 c.EmulsionDet.wall = 0
312 if with_tungsten:
313 c.EmulsionDet = AttrDict(z=0*u.cm)
314 c.EmulsionDet.PassiveOption = 1 #0 makes emulsion volumes active, 1 makes all emulsion volumes passive
315 c.EmulsionDet.row = 2
316 c.EmulsionDet.col = 2
317 c.EmulsionDet.wall= 2
318 c.EmulsionDet.target = 0 #number of neutrino target volumes
319
320 c.EmulsionDet.xdim = 42.2 *u.cm #external wall dimensions
321 c.EmulsionDet.ydim = 42.2 *u.cm
322 c.EmulsionDet.TotalWallZDim = 10.0 *u.cm #wall dimension along z, including border
323 c.EmulsionDet.WallXDim = 38.6 *u.cm #internal wall dimensions
324 c.EmulsionDet.WallYDim = 38.6 *u.cm
325 c.EmulsionDet.WallZDim = 8.15 *u.cm
326 c.EmulsionDet.WallZBorder_offset = 4.75 * u.mm
327 c.EmulsionDet.TTz = 4.0*u.cm # thickness of target station!
328 c.EmulsionDet.zdim = c.EmulsionDet.wall* c.EmulsionDet.TotalWallZDim + c.EmulsionDet.wall*c.EmulsionDet.TTz
329
330 c.EmulsionDet.n_plates = 59 # the most downstream 28 layers are W+plastic
331 c.EmulsionDet.n_tungsten_plates_tb24 = 28
332 c.EmulsionDet.EmX = 19.2 * u.cm
333 c.EmulsionDet.EmY = 19.2 * u.cm
334 c.EmulsionDet.EmTh = 0.0 * u.cm # no emulsion
335 c.EmulsionDet.PBTh = 0.0175 * u.cm + 2*0.0070 * u.cm # this is the plastic base tickness -same value as the full em+plastic basein TI18
336 c.EmulsionDet.PassiveTh = 0.1 * u.cm # this is the tungsten
337 c.EmulsionDet.EPlW = c.EmulsionDet.PBTh
338 c.EmulsionDet.AllPW = c.EmulsionDet.PassiveTh + c.EmulsionDet.EPlW # this is used to add layers of W+emulsion
339
340 c.EmulsionDet.BrPackZ = 0.*u.cm
341 c.EmulsionDet.BrPackX = 2*0.05*u.cm
342 c.EmulsionDet.BrPackY = 2*0.05*u.cm
343
344 c.EmulsionDet.BrX = c.EmulsionDet.EmX + c.EmulsionDet.BrPackX
345 c.EmulsionDet.BrY = c.EmulsionDet.EmY + c.EmulsionDet.BrPackY
346 c.EmulsionDet.BrZ = c.EmulsionDet.n_plates * c.EmulsionDet.AllPW + c.EmulsionDet.EPlW + c.EmulsionDet.BrPackZ
347
348 # do we use that? #FIXME
349 c.EmulsionDet.ShiftX = 0.#-8.0*u.cm - c.EmulsionDet.xdim/2.
350 c.EmulsionDet.ShiftY = 0.#15.5*u.cm + c.EmulsionDet.ydim/2.
351 # do we use that? #FIXME
352 c.EmulsionDet.startpos = 0.#-25.4750 * u.cm + c.EmulsionDet.z
353 c.EmulsionDet.zC = 0.# c.EmulsionDet.startpos + c.EmulsionDet.zdim/2.
354
355 # survey points in survey coordinate system!
356 c.EmulsionDet.Xpos0,c.EmulsionDet.Ypos0,c.EmulsionDet.Zpos0 = 318.6*u.mm,3218.8*u.mm,380.2*u.mm
357 c.EmulsionDet.Xpos1,c.EmulsionDet.Ypos1,c.EmulsionDet.Zpos1 = 318.6*u.mm,3383.8*u.mm,380.2*u.mm
358
359
360 # set SciFi modules
361 c.Scifi.xdim = 13.0*u.cm #sensitive only
362 c.Scifi.ydim = 13.0*u.cm
363 c.Scifi.zdim = 4.0*u.cm # z thickness of Scifi station, must match c.EmulsionDet.TTz
364 c.Scifi.nmats = 1
365 c.Scifi.nscifi = 4
366 c.Scifi.nfibers_z = 7
367 c.Scifi.scifimat_z = 0.16*u.cm # thickness of a mat of 7 fiber layers
368 c.Scifi.scifimat_length = c.Scifi.ydim
369 c.Scifi.fiber_length = c.Scifi.scifimat_length
370 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
371 c.Scifi.tedlar_to_plane = 5.42*u.mm # an air gap btw protective tedlar sheet and a sensitive plane
372 #c.Scifi.tedlar_z = 50*u.um # not included in the sw detector model
373 # offset between the edge of a baby module frame and the upstream tedlar inside it
374 c.Scifi.frame_offset = 0.8*u.cm
375 # offsets between the upstream tedlar sheet and the upstream passive block
376 c.Scifi.station_offset1 = c.Scifi.frame_offset+1.4*u.cm
377 c.Scifi.station_offset2 = c.Scifi.frame_offset+1.4*u.cm
378 c.Scifi.station_offset3 = c.Scifi.frame_offset+1.4*u.cm
379 c.Scifi.channelTimeAlignment = 1
380
381 if not with_tungsten:
382 # add variable size blocks upstream of SciFi stations
383 # complete removal of a wall goes with commenting the respective line below
384 #c.Scifi.FeTargetX3, c.Scifi.FeTargetY3, c.Scifi.FeTargetZ3 = 30.*u.cm, 10.*u.cm, 30.*u.cm
385 c.Scifi.FeTargetX4, c.Scifi.FeTargetY4, c.Scifi.FeTargetZ4 = 30.*u.cm, 10.*u.cm, 30.*u.cm
386 # flag showing whether the passive material (FeBlocks) are to be centerred around the Scifi
387 c.Scifi.PassiveBlocknotCenterred = 1
388
389# absolute edge point positions in survey coordinate system (survey is 'by eye' for now)
390 c.Scifi.Xpos0,c.Scifi.Ypos0,c.Scifi.Zpos0 = 156.6*u.mm,3186.3*u.mm,218.2*u.mm
391 c.Scifi.Xpos1,c.Scifi.Ypos1,c.Scifi.Zpos1 = 156.6*u.mm,3351.3*u.mm,221.7*u.mm # 16.5cm step
392 c.Scifi.Xpos2,c.Scifi.Ypos2,c.Scifi.Zpos2 = 156.6*u.mm,3516.3*u.mm,225.1*u.mm # 16.5cm step
393 c.Scifi.Xpos3,c.Scifi.Ypos3,c.Scifi.Zpos3 = 156.6*u.mm,3676.3*u.mm,228.5*u.mm # 16cm step
394
395 # DS1 - since there are no US planes, we take the item labelled Muon1!
396 # DS spatial alignment is included in the values
397 c.MuFilter.Muon1Dx,c.MuFilter.Muon1Dy,c.MuFilter.Muon1Dz = 318.6*u.mm, 4066.3*u.mm, 230.2*u.mm
398
399# from Scifi track alignment
400 if tb_2024_mc :
401# spatial
402 c.Scifi.LocM100, c.Scifi.LocM110 = 0.000*u.um, 0.000*u.um
403 c.Scifi.LocM200, c.Scifi.LocM210 = 0.000*u.um, 0.000*u.um
404 c.Scifi.LocM300, c.Scifi.LocM310 = 0.000*u.um, 0.000*u.um
405 c.Scifi.LocM400, c.Scifi.LocM410 = 0.000*u.um, 0.000*u.um
406 c.Scifi.RotPhiS10,c.Scifi.RotPsiS10,c.Scifi.RotThetaS10 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
407 c.Scifi.RotPhiS11,c.Scifi.RotPsiS11,c.Scifi.RotThetaS11 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
408 c.Scifi.RotPhiS20,c.Scifi.RotPsiS20,c.Scifi.RotThetaS20 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
409 c.Scifi.RotPhiS21,c.Scifi.RotPsiS21,c.Scifi.RotThetaS21 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
410 c.Scifi.RotPhiS30,c.Scifi.RotPsiS30,c.Scifi.RotThetaS30 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
411 c.Scifi.RotPhiS31,c.Scifi.RotPsiS31,c.Scifi.RotThetaS31 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
412 c.Scifi.RotPhiS40,c.Scifi.RotPsiS40,c.Scifi.RotThetaS40 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
413 c.Scifi.RotPhiS41,c.Scifi.RotPsiS41,c.Scifi.RotThetaS41 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
414# Time alignment Scifi, T0 = station 0, mat 0
415 c.Scifi.station1t,c.Scifi.station1H0t,c.Scifi.station1V0t = 0.000*u.ns, 0.000*u.ns, 0.000*u.ns
416 c.Scifi.station2t,c.Scifi.station2H0t,c.Scifi.station2V0t = 0.000*u.ns, 0.000*u.ns, 0.000*u.ns
417 c.Scifi.station3t,c.Scifi.station3H0t,c.Scifi.station3V0t = 0.000*u.ns, 0.000*u.ns, 0.000*u.ns
418 c.Scifi.station4t,c.Scifi.station4H0t,c.Scifi.station4V0t = 0.000*u.ns, 0.000*u.ns, 0.000*u.ns
419 else :
420# spatial
421 c.Scifi.LocM100, c.Scifi.LocM110 = 2800.00*u.um, 420.01*u.um
422 c.Scifi.LocM200, c.Scifi.LocM210 = 554.35*u.um, 374.96*u.um
423 c.Scifi.LocM300, c.Scifi.LocM310 = -614.41**u.um, 56.48*u.um
424 c.Scifi.LocM400, c.Scifi.LocM410 = -154.11*u.um, 472.07*u.um
425 c.Scifi.RotPhiS10,c.Scifi.RotPsiS10,c.Scifi.RotThetaS10 = 0.00*u.mrad, 0.00*u.mrad, 0.00*u.mrad
426 c.Scifi.RotPhiS11,c.Scifi.RotPsiS11,c.Scifi.RotThetaS11 = 0.00*u.mrad, 1.65*u.mrad, 0.00*u.mrad
427 c.Scifi.RotPhiS20,c.Scifi.RotPsiS20,c.Scifi.RotThetaS20 = 0.00*u.mrad, 2.23*u.mrad, 0.00*u.mrad
428 c.Scifi.RotPhiS21,c.Scifi.RotPsiS21,c.Scifi.RotThetaS21 = 0.00*u.mrad, 0.68*u.mrad, 0.00*u.mrad
429 c.Scifi.RotPhiS30,c.Scifi.RotPsiS30,c.Scifi.RotThetaS30 = 0.00*u.mrad, 1.06*u.mrad, 0.00*u.mrad
430 c.Scifi.RotPhiS31,c.Scifi.RotPsiS31,c.Scifi.RotThetaS31 = 0.00*u.mrad, 0.82*u.mrad, 0.00*u.mrad
431 c.Scifi.RotPhiS40,c.Scifi.RotPsiS40,c.Scifi.RotThetaS40 = 0.00*u.mrad, 1.61*u.mrad, 0.00*u.mrad
432 c.Scifi.RotPhiS41,c.Scifi.RotPsiS41,c.Scifi.RotThetaS41 = 0.00*u.mrad, 2.22*u.mrad, 0.00*u.mrad
433# Time alignment Scifi, T0 = station 0, mat 0
434 c.Scifi.station1t,c.Scifi.station1H0t,c.Scifi.station1V0t = 0.000*u.ns, 0.000*u.ns, -1.714*u.ns
435 c.Scifi.station2t,c.Scifi.station2H0t,c.Scifi.station2V0t = -0.510*u.ns, 0.000*u.ns, 0.003*u.ns
436 c.Scifi.station3t,c.Scifi.station3H0t,c.Scifi.station3V0t = 1.447*u.ns, 0.000*u.ns, 0.520*u.ns
437 c.Scifi.station4t,c.Scifi.station4H0t,c.Scifi.station4V0t = -0.014*u.ns, 0.000*u.ns, 2.357*u.ns
438# Additional alignment for station 2H(one mat) perhaps needed due to calibration issue
439 c.Scifi.station2000XXXt = [-1.956*u.ns,-2.231*u.ns,-1.250*u.ns,-1.474*u.ns,-2.215*u.ns,-2.027*u.ns,-1.713*u.ns,-1.899*u.ns,
440 -1.875*u.ns,-1.848*u.ns,-1.862*u.ns,-2.184*u.ns,-1.855*u.ns,-2.066*u.ns,-1.252*u.ns,-2.268*u.ns,
441 -2.370*u.ns,-1.414*u.ns,-1.577*u.ns,-1.909*u.ns,-1.952*u.ns,-1.491*u.ns,-1.269*u.ns,-1.741*u.ns,
442 -1.823*u.ns,-2.120*u.ns,-2.224*u.ns,-1.758*u.ns,-1.715*u.ns,-2.037*u.ns,-2.219*u.ns,-2.667*u.ns,
443 -2.269*u.ns,-1.526*u.ns,-1.840*u.ns,-1.258*u.ns,-2.459*u.ns,-2.631*u.ns,-1.335*u.ns,-1.448*u.ns,
444 -1.364*u.ns,-1.589*u.ns,-1.647*u.ns,-1.677*u.ns,-1.508*u.ns,-1.548*u.ns,-1.289*u.ns,-1.249*u.ns,
445 -1.388*u.ns,-1.654*u.ns,-1.945*u.ns,-2.371*u.ns,-2.308*u.ns,-1.310*u.ns,-1.497*u.ns,-1.928*u.ns,
446 -2.160*u.ns,-1.470*u.ns,-1.452*u.ns,-2.263*u.ns,-2.291*u.ns,-1.557*u.ns,-1.396*u.ns,-1.003*u.ns,
447 -0.735*u.ns,-1.354*u.ns,-1.508*u.ns,-0.978*u.ns,-1.111*u.ns,-0.750*u.ns,-1.168*u.ns,-1.891*u.ns,
448 -2.219*u.ns,-0.647*u.ns,-0.557*u.ns,-0.518*u.ns,-0.937*u.ns,-1.397*u.ns,-1.702*u.ns,-1.624*u.ns,
449 -1.055*u.ns,-1.457*u.ns,-1.114*u.ns,-1.283*u.ns,-1.562*u.ns,-1.997*u.ns,-1.717*u.ns,-0.753*u.ns,
450 -0.628*u.ns,-1.031*u.ns,-1.405*u.ns,-0.910*u.ns,-0.723*u.ns,-0.214*u.ns,-0.762*u.ns,-1.101*u.ns,
451 -0.685*u.ns,-0.799*u.ns,-1.127*u.ns,-1.164*u.ns,-1.375*u.ns,-1.535*u.ns,-1.577*u.ns,-1.604*u.ns,
452 -1.492*u.ns,-1.137*u.ns,-1.454*u.ns,-1.812*u.ns,-1.558*u.ns,-1.551*u.ns,-1.475*u.ns,-1.361*u.ns,
453 -1.071*u.ns,-1.525*u.ns,-1.471*u.ns,-0.709*u.ns,-1.032*u.ns,-0.837*u.ns,-0.945*u.ns,-0.657*u.ns,
454 -1.457*u.ns,-1.759*u.ns,-1.442*u.ns,-1.124*u.ns,-1.246*u.ns,-1.320*u.ns,-1.457*u.ns,-0.479*u.ns]
455 c.Scifi.station2001XXXt = [-1.387*u.ns,-1.486*u.ns,-1.424*u.ns,-0.859*u.ns,-0.838*u.ns,-1.253*u.ns,-1.724*u.ns,-1.970*u.ns,
456 -1.731*u.ns,-0.977*u.ns,-0.948*u.ns,-1.277*u.ns,-1.225*u.ns,-1.032*u.ns,-1.747*u.ns,-1.799*u.ns,
457 -0.858*u.ns,-1.966*u.ns,-2.336*u.ns,-1.940*u.ns,-1.809*u.ns,-1.265*u.ns,-1.688*u.ns,-1.636*u.ns,
458 -1.602*u.ns,-1.785*u.ns,-1.958*u.ns,-2.064*u.ns,-2.039*u.ns,-1.870*u.ns,-2.055*u.ns,-1.428*u.ns,
459 -1.206*u.ns,-1.516*u.ns,-1.528*u.ns,-1.546*u.ns,-1.537*u.ns,-1.860*u.ns,-1.399*u.ns,-1.403*u.ns,
460 -1.581*u.ns,-1.223*u.ns,-1.864*u.ns,-2.190*u.ns,-1.152*u.ns,-1.617*u.ns,-1.435*u.ns,-0.957*u.ns,
461 -1.370*u.ns,-1.562*u.ns,-1.780*u.ns,-1.456*u.ns,-1.493*u.ns,-1.875*u.ns,-1.713*u.ns,-1.054*u.ns,
462 -1.241*u.ns,-1.240*u.ns,-0.540*u.ns,-0.872*u.ns,-1.132*u.ns,-1.139*u.ns,-1.130*u.ns,-0.897*u.ns,
463 -1.188*u.ns,-1.626*u.ns,-1.564*u.ns,-1.865*u.ns,-2.079*u.ns,-1.943*u.ns,-0.764*u.ns,-0.936*u.ns,
464 -1.313*u.ns,-1.005*u.ns,-1.335*u.ns,-1.594*u.ns,-1.639*u.ns,-1.580*u.ns,-1.701*u.ns,-1.572*u.ns,
465 -1.927*u.ns,-1.390*u.ns,-2.407*u.ns,-2.865*u.ns,-1.308*u.ns,-1.489*u.ns,-1.921*u.ns,-2.161*u.ns,
466 -1.243*u.ns,-1.563*u.ns,-1.855*u.ns,-2.163*u.ns,-1.930*u.ns,-2.251*u.ns,-1.584*u.ns,-1.448*u.ns,
467 -2.016*u.ns,-1.721*u.ns,-1.447*u.ns,-1.220*u.ns,-1.139*u.ns,-1.932*u.ns,-1.753*u.ns,-1.781*u.ns,
468 -1.751*u.ns,-1.583*u.ns,-1.193*u.ns,-1.516*u.ns,-1.061*u.ns,-1.206*u.ns,-1.236*u.ns,-1.638*u.ns,
469 -1.777*u.ns,-1.362*u.ns,-1.385*u.ns,-1.465*u.ns,-1.185*u.ns,-0.964*u.ns,-0.869*u.ns,-1.308*u.ns,
470 -1.806*u.ns,-1.900*u.ns,-1.898*u.ns,-1.365*u.ns,-1.008*u.ns,-1.479*u.ns,-1.907*u.ns,-1.474*u.ns]
471 c.Scifi.station2002XXXt = [0.627*u.ns,0.230*u.ns,0.202*u.ns,-0.234*u.ns,-0.176*u.ns,0.004*u.ns,0.302*u.ns,0.697*u.ns,
472 0.477*u.ns,-0.210*u.ns,-0.406*u.ns,-0.086*u.ns,-0.129*u.ns,-0.625*u.ns,0.104*u.ns,-0.002*u.ns,
473 -0.683*u.ns,-0.065*u.ns,0.746*u.ns,0.408*u.ns,0.226*u.ns,-0.507*u.ns,0.324*u.ns,0.330*u.ns,
474 0.244*u.ns,0.155*u.ns,0.552*u.ns,0.442*u.ns,0.462*u.ns,0.113*u.ns,0.477*u.ns,0.086*u.ns,
475 -0.384*u.ns,0.077*u.ns,-0.004*u.ns,0.095*u.ns,-0.264*u.ns,0.418*u.ns,0.051*u.ns,-0.134*u.ns,
476 0.239*u.ns,-0.055*u.ns,-0.113*u.ns,0.542*u.ns,-0.251*u.ns,-0.027*u.ns,0.067*u.ns,-0.707*u.ns,
477 0.018*u.ns,-0.017*u.ns,0.378*u.ns,0.002*u.ns,-0.178*u.ns,0.454*u.ns,0.334*u.ns,-0.379*u.ns,
478 -0.552*u.ns,-0.050*u.ns,-0.811*u.ns,-0.457*u.ns,-0.107*u.ns,-0.078*u.ns,-0.028*u.ns,0.149*u.ns,
479 0.348*u.ns,0.301*u.ns,0.234*u.ns,0.373*u.ns,0.573*u.ns,0.047*u.ns,-0.689*u.ns,-0.639*u.ns,
480 -0.169*u.ns,-0.484*u.ns,-0.015*u.ns,0.096*u.ns,0.154*u.ns,0.057*u.ns,0.078*u.ns,-0.102*u.ns,
481 0.180*u.ns,-0.130*u.ns,0.378*u.ns,0.407*u.ns,-0.073*u.ns,0.028*u.ns,-0.049*u.ns,0.214*u.ns,
482 -0.379*u.ns,0.009*u.ns,0.044*u.ns,0.302*u.ns,-0.012*u.ns,0.533*u.ns,-0.158*u.ns,-0.383*u.ns,
483 0.548*u.ns,0.160*u.ns,-0.058*u.ns,-0.617*u.ns,-0.862*u.ns,0.477*u.ns,0.282*u.ns,0.349*u.ns,
484 0.364*u.ns,-0.080*u.ns,-0.473*u.ns,-0.022*u.ns,-0.444*u.ns,-0.181*u.ns,-0.260*u.ns,0.240*u.ns,
485 0.452*u.ns,0.114*u.ns,0.192*u.ns,0.168*u.ns,-0.085*u.ns,-0.228*u.ns,-0.370*u.ns,-0.341*u.ns,
486 0.516*u.ns,0.463*u.ns,0.323*u.ns,-0.061*u.ns,-0.245*u.ns,-0.205*u.ns,0.531*u.ns,0.772*u.ns]
487 c.Scifi.station2003XXXt = [0.800*u.ns,0.586*u.ns,0.113*u.ns,0.188*u.ns,0.760*u.ns,0.724*u.ns,0.409*u.ns,0.613*u.ns,
488 0.644*u.ns,0.538*u.ns,0.316*u.ns,0.696*u.ns,0.334*u.ns,0.490*u.ns,-0.216*u.ns,-0.351*u.ns,
489 0.776*u.ns,0.123*u.ns,0.183*u.ns,0.488*u.ns,0.406*u.ns,0.024*u.ns,-0.253*u.ns,0.376*u.ns,
490 0.319*u.ns,0.505*u.ns,0.764*u.ns,0.352*u.ns,0.220*u.ns,0.437*u.ns,0.404*u.ns,0.831*u.ns,
491 0.377*u.ns,-0.279*u.ns,-0.166*u.ns,-0.778*u.ns,0.675*u.ns,0.466*u.ns,-0.063*u.ns,0.030*u.ns,
492 0.040*u.ns,0.165*u.ns,0.257*u.ns,0.241*u.ns,0.219*u.ns,0.165*u.ns,-0.029*u.ns,-0.041*u.ns,
493 0.163*u.ns,0.187*u.ns,0.488*u.ns,0.878*u.ns,0.561*u.ns,0.010*u.ns,0.133*u.ns,0.344*u.ns,
494 0.572*u.ns,-0.019*u.ns,-0.081*u.ns,0.634*u.ns,0.706*u.ns,0.339*u.ns,0.252*u.ns,0.203*u.ns,
495 -0.029*u.ns,0.007*u.ns,-0.018*u.ns,-0.370*u.ns,-0.266*u.ns,-0.621*u.ns,-0.403*u.ns,0.073*u.ns,
496 0.009*u.ns,-0.691*u.ns,-0.736*u.ns,-0.940*u.ns,-0.656*u.ns,-0.158*u.ns,0.097*u.ns,-0.492*u.ns,
497 -0.952*u.ns,-0.185*u.ns,-0.356*u.ns,-0.156*u.ns,-0.001*u.ns,0.237*u.ns,-0.179*u.ns,-0.736*u.ns,
498 -0.879*u.ns,-0.930*u.ns,-0.139*u.ns,-0.575*u.ns,-0.959*u.ns,-1.411*u.ns,-0.991*u.ns,-0.358*u.ns,
499 -0.799*u.ns,-0.879*u.ns,-0.304*u.ns,-0.227*u.ns,-0.089*u.ns,0.099*u.ns,0.153*u.ns,0.210*u.ns,
500 0.156*u.ns,-0.186*u.ns,-0.115*u.ns,0.372*u.ns,0.183*u.ns,0.129*u.ns,0.155*u.ns,-0.129*u.ns,
501 -0.415*u.ns,0.010*u.ns,-0.138*u.ns,-0.592*u.ns,-0.460*u.ns,-0.613*u.ns,-0.628*u.ns,-0.992*u.ns,
502 0.004*u.ns,0.323*u.ns,-0.042*u.ns,-0.180*u.ns,-0.041*u.ns,-0.095*u.ns,-0.092*u.ns,-0.173*u.ns]