b661165c |
1 | // Author: Anders Vestbo <mailto:vestbo@fi.uib.no>, Uli Frankenfeld <mailto:franken@fi.uib.no> |
2 | //*-- Copyright © ASV |
108615fc |
3 | |
4 | #include "AliL3Logging.h" |
5 | #include "AliL3Transform.h" |
6 | //#include <TFile.h> |
7 | #include <math.h> |
b661165c |
8 | |
9 | //_____________________________________________________________ |
108615fc |
10 | // AliL3Transform |
11 | // |
12 | // Transformation class for ALICE TPC. |
13 | |
14 | ClassImp(AliL3Transform); |
15 | |
16 | |
17 | AliL3Transform::AliL3Transform(){ |
18 | //constructor |
19 | Init(); |
20 | } |
21 | |
22 | |
23 | AliL3Transform::~AliL3Transform(){ |
24 | } |
25 | |
a4982538 |
26 | void AliL3Transform::Init(){ |
27 | //sector: |
28 | fNTimeBins = 446; |
29 | fNRowLow = 64; |
30 | fNRowUp = 112; |
31 | fNSectorLow = 36; |
32 | fNSectorUp = 36; |
33 | fNSector = 72; |
34 | fPadPitchWidthLow = 0.400000; |
35 | fPadPitchWidthUp = 0.600000; |
36 | fZWidth = 0.56599998474121093750; |
37 | fZSigma = 0.22880849748219134199; |
38 | |
39 | //slices: |
40 | fNSlice = 36; |
41 | fNRow = 176; |
42 | fPi = 3.141592653589793; |
40f2887b |
43 | /* |
a4982538 |
44 | for(Int_t i=0;i<36;i++){ |
45 | fCos[i] = cos(fPi/9*i); |
46 | fSin[i] = sin(fPi/9*i); |
47 | } |
40f2887b |
48 | */ |
49 | for(Int_t i=0;i<36;i++){ |
50 | fCos[i] = cos( (2*fPi/18) * (i+0.5) ); |
51 | fSin[i] = sin( (2*fPi/18) * (i+0.5) ); |
52 | } |
a4982538 |
53 | |
54 | fX[0] = 84.570007324218750; |
55 | fX[1] = 85.320007324218750; |
56 | fX[2] = 86.070007324218750; |
57 | fX[3] = 86.820007324218750; |
58 | fX[4] = 87.570007324218750; |
59 | fX[5] = 88.320007324218750; |
60 | fX[6] = 89.070007324218750; |
61 | fX[7] = 89.820007324218750; |
62 | fX[8] = 90.570007324218750; |
63 | fX[9] = 91.320007324218750; |
64 | fX[10] = 92.070007324218750; |
65 | fX[11] = 92.820007324218750; |
66 | fX[12] = 93.570007324218750; |
67 | fX[13] = 94.320007324218750; |
68 | fX[14] = 95.070007324218750; |
69 | fX[15] = 95.820007324218750; |
70 | fX[16] = 96.570007324218750; |
71 | fX[17] = 97.320007324218750; |
72 | fX[18] = 98.070007324218750; |
73 | fX[19] = 98.820007324218750; |
74 | fX[20] = 99.570007324218750; |
75 | fX[21] = 100.320007324218750; |
76 | fX[22] = 101.070007324218750; |
77 | fX[23] = 101.820007324218750; |
78 | fX[24] = 102.570007324218750; |
79 | fX[25] = 103.320007324218750; |
80 | fX[26] = 104.070007324218750; |
81 | fX[27] = 104.820007324218750; |
82 | fX[28] = 105.570007324218750; |
83 | fX[29] = 106.320007324218750; |
84 | fX[30] = 107.070007324218750; |
85 | fX[31] = 107.820007324218750; |
86 | fX[32] = 108.570007324218750; |
87 | fX[33] = 109.320007324218750; |
88 | fX[34] = 110.070007324218750; |
89 | fX[35] = 110.820007324218750; |
90 | fX[36] = 111.570007324218750; |
91 | fX[37] = 112.320007324218750; |
92 | fX[38] = 113.070007324218750; |
93 | fX[39] = 113.820007324218750; |
94 | fX[40] = 114.570007324218750; |
95 | fX[41] = 115.320007324218750; |
96 | fX[42] = 116.070007324218750; |
97 | fX[43] = 116.820007324218750; |
98 | fX[44] = 117.570007324218750; |
99 | fX[45] = 118.320007324218750; |
100 | fX[46] = 119.070007324218750; |
101 | fX[47] = 119.820007324218750; |
102 | fX[48] = 120.570007324218750; |
103 | fX[49] = 121.320007324218750; |
104 | fX[50] = 122.070007324218750; |
105 | fX[51] = 122.820007324218750; |
106 | fX[52] = 123.570007324218750; |
107 | fX[53] = 124.320007324218750; |
108 | fX[54] = 125.070007324218750; |
109 | fX[55] = 125.820007324218750; |
110 | fX[56] = 126.570007324218750; |
111 | fX[57] = 127.320007324218750; |
112 | fX[58] = 128.070007324218750; |
113 | fX[59] = 128.820007324218750; |
114 | fX[60] = 129.570007324218750; |
115 | fX[61] = 130.320007324218750; |
116 | fX[62] = 131.070007324218750; |
117 | fX[63] = 131.820007324218750; |
118 | fX[64] = 135.054992675781250; |
119 | fX[65] = 136.054992675781250; |
120 | fX[66] = 137.054992675781250; |
121 | fX[67] = 138.054992675781250; |
122 | fX[68] = 139.054992675781250; |
123 | fX[69] = 140.054992675781250; |
124 | fX[70] = 141.054992675781250; |
125 | fX[71] = 142.054992675781250; |
126 | fX[72] = 143.054992675781250; |
127 | fX[73] = 144.054992675781250; |
128 | fX[74] = 145.054992675781250; |
129 | fX[75] = 146.054992675781250; |
130 | fX[76] = 147.054992675781250; |
131 | fX[77] = 148.054992675781250; |
132 | fX[78] = 149.054992675781250; |
133 | fX[79] = 150.054992675781250; |
134 | fX[80] = 151.054992675781250; |
135 | fX[81] = 152.054992675781250; |
136 | fX[82] = 153.054992675781250; |
137 | fX[83] = 154.054992675781250; |
138 | fX[84] = 155.054992675781250; |
139 | fX[85] = 156.054992675781250; |
140 | fX[86] = 157.054992675781250; |
141 | fX[87] = 158.054992675781250; |
142 | fX[88] = 159.054992675781250; |
143 | fX[89] = 160.054992675781250; |
144 | fX[90] = 161.054992675781250; |
145 | fX[91] = 162.054992675781250; |
146 | fX[92] = 163.054992675781250; |
147 | fX[93] = 164.054992675781250; |
148 | fX[94] = 165.054992675781250; |
149 | fX[95] = 166.054992675781250; |
150 | fX[96] = 167.054992675781250; |
151 | fX[97] = 168.054992675781250; |
152 | fX[98] = 169.054992675781250; |
153 | fX[99] = 170.054992675781250; |
154 | fX[100] = 171.054992675781250; |
155 | fX[101] = 172.054992675781250; |
156 | fX[102] = 173.054992675781250; |
157 | fX[103] = 174.054992675781250; |
158 | fX[104] = 175.054992675781250; |
159 | fX[105] = 176.054992675781250; |
160 | fX[106] = 177.054992675781250; |
161 | fX[107] = 178.054992675781250; |
162 | fX[108] = 179.054992675781250; |
163 | fX[109] = 180.054992675781250; |
164 | fX[110] = 181.054992675781250; |
165 | fX[111] = 182.054992675781250; |
166 | fX[112] = 183.054992675781250; |
167 | fX[113] = 184.054992675781250; |
168 | fX[114] = 185.054992675781250; |
169 | fX[115] = 186.054992675781250; |
170 | fX[116] = 187.054992675781250; |
171 | fX[117] = 188.054992675781250; |
172 | fX[118] = 189.054992675781250; |
173 | fX[119] = 190.054992675781250; |
174 | fX[120] = 191.054992675781250; |
175 | fX[121] = 192.054992675781250; |
176 | fX[122] = 193.054992675781250; |
177 | fX[123] = 194.054992675781250; |
178 | fX[124] = 195.054992675781250; |
179 | fX[125] = 196.054992675781250; |
180 | fX[126] = 197.054992675781250; |
181 | fX[127] = 198.054992675781250; |
182 | fX[128] = 199.054992675781250; |
183 | fX[129] = 200.054992675781250; |
184 | fX[130] = 201.054992675781250; |
185 | fX[131] = 202.054992675781250; |
186 | fX[132] = 203.054992675781250; |
187 | fX[133] = 204.054992675781250; |
188 | fX[134] = 205.054992675781250; |
189 | fX[135] = 206.054992675781250; |
190 | fX[136] = 207.054992675781250; |
191 | fX[137] = 208.054992675781250; |
192 | fX[138] = 209.054992675781250; |
193 | fX[139] = 210.054992675781250; |
194 | fX[140] = 211.054992675781250; |
195 | fX[141] = 212.054992675781250; |
196 | fX[142] = 213.054992675781250; |
197 | fX[143] = 214.054992675781250; |
198 | fX[144] = 215.054992675781250; |
199 | fX[145] = 216.054992675781250; |
200 | fX[146] = 217.054992675781250; |
201 | fX[147] = 218.054992675781250; |
202 | fX[148] = 219.054992675781250; |
203 | fX[149] = 220.054992675781250; |
204 | fX[150] = 221.054992675781250; |
205 | fX[151] = 222.054992675781250; |
206 | fX[152] = 223.054992675781250; |
207 | fX[153] = 224.054992675781250; |
208 | fX[154] = 225.054992675781250; |
209 | fX[155] = 226.054992675781250; |
210 | fX[156] = 227.054992675781250; |
211 | fX[157] = 228.054992675781250; |
212 | fX[158] = 229.054992675781250; |
213 | fX[159] = 230.054992675781250; |
214 | fX[160] = 231.054992675781250; |
215 | fX[161] = 232.054992675781250; |
216 | fX[162] = 233.054992675781250; |
217 | fX[163] = 234.054992675781250; |
218 | fX[164] = 235.054992675781250; |
219 | fX[165] = 236.054992675781250; |
220 | fX[166] = 237.054992675781250; |
221 | fX[167] = 238.054992675781250; |
222 | fX[168] = 239.054992675781250; |
223 | fX[169] = 240.054992675781250; |
224 | fX[170] = 241.054992675781250; |
225 | fX[171] = 242.054992675781250; |
226 | fX[172] = 243.054992675781250; |
227 | fX[173] = 244.054992675781250; |
228 | fX[174] = 245.054992675781250; |
229 | fX[175] = 246.054992675781250; |
230 | fNPads[0] = 67; |
231 | fNPads[1] = 67; |
232 | fNPads[2] = 67; |
233 | fNPads[3] = 69; |
234 | fNPads[4] = 69; |
235 | fNPads[5] = 69; |
236 | fNPads[6] = 71; |
237 | fNPads[7] = 71; |
238 | fNPads[8] = 71; |
239 | fNPads[9] = 73; |
240 | fNPads[10] = 73; |
241 | fNPads[11] = 73; |
242 | fNPads[12] = 75; |
243 | fNPads[13] = 75; |
244 | fNPads[14] = 75; |
245 | fNPads[15] = 77; |
246 | fNPads[16] = 77; |
247 | fNPads[17] = 77; |
248 | fNPads[18] = 79; |
249 | fNPads[19] = 79; |
250 | fNPads[20] = 79; |
251 | fNPads[21] = 81; |
252 | fNPads[22] = 81; |
253 | fNPads[23] = 81; |
254 | fNPads[24] = 83; |
255 | fNPads[25] = 83; |
256 | fNPads[26] = 83; |
257 | fNPads[27] = 85; |
258 | fNPads[28] = 85; |
259 | fNPads[29] = 85; |
260 | fNPads[30] = 87; |
261 | fNPads[31] = 87; |
262 | fNPads[32] = 87; |
263 | fNPads[33] = 89; |
264 | fNPads[34] = 89; |
265 | fNPads[35] = 89; |
266 | fNPads[36] = 91; |
267 | fNPads[37] = 91; |
268 | fNPads[38] = 91; |
269 | fNPads[39] = 93; |
270 | fNPads[40] = 93; |
271 | fNPads[41] = 93; |
272 | fNPads[42] = 95; |
273 | fNPads[43] = 95; |
274 | fNPads[44] = 95; |
275 | fNPads[45] = 97; |
276 | fNPads[46] = 97; |
277 | fNPads[47] = 97; |
278 | fNPads[48] = 99; |
279 | fNPads[49] = 99; |
280 | fNPads[50] = 99; |
281 | fNPads[51] = 101; |
282 | fNPads[52] = 101; |
283 | fNPads[53] = 101; |
284 | fNPads[54] = 103; |
285 | fNPads[55] = 103; |
286 | fNPads[56] = 103; |
287 | fNPads[57] = 105; |
288 | fNPads[58] = 105; |
289 | fNPads[59] = 105; |
290 | fNPads[60] = 107; |
291 | fNPads[61] = 107; |
292 | fNPads[62] = 107; |
293 | fNPads[63] = 109; |
294 | fNPads[64] = 73; |
295 | fNPads[65] = 75; |
296 | fNPads[66] = 75; |
297 | fNPads[67] = 75; |
298 | fNPads[68] = 77; |
299 | fNPads[69] = 77; |
300 | fNPads[70] = 77; |
301 | fNPads[71] = 77; |
302 | fNPads[72] = 79; |
303 | fNPads[73] = 79; |
304 | fNPads[74] = 79; |
305 | fNPads[75] = 81; |
306 | fNPads[76] = 81; |
307 | fNPads[77] = 81; |
308 | fNPads[78] = 83; |
309 | fNPads[79] = 83; |
310 | fNPads[80] = 83; |
311 | fNPads[81] = 83; |
312 | fNPads[82] = 85; |
313 | fNPads[83] = 85; |
314 | fNPads[84] = 85; |
315 | fNPads[85] = 87; |
316 | fNPads[86] = 87; |
317 | fNPads[87] = 87; |
318 | fNPads[88] = 87; |
319 | fNPads[89] = 89; |
320 | fNPads[90] = 89; |
321 | fNPads[91] = 89; |
322 | fNPads[92] = 91; |
323 | fNPads[93] = 91; |
324 | fNPads[94] = 91; |
325 | fNPads[95] = 93; |
326 | fNPads[96] = 93; |
327 | fNPads[97] = 93; |
328 | fNPads[98] = 93; |
329 | fNPads[99] = 95; |
330 | fNPads[100] = 95; |
331 | fNPads[101] = 95; |
332 | fNPads[102] = 97; |
333 | fNPads[103] = 97; |
334 | fNPads[104] = 97; |
335 | fNPads[105] = 97; |
336 | fNPads[106] = 99; |
337 | fNPads[107] = 99; |
338 | fNPads[108] = 99; |
339 | fNPads[109] = 101; |
340 | fNPads[110] = 101; |
341 | fNPads[111] = 101; |
342 | fNPads[112] = 103; |
343 | fNPads[113] = 103; |
344 | fNPads[114] = 103; |
345 | fNPads[115] = 103; |
346 | fNPads[116] = 105; |
347 | fNPads[117] = 105; |
348 | fNPads[118] = 105; |
349 | fNPads[119] = 107; |
350 | fNPads[120] = 107; |
351 | fNPads[121] = 107; |
352 | fNPads[122] = 107; |
353 | fNPads[123] = 109; |
354 | fNPads[124] = 109; |
355 | fNPads[125] = 109; |
356 | fNPads[126] = 111; |
357 | fNPads[127] = 111; |
358 | fNPads[128] = 111; |
359 | fNPads[129] = 113; |
360 | fNPads[130] = 113; |
361 | fNPads[131] = 113; |
362 | fNPads[132] = 113; |
363 | fNPads[133] = 115; |
364 | fNPads[134] = 115; |
365 | fNPads[135] = 115; |
366 | fNPads[136] = 117; |
367 | fNPads[137] = 117; |
368 | fNPads[138] = 117; |
369 | fNPads[139] = 117; |
370 | fNPads[140] = 119; |
371 | fNPads[141] = 119; |
372 | fNPads[142] = 119; |
373 | fNPads[143] = 121; |
374 | fNPads[144] = 121; |
375 | fNPads[145] = 121; |
376 | fNPads[146] = 123; |
377 | fNPads[147] = 123; |
378 | fNPads[148] = 123; |
379 | fNPads[149] = 123; |
380 | fNPads[150] = 125; |
381 | fNPads[151] = 125; |
382 | fNPads[152] = 125; |
383 | fNPads[153] = 127; |
384 | fNPads[154] = 127; |
385 | fNPads[155] = 127; |
386 | fNPads[156] = 127; |
387 | fNPads[157] = 129; |
388 | fNPads[158] = 129; |
389 | fNPads[159] = 129; |
390 | fNPads[160] = 131; |
391 | fNPads[161] = 131; |
392 | fNPads[162] = 131; |
393 | fNPads[163] = 133; |
394 | fNPads[164] = 133; |
395 | fNPads[165] = 133; |
396 | fNPads[166] = 133; |
397 | fNPads[167] = 135; |
398 | fNPads[168] = 135; |
399 | fNPads[169] = 135; |
400 | fNPads[170] = 137; |
401 | fNPads[171] = 137; |
402 | fNPads[172] = 137; |
403 | fNPads[173] = 137; |
404 | fNPads[174] = 139; |
405 | fNPads[175] = 139; |
406 | } |
108615fc |
407 | |
408 | Double_t AliL3Transform::GetEta(Float_t *xyz) |
409 | { |
410 | Double_t r3 = sqrt(xyz[0]*xyz[0]+xyz[1]*xyz[1]+xyz[2]*xyz[2]); |
411 | Double_t eta = 0.5 * log((r3+xyz[2])/(r3-xyz[2])); |
412 | return eta; |
413 | } |
414 | |
3a1c3849 |
415 | Double_t AliL3Transform::GetEta(Int_t padrow,Int_t pad,Int_t time) |
bd7e82e0 |
416 | { |
417 | Float_t xyz[3]; |
3a1c3849 |
418 | Int_t sector,row; |
419 | Slice2Sector(0,padrow,sector,row); |
420 | Raw2Local(xyz,sector,row,pad,time); |
bd7e82e0 |
421 | |
422 | return GetEta(xyz); |
423 | } |
424 | |
425 | |
108615fc |
426 | Double_t AliL3Transform::GetPhi(Float_t *xyz) |
427 | { |
428 | |
429 | Double_t phi = atan2(xyz[1],xyz[0]); |
430 | //if(phi<0) phi=phi+2*TMath::Pi(); |
431 | return phi; |
432 | } |
433 | |
434 | |
435 | Bool_t AliL3Transform::Slice2Sector(Int_t slice, Int_t slicerow, Int_t & sector, Int_t &row) const{ |
436 | if(slicerow<0&&slicerow>=fNRow) return kFALSE; |
437 | if(slice<0||slice>=fNSlice) return kFALSE; |
438 | |
439 | if(slicerow<fNRowLow){ |
440 | sector = slice; |
441 | row = slicerow; |
442 | } |
443 | else { |
444 | sector = slice+fNSlice; |
445 | row = slicerow-fNRowLow; |
446 | } |
447 | return kTRUE; |
448 | } |
449 | |
450 | Bool_t AliL3Transform::Sector2Slice(Int_t & slice, Int_t sector) const{ |
451 | if(sector<0||sector>=fNSector) return kFALSE; |
452 | if(sector<fNSectorLow) slice = sector; |
453 | else slice = sector - fNSectorLow; |
454 | return kTRUE; |
455 | } |
456 | |
457 | Bool_t AliL3Transform::Sector2Slice(Int_t & slice, Int_t & slicerow,Int_t sector, Int_t row) const{ |
458 | if(sector<0||sector>=fNSector||row<0) return kFALSE; |
459 | if(sector<fNSectorLow){ |
460 | if(row>=fNRowLow) return kFALSE; |
461 | slice = sector; |
462 | slicerow = row; |
463 | } |
464 | else{ |
465 | if(row>=fNRowUp) return kFALSE; |
466 | slice = sector - fNSectorLow; |
467 | slicerow = row + fNRowLow; |
468 | } |
469 | return kTRUE; |
470 | } |
471 | |
472 | Double_t AliL3Transform::Row2X(Int_t slicerow){ |
473 | if(slicerow<0||slicerow>=fNRow) return 0; |
474 | return fX[slicerow]; |
475 | } |
476 | |
477 | void AliL3Transform::Local2Global(Float_t *xyz,Int_t slice) |
478 | { |
479 | //Transformation to global coordinate system |
480 | Float_t x0 = xyz[0]; |
481 | Float_t y0 = xyz[1]; |
482 | Float_t cs,sn; |
483 | cs = fCos[slice]; |
484 | sn = fSin[slice]; |
485 | xyz[0]=x0*cs-y0*sn; |
486 | xyz[1]=x0*sn+y0*cs; |
487 | xyz[2]=xyz[2];//global z=local z |
488 | } |
489 | |
490 | void AliL3Transform::Local2GlobalAngle(Float_t *angle,Int_t slice){ |
40f2887b |
491 | angle[0] = fmod(angle[0]+(slice+0.5)*(2*fPi/18),2*fPi); |
108615fc |
492 | } |
493 | |
494 | void AliL3Transform::Global2LocalAngle(Float_t *angle,Int_t slice){ |
40f2887b |
495 | angle[0] = angle[0]-(slice+0.5)*(2*fPi/18); |
108615fc |
496 | if(angle[0]<0) angle[0]+=2*fPi; |
497 | } |
498 | |
499 | void AliL3Transform::Raw2Local(Float_t *xyz,Int_t sector,Int_t row,Float_t pad,Float_t time) |
500 | { |
501 | //Transformation from rawdata to local coordinate system |
502 | |
503 | Int_t slice,slicerow; |
504 | Sector2Slice(slice, slicerow, sector, row); |
505 | |
506 | xyz[0]=Row2X(slicerow); |
507 | Int_t npads= fNPads[slicerow]; |
508 | if(sector<fNSectorLow) |
509 | xyz[1]=(pad-0.5*(npads-1))*fPadPitchWidthLow; |
510 | else |
511 | xyz[1]=(pad-0.5*(npads-1))*fPadPitchWidthUp; |
512 | xyz[2]=fZWidth*time-3.*fZSigma; |
1c225117 |
513 | Int_t sign=1; |
514 | |
515 | |
0e0b1e59 |
516 | if(slice < 18) |
517 | sign = 1; |
518 | else |
519 | sign = -1; |
0e0b1e59 |
520 | |
1c225117 |
521 | xyz[2]=sign*(250.-xyz[2]); |
0e0b1e59 |
522 | |
523 | //Int_t nis=fNSectorLow; |
524 | //Int_t nos=fNSectorUp; |
525 | //if((sector<nis)/2 || ((sector-nis)<nos/2)) sign=1; |
526 | |
108615fc |
527 | } |
528 | |
529 | void AliL3Transform::Local2Global(Float_t *xyz,Int_t sector,Int_t row) |
530 | { |
531 | //Transformation to global coordinate system |
532 | Int_t slice,slicerow; |
533 | Sector2Slice(slice, slicerow, sector, row); |
534 | Float_t r=Row2X(slicerow); |
535 | Float_t cs = fCos[slice]; |
536 | Float_t sn = fSin[slice]; |
537 | |
538 | xyz[0]=r*cs-xyz[1]*sn; |
539 | xyz[1]=r*sn+xyz[1]*cs; |
540 | xyz[2]=xyz[2];//global z=local z |
541 | } |
542 | |
543 | Double_t AliL3Transform::GetMaxY(Int_t slicerow) |
544 | { |
545 | |
546 | if(slicerow < fNRowLow) |
547 | return fPadPitchWidthLow*fNPads[slicerow]/2; |
548 | |
549 | else |
550 | return fPadPitchWidthUp*fNPads[slicerow]/2; |
551 | } |
552 | |
2342214c |
553 | void AliL3Transform::Global2Local(Float_t *xyz,Int_t sector,Bool_t isSlice) |
108615fc |
554 | { |
2342214c |
555 | |
108615fc |
556 | Int_t slice; |
2342214c |
557 | if(!isSlice) |
558 | Sector2Slice(slice, sector); |
559 | else |
560 | slice = sector; |
108615fc |
561 | Float_t cs = fCos[slice]; |
562 | Float_t sn = fSin[slice]; |
563 | Float_t x1 = xyz[0]*cs + xyz[1]*sn; |
564 | Float_t y1 = -xyz[0]*sn + xyz[1]*cs; |
565 | xyz[0] = x1; |
566 | xyz[1] = y1; |
567 | } |
568 | |
569 | void AliL3Transform::Raw2Global(Float_t *xyz,Int_t sector,Int_t row,Float_t pad,Float_t time) |
570 | { |
571 | //Transformation from raw to global coordinates |
572 | |
573 | Raw2Local(xyz,sector,row,pad,time); |
574 | Local2Global(xyz,sector,row); |
575 | } |
576 | |
577 | void AliL3Transform::Local2Raw(Float_t *xyz,Int_t sector,Int_t row) |
578 | { |
579 | //Transformation from local coordinates to raw |
580 | |
581 | Int_t slice,slicerow; |
582 | Sector2Slice(slice, slicerow, sector, row); |
583 | |
584 | if(sector<fNSectorLow) |
585 | xyz[1]=xyz[1]/fPadPitchWidthLow+0.5*(fNPads[slicerow]-1); |
586 | else |
587 | xyz[1]=xyz[1]/fPadPitchWidthUp+0.5*(fNPads[slicerow]-1); |
588 | Int_t sign=-1; |
589 | Int_t nis=fNSectorLow; |
590 | Int_t nos=fNSectorUp; |
591 | |
592 | if ((sector<nis/2) || ((sector-nis)<nos/2)) sign=1; |
593 | xyz[2]=250-sign*xyz[2]; |
594 | xyz[2]=(xyz[2]+3.*fZSigma)/fZWidth; |
595 | } |
596 | |
597 | void AliL3Transform::Global2Raw(Float_t *xyz,Int_t sector,Int_t row) |
598 | { |
599 | //Transformation from global coordinates to raw. |
600 | |
601 | Global2Local(xyz,sector); |
602 | Local2Raw(xyz,sector,row); |
603 | |
604 | } |