TOF/AliTOFRawStream.cxx

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