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