]>
Commit | Line | Data |
---|---|---|
a3283a4c | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
3d1463c8 | 16 | // $Id$ |
17 | ||
18 | //----------------------------------------------------------------------------- | |
9265505b | 19 | /// \class AliMUONDigitMaker |
8c51a32f | 20 | /// MUON Digit maker from rawdata. |
a3283a4c | 21 | /// |
22 | /// Raw2Digits: | |
23 | /// Using real mapping for tracker | |
24 | /// Indranil Das (Adapted for runloader: Ch. Finck) july 05 | |
8c51a32f | 25 | /// |
26 | /// Implemented non-constant buspatch numbers for tracking | |
27 | /// with correct DDL id. | |
28 | /// (Ch. Finck, dec 05) | |
29 | /// | |
a3283a4c | 30 | /// Add reader for scaler trigger events |
31 | /// Use memcpy instead of assignment elt by elt | |
32 | /// (Ch. Finck, Jan 06) | |
8c51a32f | 33 | /// |
34 | /// Using new interface with AliMUONRawStreamTracker(Trigger) | |
35 | /// (New interface of AliMUONRawReader class) | |
36 | /// (further details could be found in Alice-note) | |
37 | /// (Ch. Finck, March 06) | |
38 | /// | |
39 | /// Add (S)Digit maker tracker (for free) | |
40 | /// and for trigger. Create trigger inverse mapping. | |
78649106 | 41 | /// |
42 | /// \author Ch. Finck, oct 06 | |
3d1463c8 | 43 | //----------------------------------------------------------------------------- |
a3283a4c | 44 | |
a3283a4c | 45 | #include "AliMUONDigitMaker.h" |
a3283a4c | 46 | |
40e382ae | 47 | #include "AliLog.h" |
a3283a4c | 48 | #include "AliMUONDDLTrigger.h" |
49 | #include "AliMUONDarcHeader.h" | |
40e382ae | 50 | #include "AliMUONVDigit.h" |
51 | #include "AliMUONVDigitStore.h" | |
52 | #include "AliMUONGlobalTrigger.h" | |
a3283a4c | 53 | #include "AliMUONLocalStruct.h" |
a3283a4c | 54 | #include "AliMUONLocalTrigger.h" |
40e382ae | 55 | #include "AliMUONRawStreamTracker.h" |
e3a2b9c9 | 56 | #include "AliMUONRawStreamTrackerHP.h" |
40e382ae | 57 | #include "AliMUONRawStreamTrigger.h" |
58 | #include "AliMUONRegHeader.h" | |
e1a10d41 | 59 | #include "AliMUONTriggerCircuit.h" |
d9b1ecb0 | 60 | #include "AliMpTriggerCrate.h" |
61 | #include "AliMpLocalBoard.h" | |
40e382ae | 62 | #include "AliMUONVTriggerStore.h" |
63 | #include "AliMpCathodType.h" | |
64 | #include "AliMpDDLStore.h" | |
65 | #include "AliMpDEManager.h" | |
66 | #include "AliMpPad.h" | |
9265505b | 67 | #include "AliMpSegmentation.h" |
a3283a4c | 68 | #include "AliMpVSegmentation.h" |
9265505b | 69 | #include "AliRawReader.h" |
7771752e | 70 | #include <TArrayS.h> |
9265505b | 71 | |
9265505b | 72 | /// \cond CLASSIMP |
a3283a4c | 73 | ClassImp(AliMUONDigitMaker) // Class implementation in ROOT context |
9265505b | 74 | /// \endcond |
75 | ||
a3283a4c | 76 | //__________________________________________________________________________ |
e3a2b9c9 | 77 | AliMUONDigitMaker::AliMUONDigitMaker(Bool_t enableErrorLogger, Bool_t useFastDecoder) |
a3283a4c | 78 | : TObject(), |
9f5dcca3 | 79 | fScalerEvent(kFALSE), |
40e382ae | 80 | fMakeTriggerDigits(kFALSE), |
e3a2b9c9 | 81 | fRawStreamTracker(NULL), |
9f5dcca3 | 82 | fRawStreamTrigger(new AliMUONRawStreamTrigger()), |
9f5dcca3 | 83 | fTrackerTimer(), |
84 | fTriggerTimer(), | |
40e382ae | 85 | fMappingTimer(), |
86 | fDigitStore(0x0), | |
87 | fTriggerStore(0x0) | |
a3283a4c | 88 | { |
40e382ae | 89 | /// ctor |
a3283a4c | 90 | |
91 | AliDebug(1,""); | |
e3a2b9c9 | 92 | |
d61f56df | 93 | CreateRawStreamTracker(useFastDecoder); |
a3283a4c | 94 | |
95 | // Standard Constructor | |
3c7f5307 | 96 | if (enableErrorLogger) { |
97 | fRawStreamTracker->EnabbleErrorLogger(); | |
98 | fRawStreamTrigger->EnabbleErrorLogger(); | |
99 | } | |
a3283a4c | 100 | |
a3283a4c | 101 | fTrackerTimer.Start(kTRUE); fTrackerTimer.Stop(); |
102 | fTriggerTimer.Start(kTRUE); fTriggerTimer.Stop(); | |
103 | fMappingTimer.Start(kTRUE); fMappingTimer.Stop(); | |
40e382ae | 104 | |
105 | SetMakeTriggerDigits(); | |
a3283a4c | 106 | |
107 | } | |
108 | ||
a3283a4c | 109 | //__________________________________________________________________________ |
9f5dcca3 | 110 | AliMUONDigitMaker::~AliMUONDigitMaker() |
a3283a4c | 111 | { |
9265505b | 112 | /// clean up |
113 | /// and time processing measure | |
114 | ||
a3283a4c | 115 | delete fRawStreamTracker; |
116 | delete fRawStreamTrigger; | |
117 | ||
a7b01aa5 | 118 | AliDebug(1, Form("Execution time for MUON tracker : R:%.2fs C:%.2fs", |
a3283a4c | 119 | fTrackerTimer.RealTime(),fTrackerTimer.CpuTime())); |
a7b01aa5 | 120 | AliDebug(1, Form(" Execution time for MUON tracker (mapping calls part) " |
a3283a4c | 121 | ": R:%.2fs C:%.2fs", |
122 | fMappingTimer.RealTime(),fMappingTimer.CpuTime())); | |
a7b01aa5 | 123 | AliDebug(1, Form("Execution time for MUON trigger : R:%.2fs C:%.2fs", |
a3283a4c | 124 | fTriggerTimer.RealTime(),fTriggerTimer.CpuTime())); |
125 | ||
a3283a4c | 126 | } |
127 | ||
d61f56df | 128 | //__________________________________________________________________________ |
129 | void AliMUONDigitMaker::CreateRawStreamTracker(Bool_t useFastDecoder) | |
130 | { | |
131 | /// Create raw stream tracker according to the passed option | |
132 | ||
133 | if (useFastDecoder) | |
134 | { | |
135 | AliInfo("Using fast decoder."); | |
136 | fRawStreamTracker = new AliMUONRawStreamTrackerHP(); | |
137 | } | |
138 | else | |
139 | fRawStreamTracker = new AliMUONRawStreamTracker(); | |
140 | } | |
141 | ||
a3283a4c | 142 | //____________________________________________________________________ |
40e382ae | 143 | Int_t AliMUONDigitMaker::Raw2Digits(AliRawReader* rawReader, |
144 | AliMUONVDigitStore* digitStore, | |
145 | AliMUONVTriggerStore* triggerStore) | |
a3283a4c | 146 | { |
9265505b | 147 | /// Main method to creates digit |
148 | /// for tracker | |
149 | /// and trigger | |
a3283a4c | 150 | |
40e382ae | 151 | AliDebug(1,Form("rawReader=%p digitStore=%p triggerStore=%p", |
152 | rawReader,digitStore,triggerStore)); | |
153 | ||
154 | fDigitStore = digitStore; | |
155 | fTriggerStore = triggerStore; | |
156 | ||
157 | if (!fDigitStore && !fTriggerStore) | |
158 | { | |
159 | AliError("No digit or trigger store given. Nothing to do..."); | |
160 | return kFALSE; | |
161 | } | |
162 | ||
163 | if ( fDigitStore ) | |
164 | { | |
165 | fDigitStore->Clear(); // insure we start with an empty container | |
166 | ReadTrackerDDL(rawReader); | |
167 | } | |
168 | ||
169 | if ( fTriggerStore || fMakeTriggerDigits ) | |
170 | { | |
171 | if ( fTriggerStore ) fTriggerStore->Clear(); | |
172 | if ( fMakeTriggerDigits && !fDigitStore ) | |
173 | { | |
174 | AliError("Asking for trigger digits but digitStore is null"); | |
175 | } | |
176 | else | |
177 | { | |
e604ac4c | 178 | ReadTriggerDDL(rawReader); |
40e382ae | 179 | } |
180 | } | |
181 | ||
a3283a4c | 182 | return kTRUE; |
a3283a4c | 183 | } |
184 | ||
185 | //____________________________________________________________________ | |
186 | Int_t AliMUONDigitMaker::ReadTrackerDDL(AliRawReader* rawReader) | |
187 | { | |
40e382ae | 188 | /// Reading tracker DDL |
189 | /// filling the fDigitStore container, which must not be null | |
a3283a4c | 190 | |
40e382ae | 191 | AliDebug(1,""); |
192 | ||
a3283a4c | 193 | fTrackerTimer.Start(kFALSE); |
194 | ||
195 | // elex info | |
196 | Int_t buspatchId; | |
197 | UChar_t channelId; | |
198 | UShort_t manuId; | |
a3283a4c | 199 | UShort_t charge; |
a3283a4c | 200 | |
a3283a4c | 201 | fRawStreamTracker->SetReader(rawReader); |
1bc885f3 | 202 | fRawStreamTracker->First(); |
203 | ||
204 | while ( fRawStreamTracker->Next(buspatchId,manuId,channelId,charge) ) | |
205 | { | |
40e382ae | 206 | // getting DE from buspatch |
207 | Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromBus(buspatchId); | |
a3283a4c | 208 | |
40e382ae | 209 | const AliMpVSegmentation* seg |
210 | = AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId, | |
211 | manuId); | |
a3283a4c | 212 | |
40e382ae | 213 | AliMp::CathodType cathodeType = AliMpDEManager::GetCathod(detElemId, |
214 | seg->PlaneType()); | |
1bc885f3 | 215 | |
40e382ae | 216 | AliMpPad pad = seg->PadByLocation(AliMpIntPair(manuId,channelId),kFALSE); |
217 | ||
218 | if (!pad.IsValid()) | |
219 | { | |
220 | AliError(Form("No pad for detElemId: %d, manuId: %d, channelId: %d", | |
221 | detElemId, manuId, channelId)); | |
222 | continue; | |
223 | } | |
224 | ||
225 | AliMUONVDigit* digit = fDigitStore->Add(detElemId,manuId,channelId,cathodeType, | |
226 | AliMUONVDigitStore::kDeny); | |
227 | if (!digit) | |
228 | { | |
229 | AliError(Form("Digit DE %04d Manu %04d Channel %02d could not be added", | |
230 | detElemId, manuId, channelId)); | |
231 | continue; | |
232 | } | |
233 | ||
234 | digit->SetPadXY(pad.GetIndices().GetFirst(), | |
235 | pad.GetIndices().GetSecond()); | |
236 | ||
237 | digit->SetADC(charge); | |
9265505b | 238 | |
40e382ae | 239 | } |
a3283a4c | 240 | |
40e382ae | 241 | fTrackerTimer.Stop(); |
a3283a4c | 242 | |
40e382ae | 243 | return kTRUE; |
a3283a4c | 244 | } |
245 | ||
246 | //____________________________________________________________________ | |
247 | Int_t AliMUONDigitMaker::ReadTriggerDDL(AliRawReader* rawReader) | |
248 | { | |
9265505b | 249 | /// reading tracker DDL |
40e382ae | 250 | /// filling the fTriggerStore container, which must not be null |
a3283a4c | 251 | |
40e382ae | 252 | AliDebug(1,""); |
253 | ||
a3283a4c | 254 | AliMUONDDLTrigger* ddlTrigger = 0x0; |
255 | AliMUONDarcHeader* darcHeader = 0x0; | |
256 | AliMUONRegHeader* regHeader = 0x0; | |
257 | AliMUONLocalStruct* localStruct = 0x0; | |
258 | ||
259 | Int_t loCircuit; | |
d2d759cf | 260 | |
a3283a4c | 261 | fTriggerTimer.Start(kFALSE); |
262 | ||
263 | fRawStreamTrigger->SetReader(rawReader); | |
264 | ||
40e382ae | 265 | while (fRawStreamTrigger->NextDDL()) |
266 | { | |
a3283a4c | 267 | ddlTrigger = fRawStreamTrigger->GetDDLTrigger(); |
268 | darcHeader = ddlTrigger->GetDarcHeader(); | |
40e382ae | 269 | |
270 | // fill global trigger information | |
271 | if (fTriggerStore) | |
272 | { | |
273 | if (darcHeader->GetGlobalFlag()) | |
274 | { | |
275 | AliMUONGlobalTrigger globalTrigger; | |
276 | globalTrigger.SetFromGlobalResponse(darcHeader->GetGlobalOutput()); | |
277 | fTriggerStore->SetGlobal(globalTrigger); | |
d2d759cf | 278 | } |
a3283a4c | 279 | } |
40e382ae | 280 | |
a3283a4c | 281 | Int_t nReg = darcHeader->GetRegHeaderEntries(); |
40e382ae | 282 | |
283 | for(Int_t iReg = 0; iReg < nReg ;iReg++) | |
284 | { //reg loop | |
285 | ||
d9b1ecb0 | 286 | |
287 | // crate info | |
288 | AliMpTriggerCrate* crate = AliMpDDLStore::Instance()-> | |
289 | GetTriggerCrate(fRawStreamTrigger->GetDDL(), iReg); | |
40e382ae | 290 | |
d53fb0de | 291 | if (!crate) |
40e382ae | 292 | AliWarning(Form("Missing crate number %d in DDL %d\n", iReg, fRawStreamTrigger->GetDDL())); |
d9b1ecb0 | 293 | |
40e382ae | 294 | |
a3283a4c | 295 | regHeader = darcHeader->GetRegHeaderEntry(iReg); |
40e382ae | 296 | |
a3283a4c | 297 | Int_t nLocal = regHeader->GetLocalEntries(); |
40e382ae | 298 | for(Int_t iLocal = 0; iLocal < nLocal; iLocal++) |
299 | { | |
300 | ||
301 | localStruct = regHeader->GetLocalEntry(iLocal); | |
302 | ||
303 | // if card exist | |
304 | if (localStruct) { | |
305 | ||
d9b1ecb0 | 306 | loCircuit = crate->GetLocalBoardId(localStruct->GetId()); |
307 | ||
308 | if ( !loCircuit ) continue; // empty slot | |
309 | ||
310 | AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(loCircuit, false); | |
311 | ||
312 | // skip copy cards | |
313 | if( !localBoard->IsNotified()) | |
314 | continue; | |
40e382ae | 315 | |
316 | if (fTriggerStore) | |
317 | { | |
318 | // fill local trigger | |
319 | AliMUONLocalTrigger localTrigger; | |
320 | localTrigger.SetLocalStruct(loCircuit, *localStruct); | |
321 | fTriggerStore->Add(localTrigger); | |
322 | } | |
323 | ||
324 | if ( fMakeTriggerDigits ) | |
325 | { | |
326 | //FIXEME should find something better than a TArray | |
327 | TArrayS xyPattern[2]; | |
40e382ae | 328 | |
d9b1ecb0 | 329 | localStruct->GetXPattern(xyPattern[0]); |
330 | localStruct->GetYPattern(xyPattern[1]); | |
40e382ae | 331 | |
332 | TriggerDigits(loCircuit, xyPattern, *fDigitStore); | |
333 | } | |
334 | } // if triggerY | |
a3283a4c | 335 | } // iLocal |
336 | } // iReg | |
337 | } // NextDDL | |
40e382ae | 338 | |
a3283a4c | 339 | fTriggerTimer.Stop(); |
340 | ||
341 | return kTRUE; | |
342 | ||
343 | } | |
7771752e | 344 | |
d2d759cf | 345 | //____________________________________________________________________ |
7771752e | 346 | Int_t AliMUONDigitMaker::TriggerDigits(Int_t nBoard, |
40e382ae | 347 | TArrayS* xyPattern, |
348 | AliMUONVDigitStore& digitStore) const | |
d2d759cf | 349 | { |
40e382ae | 350 | /// make digits for trigger from pattern, and add them to digitStore |
d2d759cf | 351 | |
d9b1ecb0 | 352 | Int_t detElemId; |
353 | ||
354 | AliMpLocalBoard* localBoard = AliMpDDLStore::Instance()->GetLocalBoard(nBoard); | |
355 | ||
8c0b5e70 | 356 | Int_t n,b; |
357 | ||
d2d759cf | 358 | // loop over x1-4 and y1-4 |
40e382ae | 359 | for (Int_t iChamber = 0; iChamber < 4; ++iChamber) |
360 | { | |
361 | for (Int_t iCath = 0; iCath < 2; ++iCath) | |
362 | { | |
7771752e | 363 | Int_t pattern = (Int_t)xyPattern[iCath].At(iChamber); |
364 | if (!pattern) continue; | |
40e382ae | 365 | |
7771752e | 366 | // get detElemId |
d9b1ecb0 | 367 | detElemId = AliMpDDLStore::Instance()->GetDEfromLocalBoard(nBoard, iChamber); |
40e382ae | 368 | |
369 | const AliMpVSegmentation* seg | |
370 | = AliMpSegmentation::Instance() | |
371 | ->GetMpSegmentation(detElemId, AliMp::GetCathodType(iCath)); | |
372 | ||
373 | // loop over the 16 bits of pattern | |
374 | for (Int_t ibitxy = 0; ibitxy < 16; ++ibitxy) | |
375 | { | |
376 | if ((pattern >> ibitxy) & 0x1) | |
377 | { | |
378 | // not quite sure about this | |
379 | Int_t offset = 0; | |
380 | if (iCath && localBoard->GetSwitch(6)) offset = -8; | |
381 | ||
382 | AliMpPad pad = seg->PadByLocation(AliMpIntPair(nBoard,ibitxy+offset),kTRUE); | |
383 | ||
384 | if (!pad.IsValid()) | |
385 | { | |
386 | AliWarning(Form("No pad for detElemId: %d, nboard %d, ibitxy: %d\n", | |
387 | detElemId, nBoard, ibitxy)); | |
388 | continue; | |
389 | } | |
8c0b5e70 | 390 | |
391 | n = pad.GetLocation(0).GetFirst(); // always take first location so that digits are not inserted several times | |
392 | b = pad.GetLocation(0).GetSecond(); | |
393 | ||
394 | AliDebug(1,Form("Using localBoard %d ixy %d instead of %d,%d", | |
395 | n,b,nBoard,ibitxy)); | |
396 | ||
397 | AliMUONVDigit* digit = digitStore.Add(detElemId,n,b,iCath,AliMUONVDigitStore::kDeny); | |
40e382ae | 398 | |
399 | if (!digit) | |
400 | { | |
8c0b5e70 | 401 | AliDebug(1, Form("Digit DE %04d LocalBoard %03d ibitxy %02d cath %d already in store", |
402 | detElemId,nBoard,ibitxy,iCath)); | |
403 | continue; | |
40e382ae | 404 | } |
405 | ||
406 | Int_t padX = pad.GetIndices().GetFirst(); | |
407 | Int_t padY = pad.GetIndices().GetSecond(); | |
408 | ||
409 | // fill digit | |
410 | digit->SetPadXY(padX,padY); | |
411 | digit->SetCharge(1.); | |
412 | }// xyPattern | |
413 | }// ibitxy | |
7771752e | 414 | }// cath |
415 | } // ichamber | |
40e382ae | 416 | |
d2d759cf | 417 | return kTRUE; |
418 | } | |
d61f56df | 419 | |
420 | //____________________________________________________________________ | |
421 | void AliMUONDigitMaker::SetFastDecoder(Bool_t useFastDecoder) | |
422 | { | |
423 | /// Set fast raw data decoder | |
424 | ||
425 | delete fRawStreamTracker; | |
426 | CreateRawStreamTracker(useFastDecoder); | |
427 | } | |
428 | ||
429 | ||
430 |