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