]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TOF/AliTOFRawStream.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Decode method updated: now it addes the CDH argument and sets the deltaBunchID per...
[u/mrichter/AliRoot.git] / TOF / AliTOFRawStream.cxx
CommitLineData
571dda3d 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
16/*
d0eb8f39 17$Log$
96f01799 18Revision 1.19 2007/05/18 13:07:53 decaro
19Error messages stored in the global raw-reader error log (Cvetan, Chiara)
20
019f6101 21Revision 1.18 2007/05/08 11:53:29 arcelli
22Improved class flexibility for further use (R.Preghenella)
23
2fd622fa 24Revision 1.17 2007/05/03 08:53:50 decaro
25Coding convention: RS3 violation -> suppression
26
19eefa2b 27Revision 1.16 2007/05/03 08:22:22 decaro
28Coding convention: RN17 violation -> suppression
29
dfcee998 30Revision 1.15 2007/04/30 15:22:06 arcelli
31Change TOF digit Time, Tot etc to int type
32
bf33f8f0 33Revision 1.14 2007/04/27 11:11:53 arcelli
34updates for the new decoder
35
dfcee998 36Revision 1.13 2007/03/16 11:46:35 decaro
37Coding convention: RN17 rule violation -> suppression
38
8994284f 39Revision 1.12 2007/02/22 09:43:45 decaro
40Added AliTOFRawStream::GetIndex method for online calibration (C.Zampolli)
41
740f0695 42Revision 1.11 2007/02/20 15:57:00 decaro
43Raw data update: to read the TOF raw data defined in UNPACKED mode
44
15ec34b9 45Revision 1.10 2006/12/15 14:01:38 cvetan
46Memory leak fixed
47
f1aa96cc 48Revision 1.9 2006/10/13 11:22:27 arcelli
49remove warnings due to uninitialized AliTOFtdcDigit data members
50
36f66fad 51Revision 1.8 2006/08/22 13:30:17 arcelli
52removal of effective c++ warnings (C.Zampolli)
53
655e379f 54Revision 1.7 2006/08/10 14:46:54 decaro
55TOF raw data format: updated version
56
d0eb8f39 57Revision 1.6.1 2006/06/28 A. De Caro, R. Preghenella:
58 Update TOF raw data format
59 according to the final version
60 (see the ALICE internal note in preparation
61 'ALICE TOF raw data format')
62 Added the methods for the correspoonding numbering
63 between the equipment IDs and the volume IDs:
64 Equip2VolNPlate(...)
65 Equip2VolNStrip(...)
66 Equip2VolNPad(...)
67
571dda3d 68Revision 0.02 2005/07/28 A. De Caro:
69 Update format TOF raw data
70 (temporary solution)
71 Correction of few wrong corrispondences
72 between 'software' and 'hardware' numberings
73
74Revision 0.01 2005/07/22 A. De Caro
75 Implement methods Next()
76 GetSector(),
77 GetPlate(),
78 GetStrip(),
79 GetPadZ(),
80 GetPadX()
81*/
82
0e46b9ae 83////////////////////////////////////////////////////////////////////////
84// //
85// This class provides access to TOF raw data in DDL files. //
86// //
87// It loops over all TOF raw data given by the AliRawReader. //
88// //
89////////////////////////////////////////////////////////////////////////
571dda3d 90
15ec34b9 91
13d6bff9 92#include "Riostream.h"
93
15ec34b9 94#include "TClonesArray.h"
95
96#include "AliDAQ.h"
d0eb8f39 97#include "AliLog.h"
571dda3d 98#include "AliRawReader.h"
99
100#include "AliTOFGeometry.h"
15ec34b9 101#include "AliTOFrawData.h"
102#include "AliTOFRawMap.h"
571dda3d 103#include "AliTOFRawStream.h"
32ead898 104//#include "AliTOFCableLengthMap.h"
571dda3d 105
5c7c93fa 106#include "AliTOFHitData.h"
d0eb8f39 107
32ead898 108#include "AliRawEventHeaderBase.h"
109
571dda3d 110ClassImp(AliTOFRawStream)
111
32ead898 112const Int_t AliTOFRawStream::fgkddlBCshift[72] =
113{
114 2, 2, -1, -1,
cfceaf15 115 2, 2, 0, 0,
32ead898 116 2, 2, 0, 0,
117 2, 2, 0, 0,
118 2, 2, 0, 0,
119 2, 2, 0, 0,
120 2, 2, 0, 0,
121 2, 2, 0, 0,
122 2, 2, 0, 0,
123 2, 2, 0, 0,
124 2, 2, -1, -1,
125 2, 2, -1, -1,
126 2, 2, -2, -2,
127 2, 2, -2, -2,
128 2, 2, -2, -2,
129 2, 2, -1, -1,
130 2, 2, -1, -1,
131 2, 2, -1, -1
132};
133
571dda3d 134
135//_____________________________________________________________________________
655e379f 136AliTOFRawStream::AliTOFRawStream(AliRawReader* rawReader):
15ec34b9 137 fRawReader(rawReader),
7813bb1a 138 fTOFrawData(0x0),
139 fDecoder(new AliTOFDecoder()),
655e379f 140 fDDL(-1),
141 fTRM(-1),
655e379f 142 fTRMchain(-1),
15ec34b9 143 fTDC(-1),
655e379f 144 fTDCchannel(-1),
15ec34b9 145 fTime(-1),
655e379f 146 fToT(-1),
15ec34b9 147 fLeadingEdge(-1),
148 fTrailingEdge(-1),
655e379f 149 fErrorFlag(-1),
655e379f 150 fSector(-1),
151 fPlate(-1),
152 fStrip(-1),
153 fPadX(-1),
154 fPadZ(-1),
15ec34b9 155 fPackedDigits(0),
655e379f 156 fWordType(-1),
157 fSlotID(-1),
158 fACQ(-1),
159 fPSbit(-1),
655e379f 160 fTDCerrorFlag(-1),
161 fInsideDRM(kFALSE),
162 fInsideTRM(kFALSE),
163 fInsideLTM(kFALSE),
164 fInsideTRMchain0(kFALSE),
38f24f3f 165 fInsideTRMchain1(kFALSE),
166 fLocalEventCounterDRM(-1),
167 fLocalEventCounterLTM(-1),
168 fLocalEventCounterTRM(0x0),
32ead898 169 fLocalEventCounterChain(0x0),
170 fCableLengthMap(0x0),
171 fEventID(0)
571dda3d 172{
173 //
174 // create an object to read TOF raw digits
175 //
176
7813bb1a 177 for (Int_t i=0;i<AliDAQ::NumberOfDdls("TOF");i++){
2fd622fa 178 fDataBuffer[i]=new AliTOFHitDataBuffer();
179 fPackedDataBuffer[i]=new AliTOFHitDataBuffer();
7813bb1a 180 }
181
182 fTOFrawData = new TClonesArray("AliTOFrawData",1000);
183 fTOFrawData->SetOwner();
184
15ec34b9 185 fRawReader->Reset();
362c9d61 186 fRawReader->Select("TOF");
38f24f3f 187
188 fLocalEventCounterTRM = new Int_t[13];//adc
189 fLocalEventCounterChain = new Int_t*[13];//adc
190 for (Int_t j=0;j<13;j++){//adc
191 fLocalEventCounterTRM[j] = -1;//adc
192 fLocalEventCounterChain[j] = new Int_t[2];//adc
193 for (Int_t k=0;k<2;k++){//adc
194 fLocalEventCounterChain[j][k] = -1;//adc
195 }//adc
196 }//adc
197
32ead898 198 fCableLengthMap = new AliTOFCableLengthMap();
199
0a8999ac 200 fEventID = fRawReader->GetBCID(); //bunch crossing
571dda3d 201}
202
203//_____________________________________________________________________________
655e379f 204AliTOFRawStream::AliTOFRawStream():
7813bb1a 205 fRawReader(0x0),
7813bb1a 206 fTOFrawData(0x0),
207 fDecoder(new AliTOFDecoder()),
655e379f 208 fDDL(-1),
209 fTRM(-1),
655e379f 210 fTRMchain(-1),
15ec34b9 211 fTDC(-1),
655e379f 212 fTDCchannel(-1),
15ec34b9 213 fTime(-1),
655e379f 214 fToT(-1),
15ec34b9 215 fLeadingEdge(-1),
216 fTrailingEdge(-1),
655e379f 217 fErrorFlag(-1),
655e379f 218 fSector(-1),
219 fPlate(-1),
220 fStrip(-1),
221 fPadX(-1),
222 fPadZ(-1),
15ec34b9 223 fPackedDigits(0),
655e379f 224 fWordType(-1),
225 fSlotID(-1),
226 fACQ(-1),
227 fPSbit(-1),
655e379f 228 fTDCerrorFlag(-1),
229 fInsideDRM(kFALSE),
230 fInsideTRM(kFALSE),
231 fInsideLTM(kFALSE),
232 fInsideTRMchain0(kFALSE),
38f24f3f 233 fInsideTRMchain1(kFALSE),
234 fLocalEventCounterDRM(-1),
235 fLocalEventCounterLTM(-1),
236 fLocalEventCounterTRM(0x0),
32ead898 237 fLocalEventCounterChain(0x0),
238 fCableLengthMap(0x0),
239 fEventID(0)
571dda3d 240{
241 //
d0eb8f39 242 // default ctr
571dda3d 243 //
7813bb1a 244 for (Int_t i=0;i<AliDAQ::NumberOfDdls("TOF");i++){
2fd622fa 245 fDataBuffer[i]=new AliTOFHitDataBuffer();
246 fPackedDataBuffer[i]=new AliTOFHitDataBuffer();
7813bb1a 247 }
15ec34b9 248
7813bb1a 249 fTOFrawData = new TClonesArray("AliTOFrawData",1000);
250 fTOFrawData->SetOwner();
38f24f3f 251
252 fLocalEventCounterTRM = new Int_t[13];//adc
253 fLocalEventCounterChain = new Int_t*[13];//adc
254 for (Int_t j=0;j<13;j++){//adc
255 fLocalEventCounterTRM[j] = -1;//adc
256 fLocalEventCounterChain[j] = new Int_t[2];//adc
257 for (Int_t k=0;k<2;k++){//adc
258 fLocalEventCounterChain[j][k] = -1;//adc
259 }//adc
260 }//adc
261
32ead898 262 fCableLengthMap = new AliTOFCableLengthMap();
263
d0eb8f39 264}
265
266//_____________________________________________________________________________
267AliTOFRawStream::AliTOFRawStream(const AliTOFRawStream& stream) :
655e379f 268 TObject(stream),
269 fRawReader(0x0),
7813bb1a 270 fTOFrawData(0x0),
271 fDecoder(new AliTOFDecoder()),
655e379f 272 fDDL(-1),
273 fTRM(-1),
655e379f 274 fTRMchain(-1),
15ec34b9 275 fTDC(-1),
655e379f 276 fTDCchannel(-1),
15ec34b9 277 fTime(-1),
655e379f 278 fToT(-1),
15ec34b9 279 fLeadingEdge(-1),
280 fTrailingEdge(-1),
655e379f 281 fErrorFlag(-1),
655e379f 282 fSector(-1),
283 fPlate(-1),
284 fStrip(-1),
285 fPadX(-1),
286 fPadZ(-1),
15ec34b9 287 fPackedDigits(0),
655e379f 288 fWordType(-1),
289 fSlotID(-1),
290 fACQ(-1),
291 fPSbit(-1),
655e379f 292 fTDCerrorFlag(-1),
293 fInsideDRM(kFALSE),
294 fInsideTRM(kFALSE),
295 fInsideLTM(kFALSE),
296 fInsideTRMchain0(kFALSE),
38f24f3f 297 fInsideTRMchain1(kFALSE),
298 fLocalEventCounterDRM(-1),
299 fLocalEventCounterLTM(-1),
300 fLocalEventCounterTRM(0x0),
32ead898 301 fLocalEventCounterChain(0x0),
302 fCableLengthMap(0x0),
303 fEventID(0)
d0eb8f39 304{
305 //
306 // copy constructor
307 //
308
309 fRawReader = stream.fRawReader;
15ec34b9 310
311 fTOFrawData = stream.fTOFrawData;
312
d0eb8f39 313 fDDL = stream.fDDL;
314 fTRM = stream.fTRM;
d0eb8f39 315 fTRMchain = stream.fTRMchain;
15ec34b9 316 fTDC = stream.fTDC;
d0eb8f39 317 fTDCchannel = stream.fTDCchannel;
15ec34b9 318 fTime = stream.fTime;
d0eb8f39 319 fToT = stream.fToT;
15ec34b9 320 fLeadingEdge = stream.fLeadingEdge;
321 fTrailingEdge = stream.fTrailingEdge;
322
d0eb8f39 323 fErrorFlag = stream.fErrorFlag;
324
325 fSector = stream.fSector;
326 fPlate = stream.fPlate;
327 fStrip = stream.fStrip;
328 fPadX = stream.fPadX;
329 fPadZ = stream.fPadZ;
330
15ec34b9 331 fPackedDigits = stream.fPackedDigits;
332
d0eb8f39 333 fWordType = stream.fWordType;
334 fSlotID = stream.fSlotID;
335 fACQ = stream.fACQ;
336 fPSbit = stream.fPSbit;
d0eb8f39 337 fTDCerrorFlag = stream.fTDCerrorFlag;
338 fInsideDRM = stream.fInsideDRM;
339 fInsideTRM = stream.fInsideTRM;
340 fInsideLTM = stream.fInsideLTM;
341 fInsideTRMchain0 = stream.fInsideTRMchain0;
342 fInsideTRMchain1 = stream.fInsideTRMchain1;
d0eb8f39 343
7813bb1a 344 for (Int_t i=0;i<AliDAQ::NumberOfDdls("TOF");i++){
345 fDataBuffer[i]= new AliTOFHitDataBuffer(*stream.fDataBuffer[i]);
346 fPackedDataBuffer[i]= new AliTOFHitDataBuffer(*stream.fPackedDataBuffer[i]);
347 }
348
349 fTOFrawData = new TClonesArray(*stream.fTOFrawData);
38f24f3f 350
351 fLocalEventCounterDRM = stream.fLocalEventCounterDRM;//adc
352 fLocalEventCounterLTM = stream.fLocalEventCounterLTM;//adc
353 for (Int_t j=0;j<13;j++){//adc
354 fLocalEventCounterTRM[j] = stream.fLocalEventCounterTRM[j];//adc
355 for (Int_t k=0;k<2;k++){//adc
356 fLocalEventCounterChain[j][k] = stream.fLocalEventCounterChain[j][k];//adc
357 }//adc
358 }//adc
359
32ead898 360 fCableLengthMap = stream.fCableLengthMap;
361
362 fEventID = stream.fEventID;
363
571dda3d 364}
365
366//_____________________________________________________________________________
367AliTOFRawStream& AliTOFRawStream::operator = (const AliTOFRawStream& stream)
368{
369 //
370 // assignment operator
371 //
372
373 fRawReader = stream.fRawReader;
15ec34b9 374
375 fTOFrawData = stream.fTOFrawData;
376
571dda3d 377 fDDL = stream.fDDL;
378 fTRM = stream.fTRM;
d0eb8f39 379 fTRMchain = stream.fTRMchain;
15ec34b9 380 fTDC = stream.fTDC;
d0eb8f39 381 fTDCchannel = stream.fTDCchannel;
15ec34b9 382 fTime = stream.fTime;
d0eb8f39 383 fToT = stream.fToT;
15ec34b9 384 fLeadingEdge = stream.fLeadingEdge;
385 fTrailingEdge = stream.fTrailingEdge;
571dda3d 386 fErrorFlag = stream.fErrorFlag;
d0eb8f39 387
388 fSector = stream.fSector;
389 fPlate = stream.fPlate;
390 fStrip = stream.fStrip;
391 fPadX = stream.fPadX;
392 fPadZ = stream.fPadZ;
571dda3d 393
15ec34b9 394 fPackedDigits = stream.fPackedDigits;
395
d0eb8f39 396 fWordType = stream.fWordType;
397 fSlotID = stream.fSlotID;
398 fACQ = stream.fACQ;
399 fPSbit = stream.fPSbit;
d0eb8f39 400 fTDCerrorFlag = stream.fTDCerrorFlag;
401 fInsideDRM = stream.fInsideDRM;
402 fInsideTRM = stream.fInsideTRM;
403 fInsideLTM = stream.fInsideLTM;
404 fInsideTRMchain0 = stream.fInsideTRMchain0;
405 fInsideTRMchain1 = stream.fInsideTRMchain1;
7813bb1a 406
407 for (Int_t i=0;i<AliDAQ::NumberOfDdls("TOF");i++){
408 fDataBuffer[i] = stream.fDataBuffer[i];
409 fPackedDataBuffer[i] = stream.fPackedDataBuffer[i];
410 }
411
412 fTOFrawData = stream.fTOFrawData;
38f24f3f 413
414 fLocalEventCounterDRM = stream.fLocalEventCounterDRM;//adc
415 fLocalEventCounterLTM = stream.fLocalEventCounterLTM;//adc
416 for (Int_t j=0;j<13;j++){//adc
417 fLocalEventCounterTRM[j] = stream.fLocalEventCounterTRM[j];//adc
418 for (Int_t k=0;k<2;k++){//adc
419 fLocalEventCounterChain[j][k] = stream.fLocalEventCounterChain[j][k];//adc
420 }//adc
421 }//adc
422
32ead898 423 fCableLengthMap = stream.fCableLengthMap;
424
425 fEventID = stream.fEventID;
426
571dda3d 427 return *this;
428
429}
430
431//_____________________________________________________________________________
432AliTOFRawStream::~AliTOFRawStream()
433{
434// destructor
435
15ec34b9 436 fPackedDigits = 0;
437
7813bb1a 438 for (Int_t i=0;i<72;i++){
439 delete fDataBuffer[i];
440 delete fPackedDataBuffer[i];
441 }
d3c7bfac 442
7813bb1a 443 delete fDecoder;
15ec34b9 444
7813bb1a 445 fTOFrawData->Clear();
446 delete fTOFrawData;
38f24f3f 447
415c84fa 448 delete [] fLocalEventCounterTRM;
32ead898 449 for (Int_t ii=0; ii<2; ii++)
450 delete [] fLocalEventCounterChain[ii];
451
452 delete fCableLengthMap;
38f24f3f 453
571dda3d 454}
455
456
15ec34b9 457//_____________________________________________________________________________
458
459void AliTOFRawStream::LoadRawData(Int_t indexDDL)
460{
0192b95f 461 //
462 // To load raw data
463 //
15ec34b9 464
74cd22f8 465 fTOFrawData->Clear();
0192b95f 466
74cd22f8 467 TClonesArray &arrayTofRawData = *fTOFrawData;
468
15ec34b9 469 fPackedDigits = 0;
470
471 // create raw data map
472 AliTOFRawMap *rawMap = new AliTOFRawMap(fTOFrawData);
473 rawMap->Clear();
474
475 Int_t slot[4] = {-1, -1, -1, -1};
476
415c84fa 477 fLocalEventCounterDRM = -1;
478 fLocalEventCounterLTM = -1;
479 for (Int_t ii=0; ii<13; ii++)
480 fLocalEventCounterTRM[ii] = -1;
481 for (Int_t ii=0; ii<13; ii++)
482 for (Int_t jj=0; jj<2; jj++)
483 fLocalEventCounterChain[ii][jj] = -1;
484
15ec34b9 485 fRawReader->Reset();
486 fRawReader->Select("TOF", indexDDL, indexDDL);
487
488 Bool_t signal = kFALSE;
489
74cd22f8 490 AliTOFrawData *rawDigit = NULL;
7813bb1a 491
15ec34b9 492 while(Next()) {
493
494 signal = (fSector!=-1 && fPlate!=-1 && fStrip!=-1 && fPadZ!=-1 && fPadX!=-1);
495 if (signal) {
38f24f3f 496 AliDebug(2,Form(" %2i %1i %2i %1i %2i", fSector, fPlate, fStrip, fPadZ, fPadX));
15ec34b9 497
498 slot[0] = fTRM;
499 slot[1] = fTRMchain;
500 slot[2] = fTDC;
501 slot[3] = fTDCchannel;
502
503 if (rawMap->TestHit(slot) != kEmpty) {
504
74cd22f8 505 rawDigit = static_cast<AliTOFrawData*>(rawMap->GetHit(slot));
15ec34b9 506
507 if (rawDigit->GetLeading()!=-1 && rawDigit->GetTrailing()==-1 &&
508 fLeadingEdge==-1 && fTrailingEdge!=-1) {
509
510 rawDigit->Update(fTime, fToT, fLeadingEdge, fTrailingEdge, fPSbit, fACQ, fErrorFlag);
511 }
512 else if ( ((rawDigit->GetTOF()!=-1 || rawDigit->GetLeading()!=-1 || rawDigit->GetTrailing()!=-1) &&
513 (fLeadingEdge!=-1 || fTrailingEdge!=-1 || fTime!=-1) )
514
0192b95f 515 )
15ec34b9 516 {
517
74cd22f8 518 new (arrayTofRawData[fPackedDigits++]) AliTOFrawData(fTRM, fTRMchain, fTDC, fTDCchannel, fTime, fToT, fLeadingEdge, fTrailingEdge, fPSbit, fACQ, fErrorFlag);
15ec34b9 519
74cd22f8 520 rawMap->SetHit(slot);
15ec34b9 521
74cd22f8 522 }
15ec34b9 523
524
525 }
526 else {
527
15ec34b9 528 new (arrayTofRawData[fPackedDigits++]) AliTOFrawData(fTRM, fTRMchain, fTDC, fTDCchannel, fTime, fToT, fLeadingEdge, fTrailingEdge, fPSbit, fACQ, fErrorFlag);
529
530 rawMap->SetHit(slot);
531
532 } // else if (rawMap->TestHit(slot) == kEmpty)
533
534 } // if (signal)
535
536 } // closed -> while (Next())
537
74cd22f8 538 rawMap->Delete();
7813bb1a 539
15ec34b9 540}
541
571dda3d 542//_____________________________________________________________________________
543Bool_t AliTOFRawStream::Next()
544{
d0eb8f39 545 //
546 // Read next 32-bit word in TOF raw data files
547 // returns kFALSE if there is no word left
548 //
571dda3d 549
550 UInt_t data;
551
38f24f3f 552 Int_t dummy = 0;
553
571dda3d 554 if (!fRawReader->ReadNextInt(data)) return kFALSE;
555
d0eb8f39 556 if (fSector!=-1 && fPlate!=-1 && fStrip!=-1 && fPadZ!=-1 && fPadX!=-1) {
557 fSector = -1;
15ec34b9 558 fPlate = -1;
559 fStrip = -1;
560 fPadZ = -1;
561 fPadX = -1;
562 fTime = -1;
563 fToT = -1;
564 fLeadingEdge = -1;
565 fTrailingEdge = -1;
d0eb8f39 566 }
571dda3d 567
571dda3d 568 fDDL = fRawReader->GetDDLID();
569
d0eb8f39 570 fWordType = GetField(data,WORD_TYPE_MASK,WORD_TYPE_POSITION);
571
572 switch (fWordType) { // switch word type
573
574 case GLOBAL_HEADER_TYPE: // global header
575 fSlotID = GetField(data, HEADER_SLOT_ID_MASK, HEADER_SLOT_ID_POSITION);
576 fTRM = fSlotID;
577
578
579 switch (fSlotID) { // switch global header slot ID
580
581 case DRM_ID_NUMBER: //DRM global header
582 if (fInsideDRM) { // unexpected DRM global headers -> exit
583 break;
584 }
585 fInsideDRM = kTRUE; // DRM global header accepted
586 break;
587
588 case LTM_ID_NUMBER: // LTM global header
589 if (fInsideLTM) { // unexpected LTM global headers -> exit
590 break;
591 }
592 fInsideLTM = kTRUE; // LTM global header accepted
593 break;
594
595 case 3: //TRM header
596 case 4: //TRM header
597 case 5: //TRM header
598 case 6: //TRM header
599 case 7: //TRM header
600 case 8: //TRM header
601 case 9: //TRM header
602 case 10: //TRM header
603 case 11: //TRM header
604 case 12: //TRM header
605 if (fInsideTRM) { // unexpected TRM global headers -> exit
606 break;
607 }
608 fInsideTRM = kTRUE; // TRM global header accepted
609 fACQ = GetField(data,TRM_ACQ_BITS_MASK,TRM_ACQ_BITS_POSITION);
610 break;
611
612 default: // unexpected global header slot ID
613 break;
614
615 } //end switch global header slot id
616
617 break;
618
619
620 case GLOBAL_TRAILER_TYPE: // global trailer
621 fSlotID = GetField(data,HEADER_SLOT_ID_MASK,HEADER_SLOT_ID_POSITION);
d0eb8f39 622
623 switch (fSlotID) { // switch global trailer slot ID
624
625 case DRM_ID_NUMBER: // DRM global trailer
626 if (!fInsideDRM) { // unexpected DRM global trailers -> exit
627 break;
628 }
38f24f3f 629 dummy = 0x0000fff0;
630 //AliInfo(Form(" DRM local event counter = %i", GetField(data,dummy,4)));
631 fLocalEventCounterDRM = GetField(data,dummy,4);//adc
d0eb8f39 632 fInsideDRM = kFALSE; // DRM global trailer accepted
633 fInsideTRM = kFALSE;
634 fInsideLTM = kFALSE;
635 fInsideTRMchain0 = kFALSE;
636 fInsideTRMchain1 = kFALSE;
d0eb8f39 637 fSector = -1;
15ec34b9 638 fPlate = -1;
639 fStrip = -1;
640 fPadZ = -1;
641 fPadX = -1;
642 fDDL = -1;
643 fTRM = -1;
644 fTDC = -1;
645 fTRMchain = -1;
d0eb8f39 646 fTDCchannel = -1;
15ec34b9 647 fTime = -1;
648 fToT = -1;
649 fLeadingEdge = -1;
650 fTrailingEdge = -1;
d0eb8f39 651 fErrorFlag = -1;
15ec34b9 652 fACQ = -1;
d0eb8f39 653 fPSbit = -1;
d0eb8f39 654 fTDCerrorFlag = -1;
655 break;
656 case LTM_ID_NUMBER: // LTM global trailer
657 if (!fInsideLTM) { // unexpected LTM global trailer -> exit
658 break;
659 }
38f24f3f 660 dummy = 0x0fff0000;
661 //AliInfo(Form(" LTM local event counter = %i", GetField(data,dummy,16)));
662 fLocalEventCounterLTM = GetField(data,dummy,16);//adc
d0eb8f39 663 fInsideLTM = kFALSE; // LTM global trailer accepted
664 break;
665 case 15: //TRM global trailer
666 if (!fInsideTRM) { // unexpected TRM global trailers -> exit
667 break;
668 }
38f24f3f 669 dummy = 0x0fff0000;
670 //AliInfo(Form(" TRM local event counter = %i", GetField(data,dummy,16)));
671 fLocalEventCounterTRM[fTRM] = GetField(data,dummy,16);//adc
d0eb8f39 672 fInsideTRM = kFALSE; // TRM global trailer accepted
673 break;
674 default: // unexpected global trailer slot ID
675 break;
676 } //end switch global trailer slot id
677
678
679 break;
680
681
682 case ERROR_TYPE: // TDC error
15ec34b9 683 fTDC = GetField(data,TRM_TDC_ERROR_TDC_ID_MASK,TRM_TDC_ERROR_TDC_ID_POSITION);
d0eb8f39 684 fTDCerrorFlag = GetField(data,TRM_TDC_ERROR_FLAGS_MASK,TRM_TDC_ERROR_FLAGS_POSITION);
685 break;
686
687
688 case FILLER_TYPE: // filler
689 break;
690
691
692 default: // other word types
693
694 if (fInsideTRM) { // inside TRM
695
696 switch (fWordType) { // switch word type inside TRM
697 case TRM_CHAIN0_HEADER_TYPE: // TRM chain0 header
698 if (fInsideTRMchain0) { // unexpected TRM chain0 header
699 break;
700 }
701 fInsideTRMchain0 = kTRUE;
702 fTRMchain = 0;
703 break;
704 case TRM_CHAIN0_TRAILER_TYPE: // TRM chain0 trailer
705 if (!fInsideTRMchain0) { // unexpected TRM chain0 trailer
706 break;
707 }
38f24f3f 708 dummy = 0x0fff0000;
709 //AliInfo(Form(" chain local event counter = %i", GetField(data,dummy,16)));
710 fLocalEventCounterChain[fTRM][fTRMchain] = GetField(data,dummy,16);//adc
d0eb8f39 711 fInsideTRMchain0 = kFALSE;
712 fTRMchain = -1;
713 break;
714 case TRM_CHAIN1_HEADER_TYPE: // TRM chain1 header
715 if (fInsideTRMchain1) { // unexpected TRM chain1 header
716 break;
717 }
718 fInsideTRMchain1 = kTRUE;
719 fTRMchain = 1;
720 break;
721 case TRM_CHAIN1_TRAILER_TYPE: // TRM chain1 trailer
722 if (!fInsideTRMchain1) { // unexpected TRM chain1 trailer
723 break;
724 }
38f24f3f 725 dummy = 0x0fff0000;
726 //AliInfo(Form(" chain local event counter = %i", GetField(data,dummy,16)));
727 fLocalEventCounterChain[fTRM][fTRMchain] = GetField(data,dummy,16);//adc
d0eb8f39 728 fInsideTRMchain1 = kFALSE;
729 fTRMchain = -1;
730 break;
731 } // end switch word type inside TRM
732
733 } // end if (fInsideTRM)
734
735
736 if (
737 ((fInsideTRMchain0&&!fInsideTRMchain1) || (!fInsideTRMchain0&&fInsideTRMchain1))
738 && fWordType!=TRM_CHAIN0_HEADER_TYPE && fWordType!=TRM_CHAIN0_TRAILER_TYPE
739 && fWordType!=TRM_CHAIN1_HEADER_TYPE && fWordType!=TRM_CHAIN1_TRAILER_TYPE
740 ){ // inside TRM chains
15ec34b9 741
742 fPSbit = GetField(data,TRM_PS_BITS_MASK,TRM_PS_BITS_POSITION);
743 fTDC = GetField(data,TRM_TDC_ID_MASK,TRM_TDC_ID_POSITION);
d0eb8f39 744 fTDCchannel = GetField(data,TRM_CHAN_MASK,TRM_CHAN_POSITION);
15ec34b9 745 fErrorFlag = GetField(data,TRM_E_BIT_MASK,TRM_E_BIT_POSITION);
746
747 SetSector();
748 SetPlate();
749 SetStrip();
750 SetPadZ();
751 SetPadX();
752
d0eb8f39 753
754 switch (fPSbit) { // switch fPSbit bits inside TRM chains
15ec34b9 755
756 case 0: // packing ok, digit time and TOT
757 fToT = GetField(data,TRM_TOT_WIDTH_MASK, TRM_TOT_WIDTH_POSITION);
32ead898 758 fTime = GetField(data,TRM_DIGIT_TIME_MASK,TRM_DIGIT_TIME_POSITION)
cfceaf15 759 /*-
760 fCableLengthMap->GetCableTimeShiftBin(fDDL, fTRM, fTRMchain, fTDC)*/
32ead898 761 /*+
762 (Int_t)(fgkddlBCshift[fDDL]*25.*1000./AliTOFGeometry::TdcBinWidth())
763 +
764 (Int_t)((fLocalEventCounterChain[fTRM][fTRMchain]-fEventID)*25.*1000./AliTOFGeometry::TdcBinWidth())*/
765 ;
d0eb8f39 766 break;
767
768 case 1: // leading edge digit, long digit time, no TOT
15ec34b9 769 //fToT = -1;
770 //fTime = -1;
32ead898 771 fLeadingEdge = GetField(data,TRM_LONG_DIGIT_TIME_MASK,TRM_LONG_DIGIT_TIME_POSITION)
cfceaf15 772 /*-
773 fCableLengthMap->GetCableTimeShiftBin(fDDL, fTRM, fTRMchain, fTDC)*/
32ead898 774 /*+
775 (Int_t)(fgkddlBCshift[fDDL]*25.*1000./AliTOFGeometry::TdcBinWidth())
776 +
777 (Int_t)((fLocalEventCounterChain[fTRM][fTRMchain]-fEventID)*25.*1000./AliTOFGeometry::TdcBinWidth())*/
778 ;
d0eb8f39 779 break;
780
781 case 2: // trailing edge digit, long digit time, no TOT
15ec34b9 782 //fToT = -1;
783 //fTime = -1;
32ead898 784 fTrailingEdge = GetField(data,TRM_LONG_DIGIT_TIME_MASK,TRM_LONG_DIGIT_TIME_POSITION)
cfceaf15 785 /*-
786 fCableLengthMap->GetCableTimeShiftBin(fDDL, fTRM, fTRMchain, fTDC)*/
32ead898 787 /*+
788 (Int_t)(fgkddlBCshift[fDDL]*25.*1000./AliTOFGeometry::TdcBinWidth())
789 +
790 (Int_t)((fLocalEventCounterChain[fTRM][fTRMchain]-fEventID)*25.*1000./AliTOFGeometry::TdcBinWidth())*/
791 ;
d0eb8f39 792 break;
15ec34b9 793
d0eb8f39 794 case 3: // TOT overflow
15ec34b9 795 fToT = GetField(data,TRM_TOT_WIDTH_MASK, TRM_TOT_WIDTH_POSITION);
32ead898 796 fTime = GetField(data,TRM_DIGIT_TIME_MASK,TRM_DIGIT_TIME_POSITION)
cfceaf15 797 /*-
798 fCableLengthMap->GetCableTimeShiftBin(fDDL, fTRM, fTRMchain, fTDC)*/
32ead898 799 /*+
800 (Int_t)(fgkddlBCshift[fDDL]*25.*1000./AliTOFGeometry::TdcBinWidth())
801 +
802 (Int_t)((fLocalEventCounterChain[fTRM][fTRMchain]-fEventID)*25.*1000./AliTOFGeometry::TdcBinWidth())*/
803 ;
d0eb8f39 804 break;
d0eb8f39 805
15ec34b9 806 } // end switch PS bits inside TRM chains
d0eb8f39 807
808 } // end if is inside TRM chains
809
810 } // end switch on fWordType
571dda3d 811
15ec34b9 812
571dda3d 813 return kTRUE;
7813bb1a 814
571dda3d 815}
816//_____________________________________________________________________________
817
d0eb8f39 818void AliTOFRawStream::SetSector()
571dda3d 819{
820 //
d0eb8f39 821 // Evaluate the TOF sector number -> [ 0;17]
822 // corresponding to the TOF equipment IDs:
823 // fDDL -> [ 0;71]
824 // fTRM -> [ 3;12]
15ec34b9 825 // fTRMchain -> [ 0; 1]
d0eb8f39 826 // fTDC -> [ 0;14]
827 // fTDCchannel -> [ 0; 7]
571dda3d 828 //
829
830 Int_t iSector = -1;
571dda3d 831
d0eb8f39 832 if (!(fDDL==-1)) iSector = Int_t((Float_t)(fDDL)/AliTOFGeometry::NDDL());
833
834 fSector = iSector;
571dda3d 835
836}
837//_____________________________________________________________________________
838
d0eb8f39 839
840void AliTOFRawStream::SetPlate()
571dda3d 841{
842 //
d0eb8f39 843 // Evaluate the TOF plate number ->[ 0; 4]
844 // corresponding to the TOF equipment IDs:
845 // fDDL -> [ 0;71]
846 // fTRM -> [ 3;12]
15ec34b9 847 // fTRMchain -> [ 0; 1]
d0eb8f39 848 // fTDC -> [ 0;14]
849 // fTDCchannel -> [ 0; 7]
571dda3d 850 //
851
852 Int_t iPlate = -1;
d0eb8f39 853 if (!(fDDL==-1 || fTRM==-1 || fTDC==-1
854 || fSector==-1))
855 iPlate = Equip2VolNplate(GetDDLnumberPerSector(fDDL), fTRM, fTDC);
571dda3d 856
d0eb8f39 857 fPlate = iPlate;
571dda3d 858
859}
860//_____________________________________________________________________________
861
d0eb8f39 862void AliTOFRawStream::SetStrip()
571dda3d 863{
864 //
d0eb8f39 865 // Evaluate the TOF strip number per module -> [ 0; 14/18]
866 // corresponding to the TOF equipment IDs:
867 // fDDL -> [ 0;71]
868 // fTRM -> [ 3;12]
15ec34b9 869 // fTRMchain -> [ 0; 1]
d0eb8f39 870 // fTDC -> [ 0;14]
871 // fTDCchannel -> [ 0; 7]
571dda3d 872 //
873
874 Int_t iStrip = -1;
875
d0eb8f39 876 if (!(fDDL==-1 || fTRM==-1 || fTDC==-1
877 || fSector==-1 || fPlate==-1))
878 iStrip = Equip2VolNstrip(GetDDLnumberPerSector(fDDL), fTRM, fTDC);
571dda3d 879
d0eb8f39 880 fStrip = iStrip;
571dda3d 881
882}
883//_____________________________________________________________________________
884
d0eb8f39 885void AliTOFRawStream::SetPadZ()
571dda3d 886{
887 //
d0eb8f39 888 // Evaluate the TOF padRow number per strip -> [ 0; 1]
889 // corresponding to the TOF equipment IDs:
890 // fDDL -> [ 0;71]
891 // fTRM -> [ 3;12]
15ec34b9 892 // fTRMchain -> [ 0; 1]
d0eb8f39 893 // fTDC -> [ 0;14]
894 // fTDCchannel -> [ 0; 7]
571dda3d 895 //
896
897 Int_t iPadZ = -1;
898
d0eb8f39 899 if (!(fDDL==-1 || fTRM==-1 || fTRMchain==-1 || fTDC==-1 || fTDCchannel==-1
900 || fSector==-1 || fPlate==-1 || fStrip==-1))
901 {
902 Int_t iPadAlongTheStrip = Equip2VolNpad(GetDDLnumberPerSector(fDDL), fTRMchain, fTDC, fTDCchannel);
903 if (iPadAlongTheStrip!=-1)
904 iPadZ = iPadAlongTheStrip%AliTOFGeometry::NpadZ();
905 }
571dda3d 906
38f24f3f 907 //iPadZ = Equip2VolNpad(GetDDLnumberPerSector(fDDL), fTRMchain, fTDC, fTDCchannel)%AliTOFGeometry::NpadZ();
908 //iPadZ = Equip2VolNpadZ(GetDDLnumberPerSector(fDDL), fTRMchain, fTDC, fTDCchannel);
909
d0eb8f39 910 fPadZ = iPadZ;
571dda3d 911
d0eb8f39 912}
913//_____________________________________________________________________________
914
915void AliTOFRawStream::SetPadX()
916{
917 //
918 // Evaluate the TOF pad number per strip padRow -> [ 0;47]
919 // corresponding to the TOF equipment IDs:
920 // fDDL -> [ 0;71]
921 // fTRM -> [ 3;12]
15ec34b9 922 // fTRMchain -> [ 0; 1]
d0eb8f39 923 // fTDC -> [ 0;14]
924 // fTDCchannel -> [ 0; 7]
925 //
571dda3d 926
d0eb8f39 927 Int_t iPadX = -1;
571dda3d 928
d0eb8f39 929 if (!(fDDL==-1 || fTRM==-1 || fTRMchain==-1 || fTDC==-1 || fTDCchannel==-1
930 || fSector==-1 || fPlate==-1 || fStrip==-1))
931 {
932 Int_t iPadAlongTheStrip = Equip2VolNpad(GetDDLnumberPerSector(fDDL), fTRMchain, fTDC, fTDCchannel);
933 if (iPadAlongTheStrip!=-1)
934 iPadX = (Int_t)(iPadAlongTheStrip/(Float_t(AliTOFGeometry::NpadZ())));
935 }
936
38f24f3f 937 //iPadX = (Int_t)(Equip2VolNpad(GetDDLnumberPerSector(fDDL), fTRMchain, fTDC, fTDCchannel)/(Float_t(AliTOFGeometry::NpadZ())));
938 //iPadX = Equip2VolNpadX(GetDDLnumberPerSector(fDDL), fTRMchain, fTDC, fTDCchannel);
939
d0eb8f39 940 fPadX = iPadX;
571dda3d 941
942}
571dda3d 943
d0eb8f39 944//----------------------------------------------------------------------------
945Int_t AliTOFRawStream::GetField(UInt_t word, Int_t fieldMask, Int_t fieldPosition) const
946{
947 //
948 //
949 //
950
951 return ((word & fieldMask) >> fieldPosition);
952}
953
954//----------------------------------------------------------------------------
13d6bff9 955Int_t AliTOFRawStream::Equip2VolNplate(Int_t iDDL, Int_t nTRM, Int_t nTDC)
571dda3d 956{
957 //
d0eb8f39 958 // Returns the TOF plate number [0;4]
959 // corresponding to the TOF equipment ID numbers:
960 // iDDL -> DDL number per sector [0;3]
961 // nTRM -> TRM number [3;12]
962 // nTDC -> TDC number [0;14]
571dda3d 963 //
964
d0eb8f39 965 Int_t iPlate = -1;
966 if (iDDL==0) {
967
968 if (nTRM>=4 && nTRM<7) {
969 iPlate = 0;
970 } else if (nTRM==7) {
971 if (nTDC<12) iPlate = 0;
972 else iPlate = 1;
973 } else if (nTRM>=8 && nTRM<11) {
974 iPlate = 1;
975 } else if (nTRM==11) {
976 if (nTDC<9) iPlate = 1;
977 else iPlate = 2;
978 }else if (nTRM==12) {
979 iPlate = 2;
980 }
981
982 } else if (iDDL==1) {
983
984 if (nTRM==3) {
985 if (nTDC<3) iPlate = 0;
986 } else if (nTRM>=4 && nTRM<7) {
987 iPlate = 0;
988 } else if (nTRM==7) {
989 if (nTDC<6) iPlate = 1;
990 else iPlate = 0;
991 } else if (nTRM>=8 && nTRM<11) {
992 iPlate = 1;
993 } else if (nTRM==11) {
994 if (nTDC<9) iPlate = 2;
995 else iPlate = 1;
996 } else if (nTRM==12) {
997 iPlate = 2;
998 }
999
1000 } else if (iDDL==2) {
1001
1002 if (nTRM>=4 && nTRM<7) {
1003 iPlate = 4;
1004 } else if (nTRM==7) {
1005 if (nTDC<12) iPlate = 4;
1006 else iPlate = 3;
1007 } else if (nTRM>=8 && nTRM<11) {
1008 iPlate = 3;
1009 } else if (nTRM==11) {
1010 if (nTDC<9) iPlate = 3;
1011 else iPlate = 2;
1012 }else if (nTRM==12) {
1013 iPlate = 2;
1014 }
1015
1016 } else if (iDDL==3) {
1017
1018 if (nTRM==3) {
1019 if (nTDC<3) iPlate = 4;
1020 } else if (nTRM>=4 && nTRM<7) {
1021 iPlate = 4;
1022 } else if (nTRM==7) {
1023 if (nTDC<6) iPlate = 3;
1024 else iPlate = 4;
1025 } else if (nTRM>=8 && nTRM<11) {
1026 iPlate = 3;
1027 } else if (nTRM==11) {
1028 if (nTDC<9) iPlate = 2;
1029 else iPlate = 3;
1030 } else if (nTRM==12) {
1031 iPlate = 2;
1032 }
571dda3d 1033
d0eb8f39 1034 }
571dda3d 1035
d0eb8f39 1036 return iPlate;
571dda3d 1037
d0eb8f39 1038}
571dda3d 1039
d0eb8f39 1040//----------------------------------------------------------------------------
13d6bff9 1041Int_t AliTOFRawStream::Equip2VolNstrip(Int_t iDDL, Int_t nTRM, Int_t nTDC)
d0eb8f39 1042{
1043 //
1044 // Returns the TOF strip number per module:
1045 // [0;14], in the central plates,
1046 // [0;18], in the intermediate and external plates
1047 // corresponding to the TOF equipment ID numbers:
1048 // iDDL -> DDL number per sector [0;3]
1049 // nTRM -> TRM number [3;12]
1050 // nTDC -> TDC number [0;14]
1051 //
571dda3d 1052
d0eb8f39 1053 Int_t iStrip = -1;
571dda3d 1054
d0eb8f39 1055 if (iDDL==0) {
1056
1057 if (nTRM== 4) iStrip = (Int_t)(nTDC/3.);
1058 else if (nTRM== 5) iStrip = 5 + (Int_t)(nTDC/3.);
1059 else if (nTRM== 6) iStrip = 10 + (Int_t)(nTDC/3.);
1060 else if (nTRM== 7) {
1061 if (nTDC<12) iStrip = 15 + (Int_t)(nTDC/3.);
1062 else iStrip = (Int_t)(nTDC/3.) - 4;
1063 }
1064 else if (nTRM== 8) iStrip = 1 + (Int_t)(nTDC/3.);
1065 else if (nTRM== 9) iStrip = 6 + (Int_t)(nTDC/3.);
1066 else if (nTRM==10) iStrip = 11 + (Int_t)(nTDC/3.);
1067 else if (nTRM==11) {
1068 if (nTDC<9) iStrip = 16 + (Int_t)(nTDC/3.);
1069 else iStrip = (Int_t)(nTDC/3.) - 3;
1070 }
1071 else if (nTRM==12) iStrip = 2 + (Int_t)(nTDC/3.);
1072
1073 } else if (iDDL==1) {
1074
1075 if (nTRM==3 && nTDC<3) iStrip = (Int_t)(nTDC/3.);
1076 else if (nTRM== 4) iStrip = 5 - (Int_t)(nTDC/3.);
1077 else if (nTRM== 5) iStrip = 10 - (Int_t)(nTDC/3.);
1078 else if (nTRM== 6) iStrip = 15 - (Int_t)(nTDC/3.);
1079 else if (nTRM== 7) {
1080 if (nTDC<6) iStrip = 1 - (Int_t)(nTDC/3.);
1081 else iStrip = 20 - (Int_t)(nTDC/3.);
1082 }
1083 else if (nTRM== 8) iStrip = 6 - (Int_t)(nTDC/3.);
1084 else if (nTRM== 9) iStrip = 11 - (Int_t)(nTDC/3.);
1085 else if (nTRM==10) iStrip = 16 - (Int_t)(nTDC/3.);
1086 else if (nTRM==11) {
1087 if (nTDC<9) iStrip = 2 - (Int_t)(nTDC/3.);
1088 else iStrip = 21 - (Int_t)(nTDC/3.);
1089 }
1090 else if (nTRM==12) iStrip = 7 - (Int_t)(nTDC/3.);
1091
1092 } else if (iDDL==2) {
1093
1094 if (nTRM== 4) iStrip = 18 - (Int_t)(nTDC/3.);
1095 else if (nTRM== 5) iStrip = 18 - ( 5 + (Int_t)(nTDC/3.));
1096 else if (nTRM== 6) iStrip = 18 - (10 + (Int_t)(nTDC/3.));
1097 else if (nTRM== 7) {
1098 if (nTDC<12) iStrip = 18 - (15 + (Int_t)(nTDC/3.));
1099 else iStrip = 18 - ((Int_t)(nTDC/3.) - 4);
1100 }
1101 else if (nTRM== 8) iStrip = 18 - ( 1 + (Int_t)(nTDC/3.));
1102 else if (nTRM== 9) iStrip = 18 - ( 6 + (Int_t)(nTDC/3.));
1103 else if (nTRM==10) iStrip = 18 - (11 + (Int_t)(nTDC/3.));
1104 else if (nTRM==11) {
1105 if (nTDC<9) iStrip = 18 - (16 + (Int_t)(nTDC/3.));
1106 else iStrip = 14 - ((Int_t)(nTDC/3.) - 3);
1107 }
1108 else if (nTRM==12) iStrip = 14 - ( 2 + (Int_t)(nTDC/3.));
1109
1110 } else if (iDDL==3) {
1111
1112 if (nTRM==3 && nTDC<3) iStrip = 18 - (Int_t)(nTDC/3.);
1113 else if (nTRM== 4) iStrip = 18 - ( 5 - (Int_t)(nTDC/3.));
1114 else if (nTRM== 5) iStrip = 18 - (10 - (Int_t)(nTDC/3.));
1115 else if (nTRM== 6) iStrip = 18 - (15 - (Int_t)(nTDC/3.));
1116 else if (nTRM== 7) {
1117 if (nTDC<6) iStrip = 18 - (1 - (Int_t)(nTDC/3.));
1118 else iStrip = 18 - (20 - (Int_t)(nTDC/3.));
1119 }
1120 else if (nTRM== 8) iStrip = 18 - ( 6 - (Int_t)(nTDC/3.));
1121 else if (nTRM== 9) iStrip = 18 - (11 - (Int_t)(nTDC/3.));
1122 else if (nTRM==10) iStrip = 18 - (16 - (Int_t)(nTDC/3.));
1123 else if (nTRM==11) {
1124 if (nTDC<9) iStrip = 14 - ( 2 - (Int_t)(nTDC/3.));
1125 else iStrip = 18 - (21 - (Int_t)(nTDC/3.));
1126 }
1127 else if (nTRM==12) iStrip = 14 - ( 7 - (Int_t)(nTDC/3.));
1128
1129 }
1130
1131 return iStrip;
1132
1133}
1134
1135//----------------------------------------------------------------------------
1136Int_t AliTOFRawStream::Equip2VolNpad(Int_t iDDL, Int_t iChain, Int_t nTDC,
13d6bff9 1137 Int_t iCH)
d0eb8f39 1138{
1139 //
1140 // Returns the TOF pad number per strip [0;95]
1141 // corresponding to the TOF equipment ID numbers:
1142 // iDDL -> DDL number per sector [0;3]
1143 // iChain -> TRM chain number [0;1]
1144 // nTDC -> TDC number [0;14]
1145 // iCH -> TDC channel number [0;7]
1146 //
1147
1148 Int_t iPadAlongTheStrip = -1;
1149
38f24f3f 1150 // wrong
1151 //Int_t iTDClocal = nTDC%3 + (1-iChain)*3;
1152 //if (iDDL==0 || iDDL==3) iTDClocal = 5 - iTDClocal;
1153 //else if (iDDL==1 || iDDL==2) iTDClocal = 6 + (5 - iTDClocal);
d0eb8f39 1154
38f24f3f 1155 // right
1156 Int_t iTDClocal = -1;
1157 Int_t iTDClocal03 = nTDC%3 + (1-iChain)*3;
1158 Int_t iTDClocal12 = 2-nTDC%3 + iChain*3;
1159 if (iDDL==0 || iDDL==3) iTDClocal = 5 - iTDClocal03;
1160 else if (iDDL==1 || iDDL==2) iTDClocal = 6 + (5 - iTDClocal12);
d0eb8f39 1161
1162 Int_t iCHlocal = iCH;
1163 if (iDDL==0 || iDDL==3) iCHlocal = 7 - iCH;
1164
1165 iPadAlongTheStrip = iTDClocal*AliTOFGeometry::NCh() + iCHlocal;
1166
1167 if (((iDDL==1 || iDDL==2) && iPadAlongTheStrip< AliTOFGeometry::NpadX()) ||
019f6101 1168 ((iDDL==0 || iDDL==3) && iPadAlongTheStrip>=AliTOFGeometry::NpadX())) {
13d6bff9 1169 std::cerr << "Problems with the padX number!" << endl;
1170 //AliWarning("Problems with the padX number!");
019f6101 1171 }
d0eb8f39 1172 return iPadAlongTheStrip;
1173
1174}
1175
38f24f3f 1176//----------------------------------------------------------------------------
1177Int_t AliTOFRawStream::Equip2VolNpadX(Int_t iDDL, Int_t iChain, Int_t nTDC,
1178 Int_t iCH)
1179{
1180 //
1181 // Returns the TOF padX number [0;47]
1182 // corresponding to the TOF equipment ID numbers:
1183 // iDDL -> DDL number per sector [0;3]
1184 // iChain -> TRM chain number [0;1]
1185 // nTDC -> TDC number [0;14]
1186 // iCH -> TDC channel number [0;7]
1187 //
1188
1189 Int_t iPadX = (Int_t)(AliTOFRawStream::Equip2VolNpad(iDDL, iChain, nTDC, iCH)/
1190 (Float_t(AliTOFGeometry::NpadZ())));
1191
1192 return iPadX;
1193
1194}
1195
1196//----------------------------------------------------------------------------
1197Int_t AliTOFRawStream::Equip2VolNpadZ(Int_t iDDL, Int_t iChain, Int_t nTDC,
1198 Int_t iCH)
1199{
1200 //
1201 // Returns the TOF padZ number [0;1]
1202 // corresponding to the TOF equipment ID numbers:
1203 // iDDL -> DDL number per sector [0;3]
1204 // iChain -> TRM chain number [0;1]
1205 // nTDC -> TDC number [0;14]
1206 // iCH -> TDC channel number [0;7]
1207 //
1208
1209 Int_t iPadZ = AliTOFRawStream::Equip2VolNpad(iDDL, iChain, nTDC, iCH)%AliTOFGeometry::NpadZ();
1210
1211 return iPadZ;
1212
1213}
1214
d0eb8f39 1215//----------------------------------------------------------------------------
1216Int_t AliTOFRawStream::GetSectorNumber(Int_t nDDL) const
1217{
1218 //
1219 // Returns the sector number [0;17]
1220 // corresponing to the assigned DRM/DDL number [0;71]
1221 //
1222
1223 Int_t iSector = Int_t((Float_t)(nDDL)/AliTOFGeometry::NDDL());
1224
1225 return iSector;
1226
1227}
1228//----------------------------------------------------------------------------
1229Int_t AliTOFRawStream::GetDDLnumberPerSector(Int_t nDDL) const
1230{
1231 //
1232 // Return the DRM/DDL number per sector [0;3]
1233 // corresponing to the assigned DRM/DDL number [0;71]
1234 //
1235
1236 Int_t iDDL = nDDL%AliTOFGeometry::NDDL();
1237
1238 return iDDL;
1239
1240}
1241
7813bb1a 1242//----------------------------------------------------------------------------
1243void AliTOFRawStream::EquipmentId2VolumeId(AliTOFHitData *hitData, Int_t *volume) const
1244{
1245 EquipmentId2VolumeId(hitData->GetDDLID(),hitData->GetSlotID(),hitData->GetChain(),hitData->GetTDC(),hitData->GetChan(),volume);
1246}
d0eb8f39 1247//----------------------------------------------------------------------------
1248void AliTOFRawStream::EquipmentId2VolumeId(Int_t nDDL, Int_t nTRM, Int_t iChain,
1249 Int_t nTDC, Int_t iCH,
1250 Int_t *volume) const
1251{
1252 //
1253 // To convert:
1254 // nDDL (variable in [0;71]) -> number of the DDL file
1255 // nTRM (variable in [3;12]) -> number of the TRM slot
1256 // iChain (variable in [0; 1]) -> number of the TRM chain
1257 // nTDC (variable in [0;14]) -> number of the TDC
1258 // iCH (variable in [0; 7]) -> number of the TDC channel
1259 //
1260 // in:
1261 // sector number, i.e. volume[0] (variable in [0,17])
1262 // plate number, i.e. volume[1] (variable in [0, 5])
1263 // strip number, i.e. volume[2] (variable in [0,14/18])
1264 // padX number, i.e. volume[3] (variable in [0,47])
1265 // padZ number, i.e. volume[4] (variable in [0, 1])
1266 //
1267
1268 Int_t iDDL = GetDDLnumberPerSector(nDDL);
1269
1270 Int_t iSector = GetSectorNumber(nDDL);
1271
1272 Int_t iPlate = Equip2VolNplate(iDDL, nTRM, nTDC);
019f6101 1273 if (iPlate==-1) {
13d6bff9 1274 if (fRawReader)
1275 fRawReader->AddMajorErrorLog(kPlateError,"plate = -1");
019f6101 1276 AliWarning("Problems with the plate number!");
1277 }
d0eb8f39 1278
1279 Int_t iStrip = Equip2VolNstrip(iDDL, nTRM, nTDC);
019f6101 1280 if (iStrip==-1) {
13d6bff9 1281 if (fRawReader)
1282 fRawReader->AddMajorErrorLog(kStripError,"strip = -1");
019f6101 1283 AliWarning("Problems with the strip number!");
1284 }
d0eb8f39 1285
1286 Int_t iPadAlongTheStrip = Equip2VolNpad(iDDL, iChain, nTDC, iCH);
019f6101 1287 if (iPadAlongTheStrip==-1){
13d6bff9 1288 if (fRawReader)
1289 fRawReader->AddMajorErrorLog(kPadAlongStripError,"pad = -1");
019f6101 1290 AliWarning("Problems with the pad number along the strip!");
1291 }
13d6bff9 1292
d0eb8f39 1293 Int_t iPadX = (Int_t)(iPadAlongTheStrip/(Float_t(AliTOFGeometry::NpadZ())));
1294 Int_t iPadZ = iPadAlongTheStrip%AliTOFGeometry::NpadZ();
38f24f3f 1295
1296 //Int_t iPadX = (Int_t)(Equip2VolNpad(iDDL, iChain, nTDC, iCH)/(Float_t(AliTOFGeometry::NpadZ())));
1297 //Int_t iPadZ = Equip2VolNpad(iDDL, iChain, nTDC, iCH)%AliTOFGeometry::NpadZ();
1298
1299 //Int_t iPadX = Equip2VolNpadX(iDDL, iChain, nTDC, iCH);
1300 //Int_t iPadZ = Equip2VolNpadZ(iDDL, iChain, nTDC, iCH);
1301
d0eb8f39 1302 volume[0] = iSector;
1303 volume[1] = iPlate;
1304 volume[2] = iStrip;
1305 volume[3] = iPadX;
1306 volume[4] = iPadZ;
571dda3d 1307
1308}
7813bb1a 1309//-----------------------------------------------------------------------------
019f6101 1310Bool_t AliTOFRawStream::DecodeDDL(Int_t nDDLMin, Int_t nDDLMax, Int_t verbose = 0) {
1311 //
1312 // To decode raw data for DDL number in [nDDLmin; nDDLmax]
1313 //
1314
7813bb1a 1315 //check and fix valid DDL range
019f6101 1316 if (nDDLMin < 0){
1317 nDDLMin = 0;
1318 fRawReader->AddMinorErrorLog(kDDLMinError);
1319 AliWarning("Wrong DDL range: setting first DDL ID to 0");
7813bb1a 1320 }
019f6101 1321 if (nDDLMax > 71){
1322 nDDLMax = 71;
1323 fRawReader->AddMinorErrorLog(kDDLMaxError);
1324 AliWarning("Wrong DDL range: setting last DDL ID to 71");
7813bb1a 1325 }
1326
1327 //select required DDLs
019f6101 1328 fRawReader->Select("TOF", nDDLMin, nDDLMax);
7813bb1a 1329
1330 if (verbose)
019f6101 1331 AliInfo(Form("Selected TOF DDL range: %d-%d", nDDLMin, nDDLMax));
7813bb1a 1332
2fd622fa 1333 return(Decode(verbose));
1334}
1335//-----------------------------------------------------------------------------
019f6101 1336Bool_t AliTOFRawStream::Decode(Int_t verbose = 0) {
1337 //
1338 // New decoder method
1339 //
1340
2fd622fa 1341 Int_t currentEquipment;
1342 Int_t currentDDL;
1343
1344 //pointers
1345 UChar_t *data = 0x0;
1346
7813bb1a 1347 //loop and read DDL headers
1348 while(fRawReader->ReadHeader()){
1349
1350 //memory leak prevention (actually data should be always 0x0 here)
1351 if (data != 0x0)
1352 delete [] data;
1353
1354 //get equipment infos
1355 currentEquipment = fRawReader->GetEquipmentId();
1356 currentDDL = fRawReader->GetDDLID();
1357 const Int_t kDataSize = fRawReader->GetDataSize();
1358 const Int_t kDataWords = kDataSize / 4;
1359 data = new UChar_t[kDataSize];
1360
1361 if (verbose)
1362 AliInfo(Form("Found equipment # %d header (DDL # %d): %d bytes (%d words)", currentEquipment, currentDDL, kDataSize, kDataWords));
1363
1364 if (verbose)
1365 AliInfo(Form("Reading equipment #%d (DDL # %d) data...", currentEquipment, currentDDL));
1366
1367 //read equipment payload
1368 if (!fRawReader->ReadNext(data, kDataSize))
1369 {
019f6101 1370 fRawReader->AddMajorErrorLog(kDDLdataReading);
7813bb1a 1371 if (verbose)
019f6101 1372 AliWarning("Error while reading DDL data. Go to next equipment");
7813bb1a 1373 delete [] data;
1374 data = 0x0;
1375 continue;
1376 }
1377
1378 if (verbose)
1379 AliInfo(Form("Equipment # %d (DDL # %d) data has been readed", currentEquipment, currentDDL));
1380
1381
1382 //set up the decoder
1383 fDecoder->SetVerbose(verbose);
1384 fDecoder->SetDataBuffer(fDataBuffer[currentDDL]);
1385 fDecoder->SetPackedDataBuffer(fPackedDataBuffer[currentDDL]);
1386
1387 //start decoding
019f6101 1388 if (fDecoder->Decode((UInt_t *)data, kDataWords) == kTRUE) {
1389 fRawReader->AddMajorErrorLog(kDDLDecoder,Form("DDL # = %d",currentDDL));
1390 AliWarning(Form("Error while decoding DDL # %d: decoder returned with errors", currentDDL));
1391 }
7813bb1a 1392
1393 delete [] data;
1394 data = 0x0;
1395 }
1396
1397 //reset reader
1398 fRawReader->Reset();
1399
1400 if (verbose)
1401 AliInfo("All done");
1402
1403 return kFALSE;
1404
1405}
1406//---------------------------------------------------------------------------
1407void
1408AliTOFRawStream::ResetBuffers()
1409{
19eefa2b 1410 //
1411 // To reset the buffers
1412 //
1413
7813bb1a 1414 for (Int_t iDDL = 0; iDDL < AliDAQ::NumberOfDdls("TOF"); iDDL++){
1415 ResetDataBuffer(iDDL);
1416 ResetPackedDataBuffer(iDDL);
1417 }
1418}
1419
1420//---------------------------------------------------------------------------
1421Bool_t
1422AliTOFRawStream::LoadRawDataBuffers(Int_t indexDDL, Int_t verbose)
1423{
19eefa2b 1424 //
1425 // To load the buffers
1426 //
1427
7813bb1a 1428 fTOFrawData->Clear();
1429 fPackedDigits = 0;
1430
1431 if (verbose > 0)
1432 AliInfo(Form("Decoding raw data for DDL # %d ...", indexDDL));
1433
1434 if (DecodeDDL(indexDDL, indexDDL, verbose) != 0){ //decode DDL
019f6101 1435 fRawReader->AddMajorErrorLog(kDDLDecoder,Form("DDL # = %d",indexDDL));
1436 AliWarning(Form("Error while decoding DDL # %d", indexDDL));
7813bb1a 1437 return kTRUE;
1438 }
1439
1440 if (verbose > 0)
1441 AliInfo(Form("Done. %d packed %s been found.", fPackedDataBuffer[indexDDL]->GetEntries(), fPackedDataBuffer[indexDDL]->GetEntries() > 1 ? "hits have" : "hit has"));
1442
1443 AliTOFHitData *hitData; //hit data pointer
1444
1445 if (verbose > 0)
1446 AliInfo("Filling TClonesArray ...");
1447
1448 //loop over DDL packed hits
1449 for (Int_t iHit = 0; iHit < fPackedDataBuffer[indexDDL]->GetEntries(); iHit++){
1450 hitData = fPackedDataBuffer[indexDDL]->GetHit(iHit); //get hit data
1451 Int_t hitACQ = hitData->GetACQ();
1452 Int_t hitPS = hitData->GetPS();
1453 Int_t hitSlotID = hitData->GetSlotID();
1454 Int_t hitChain = hitData->GetChain();
1455 Int_t hitTDC = hitData->GetTDC();
1456 Int_t hitChan = hitData->GetChan();
7813bb1a 1457 Int_t hitTimeBin = hitData->GetTimeBin();
7813bb1a 1458 Int_t hitTOTBin = hitData->GetTOTBin();
1459
32ead898 1460 Int_t hitLeading = hitData->GetTimeBin()
cfceaf15 1461 /*-
1462 fCableLengthMap->GetCableTimeShiftBin(indexDDL, hitSlotID, hitChain, hitTDC)*/;//-1; // adc
bf33f8f0 1463 Int_t hitTrailing = -1;
1464 Int_t hitError = -1;
7813bb1a 1465
1466 TClonesArray &arrayTofRawData = *fTOFrawData;
bf33f8f0 1467 new (arrayTofRawData[fPackedDigits++]) AliTOFrawData(hitSlotID, hitChain, hitTDC, hitChan, hitTimeBin, hitTOTBin, hitLeading, hitTrailing, hitPS, hitACQ, hitError);
7813bb1a 1468 }
1469
1470 if (verbose > 0)
1471 AliInfo("Done.");
1472
1473 if (verbose > 0)
1474 AliInfo("Resetting buffers ...");
1475
1476 fDataBuffer[indexDDL]->Reset();
1477 fPackedDataBuffer[indexDDL]->Reset();
1478
1479 if (verbose > 0)
1480 AliInfo("Done.");
1481
1482 return kFALSE;
1483}
1484