]>
Commit | Line | Data |
---|---|---|
e0679afb | 1 | /************************************************************************** |
b92524d0 | 2 | * This file is property of and copyright by the ALICE HLT Project * |
3 | * All rights reserved. * | |
e0679afb | 4 | * * |
b92524d0 | 5 | * Primary Authors: * |
6 | * Artur Szostak <artursz@iafrica.com> * | |
e0679afb | 7 | * * |
8 | * Permission to use, copy, modify and distribute this software and its * | |
9 | * documentation strictly for non-commercial purposes is hereby granted * | |
10 | * without fee, provided that the above copyright notice appears in all * | |
11 | * copies and that both the copyright notice and this permission notice * | |
12 | * appear in the supporting documentation. The authors make no claims * | |
b92524d0 | 13 | * about the suitability of this software for any purpose. It is * |
e0679afb | 14 | * provided "as is" without express or implied warranty. * |
15 | **************************************************************************/ | |
16 | ||
17 | /* $Id$ */ | |
18 | ||
b92524d0 | 19 | /** |
20 | * @file AliHLTMUONMansoTrackerFSM.cxx | |
21 | * @author Artur Szostak <artursz@iafrica.com> | |
22 | * @date | |
23 | * @brief Implementation of AliHLTMUONMansoTrackerFSM class. | |
24 | */ | |
e0679afb | 25 | |
26 | #include "AliHLTMUONMansoTrackerFSM.h" | |
27 | #include "AliHLTMUONCalculations.h" | |
b92524d0 | 28 | #include "AliHLTMUONConstants.h" |
e0679afb | 29 | #include "AliHLTMUONUtils.h" |
30 | #include <cmath> | |
31 | ||
32 | ||
33 | #ifdef DEBUG | |
34 | #include <ostream> | |
35 | namespace | |
36 | { | |
37 | ||
38 | std::ostream& operator << (std::ostream& os, AliHLTMUONMansoTrackerFSM::StatesSM4 state) | |
39 | { | |
40 | switch (state) | |
41 | { | |
42 | case AliHLTMUONMansoTrackerFSM::kSM4Idle: os << "kSM4Idle"; break; | |
43 | case AliHLTMUONMansoTrackerFSM::kWaitChamber8: os << "kWaitChamber8"; break; | |
44 | case AliHLTMUONMansoTrackerFSM::kWaitMoreChamber8: os << "kWaitMoreChamber8"; break; | |
45 | case AliHLTMUONMansoTrackerFSM::kWaitChamber7: os << "kWaitChamber7"; break; | |
46 | case AliHLTMUONMansoTrackerFSM::kWaitMoreChamber7: os << "kWaitMoreChamber7"; break; | |
b92524d0 | 47 | default: os << "FAULT!!"; |
e0679afb | 48 | } |
49 | return os; | |
50 | } | |
51 | ||
52 | std::ostream& operator << (std::ostream& os, AliHLTMUONMansoTrackerFSM::StatesSM5 state) | |
53 | { | |
54 | switch (state) | |
55 | { | |
56 | case AliHLTMUONMansoTrackerFSM::kSM5Idle: os << "kSM5Idle"; break; | |
57 | case AliHLTMUONMansoTrackerFSM::kWaitChamber10: os << "kWaitChamber10"; break; | |
58 | case AliHLTMUONMansoTrackerFSM::kWaitMoreChamber10: os << "kWaitMoreChamber10"; break; | |
59 | case AliHLTMUONMansoTrackerFSM::kWaitChamber9: os << "kWaitChamber9"; break; | |
60 | case AliHLTMUONMansoTrackerFSM::kWaitMoreChamber9: os << "kWaitMoreChamber9"; break; | |
61 | case AliHLTMUONMansoTrackerFSM::kSM5Done: os << "kSM5Done"; break; | |
b92524d0 | 62 | default: os << "FAULT!!"; |
e0679afb | 63 | } |
64 | return os; | |
65 | } | |
66 | ||
67 | } // end of namespace | |
68 | #endif // DEBUG | |
69 | ||
70 | ||
71 | // Deviate from the Manso implementation by allowing a and b | |
72 | // parameters per chamber and not just per station. | |
73 | // The default values are derived from the work done in | |
74 | // "A first algorithm for dimuon High Level Trigger" | |
75 | // Ref ID: ALICE-INT-2002-04 version 1.0 | |
76 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgA7 = 0.016f; | |
77 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgB7 = 2.0f; | |
78 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgA8 = 0.016f; | |
79 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgB8 = 2.0f; | |
80 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgA9 = 0.020f; | |
81 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgB9 = 3.0f; | |
82 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgA10 = 0.020f; | |
83 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgB10 = 3.0f; | |
b92524d0 | 84 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgZ7 = -1274.5f; |
85 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgZ8 = -1305.5f; | |
86 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgZ9 = -1408.6f; | |
87 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgZ10 = -1439.6f; | |
88 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgZ11 = -1603.5f; | |
89 | AliHLTFloat32_t AliHLTMUONMansoTrackerFSM::fgZ13 = -1703.5f; | |
e0679afb | 90 | |
91 | ||
b92524d0 | 92 | void AliHLTMUONMansoTrackerFSM::AliRegionOfInterest::Create( |
93 | AliHLTMUONRecHitStruct p, AliHLTFloat32_t a, AliHLTFloat32_t b | |
94 | ) | |
e0679afb | 95 | { |
96 | // Creates a region of interest specific to the Manso algorithm from a point and | |
97 | // two Manso specific parameters. | |
98 | ||
99 | fCentre = p; | |
100 | // Compute the radius Rp | |
b92524d0 | 101 | AliHLTFloat32_t rp = (AliHLTFloat32_t) sqrt( p.fX * p.fX + p.fY * p.fY ); |
e0679afb | 102 | |
103 | // The radius Rs for the region of interest is computed from the | |
104 | // specification given in the document: | |
105 | // "A first algorithm for dimuon High Level Trigger" | |
106 | // Ref ID: ALICE-INT-2002-04 version 1.0 | |
107 | // equation: | |
108 | // Rs = a * Rp + b | |
109 | // given on page 3 section 4. | |
110 | fRs = a * rp + b; | |
111 | } | |
112 | ||
113 | ||
b92524d0 | 114 | bool AliHLTMUONMansoTrackerFSM::AliRegionOfInterest::Contains(AliHLTMUONRecHitStruct p) const |
e0679afb | 115 | { |
116 | // Compute the distance between the centre of the region of interest and | |
117 | // the point p. This distance must be less than the radius of the region | |
118 | // of interest for p to be contained in the region of interest. | |
b92524d0 | 119 | register AliHLTFloat32_t lx = fCentre.fX - p.fX; |
120 | register AliHLTFloat32_t ly = fCentre.fY - p.fY; | |
e0679afb | 121 | register AliHLTFloat32_t r = (AliHLTFloat32_t) sqrt( lx * lx + ly * ly ); |
122 | DebugTrace("\tAliRegionOfInterest::Contains : p = " << p | |
123 | << " , centre = " << fCentre << " , r = " << r << " , Rs = " << fRs | |
124 | ); | |
125 | return r <= fRs; | |
126 | } | |
127 | ||
128 | ||
129 | void AliHLTMUONMansoTrackerFSM::AliRegionOfInterest::GetBoundaryBox( | |
b92524d0 | 130 | AliHLTFloat32_t& left, AliHLTFloat32_t& right, |
131 | AliHLTFloat32_t& bottom, AliHLTFloat32_t& top | |
e0679afb | 132 | ) const |
133 | { | |
134 | // Works out the smallest boundary box that will contain the region of interest. | |
135 | ||
b92524d0 | 136 | left = fCentre.fX - fRs; |
137 | right = fCentre.fX + fRs; | |
138 | bottom = fCentre.fY - fRs; | |
139 | top = fCentre.fY + fRs; | |
e0679afb | 140 | } |
141 | ||
142 | ||
b92524d0 | 143 | AliHLTMUONMansoTrackerFSM::AliVertex::AliVertex( |
144 | AliHLTFloat32_t x, AliHLTFloat32_t y, AliHLTFloat32_t z | |
145 | ) | |
e0679afb | 146 | : fX(x), fY(y), fZ(z) |
147 | { | |
148 | // Constructor for vertex. | |
e0679afb | 149 | } |
150 | ||
151 | ||
b92524d0 | 152 | AliHLTMUONMansoTrackerFSM::AliVertex::AliVertex(AliHLTMUONRecHitStruct xy, AliHLTFloat32_t z) |
153 | : fX(xy.fX), fY(xy.fY), fZ(z) | |
e0679afb | 154 | { |
155 | // Construct vertex from a point on the XY plane and z coordinate. | |
e0679afb | 156 | } |
157 | ||
158 | ||
159 | AliHLTMUONMansoTrackerFSM::AliLine::AliLine( | |
b92524d0 | 160 | AliHLTFloat32_t ax, AliHLTFloat32_t ay, AliHLTFloat32_t az, |
161 | AliHLTFloat32_t bx, AliHLTFloat32_t by, AliHLTFloat32_t bz | |
162 | ) : | |
163 | fMx(ax - bx), | |
164 | fMy(ay - by), | |
165 | fMz(az - bz), | |
166 | fCx(bx), | |
167 | fCy(by), | |
168 | fCz(bz) | |
e0679afb | 169 | { |
170 | // Construct a line defined by L = M*t + C = (A-B)*t + B | |
171 | // where M and C are 3D vectors and t is a free parameter. | |
172 | // A = (ax, ay, az) and B = (bx, by, bz) | |
e0679afb | 173 | } |
174 | ||
175 | ||
176 | AliHLTMUONMansoTrackerFSM::AliLine::AliLine(AliVertex a, AliVertex b) : | |
b92524d0 | 177 | fMx(a.X() - b.X()), |
178 | fMy(a.Y() - b.Y()), | |
179 | fMz(a.Z() - b.Z()), | |
180 | fCx(b.X()), | |
181 | fCy(b.Y()), | |
182 | fCz(b.Z()) | |
e0679afb | 183 | { |
184 | // Contruct a line to go through two vertices a and b. | |
e0679afb | 185 | } |
186 | ||
187 | ||
b92524d0 | 188 | AliHLTMUONRecHitStruct AliHLTMUONMansoTrackerFSM::AliLine::FindIntersectWithXYPlain( |
189 | AliHLTFloat32_t z | |
190 | ) const | |
e0679afb | 191 | { |
192 | // Find the point of intersection of the line and the XY plane at z. | |
193 | ||
b92524d0 | 194 | assert( fMz != 0.0 ); // Should not have a ray perpendicular to the beam axis. |
e0679afb | 195 | AliHLTFloat32_t t = (z - fCz) / fMz; |
b92524d0 | 196 | AliHLTMUONRecHitStruct p; |
197 | p.fX = fMx*t + fCx; | |
198 | p.fY = fMy*t + fCy; | |
199 | p.fZ = z; | |
200 | return p; | |
e0679afb | 201 | } |
202 | ||
203 | ||
204 | AliHLTMUONMansoTrackerFSM::AliHLTMUONMansoTrackerFSM() : | |
205 | fCallback(NULL), | |
206 | fSm4state(kSM4Idle), | |
207 | fSm5state(kSM5Idle), | |
208 | fRequestsCompleted(0), | |
209 | fSt4chamber(kChamber1), | |
210 | fV1(), | |
211 | fMc1(), | |
212 | fSt5z(0), | |
213 | fSt5data(), | |
214 | fSt4z(0), | |
215 | fSt4points(), | |
216 | fSt5rec(), | |
b92524d0 | 217 | fFoundPoint(), |
218 | fTriggerId(-1), | |
219 | fTrackId(0) | |
e0679afb | 220 | { |
221 | // Default constructor | |
222 | } | |
223 | ||
224 | ||
b92524d0 | 225 | void AliHLTMUONMansoTrackerFSM::FindTrack(const AliHLTMUONTriggerRecordStruct& trigger) |
e0679afb | 226 | { |
227 | // Tries to find the track from the trigger seed. | |
228 | ||
229 | DebugTrace("SM5 state = " << fSm5state << " , SM4 state = " << fSm4state); | |
b92524d0 | 230 | DebugTrace("Processing trigger with ID = " << trigger.fId); |
231 | ||
232 | fTriggerId = trigger.fId; | |
233 | ||
234 | AliHLTMUONParticleSign sign; | |
235 | bool hitset[4]; | |
236 | AliHLTMUONUtils::UnpackTriggerRecordFlags(trigger.fFlags, sign, hitset); | |
237 | DebugTrace("hitset = {" << hitset[0] << ", " << hitset[1] << ", " | |
238 | << hitset[2] << ", " << hitset[3] << "}" | |
239 | ); | |
240 | ||
241 | // Find first valid hit in the trigger record from chambers 11 and 12. | |
242 | int firstHit = -1; | |
243 | for (int i = 0; i < 2; i++) | |
244 | { | |
245 | if (hitset[i]) | |
246 | { | |
247 | firstHit = i; | |
248 | break; | |
249 | } | |
250 | } | |
251 | // Find the second valid hit in the trigger record from chambers 13 and 14. | |
252 | int secondHit = -1; | |
253 | for (int i = 2; i < 4; i++) | |
254 | { | |
255 | if (hitset[i]) | |
256 | { | |
257 | secondHit = i; | |
258 | break; | |
259 | } | |
260 | } | |
261 | if (firstHit == -1 or secondHit == -1) | |
262 | { | |
263 | NoTrackFound(); | |
264 | return; | |
265 | } | |
266 | ||
267 | fV1 = AliVertex( | |
268 | trigger.fHit[firstHit].fX, | |
269 | trigger.fHit[firstHit].fY, | |
270 | trigger.fHit[firstHit].fZ | |
271 | ); | |
272 | AliVertex v2 = AliVertex( | |
273 | trigger.fHit[secondHit].fX, | |
274 | trigger.fHit[secondHit].fY, | |
275 | trigger.fHit[secondHit].fZ | |
276 | ); | |
277 | ||
278 | DebugTrace("Using fV1 = {x = " << fV1.X() << ", y = " << fV1.Y() << ", " | |
279 | << fV1.Z() << "}, with firstHit = " << firstHit | |
280 | ); | |
281 | DebugTrace("Using v2 = {x = " << v2.X() << ", y = " << v2.Y() << ", " | |
282 | << v2.Z() << "}, with secondHit = " << secondHit | |
283 | ); | |
e0679afb | 284 | |
285 | // Form the vector line between the above two impact points and | |
286 | // find the crossing point of the line with chamber 10 (i.e. station 5). | |
287 | fMc1.fLine = AliLine(fV1, v2); | |
b92524d0 | 288 | AliHLTMUONRecHitStruct p10 = fMc1.fLine.FindIntersectWithXYPlain( fgZ10 ); |
e0679afb | 289 | |
290 | // Build a region of interest for tracking station 5 (chamber 10). | |
291 | // Remember the parameters a and b are station specific. | |
292 | fMc1.fChamber = kChamber10; | |
293 | fMc1.fRoi.Create(p10, fgA10, fgB10); | |
294 | ||
295 | // Make SM5 state transition before the call to RequestClusters since | |
296 | // that method could call one of our methods again, so we need to be | |
297 | // in a consistant internal state. | |
298 | fSm5state = kWaitChamber10; | |
299 | ||
300 | AliHLTFloat32_t left, right, bottom, top; | |
301 | fMc1.fRoi.GetBoundaryBox(left, right, bottom, top); | |
302 | RequestClusters(left, right, bottom, top, kChamber10, &fMc1); | |
303 | } | |
304 | ||
305 | ||
b92524d0 | 306 | void AliHLTMUONMansoTrackerFSM::ReturnClusters( |
307 | void* tag, const AliHLTMUONRecHitStruct* clusters, | |
308 | AliHLTUInt32_t count | |
309 | ) | |
e0679afb | 310 | { |
b92524d0 | 311 | // Implementation of AliHLTMUONMansoTrackerFSM::ReturnClusters. |
e0679afb | 312 | |
b92524d0 | 313 | assert( count > 0 ); |
314 | assert( clusters != NULL ); | |
e0679afb | 315 | |
316 | AliTagData* data = (AliTagData*)tag; | |
317 | DebugTrace("Got AliHLTMUONMansoTrackerFSM::ReturnClusters(tag = " << tag | |
318 | << ", chamber = " << data->fChamber | |
319 | << ", clusters = " << clusters << ", count = " << count << ")" | |
320 | ); | |
321 | DebugTrace("SM5 state = " << fSm5state << " , SM4 state = " << fSm4state); | |
322 | ||
323 | switch (data->fChamber) | |
324 | { | |
325 | case kChamber7: ReceiveClustersChamber7(clusters, count, data); break; | |
326 | case kChamber8: ReceiveClustersChamber8(clusters, count, data); break; | |
327 | case kChamber9: ReceiveClustersChamber9(clusters, count); break; | |
328 | case kChamber10: ReceiveClustersChamber10(clusters, count); break; | |
329 | default: | |
330 | // Error | |
331 | DebugTrace("ERROR: Got tag with an invalid value: " << data->fChamber); | |
332 | } | |
333 | } | |
334 | ||
335 | ||
336 | void AliHLTMUONMansoTrackerFSM::EndOfClusters(void* tag) | |
337 | { | |
b92524d0 | 338 | // Implementation of AliHLTMUONMansoTrackerFSM::EndOfClusters. |
e0679afb | 339 | |
340 | AliTagData* data = (AliTagData*)tag; | |
341 | DebugTrace("Got AliHLTMUONMansoTrackerFSM::EndOfClusters(chamber = " << data->fChamber << ")"); | |
342 | DebugTrace("SM5 state = " << fSm5state << " , SM4 state = " << fSm4state); | |
343 | ||
344 | switch (data->fChamber) | |
345 | { | |
346 | case kChamber7: EndOfClustersChamber7(); break; | |
347 | case kChamber8: EndOfClustersChamber8(); break; | |
348 | case kChamber9: EndOfClustersChamber9(); break; | |
349 | case kChamber10: EndOfClustersChamber10(); break; | |
350 | default: | |
351 | // Error | |
352 | DebugTrace("ERROR: Got tag with an invalid value: " << data->fChamber); | |
353 | } | |
354 | } | |
355 | ||
356 | ||
b92524d0 | 357 | bool AliHLTMUONMansoTrackerFSM::FillTrackData(AliHLTMUONMansoTrackStruct& track) |
e0679afb | 358 | { |
b92524d0 | 359 | // Implementation of AliHLTMUONMansoTrackerFSM::FillTrackData |
e0679afb | 360 | |
361 | DebugTrace("FillTrack: st5 = " << fSt5rec->fClusterPoint << ", st4 = " << fFoundPoint->fClusterPoint); | |
362 | ||
b92524d0 | 363 | track.fId = fTrackId; |
364 | // Increment track Id and keep it positive. | |
365 | //TODO: handle the wrapparound better. | |
366 | if (fTrackId < 0x7FFFFFFF) | |
367 | fTrackId++; | |
368 | else | |
369 | fTrackId = 0; | |
370 | ||
371 | track.fTrigRec = fTriggerId; | |
372 | ||
373 | AliHLTFloat32_t x1 = fFoundPoint->fClusterPoint.fX; | |
374 | AliHLTFloat32_t y1 = fFoundPoint->fClusterPoint.fY; | |
375 | AliHLTFloat32_t z1 = fFoundPoint->fClusterPoint.fZ; | |
376 | AliHLTFloat32_t y2 = fSt5rec->fClusterPoint.fY; | |
377 | AliHLTFloat32_t z2 = fSt5rec->fClusterPoint.fZ; | |
378 | ||
379 | bool calculated = AliHLTMUONCalculations::ComputeMomentum(x1, y1, y2, z1, z2); | |
380 | ||
381 | track.fPx = AliHLTMUONCalculations::Px(); | |
382 | track.fPy = AliHLTMUONCalculations::Py(); | |
383 | track.fPz = AliHLTMUONCalculations::Pz(); | |
384 | DebugTrace("Calculated Px = " << track.fPx << ", Py = " << track.fPy | |
385 | << ", Pz = " << track.fPx | |
386 | ); | |
e0679afb | 387 | DebugTrace("\tusing x1 = " << x1 << " , y1 = " << y1 << " , y2 = " << y2 |
b92524d0 | 388 | << " , z1 = " << z1 << " , z2 = " << z2 |
e0679afb | 389 | ); |
e0679afb | 390 | |
b92524d0 | 391 | track.fChi2 = 0; |
e0679afb | 392 | |
b92524d0 | 393 | bool hitset[4]; |
394 | // Depending on which chamber we found reconstructed hits, fill the hit | |
395 | // in the appropriate location. This is done for station 4 then 5. | |
e0679afb | 396 | if (fSt4chamber == kChamber8) |
397 | { | |
b92524d0 | 398 | track.fHit[0] = AliHLTMUONConstants::NilRecHitStruct(); |
399 | hitset[0] = false; | |
400 | track.fHit[1] = fFoundPoint->fClusterPoint; | |
401 | hitset[1] = true; | |
e0679afb | 402 | } |
403 | else | |
404 | { | |
b92524d0 | 405 | track.fHit[0] = fFoundPoint->fClusterPoint; |
406 | hitset[0] = true; | |
407 | track.fHit[1] = AliHLTMUONConstants::NilRecHitStruct(); | |
408 | hitset[1] = false; | |
e0679afb | 409 | } |
410 | if (fMc1.fChamber == kChamber10) | |
411 | { | |
b92524d0 | 412 | track.fHit[2] = AliHLTMUONConstants::NilRecHitStruct(); |
413 | hitset[2] = false; | |
414 | track.fHit[3] = fSt5rec->fClusterPoint; | |
415 | hitset[3] = true; | |
e0679afb | 416 | } |
417 | else | |
418 | { | |
b92524d0 | 419 | track.fHit[2] = fSt5rec->fClusterPoint; |
420 | hitset[2] = true; | |
421 | track.fHit[3] = AliHLTMUONConstants::NilRecHitStruct(); | |
422 | hitset[3] = false; | |
e0679afb | 423 | } |
b92524d0 | 424 | |
425 | track.fFlags = AliHLTMUONUtils::PackMansoTrackFlags( | |
426 | AliHLTMUONCalculations::Sign(), hitset | |
427 | ); | |
428 | return calculated; | |
e0679afb | 429 | } |
430 | ||
431 | ||
432 | void AliHLTMUONMansoTrackerFSM::Reset() | |
433 | { | |
b92524d0 | 434 | // Implementation of AliHLTMUONMansoTrackerFSM::Reset |
e0679afb | 435 | |
436 | DebugTrace("SM5 state = " << fSm5state << " , SM4 state = " << fSm4state); | |
437 | fSt5data.Clear(); | |
438 | fSt4points.Clear(); | |
439 | fSm4state = kSM4Idle; | |
440 | fSm5state = kSM5Idle; | |
441 | fRequestsCompleted = 0; | |
b92524d0 | 442 | fTriggerId = -1; |
e0679afb | 443 | } |
444 | ||
445 | ||
446 | // Note: In the following ReceiveClustersXXX and EndOfClustersXXX methods we make | |
447 | // the state machine transitions before calls to RequestClusters, FoundTrack, | |
448 | // NoTrackFound or EndOfClusterRequests. This is important since the callback | |
449 | // object will make recursive calls to the tracker's methods so we need to maintain | |
450 | // a consistant internal state. | |
451 | // The same would go for updating internal variables. | |
452 | // In general one should only call the callback methods at the end of any of the | |
453 | // following routines. | |
454 | ||
455 | void AliHLTMUONMansoTrackerFSM::ReceiveClustersChamber7( | |
b92524d0 | 456 | const AliHLTMUONRecHitStruct* clusters, AliHLTUInt32_t count, |
457 | const AliTagData* data | |
e0679afb | 458 | ) |
459 | { | |
460 | // State change method for Station 4 state machine. | |
461 | ||
462 | switch (fSm4state) | |
463 | { | |
464 | case kWaitChamber7: | |
465 | fSm4state = kWaitMoreChamber7; | |
466 | ||
467 | case kWaitMoreChamber7: | |
b92524d0 | 468 | for (AliHLTUInt32_t j = 0; j < count; j++) |
e0679afb | 469 | { |
b92524d0 | 470 | AliHLTMUONRecHitStruct cluster = clusters[j]; |
e0679afb | 471 | // Check that the cluster actually is in our region of interest on station 4. |
472 | if ( data->fRoi.Contains(cluster) ) | |
473 | { | |
b92524d0 | 474 | DebugTrace("Adding cluster [" << cluster.fX << ", " << cluster.fY << "] from chamber 7."); |
475 | AliStation4Data* newdata = fSt4points.Add(); | |
e0679afb | 476 | newdata->fClusterPoint = cluster; |
477 | newdata->fSt5tag = data; | |
478 | } | |
479 | } | |
480 | break; | |
481 | ||
482 | default: | |
483 | DebugTrace("ERROR: Unexpected state for SM4 in AliHLTMUONMansoTrackerFSM::ReceiveClustersChamber7!"); | |
484 | } | |
485 | } | |
486 | ||
487 | ||
488 | void AliHLTMUONMansoTrackerFSM::ReceiveClustersChamber8( | |
b92524d0 | 489 | const AliHLTMUONRecHitStruct* clusters, AliHLTUInt32_t count, |
490 | const AliTagData* data | |
e0679afb | 491 | ) |
492 | { | |
493 | // State change method for Station 4 state machine. | |
494 | ||
495 | switch (fSm4state) | |
496 | { | |
497 | case kWaitChamber8: | |
498 | fSm4state = kWaitMoreChamber8; | |
499 | fSt4z = fgZ8; | |
500 | fSt4chamber = kChamber8; | |
501 | ||
502 | case kWaitMoreChamber8: | |
b92524d0 | 503 | for (AliHLTUInt32_t j = 0; j < count; j++) |
e0679afb | 504 | { |
b92524d0 | 505 | AliHLTMUONRecHitStruct cluster = clusters[j]; |
e0679afb | 506 | // Check that the cluster actually is in our region of interest on station 4. |
507 | if ( data->fRoi.Contains(cluster) ) | |
508 | { | |
b92524d0 | 509 | DebugTrace("Adding cluster [" << cluster.fX << ", " << cluster.fY << "] from chamber 8."); |
510 | AliStation4Data* newdata = fSt4points.Add(); | |
e0679afb | 511 | newdata->fClusterPoint = cluster; |
512 | newdata->fSt5tag = data; | |
513 | } | |
514 | } | |
515 | break; | |
516 | ||
517 | default: | |
518 | DebugTrace("ERROR: Unexpected state for SM4 in AliHLTMUONMansoTrackerFSM::ReceiveClustersChamber8!"); | |
519 | } | |
520 | } | |
521 | ||
522 | ||
b92524d0 | 523 | void AliHLTMUONMansoTrackerFSM::ReceiveClustersChamber9( |
524 | const AliHLTMUONRecHitStruct* clusters, AliHLTUInt32_t count | |
525 | ) | |
e0679afb | 526 | { |
527 | // State change method for Station 5 state machine. | |
528 | ||
529 | switch (fSm5state) | |
530 | { | |
531 | case kWaitChamber9: | |
532 | fSm5state = kWaitMoreChamber9; | |
533 | fSm4state = kWaitChamber8; // Start SM4. | |
534 | ||
535 | case kWaitMoreChamber9: | |
b92524d0 | 536 | for (AliHLTUInt32_t j = 0; j < count; j++) |
e0679afb | 537 | { |
b92524d0 | 538 | AliHLTMUONRecHitStruct cluster = clusters[j]; |
e0679afb | 539 | // Check that the cluster actually is in our region of interest on station 5. |
540 | if ( fMc1.fRoi.Contains(cluster) ) | |
541 | { | |
b92524d0 | 542 | DebugTrace("Adding cluster [" << cluster.fX << ", " << cluster.fY << "] from chamber 9."); |
543 | AliStation5Data* data = fSt5data.Add(); | |
e0679afb | 544 | data->fClusterPoint = cluster; |
545 | ProjectToStation4(data, fgZ9); // This adds a new request for station 4. | |
546 | } | |
547 | } | |
548 | break; | |
549 | ||
550 | default: | |
551 | DebugTrace("ERROR: Unexpected state for SM5 in AliHLTMUONMansoTrackerFSM::ReceiveClustersChamber9!"); | |
552 | } | |
553 | } | |
554 | ||
555 | ||
b92524d0 | 556 | void AliHLTMUONMansoTrackerFSM::ReceiveClustersChamber10( |
557 | const AliHLTMUONRecHitStruct* clusters, AliHLTUInt32_t count | |
558 | ) | |
e0679afb | 559 | { |
560 | // State change method for Station 5 state machine. | |
561 | ||
562 | switch (fSm5state) | |
563 | { | |
564 | case kWaitChamber10: | |
565 | fSm5state = kWaitMoreChamber10; | |
566 | fSt5z = fgZ10; | |
567 | fSm4state = kWaitChamber8; // Start SM4. | |
568 | ||
569 | case kWaitMoreChamber10: | |
b92524d0 | 570 | for (AliHLTUInt32_t j = 0; j < count; j++) |
e0679afb | 571 | { |
b92524d0 | 572 | AliHLTMUONRecHitStruct cluster = clusters[j]; |
e0679afb | 573 | // Check that the cluster actually is in our region of interest on station 5. |
574 | if ( fMc1.fRoi.Contains(cluster) ) | |
575 | { | |
b92524d0 | 576 | DebugTrace("Adding cluster [" << cluster.fX << ", " << cluster.fY << "] from chamber 10."); |
577 | AliStation5Data* data = fSt5data.Add(); | |
e0679afb | 578 | data->fClusterPoint = cluster; |
579 | ProjectToStation4(data, fgZ10); // This adds a new request for station 4. | |
580 | } | |
581 | } | |
582 | break; | |
583 | ||
584 | default: | |
585 | DebugTrace("ERROR: Unexpected state for SM5 in AliHLTMUONMansoTrackerFSM::ReceiveClustersChamber10!"); | |
586 | } | |
587 | } | |
588 | ||
589 | ||
590 | void AliHLTMUONMansoTrackerFSM::EndOfClustersChamber7() | |
591 | { | |
592 | // State change method for Station 4 state machine. | |
593 | ||
594 | fRequestsCompleted++; // Increment the number of requests completed for station 4. | |
595 | DebugTrace("fRequestsCompleted = " << fRequestsCompleted ); | |
596 | ||
597 | switch (fSm4state) | |
598 | { | |
599 | case kWaitChamber7: | |
600 | // If all data from station 5 is received and no data found on | |
601 | // chambers 7 or 8 then we can not find a track. | |
602 | if (fSm5state == kSM5Done) NoTrackFound(); | |
603 | break; | |
604 | ||
605 | case kWaitMoreChamber7: | |
606 | if (fRequestsCompleted == fSt5data.Count() && fSm5state == kSM5Done) | |
607 | ProcessClusters(); | |
608 | break; | |
609 | ||
610 | default: | |
611 | DebugTrace("ERROR: Unexpected state for SM4 in AliHLTMUONMansoTrackerFSM::EndOfClustersChamber7!"); | |
612 | } | |
613 | } | |
614 | ||
615 | ||
616 | void AliHLTMUONMansoTrackerFSM::EndOfClustersChamber8() | |
617 | { | |
618 | // State change method for Station 4 state machine. | |
619 | ||
620 | fRequestsCompleted++; // Increment the number of requests completed for station 4. | |
621 | DebugTrace("fRequestsCompleted = " << fRequestsCompleted ); | |
622 | ||
623 | switch (fSm4state) | |
624 | { | |
625 | case kWaitChamber7: | |
626 | // Ignore. The requests for chamber 8 are already re-requested below. | |
627 | break; | |
628 | ||
629 | case kWaitChamber8: | |
630 | { | |
631 | fSm4state = kWaitChamber7; | |
632 | fSt4z = fgZ7; | |
633 | fSt4chamber = kChamber7; | |
634 | ||
635 | // We need to resend the requests for chamber 8, but change the request | |
636 | // to get data for chamber 7 instead: | |
b92524d0 | 637 | AliHLTUInt32_t reqlistsize = fSt5data.Count(); |
e0679afb | 638 | DebugTrace("Re-requesting clusters from chamber 7... reqlistsize = " << reqlistsize); |
639 | ||
640 | Station5List::Iterator rec = fSt5data.First(); | |
b92524d0 | 641 | for (AliHLTUInt32_t i = 0; i < reqlistsize; i++, rec++) |
e0679afb | 642 | { |
643 | // Need to create a new st5 data block for the request. | |
b92524d0 | 644 | AliStation5Data* data = fSt5data.Add(); |
e0679afb | 645 | data->fClusterPoint = rec->fClusterPoint; |
646 | data->fTag.fLine = rec->fTag.fLine; | |
647 | ||
648 | // Rebuild a region of interest for chamber 7. | |
649 | // Remember the parameters a and b are station specific. | |
b92524d0 | 650 | AliHLTMUONRecHitStruct p7 = data->fTag.fLine.FindIntersectWithXYPlain( fgZ7 ); |
e0679afb | 651 | data->fTag.fChamber = kChamber7; |
652 | data->fTag.fRoi.Create(p7, fgA7, fgB7); | |
653 | ||
654 | AliHLTFloat32_t left, right, bottom, top; | |
655 | data->fTag.fRoi.GetBoundaryBox(left, right, bottom, top); | |
656 | // Make request for chamber 7 data. | |
657 | RequestClusters(left, right, bottom, top, kChamber7, &data->fTag); | |
658 | } | |
659 | } | |
660 | break; | |
661 | ||
662 | case kWaitMoreChamber8: | |
663 | if (fRequestsCompleted == fSt5data.Count() && fSm5state == kSM5Done) | |
664 | ProcessClusters(); | |
665 | break; | |
666 | ||
667 | default: | |
668 | DebugTrace("ERROR: Unexpected state for SM4 in AliHLTMUONMansoTrackerFSM::EndOfClustersChamber8!"); | |
669 | } | |
670 | } | |
671 | ||
672 | ||
673 | void AliHLTMUONMansoTrackerFSM::EndOfClustersChamber9() | |
674 | { | |
675 | // State change method for Station 5 state machine. | |
676 | ||
677 | switch (fSm5state) | |
678 | { | |
679 | case kWaitChamber9: | |
680 | fSm5state = kSM5Done; | |
681 | EndOfClusterRequests(); | |
682 | NoTrackFound(); | |
683 | break; | |
684 | ||
685 | case kWaitMoreChamber9: | |
686 | fSm5state = kSM5Done; | |
687 | EndOfClusterRequests(); | |
688 | if (fRequestsCompleted == fSt5data.Count()) | |
689 | ProcessClusters(); | |
690 | break; | |
691 | ||
692 | default: | |
693 | DebugTrace("ERROR: Unexpected state for SM5 in AliHLTMUONMansoTrackerFSM::EndOfClustersChamber9!"); | |
694 | } | |
695 | } | |
696 | ||
697 | ||
698 | void AliHLTMUONMansoTrackerFSM::EndOfClustersChamber10() | |
699 | { | |
700 | // State change method for Station 5 state machine. | |
701 | ||
702 | switch (fSm5state) | |
703 | { | |
704 | case kWaitChamber10: | |
705 | { | |
706 | fSm5state = kWaitChamber9; | |
707 | fSt5z = fgZ9; | |
708 | ||
709 | // No clusters found on chamber 10 so we need to make a request for | |
710 | // clusters from chamber 9: | |
b92524d0 | 711 | AliHLTMUONRecHitStruct p9 = fMc1.fLine.FindIntersectWithXYPlain( fgZ9 ); |
e0679afb | 712 | |
713 | // Build a region of interest for tracking station 5 (chamber 9). | |
714 | // Remember the parameters a and b are station specific. | |
715 | fMc1.fChamber = kChamber9; | |
716 | fMc1.fRoi.Create(p9, fgA9, fgB9); | |
717 | ||
718 | AliHLTFloat32_t left, right, bottom, top; | |
719 | fMc1.fRoi.GetBoundaryBox(left, right, bottom, top); | |
720 | RequestClusters(left, right, bottom, top, kChamber9, &fMc1); | |
721 | } | |
722 | break; | |
723 | ||
724 | case kWaitMoreChamber10: | |
725 | fSm5state = kSM5Done; | |
726 | EndOfClusterRequests(); | |
727 | if (fRequestsCompleted == fSt5data.Count()) | |
728 | ProcessClusters(); | |
729 | break; | |
730 | ||
731 | default: | |
732 | DebugTrace("ERROR: Unexpected state for SM5 in AliHLTMUONMansoTrackerFSM::EndOfClustersChamber10!"); | |
733 | } | |
734 | } | |
735 | ||
736 | ||
b92524d0 | 737 | void AliHLTMUONMansoTrackerFSM::ProjectToStation4( |
738 | AliStation5Data* data, register AliHLTFloat32_t station5z | |
739 | ) | |
e0679afb | 740 | { |
741 | // Perform chamber specific operations: | |
742 | // Since certain states of SM4 means that it is fetching for Chamber8 | |
743 | // and other states are for fetching from Chamber7. We need to make | |
744 | // requests for the correct chamber. | |
b92524d0 | 745 | assert( fSm4state == kWaitChamber8 |
e0679afb | 746 | || fSm4state == kWaitMoreChamber8 |
747 | || fSm4state == kWaitChamber7 | |
748 | || fSm4state == kWaitMoreChamber7 | |
749 | ); | |
750 | AliTagData* tag = &data->fTag; | |
751 | if (fSm4state == kWaitChamber8 || fSm4state == kWaitMoreChamber8) | |
752 | { | |
753 | // Form the vector line between trigger station 1 and tracking station 5, | |
754 | // and find the intersection point of the line with station 4 (chamber8). | |
755 | AliLine line51( AliVertex(data->fClusterPoint, station5z), fV1 ); | |
b92524d0 | 756 | AliHLTMUONRecHitStruct intercept = line51.FindIntersectWithXYPlain( fgZ8 ); |
e0679afb | 757 | tag->fLine = line51; |
758 | ||
759 | // Build a region of interest for tracking station 4. | |
760 | tag->fChamber = kChamber8; | |
761 | tag->fRoi.Create(intercept, fgA8, fgB8); | |
762 | } | |
763 | else | |
764 | { | |
765 | // Form the vector line between trigger station 1 and tracking station 5, | |
766 | // and find the intersection point of the line with station 4 (chamber7). | |
767 | AliLine line51( AliVertex(data->fClusterPoint, station5z), fV1 ); | |
b92524d0 | 768 | AliHLTMUONRecHitStruct intercept = line51.FindIntersectWithXYPlain( fgZ7 ); |
e0679afb | 769 | tag->fLine = line51; |
770 | ||
771 | // Build a region of interest for tracking station 4. | |
772 | tag->fChamber = kChamber7; | |
773 | tag->fRoi.Create(intercept, fgA7, fgB7); | |
774 | } | |
775 | ||
776 | // Make the request for clusters from station 4. | |
777 | AliHLTFloat32_t left, right, bottom, top; | |
778 | tag->fRoi.GetBoundaryBox(left, right, bottom, top); | |
779 | RequestClusters(left, right, bottom, top, tag->fChamber, tag); | |
780 | } | |
781 | ||
782 | ||
783 | void AliHLTMUONMansoTrackerFSM::ProcessClusters() | |
784 | { | |
785 | // Process clusters that have been received. | |
786 | // This is called once all clusters have been found. | |
787 | ||
788 | DebugTrace("ProcessClusters..."); | |
789 | ||
790 | // Check if the cluster point list on station 4 is empty. | |
791 | // If it is then we have not found any tracks. | |
792 | fFoundPoint = fSt4points.First(); | |
793 | if (fFoundPoint == fSt4points.End()) | |
794 | { | |
795 | NoTrackFound(); | |
796 | return; | |
797 | } | |
798 | ||
799 | fSt5rec = fSt5data.First(); | |
800 | if (fSt5rec != fSt5data.End()) | |
801 | { | |
802 | // Only look at station 5 data records that are for the found chamber number. | |
803 | // Note: either we only have chamber 8 data or we have chamber 7 data followed | |
804 | // by chamber 8 data. | |
805 | // Thus if we hit records that we are not interested in already then the list | |
806 | // contains no interesting data and we can signal no track found. | |
807 | if (fSt5rec->fTag.fChamber != fSt4chamber) | |
808 | { | |
809 | NoTrackFound(); | |
810 | return; | |
811 | } | |
812 | ||
813 | // For all combinations of cluster point pairs from station 4 and 5 | |
814 | // signal a found track: | |
815 | do | |
816 | { | |
817 | DebugTrace("\tfSt5rec->fTag.chamber = " << fSt5rec->fTag.fChamber | |
818 | << " , fSt4chamber = " << fSt4chamber | |
819 | ); | |
820 | ||
821 | for (fFoundPoint = fSt4points.First(); fFoundPoint != fSt4points.End(); fFoundPoint++) | |
822 | { | |
823 | if (fFoundPoint->fSt5tag == &fSt5rec->fTag) | |
824 | FoundTrack(); | |
825 | } | |
826 | ||
827 | fSt5rec++; // Get next station 5 cluster point. | |
828 | } while (fSt5rec != fSt5data.End() && fSt5rec->fTag.fChamber == fSt4chamber); | |
829 | } | |
830 | else | |
831 | NoTrackFound(); | |
832 | } | |
833 |