HMPID/AliHMPIDRawStream.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 ///
  18 /// This is a class for reading the HMPID raw data
  19 /// The format of the raw data corresponds to the one
  20 /// which was documented by Paolo Martinengo.
  21 ///
  22 ///////////////////////////////////////////////////////////////////////////////
  23
  24 #include "AliHMPIDRawStream.h"
  25 #include "AliRawReader.h"
  26 #include "AliLog.h"
  27
  28 Int_t fPos[170000];
  29 Int_t iPos = 0;
  30 //static Bool_t stDeb = kTRUE;
  31
  32 ClassImp(AliHMPIDRawStream)
  33
  34 //_____________________________________________________________________________
  35 AliHMPIDRawStream::AliHMPIDRawStream(AliRawReader* rawReader) :
  36   fDDLNumber(-1),
  37   fRawReader(rawReader),
  38   fData(NULL),
  39   fPosition(-1)
  40 {
  41   // Constructor
  42   fNPads = 0;
  43   fCharge = 0x0;
  44   fPad =0x0;
  45
  46   for(Int_t l=1; l < kSumErr; l++) fNumOfErr[l]=0;               //reset errors
  47
  48   fRawReader->Reset();
  49   fRawReader->Select("HMPID");
  50 }
  51 //-----------------------------------------------------------------------------
  52 AliHMPIDRawStream::AliHMPIDRawStream() :
  53   fDDLNumber(-1),
  54   fRawReader(0x0),
  55   fData(NULL),
  56   fPosition(-1)
  57 {
  58   // Constructor
  59   for(Int_t l=1; l < kSumErr; l++) fNumOfErr[l]=0;               //reset errors
  60 }
  61 //_____________________________________________________________________________
  62 AliHMPIDRawStream::~AliHMPIDRawStream()
  63 {
  64   // destructor
  65     DelVars();
  66 }
  67 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  68 void AliHMPIDRawStream::Reset()
  69 {
  70   // reset raw stream params
  71   // Reinitalize the containers
  72   fDDLNumber = -1;
  73   fPosition = -1;
  74   fData = NULL;
  75   if (fRawReader) fRawReader->Reset();
  76 }
  77 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  78 Bool_t AliHMPIDRawStream::Next()
  79 {
  80   // read next DDL raw data from the HMPID raw data stream
  81   // return kFALSE in case of error or no data left
  82   AliDebug(1,"Start.");
  83   do {
  84     if (!fRawReader->ReadNextData(fData)) return kFALSE;
  85   } while (fRawReader->GetDataSize() == 0);
  86
  87
  88   /*
  89   Event type is selected as in $ALICE_ROOT/RAW/event.h
  90   #define START_OF_RUN                    ((eventTypeType) 1)
  91   #define END_OF_RUN                      ((eventTypeType) 2)
  92   #define START_OF_RUN_FILES              ((eventTypeType) 3)
  93   #define END_OF_RUN_FILES                ((eventTypeType) 4)
  94   #define START_OF_BURST                  ((eventTypeType) 5)
  95   #define END_OF_BURST                    ((eventTypeType) 6)
  96   #define PHYSICS_EVENT                   ((eventTypeType) 7) <<---------------
  97   #define CALIBRATION_EVENT               ((eventTypeType) 8)
  98   #define EVENT_FORMAT_ERROR              ((eventTypeType) 9)
  99   #define START_OF_DATA                   ((eventTypeType)10)
 100   #define END_OF_DATA                     ((eventTypeType)11)
 101   #define SYSTEM_SOFTWARE_TRIGGER_EVENT   ((eventTypeType)12)
 102   #define DETECTOR_SOFTWARE_TRIGGER_EVENT ((eventTypeType)13)
 103   #define EVENT_TYPE_MIN                  1
 104   #define EVENT_TYPE_MAX                  13
 105   */
 106
 107   fPosition = 0;
 108   Bool_t status;
 109
 110   if(fRawReader->GetType() == 7)  {                             //New: Select Physics events, Old: Raw data size is not 0 and not 47148 (pedestal)
 111     fDDLNumber = fRawReader->GetDDLID();
 112 //    Printf("DDL %i started to be decoded!.",fDDLNumber);
 113     InitVars(fRawReader->GetDataSize()/4);
 114     status = ReadHMPIDRawData();
 115 //    if(status) Printf("Event DDL %i successfully decoded!.",fDDLNumber);
 116 //    else Printf("Event DDL %i ERROR in decoding!.",fDDLNumber);
 117 //    if(stDeb) DumpData(fRawReader->GetDataSize());
 118 //    stDeb=kFALSE;
 119   }
 120 //  return status;
 121   return kTRUE;
 122 }
 123 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 124 void AliHMPIDRawStream::InitVars(Int_t n)
 125 {
 126   fNPads = 0;
 127   fCharge = new Int_t[n];
 128   fPad = new Int_t[n];
 129   //for debug purpose
 130   for(Int_t i=0; i < 170000; i++) fPos[i]=0;                     //reset debug
 131   iPos = 0;
 132 }
 133 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 134 void AliHMPIDRawStream::DelVars()
 135 {
 136   fNPads = 0;
 137   delete fCharge;
 138   delete fPad;
 139 }
 140 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 141 Bool_t AliHMPIDRawStream::ReadHMPIDRawData()
 142 {
 143   Int_t cntGlob = fRawReader->GetDataSize()/4;
 144   Int_t cnt = cntGlob;
 145   Int_t nwSeg;
 146   Int_t cntSegment;
 147
 148   fWord = GetWord(cnt);cnt--;
 149
 150
 151   while (cnt>0) {
 152
 153     nwSeg = (fWord >> kbit8) & 0xfff;
 154     if(!CheckSegment()) return kFALSE;
 155     if(!ReadSegment(cntSegment)) return kFALSE;
 156
 157     if(nwSeg != cntSegment) {Printf("Error in Segment counters: %i different wrt %i",nwSeg,cntSegment);return kFALSE;}
 158 //    Printf(" cnt %i cntSegment %i",cnt,cntSegment);
 159     fWord = GetWord(cntSegment+1,kBwd);
 160     cnt-=cntSegment+1;
 161   }
 162
 163   return kTRUE;
 164
 165 }
 166 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 167 Bool_t AliHMPIDRawStream::ReadSegment(Int_t &cntSegment)
 168 {
 169   cntSegment = (fWord >> kbit8) & 0xfff;
 170   Int_t cnt = cntSegment;
 171   Int_t cntRow;
 172   Int_t nwRow;
 173
 174   fWord = GetWord(cnt,kBwd);
 175
 176   while (cnt>0) {
 177
 178     cntRow  = (fWord >> kbit16) & 0xfff;
 179     if(!CheckRowMarker()) return kFALSE;
 180     if(!ReadRow(nwRow)) return kFALSE;
 181
 182     if(nwRow != cntRow) {Printf("Error in Row counters: %i different wrt %i",nwRow,cntRow);return kFALSE;}
 183     fWord = GetWord(cntRow+1);
 184     cnt -= cntRow+1;
 185
 186   }
 187
 188   cntSegment -= cnt;
 189
 190   return kTRUE;
 191
 192 }
 193 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 194 Bool_t AliHMPIDRawStream::ReadRow(Int_t &cntRow)
 195 {
 196   Int_t cnt;
 197   Int_t cntDilogic;
 198   Int_t nwDil;
 199
 200   cntRow  = (fWord >> kbit16) & 0xfff;
 201   cnt = cntRow;
 202
 203   fWord = GetWord(cntRow);
 204
 205   while (cnt>0) {
 206
 207     if(!CheckEoE(nwDil)) return kFALSE;
 208     if(!ReadDilogic(cntDilogic)) return kFALSE;
 209
 210     if(nwDil != cntDilogic) {Printf("Error in Dilogic counters: %i different wrt %i",nwDil,cntDilogic);return kFALSE;}
 211     cnt -= cntDilogic;
 212     fWord = GetWord(1,kBwd); // go to next Dilogic bank...
 213     cnt--;
 214 //    Printf(" cnt %i cntDilogic %i ",cnt,cntDilogic);
 215   }
 216
 217   cntRow -= cnt;
 218
 219   return kTRUE;
 220
 221 }
 222 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 223 Bool_t AliHMPIDRawStream::ReadDilogic(Int_t &cntDilogic)
 224 {
 225   cntDilogic = fWord & 0x7f;
 226
 227   Int_t cnt = cntDilogic;
 228
 229 //  Printf(" cnt DILOGIC %i at %i word %08X",cnt,fPosition,fWord);
 230
 231   for(Int_t iDil=0;iDil<cntDilogic;iDil++) {
 232     UInt_t dilogic = 0, row = 0;
 233     fWord = GetWord(1,kBwd);
 234 //check on row number
 235     cnt--;
 236     row = (fWord >> kbit22) & 0x1f;
 237     if(!CheckRow(row)) continue;
 238 //check dilogic number
 239     dilogic = (fWord >> kbit18) & 0xf;                                              //dilogic info in raw word is between bits: 18...21
 240     if(!CheckDilogic(dilogic)) continue;
 241 //check pad number
 242     UInt_t pad = (fWord >> kbit12) & 0x3f;                                          //pad info in raw word is between bits: 12...17
 243     if(!CheckPad(pad)) continue;
 244     Int_t charge = fWord & 0xfff;
 245 //    if(charge==0) Printf(" CARICA ZERO!!! at %i",fPosition);
 246 //    if(charge==0) return kFALSE;
 247     fPad[fNPads] = GetPad(fDDLNumber,row,dilogic,pad);
 248     fCharge[fNPads] = charge;
 249     fNPads++;
 250   }
 251
 252   cntDilogic -= cnt;
 253   return kTRUE;
 254 }
 255 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 256 Bool_t AliHMPIDRawStream::CheckSegment()
 257 {
 258   UInt_t markSegment = 0xAB0;
 259   /*
 260   if (iRow%8 == 0) {
 261     UInt_t segWord = GetWord();
 262     if ((segWord >> 20) != statusSegWord) {
 263       fRawReader->AddMajorErrorLog(kBadSegWordErr);
 264       AliWarning(Form("Wrong segment word signature: %x, expected 0xab0!",(segWord >> 20)));
 265       fNumOfErr[kBadSegWordErr]++;
 266       return kTRUE;
 267     }
 268 */
 269   UInt_t segMarker = (fWord >> kbit20) & 0xfff;
 270   if (segMarker != markSegment ) {
 271     fRawReader->AddMajorErrorLog(kWrongSegErr,Form("Segment marker %0X wrong (expected %0X) at %i in word %0X!",segMarker,markSegment,fPosition,fWord));
 272     AliWarning(Form("Segment marker %X wrong (expected %0X)! at %i in word %0X!",segMarker,markSegment,fPosition,fWord));
 273     fNumOfErr[kWrongSegErr]++;
 274     return kFALSE;
 275   }
 276
 277   UInt_t segAddress = fWord & 0xff;
 278   if (segAddress<1 ||segAddress>3) {
 279     fRawReader->AddMajorErrorLog(kWrongSegErr,Form("Segment address %d not in the valid range [1-3] at %i in word %0X",segAddress,fPosition,fWord));
 280     AliWarning(Form("Segment address %d not in the valid range [1-3]",segAddress));
 281     fNumOfErr[kWrongSegErr]++;
 282     return kFALSE;
 283   }
 284 //  Printf("Segment Marker found at %i! Number of segment is %i",fPosition,segAddress);
 285   return kTRUE;
 286 }
 287 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 288 Bool_t AliHMPIDRawStream::CheckRow(UInt_t row)
 289 {
 290 //check on row number
 291 //  Printf("ROW %i word %0X",row,fWord);
 292   if(row>=1 && row <=kNRows) return kTRUE;
 293
 294   fRawReader->AddMajorErrorLog(kWrongRowErr,Form("row %d",row));
 295   AliWarning(Form("Wrong row index: %d, expected (1 -> %d) word %0X at %i...",row,kNRows,fWord,fPosition));
 296   fNumOfErr[kWrongRowErr]++;
 297   return kFALSE;
 298 }
 299 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 300 Bool_t AliHMPIDRawStream::CheckDilogic(UInt_t dilogic)
 301 {
 302 //check dilogic number
 303   if (dilogic>= 1 && dilogic <=kNDILOGICAdd) return kTRUE;
 304
 305   fRawReader->AddMajorErrorLog(kWrongDilogicErr,Form("dil %d",dilogic));
 306   AliWarning(Form("Wrong DILOGIC index: %d, expected (1 -> %d)!",dilogic,kNDILOGICAdd));
 307   fNumOfErr[kWrongDilogicErr]++;
 308   //dilogic = iDILOGIC;
 309   return kFALSE;
 310 }
 311 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 312 Bool_t AliHMPIDRawStream::CheckPad(UInt_t pad)
 313 {
 314 //check pad number
 315   if (pad < kNPadAdd) return kTRUE;
 316
 317   fRawReader->AddMajorErrorLog(kWrongPadErr,Form("pad %d",pad));
 318   AliWarning(Form("Wrong pad index: %d, expected (0 -> %d)!",pad,kNPadAdd));
 319   fNumOfErr[kWrongPadErr]++;
 320   return kFALSE;
 321 }
 322 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 323 Bool_t AliHMPIDRawStream::CheckEoE(Int_t &nDil)
 324 {
 325   if (!((fWord >> kbit27) & 0x1)) {                                                //check 27th bit in EoE. It must be 1!
 326     fRawReader->AddMajorErrorLog(kEoEFlagErr);
 327     AliWarning(Form("Missing end-of-event flag! (%08X) at %i",fWord,fPosition));
 328     fNumOfErr[kEoEFlagErr]++;
 329     return kFALSE;
 330   }
 331   nDil = fWord & 0x7f;
 332   if(nDil < 0 || nDil > 48 ) {
 333     fRawReader->AddMajorErrorLog(kEoESizeErr,Form("EoE size=%d",nDil));
 334     AliWarning(Form("Wrong end-of-event word-count: %08X",fWord));
 335     fNumOfErr[kEoESizeErr]++;
 336     return kFALSE;
 337   }
 338 //  UInt_t da = (eOfEvent >> 18) & 0xf;
 339 //  if (cntData!=0 && da != dilogic) {
 340 //    fRawReader->AddMajorErrorLog(kEoEDILOGICErr,Form("eoe dil %d != %d",da,dilogic));
 341 //    AliWarning(Form("Wrong DILOGIC address found in end-of-event: %d, expected %d!",da,dilogic));
 342 //    fNumOfErr[kEoEDILOGICErr]++;
 343 //    return kFALSE;  AliQAChecker::Instance()->Run(AliQA::kHMPID, task, obj) ;
 344
 345 //  }
 346 //  UInt_t ca = (eOfEvent >> 22) & 0x1f;
 347 //  if (cntData!=0 &&  ca != row) {
 348 //    fRawReader->AddMajorErrorLog(kEoERowErr,Form("eoe row %d != %d",ca,row));
 349 //    AliWarning(Form("Wrong row index found in end-of-event: %d, expected %d!",ca,row));
 350 //    fNumOfErr[kEoERowErr]++;
 351 //    return kFALSE;
 352 //  }
 353   return kTRUE;
 354 }
 355 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 356 Bool_t AliHMPIDRawStream::CheckRowMarker()
 357 {
 358   UInt_t nMAXwordsInRow = 0x1EA;
 359   UInt_t statusControlRow = 0x32a8; // 0x36a8 for zero suppression
 360 //First check on row marker
 361   UInt_t rowControlWord = fWord >> kbit0 & 0xfbff;
 362
 363   if(rowControlWord != statusControlRow) {
 364     fRawReader->AddMajorErrorLog(kRowMarkerErr);
 365     AliWarning(Form("Wrong row marker %x expected 0x32a8!",rowControlWord));
 366     fNumOfErr[kRowMarkerErr]++;
 367     return kFALSE;
 368   }
 369 //Second check on row marker
 370    UInt_t wordsInRow = fWord >> kbit16 & 0x0fff;                // Number of words after the row marker
 371
 372   if (wordsInRow > nMAXwordsInRow) {
 373     fRawReader->AddMajorErrorLog(kRowMarkerSizeErr);
 374     AliWarning(Form(" FATAL: Number of words %x in a row exceeds the expected value: 0x1EA !",wordsInRow));
 375     fNumOfErr[kRowMarkerSizeErr]++;
 376     return kFALSE;
 377   }
 378
 379   return kTRUE;
 380 }
 381 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 382 UInt_t AliHMPIDRawStream::GetWord(Int_t n,EDirection dir)
 383 {
 384   // This method returns the n-th 32 bit word
 385   // inside the raw data payload.
 386   // The method is supposed to be endian (platform)
 387   // independent.
 388   if (!fData || fPosition < 0) AliFatal("Raw data payload buffer is not yet initialized !");
 389
 390   UInt_t word = 0;
 391
 392   if(dir==kBwd) n = -n;
 393   fPosition+=4*n-4;
 394
 395   StorePosition();
 396
 397   word |= fData[fPosition++];
 398   word |= fData[fPosition++] << 8;
 399   word |= fData[fPosition++] << 16;
 400   word |= fData[fPosition++] << 24;
 401
 402   return word;
 403 }
 404 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 405 void AliHMPIDRawStream::DumpData(Int_t nw)
 406 {
 407    for(Int_t i=0;i<nw;i+=4) {
 408      if(!(i%16)) printf(" \n %8i) ",i);
 409      printf("%02X%02X%02X%02X [ %06i ] ",fData[i+3],fData[i+2],fData[i+1],fData[i+0],fPos[i]);
 410    }
 411    Printf(" \n -----end of dump ----------- ");
 412 }
 413 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 414 void AliHMPIDRawStream::StorePosition()
 415 {
 416   if(fPos[fPosition]!=0) {
 417 //    Printf("Position already stored!!! Value %i at address %i",fPos[fPosition],fPosition);
 418     return;
 419   }
 420   iPos++;
 421   fPos[fPosition] = iPos;
 422 //  if(stDeb)Printf("%i - Actual position %i",iPos,fPosition);
 423 }