]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - HMPID/AliHMPIDRawStream.cxx
Modifications in CMake* files to export all required headers.
[u/mrichter/AliRoot.git] / HMPID / AliHMPIDRawStream.cxx
index ec7bfff29a7f9ec4d125aac9c2ce4d71ecf173c1..e9d761d56cb4ba2c96902474ffeff3cc9b291025 100644 (file)
 #include "AliRawReader.h"
 #include "AliLog.h"
 
-
 ClassImp(AliHMPIDRawStream)
 
-//_____________________________________________________________________________
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 AliHMPIDRawStream::AliHMPIDRawStream(AliRawReader* rawReader) :
+  fNPads(0),
+  fCharge(0x0),
+  fPad(0x0),
   fDDLNumber(-1),
+  fnDDLInStream(0x0),
+  fnDDLOutStream(0x0),
+  fLDCNumber( 0),
+  fTimeStamp( 0),
   fRawReader(rawReader),
-  fData(NULL),
-  fPosition(-1)
+  fData(0x0),
+  fNumOfErr(0x0),
+  fPosition(-1),
+  fWord(0),
+  fZeroSup(kTRUE),
+  fPos(0x0),
+  fiPos(0),
+  fTurbo(kFALSE),
+  fRawDataSize(0)
 {
+  //
   // Constructor
-  Init();
-
+  //
+  fNumOfErr = new Int_t*[kNDDL];                                 // Store the numner of errors for a given error type and a given DD
+  for(Int_t i=0;i<kNDDL;i++) {
+    fNumOfErr[i] = new Int_t [kSumErr];
+  }
+  
+  fnDDLInStream=new Int_t[kNDDL];
+  fnDDLOutStream=new Int_t[kNDDL];
+  for(Int_t iddl=0;iddl<kNDDL;iddl++) { fnDDLInStream[iddl]=-1;fnDDLOutStream[iddl]=-1;}
+  for(Int_t iddl=0;iddl<kNDDL;iddl++) 
+    for(Int_t ierr=0; ierr < kSumErr; ierr++) fNumOfErr[iddl][ierr]=0;               //reset errors  
   fRawReader->Reset();
   fRawReader->Select("HMPID");
+  
+  
+  
 }
-//-----------------------------------------------------------------------------
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 AliHMPIDRawStream::AliHMPIDRawStream() :
+  fNPads(0),
+  fCharge(0x0),
+  fPad(0x0),
   fDDLNumber(-1),
+  fnDDLInStream(0x0),
+  fnDDLOutStream(0x0),
+  fLDCNumber( 0),
+  fTimeStamp( 0),
   fRawReader(0x0),
-  fData(NULL),
-  fPosition(-1)
+  fData(0x0),
+  fNumOfErr(0x0),  
+  fPosition(-1),
+  fWord(0),
+  fZeroSup(kTRUE),
+  fPos(0x0),
+  fiPos(0),
+  fTurbo(kFALSE) ,
+  fRawDataSize(0)
 {
+  //
   // Constructor
-  Init();
+  //
+  fNumOfErr = new Int_t*[kNDDL];                                 // Store the numner of errors for a given error type and a given DD
+  for(Int_t i=0;i<kNDDL;i++) {
+    fNumOfErr[i] = new Int_t [kSumErr];
+  }
+  fnDDLInStream=new Int_t[kNDDL];
+  fnDDLOutStream=new Int_t[kNDDL];
+  for(Int_t iddl=0;iddl<kNDDL;iddl++) { fnDDLInStream[iddl]=-1;fnDDLOutStream[iddl]=-1;}
+
+  
+  for(Int_t iddl=0;iddl<kNDDL;iddl++) 
+    for(Int_t ierr=0; ierr < kSumErr; ierr++) fNumOfErr[iddl][ierr]=0;               //reset errors  
+
+  
 }
-//_____________________________________________________________________________
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 AliHMPIDRawStream::~AliHMPIDRawStream()
 {
+  //
   // destructor
-}
+  //
+  DelVars();
+  
+  
+  fDDLNumber=0;
+  fLDCNumber=0;
+  fTimeStamp=0;
+  fPosition=0;
+  fWord=0;
+  fZeroSup=0;
+  fTurbo=0;
+  fRawDataSize=0;
+  for(Int_t i=0;i<kNDDL;i++) delete [] fNumOfErr[i]; 
+  delete [] fNumOfErr; 
 
-//_____________________________________________________________________________
-void AliHMPIDRawStream::Init()
-{
-  // Initalize the container
-  // with the pad charges
-  for(Int_t h = 0; h < kNDDL; h++) {  
-    for(Int_t i = 0; i < kNRows; i++){
-      for(Int_t j = 0; j < kNDILOGICAdd; j++){
-       for(Int_t k = 0; k < kNPadAdd; k++){
-         fCharge[h][i][j][k] = -1;
-         fPad[h][i][j][k]=-1;
-       }
-      }
-    }
-  }
-  fZeroSup=kTRUE;
-}//Init()
-//_____________________________________________________________________________
+  delete [] fnDDLInStream;
+  delete [] fnDDLOutStream;
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 void AliHMPIDRawStream::Reset()
 {
   // reset raw stream params
   // Reinitalize the containers
-  Init();
   fDDLNumber = -1;
+  fLDCNumber =  0;
+  fTimeStamp =  0;
   fPosition = -1;
   fData = NULL;
   if (fRawReader) fRawReader->Reset();
 }
-//_____________________________________________________________________________
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::Turbo()
+{
+  
+  Int_t row,dilogic;UInt_t pad;
+  Int_t cntGlob = fRawReader->GetDataSize()/4;
+  fPosition=0;
+  fNPads=0;
+//  Int_t gw=0;
+  for(Int_t i=1;i<cntGlob;i++) {
+    if(!GetWord(1)) return kFALSE;
+    if (((fWord >> kbit27) & 1)) continue;
+    UInt_t statusControlRow = 0x32a8; 
+    UInt_t rowControlWord = fWord >> kbit0 & 0xfbff;
+    if(rowControlWord == statusControlRow) continue;
+
+    row = (fWord >> kbit22) & 0x1f;
+    dilogic = (fWord >> kbit18) & 0xf;                                              //dilogic info in raw word is between bits: 18...21
+    
+    pad = (fWord >> kbit12) & 0x3f;                                          //pad info in raw word is between bits: 12...17
+    if(!CheckPad(pad)) continue;
+    Int_t charge = fWord & 0xfff;
+    if(GetPad(fDDLNumber,row,dilogic,pad)<0) continue;
+    fPad[fNPads] = GetPad(fDDLNumber,row,dilogic,pad);
+    fCharge[fNPads] = charge; 
+    fNPads++;
+    if(charge==0) fNumOfErr[fDDLNumber][kPedQZero]++;
+  }//word loop
+  //Printf("Size: %i  DDL %i row %i dilogic %i pad %i fPos %i fNPads: %i Charge: %d Word %4.4x GoodW: %i",cntGlob,fDDLNumber,row,dilogic,pad,fPosition,fNPads,fCharge[fNPads-1],fWord,gw++);
+  return kTRUE;
+}//Turbo()
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 Bool_t AliHMPIDRawStream::Next()
 {
   // read next DDL raw data from the HMPID raw data stream
@@ -94,220 +181,441 @@ Bool_t AliHMPIDRawStream::Next()
   do {
     if (!fRawReader->ReadNextData(fData)) return kFALSE;
   } while (fRawReader->GetDataSize() == 0);
+  Int_t runNumber = fRawReader->GetRunNumber();
+  
+  Int_t ddlArray[] = {0,1,2,3,4,5,6,7,8,9,10,11,12,13}; 
+
+  // next line: fix for link chamber2 cables inversion for period LHC11d, LHC11e, LHC11f and LHC11h. 
   
-  if ( fRawReader->GetDataSize() > 47148) {
-    AliWarning(Form("Raw data event size is larger (%d) than possible for HMPID!!!! ",fRawReader->GetDataSize()));
+  if(runNumber>=156620 && runNumber<=170593) {ddlArray[4] = 5; ddlArray[5] = 4;}  
+  else {ddlArray[4] = 4; ddlArray[5] = 5;}
+     
+  /*
+  Event type is selected as in $ALICE_ROOT/RAW/event.h  
+  #define START_OF_RUN                    ((eventTypeType) 1)
+  #define END_OF_RUN                      ((eventTypeType) 2)
+  #define START_OF_RUN_FILES              ((eventTypeType) 3)
+  #define END_OF_RUN_FILES                ((eventTypeType) 4)
+  #define START_OF_BURST                  ((eventTypeType) 5)
+  #define END_OF_BURST                    ((eventTypeType) 6)
+  #define PHYSICS_EVENT                   ((eventTypeType) 7) <<---------------  
+  #define CALIBRATION_EVENT               ((eventTypeType) 8)
+  #define EVENT_FORMAT_ERROR              ((eventTypeType) 9)
+  #define START_OF_DATA                   ((eventTypeType)10)
+  #define END_OF_DATA                     ((eventTypeType)11)
+  #define SYSTEM_SOFTWARE_TRIGGER_EVENT   ((eventTypeType)12)
+  #define DETECTOR_SOFTWARE_TRIGGER_EVENT ((eventTypeType)13)
+  #define EVENT_TYPE_MIN                  1
+  #define EVENT_TYPE_MAX                  13 
+  */
+
+  fPosition = 0;
+  Bool_t status=kFALSE;
+  fRawDataSize=0;        
+  fDDLNumber = -1;
+  if (fRawReader->GetDDLID()>=0)
+    fDDLNumber = ddlArray[fRawReader->GetDDLID()];
+  
+  if(fDDLNumber<0) {
+    AliWarning(Form("fDDLNumber not a acceptable value %i",fDDLNumber));
+    return kFALSE;
   }
+    
+  if(fRawReader->GetType() == 7 || fRawReader->GetType() == 8 )  {           //New: Select Physics events, Old: Raw data size is not 0 and not 47148 (pedestal)
+    fnDDLInStream[fDDLNumber]=1; fnDDLOutStream[fDDLNumber]=0;
+    
+    fLDCNumber = fRawReader->GetLDCId();
+    fTimeStamp = fRawReader->GetTimestamp();
+    
+    AliDebug(1,Form("DDL %i started to be decoded!",fDDLNumber));
+    fRawDataSize=fRawReader->GetDataSize()/4;
+    DelVars();                                         //We have to delete the variables initialized in the InitVars before recall IntiVars!!!
+    InitVars(fRawDataSize);                            //To read the charge and pads we cannot delete before the status return
+    
+    if(fTurbo==kTRUE) status=Turbo();
+    else status = ReadHMPIDRawData();
+   
+    if(status) AliDebug(1,Form("Event DDL %i successfully decoded!.",fDDLNumber));
+    else AliDebug(1,Form("Event DDL %i ERROR in decoding!.",fDDLNumber));
+    //DumpData(fRawReader->GetDataSize());
+  
+   
+    }
+    if(status==kTRUE) {fnDDLOutStream[fDDLNumber]++; }//Printf("fnDDLOutStream[%d]=%d",fDDLNumber,fnDDLOutStream[fDDLNumber]); } //Count the number of events when the DDL was succesfully decoded
+   
+//    return status;  // temporary solution...
+   return kTRUE;
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+void AliHMPIDRawStream::InitVars(Int_t n)
+{
+  //
+  //
+  //
+  fNPads = 0;
+  fCharge = new Int_t[n]; 
+  fPad = new Int_t[n];
+  //for debug purpose
+  fPos = new Int_t[4*n+4];                     //reset debug
+  for(Int_t ie = 0 ; ie < 4*n+4; ie++) fPos[ie] = 0; //initialize for 0, otherwise the position is considered filled and will not be updated for the dump
+  fiPos = 0;
+}    
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+void AliHMPIDRawStream::DelVars()
+{
+  //
+  //
+  // 
+  //Clean the initvars!!!!!!!!
+  fNPads = 0; 
+  if (fCharge)     { delete [] fCharge;    fCharge = 0x0; }
+  if (fPad)        { delete [] fPad;       fPad = 0x0;       }   
+  if (fPos)        { delete [] fPos;       fPos = 0x0;       }     
+  fiPos=0;
+  
+     
+}    
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::ReadHMPIDRawData()
+{
+  //Here the loop on the decoding the raw bank 
+  //for one ddl starts.
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
+  Int_t cntGlob = fRawReader->GetDataSize()/4;
+  if(cntGlob==0) {fNumOfErr[fDDLNumber][kRawDataSizeErr]++; return kFALSE; }
+  
+  Int_t cnt = cntGlob;
+  Int_t nwSeg;
+  Int_t cntSegment;
+
+  if(!GetWord(cnt)) return kFALSE;
+  cnt--;
+
+  
+  while (cnt>20) { //counter limit is changed from 0 to 20 to take into account (skipp) the 5 extra words in the equipment header
+    nwSeg = (fWord >> kbit8) & 0xfff;
+    if(!CheckSegment()) return kFALSE;
+    if(!ReadSegment(cntSegment)) return kFALSE;
+
+    if(nwSeg != cntSegment) {AliDebug(1,Form("Error in Segment counters: %i different wrt %i",nwSeg,cntSegment)); return kFALSE;}
+    if(!GetWord(cntSegment+1,kBwd)) return kFALSE;
+    cnt-=cntSegment+1;
+  }
+  
+  return kTRUE;
+  
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::ReadSegment(Int_t &cntSegment)
+{
+  //Read the segment
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
+  cntSegment = (fWord >> kbit8) & 0xfff;
+  Int_t cnt = cntSegment;
+  Int_t cntRow;
+  Int_t nwRow;
+
+  if(!GetWord(cnt,kBwd)) return kFALSE;
+  
+  while (cnt>0) {
+
+    cntRow  = (fWord >> kbit16) & 0xfff;
+    if(!CheckRowMarker()) return kFALSE;
+    if(!ReadRow(nwRow)) return kFALSE;
+
+    if(nwRow != cntRow) {AliDebug(1,Form("Error in Row counters: %i different wrt %i",nwRow,cntRow)); return kFALSE;}
+    if(!GetWord(cntRow+1)) return kFALSE;
+    cnt -= cntRow+1;
+    
+  }
+
+  cntSegment -= cnt;
+  
+  return kTRUE;
+    
+}    
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::ReadRow(Int_t &cntRow)
+{
+  // Read the row
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
+
+  Int_t cnt;
+  Int_t cntDilogic;
+  Int_t nwDil;
   
-  else if (fRawReader->GetDataSize() == 47148) {
-    fDDLNumber = fRawReader->GetDDLID();
-    fPosition = 0;
+  cntRow  = (fWord >> kbit16) & 0xfff;
+  cnt = cntRow;  
+  
+  if(!GetWord(cntRow)) return kFALSE;
+  
+  while (cnt>0) {
     
-    Init();
+    if(!CheckEoE(nwDil)) return kFALSE;
+    if(!ReadDilogic(cntDilogic)) return kFALSE;
     
-    // Look over rows
-    for(Int_t iRow = 1; iRow <= kNRows; iRow++) {
-      // Read row marker
-      UInt_t rowMarker = GetNextWord() & 0x1ffffff;
-      if (rowMarker != 0x1ea32a8) {
-       fRawReader->AddMajorErrorLog(kRowMarkerErr);
-       AliWarning(Form("Wrong row marker %x for row %d, expected 0x1ea32a8!",rowMarker,iRow));
-       return kTRUE;
-      }//check for row marker
-      UInt_t dilogic = 0, row = 0;
-      for(Int_t iDILOGIC = 1; iDILOGIC <= kNDILOGICAdd; iDILOGIC++) {
-       // Read pad charges
-       for(Int_t iPad = 0; iPad < kNPadAdd; iPad++) {
-         UInt_t data = GetNextWord();
-         row = (data >> 22) & 0x1f;                                                      //row information in raw word is between bits: 22...26
-         if (row < 1 || row > kNRows) {
-           fRawReader->AddMajorErrorLog(kWrongRowErr,Form("row %d",row));
-           AliWarning(Form("Wrong row index: %d, expected (1 -> %d)!",row,kNRows));
-           row = iRow;
-         }
-         dilogic = (data >> 18) & 0xf;                                                   //dilogic info in raw word is between bits: 18...21
-         if (dilogic < 1 || dilogic > kNDILOGICAdd) {
-           fRawReader->AddMajorErrorLog(kWrongDilogicErr,Form("dil %d",dilogic));
-           AliWarning(Form("Wrong DILOGIC index: %d, expected (1 -> %d)!",dilogic,kNDILOGICAdd));
-           dilogic = iDILOGIC;
-         }
-         UInt_t pad = (data >> 12) & 0x3f;                                               //pad info in raw word is between bits: 12...17
-         if (pad >= kNPadAdd) {
-           fRawReader->AddMajorErrorLog(kWrongPadErr,Form("pad %d",pad));
-           AliWarning(Form("Wrong pad index: %d, expected (0 -> %d)!",pad,kNPadAdd));
-           pad = iPad;
-         }
-         fCharge[fDDLNumber][row][dilogic][pad] = data & 0xfff;
-       }
-       // Now read the end-of-event word
-       UInt_t eOfEvent = GetNextWord() & 0xfffffff;
-       if (!((eOfEvent >> 27) & 0x1)) {                                                  //check 27th bit in EoE. It must be 1!
-         fRawReader->AddMajorErrorLog(kEoEFlagErr);
-         AliWarning(Form("Missing end-of-event flag! (%x)",eOfEvent));
-         return kTRUE;
-       }
-       UInt_t wc = eOfEvent & 0x7f;
-       if (wc != 48) {
-         fRawReader->AddMajorErrorLog(kEoESizeErr,Form("eoe size=%d",wc));
-         AliWarning(Form("Wrong end-of-event word-count:%d, expected 48!",wc));
-         return kTRUE;
-       }
-       UInt_t da = (eOfEvent >> 18) & 0xf;
-       if (da != dilogic) {
-         fRawReader->AddMajorErrorLog(kEoEDILOGICErr,Form("eoe dil %d != %d",da,dilogic));
-         AliWarning(Form("Wrong DILOGIC address found in end-of-event: %d, expected %d!",da,dilogic));
-         return kTRUE;
-       }
-       UInt_t ca = (eOfEvent >> 22) & 0x1f;
-       if (ca != row) {
-         fRawReader->AddMajorErrorLog(kEoERowErr,Form("eoe row %d != %d",ca,row));
-         AliWarning(Form("Wrong row index found in end-of-event: %d, expected %d!",ca,row));
-         return kTRUE;
-       }
-      }//DILOGIC loop
-      
-      // Read the segment marker
-      // One maker per 8 rows
-      
-      if (iRow%8 == 0) {
-       UInt_t segWord = GetNextWord();
-       if ((segWord >> 8) != 0xab0f59) {
-         fRawReader->AddMajorErrorLog(kBadSegWordErr);
-         AliWarning(Form("Wrong segment word signature: %x, expected 0xab0f59!",(segWord >> 8)));
-         return kTRUE;
-       }
-       
-       if ((segWord & 0xff) != (((UInt_t)iRow + 7) / 8)) {
-         fRawReader->AddMajorErrorLog(kWrongSegErr,Form("seg %d != %d",segWord & 0xff,(iRow + 7) / 8));
-         AliWarning(Form("Segment index (%d) does not correspond to the one expected from row index (%d)!",segWord & 0xff,(iRow + 7) / 8));
-         return kTRUE;
-       }
-      }
-    }//loop of Row   
-  }//Pedestal files selected by data lenght
+    if(nwDil != cntDilogic) {AliDebug(1,Form("Error in Dilogic counters: %i different wrt %i",nwDil,cntDilogic));return kFALSE;}
+    cnt -= cntDilogic;
+    if(!GetWord(1,kBwd)) return kFALSE; // go to next Dilogic bank...
+    cnt--;
+//    Printf(" cnt %i cntDilogic %i ",cnt,cntDilogic);
+  }
+  
+  cntRow -= cnt;  
+  
+  return kTRUE;
+  
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::ReadDilogic(Int_t &cntDilogic)
+{
+  // Read the dilogic bank
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
 
-  else {                                                                  //Raw data size is not 0 and not 47148 (pedestal)
-    fDDLNumber = fRawReader->GetDDLID();
-    fPosition = 0;
+  cntDilogic = fWord & 0x7f;
+  
+  Int_t cnt = cntDilogic;
   
-    Init();
+//  Printf(" cnt DILOGIC %i at %i word %08X",cnt,fPosition,fWord);
+
+  for(Int_t iDil=0;iDil<cntDilogic;iDil++) {
+    UInt_t dilogic = 0, row = 0;
+    if(!GetWord(1,kBwd)) return kFALSE;
+//check on row number      
+    cnt--;
+    row = (fWord >> kbit22) & 0x1f;
+    if(!CheckRow(row)) continue;
+//check dilogic number     
+    dilogic = (fWord >> kbit18) & 0xf;                                              //dilogic info in raw word is between bits: 18...21
+    if(!CheckDilogic(dilogic)) continue;
+//check pad number
+    UInt_t pad = (fWord >> kbit12) & 0x3f;                                          //pad info in raw word is between bits: 12...17
+    if(!CheckPad(pad)) continue;
+    Int_t charge = fWord & 0xfff;
+    if(GetPad(fDDLNumber,row,dilogic,pad)<0) continue;
+    fPad[fNPads] = GetPad(fDDLNumber,row,dilogic,pad);
+    fCharge[fNPads] = charge; 
+    fNPads++;
     
-    // Look over rows
-    for(Int_t iRow = 1; iRow <= kNRows; iRow++) {
-      
-      UInt_t rowMarker = GetNextWord();                                 // Read row marker
-      
-      Int_t numRows= rowMarker >> 16 & 0xffffff;
-      
-      if ((rowMarker >> 0 & 0xffff) != 0x32a8) {
-       fRawReader->AddMajorErrorLog(kRowMarkerErr);
-        AliWarning(Form("Wrong row marker %x for row %d, expected 0x32a8!",rowMarker,iRow));
-       return kTRUE;
-      }//check for row marker
-      UInt_t dilogic = 0, row = 0;
-      UInt_t cntData=0;
-      UInt_t cntEoE=0;
-      //molnarl: lets read how many rows we have from the marker; 10 dilogic EoE so we can get the number of data words they should be in order
-      for(Int_t iWordInRow=0; iWordInRow<numRows;iWordInRow++) //loop over words 
-       {     
-         UInt_t tmpword=GetNextWord();               
-         UInt_t eOfEvent = tmpword;                     // always assume that it is an EoE. If bit 
-         if (!((eOfEvent >> 27) & 0x1)) {               // if it is not EoE then data!
-           UInt_t data=tmpword;
-           row = (data >> 22) & 0x1f;                                                  //row information in raw word is between bits: 22...26
-           if (row < 1 || row > kNRows) {                                              //select bits from 22 and with 0x1f ask for the next 5 bits 
-             fRawReader->AddMajorErrorLog(kWrongRowErr,Form("row %d",row));
-             AliWarning(Form("Wrong row index: %d, expected (1 -> %d)!",row,kNRows));
-             // row = iRow;
-           }
-           dilogic = (data >> 18) & 0xf;                                              //dilogic info in raw word is between bits: 18...21
-           if (dilogic < 1 || dilogic > kNDILOGICAdd) {
-             fRawReader->AddMajorErrorLog(kWrongDilogicErr,Form("dil %d",dilogic));
-             AliWarning(Form("Wrong DILOGIC index: %d, expected (1 -> %d)!",dilogic,kNDILOGICAdd));
-             //dilogic = iDILOGIC;
-           }
-           UInt_t pad = (data >> 12) & 0x3f;                                                    //pad info in raw word is between bits: 12...17
-           if (pad >= kNPadAdd) {
-             fRawReader->AddMajorErrorLog(kWrongPadErr,Form("pad %d",pad));
-             AliWarning(Form("Wrong pad index: %d, expected (0 -> %d)!",pad,kNPadAdd));
-             //pad = iPad;
-           }
-           fCharge[fDDLNumber][row][dilogic][pad] = data & 0xfff;  cntData++;
-           
-         }//not EoE but data!
-         //if it is EoE
-         else{
-           //Printf("EoE word");
-           if (!((eOfEvent >> 27) & 0x1)) {                                                  //check 27th bit in EoE. It must be 1!
-             fRawReader->AddMajorErrorLog(kEoEFlagErr);
-             AliWarning(Form("Missing end-of-event flag! (%x)",eOfEvent));
-             return kTRUE;
-           }
-           UInt_t wc = eOfEvent & 0x7f;
-           if (wc != cntData) {
-             fRawReader->AddMajorErrorLog(kEoESizeErr,Form("eoe size=%d",wc));
-             AliWarning(Form("Wrong end-of-event word-count:%d, expected 48!",wc));
-             return kTRUE;
-           }
-           UInt_t da = (eOfEvent >> 18) & 0xf;
-           if (cntData!=0 && da != dilogic) {
-             fRawReader->AddMajorErrorLog(kEoEDILOGICErr,Form("eoe dil %d != %d",da,dilogic));
-             AliWarning(Form("Wrong DILOGIC address found in end-of-event: %d, expected %d!",da,dilogic));
-             return kTRUE;
-           }
-           UInt_t ca = (eOfEvent >> 22) & 0x1f;
-           if (cntData!=0 &&  ca != row) {
-             fRawReader->AddMajorErrorLog(kEoERowErr,Form("eoe row %d != %d",ca,row));
-             AliWarning(Form("Wrong row index found in end-of-event: %d, expected %d!",ca,row));
-             return kTRUE;
-           }
-           cntData=0;  cntEoE++;//zero it and wait for new data words
-         } //EoE word
-         
-       } //loop over rows within two markers
-      
-      
-      
-      // Read the segment marker
-      // One maker per 8 rows
-      
-      if (iRow%8 == 0) {
-       UInt_t segWord = GetNextWord();
-       if ((segWord >> 8) != 0xab0f59) {
-         fRawReader->AddMajorErrorLog(kBadSegWordErr);
-         AliWarning(Form("Wrong segment word signature: %x, expected 0xab0f59!",(segWord >> 8)));
-         return kTRUE;
-       }
-       
-       if ((segWord & 0xff) != (((UInt_t)iRow + 7) / 8)) {
-         fRawReader->AddMajorErrorLog(kWrongSegErr,Form("seg %d != %d",segWord & 0xff,(iRow + 7) / 8));
-         AliWarning(Form("Segment index (%d) does not correspond to the one expected from row index (%d)!",segWord & 0xff,(iRow + 7) / 8));
-         return kTRUE;
-       }
-      }
-    }//loop of Row   
-  }//REal data files selected by data lenght
-  
-  AliDebug(1,"Stop."); 
+    if(charge==0) 
+    {
+      AliDebug(1,Form("If PEDESTAL run -> WARNING: ZERO charge is read from DDL: %d row: %d dil: %d pad: %d",fDDLNumber,row,dilogic,pad));
+      fNumOfErr[fDDLNumber][kPedQZero]++;
+    }
+    
+  }//iDil
+
+  cntDilogic -= cnt;  
+  return kTRUE;  
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::CheckSegment()
+{
+  // Check the segment marker
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
+
+  UInt_t markSegment = 0xAB0;
+  /*
+  if (iRow%8 == 0) {
+    UInt_t segWord = GetWord();          
+    if ((segWord >> 20) != statusSegWord) {
+      fRawReader->AddMajorErrorLog(kBadSegWordErr);
+      AliDebug(1,Form("Wrong segment word signature: %x, expected 0xab0!",(segWord >> 20)));
+      fNumOfErr[kBadSegWordErr]++;
+      return kTRUE;
+    }
+*/
+  UInt_t segMarker = (fWord >> kbit20) & 0xfff;
+  if (segMarker != markSegment ) {
+    //fRawReader->AddMajorErrorLog(kWrongSegErr,Form("Segment marker %0X wrong (expected %0X) at %i in word %0X!",segMarker,markSegment,fPosition,fWord));
+    AliDebug(1,Form("Segment marker %X wrong (expected %0X)! at %i in word %0X!",segMarker,markSegment,fPosition,fWord));
+    fNumOfErr[fDDLNumber][kWrongSegErr]++;
+    return kFALSE;
+  }
+  
+  UInt_t segAddress = fWord & 0xff;
+  if (segAddress<1 ||segAddress>3) {
+    //fRawReader->AddMajorErrorLog(kWrongSegErr,Form("Segment address %d not in the valid range [1-3] at %i in word %0X",segAddress,fPosition,fWord));
+    AliDebug(1,Form("Segment address %d not in the valid range [1-3]",segAddress));
+    fNumOfErr[fDDLNumber][kWrongSegErr]++;
+    return kFALSE;
+  }
+//  Printf("Segment Marker found at %i! Number of segment is %i",fPosition,segAddress);
   return kTRUE;
 }
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::CheckRow(UInt_t row)
+{
+  //check on row number      
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
+//  Printf("ROW %i word %0X",row,fWord);
+  if(row>=1 && row <=kNRows) return kTRUE;
+  
+  //fRawReader->AddMajorErrorLog(kWrongRowErr,Form("row %d",row));
+  AliDebug(1,Form("Wrong row index: %d, expected (1 -> %d) word %0X at %i...",row,kNRows,fWord,fPosition));
+  fNumOfErr[fDDLNumber][kWrongRowErr]++;
+  return kFALSE;
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::CheckDilogic(UInt_t dilogic)
+{
+//check dilogic number     
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
+  if (dilogic>= 1 && dilogic <=kNDILOGICAdd) return kTRUE;
 
-//_____________________________________________________________________________
-UInt_t AliHMPIDRawStream::GetNextWord()
+  //fRawReader->AddMajorErrorLog(kWrongDilogicErr,Form("dil %d",dilogic));
+  AliDebug(1,Form("Wrong DILOGIC index: %d, expected (1 -> %d)!",dilogic,kNDILOGICAdd));
+  fNumOfErr[fDDLNumber][kWrongDilogicErr]++;
+  //dilogic = iDILOGIC;
+  return kFALSE;
+}      
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::CheckPad(UInt_t pad)
 {
-  // This method returns the next 32 bit word
+//check pad number     
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
+  if (pad < kNPadAdd) return kTRUE;
+  
+  //fRawReader->AddMajorErrorLog(kWrongPadErr,Form("pad %d",pad));
+  AliDebug(1,Form("Wrong pad index: %d, expected (0 -> %d)!",pad,kNPadAdd));
+  fNumOfErr[fDDLNumber][kWrongPadErr]++;
+  return kFALSE;
+}    
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::CheckEoE(Int_t &nDil)
+{
+  //check the End of Event
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
+  if (!((fWord >> kbit27) & 0x1)) {                                                //check 27th bit in EoE. It must be 1!
+    //fRawReader->AddMajorErrorLog(kEoEFlagErr);
+    AliDebug(1,Form("Missing end-of-event flag! (%08X) at %i",fWord,fPosition));
+    fNumOfErr[fDDLNumber][kEoEFlagErr]++;
+    return kFALSE;
+  }
+  nDil = fWord & 0x7f;                                               //nDil=EoE word count
+  if(nDil < 0 || nDil > 48 ) { 
+
+    //fRawReader->AddMajorErrorLog(kEoESizeErr,Form("EoE size=%d",nDil));
+    AliDebug(1,Form("Wrong end-of-event word-count: %08X",fWord));
+    fNumOfErr[fDDLNumber][kEoESizeErr]++;
+    return kFALSE;
+  }
+//  UInt_t da = (eOfEvent >> 18) & 0xf;
+//  if (cntData!=0 && da != dilogic) {
+//    fRawReader->AddMajorErrorLog(kEoEDILOGICErr,Form("eoe dil %d != %d",da,dilogic));
+//    AliDebug(1,Form("Wrong DILOGIC address found in end-of-event: %d, expected %d!",da,dilogic));
+//    fNumOfErr[kEoEDILOGICErr]++;
+//    return kFALSE;  AliQAChecker::Instance()->Run(AliQAv1::kHMPID, task, obj) ;  
+
+//  }
+//  UInt_t ca = (eOfEvent >> 22) & 0x1f;
+//  if (cntData!=0 &&  ca != row) {
+//    fRawReader->AddMajorErrorLog(kEoERowErr,Form("eoe row %d != %d",ca,row));
+//    AliDebug(1,Form("Wrong row index found in end-of-event: %d, expected %d!",ca,row));
+//    fNumOfErr[kEoERowErr]++;
+//    return kFALSE;
+//  }
+  return kTRUE;
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::CheckRowMarker()
+{
+  //check the row marker
+  //It returns: kFALSE if any error occurs
+  //            kTRUE  if all OK
+  UInt_t nMAXwordsInRow = 0x1EA;
+  UInt_t statusControlRow = 0x32a8; // 0x36a8 for zero suppression
+//First check on row marker    
+  UInt_t rowControlWord = fWord >> kbit0 & 0xfbff;
+  
+  if(rowControlWord != statusControlRow) {
+    //fRawReader->AddMajorErrorLog(kRowMarkerErr);
+    AliDebug(1,Form("Wrong row marker %x expected 0x32a8!",rowControlWord));
+              fNumOfErr[fDDLNumber][kRowMarkerErr]++;
+              return kFALSE;
+   }
+//Second check on row marker    
+   UInt_t wordsInRow = fWord >> kbit16 & 0x0fff;                // Number of words after the row marker, bit 10 is skipped in this check
+  
+  if (wordsInRow > nMAXwordsInRow) {
+    //fRawReader->AddMajorErrorLog(kRowMarkerSizeErr);
+    AliDebug(1,Form(" FATAL: Number of words %x in a row exceeds the expected value: 0x1EA !",wordsInRow));
+    fNumOfErr[fDDLNumber][kRowMarkerSizeErr]++;
+    return kFALSE;
+  }
+  
+  return kTRUE;
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Bool_t AliHMPIDRawStream::GetWord(Int_t n,EDirection dir)
+{
+  // This method returns the n-th 32 bit word
   // inside the raw data payload.
   // The method is supposed to be endian (platform)
   // independent.
-  if (!fData || fPosition < 0) AliFatal("Raw data payload buffer is not yet initialized !");
   
-  UInt_t word = 0;
-  word |= fData[fPosition++];
-  word |= fData[fPosition++] << 8;
-  word |= fData[fPosition++] << 16;
-  word |= fData[fPosition++] << 24;
+  fWord = 0;
+  if (fPosition < 0) {
+    AliError("fPosition < 0 !!! Event skipped.");
+    fRawReader->AddMajorErrorLog(kRawDataSizeErr,"fPosition<0");
+    return kFALSE;
+  }
+
+  if(dir==kBwd) n = -n; 
+  fPosition+=4*n-4;
+
+  if(fPosition==-4) return kTRUE;
   
-  return word;
-}
+  if(fPosition<0 || fPosition > fRawReader->GetDataSize()) {
+    AliWarning(Form("fPosition out of boundaries %i",fPosition));
+    return kFALSE;
+  }
+    
+  StorePosition();
+  
+  fWord |= fData[fPosition++];
+  fWord |= fData[fPosition++] << 8;
+  fWord |= fData[fPosition++] << 16;
+  fWord |= fData[fPosition++] << 24;
 
-//_____________________________________________________________________________
+  return kTRUE;
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+void AliHMPIDRawStream::DumpData(Int_t nw)
+{
+  //just a simple raw data dump
+  // in () is the position in bytes
+  //--
+   for(Int_t i=0;i<nw;i+=4) {
+     if(!(i%16)) printf(" \n %8i) ",i);
+     printf("%02X%02X%02X%02X [ %06i ] ",fData[i+3],fData[i+2],fData[i+1],fData[i+0],fPos[i]);
+   }
+   Printf(" \n -----end of dump ----------- ");
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+void AliHMPIDRawStream::StorePosition()
+{
+  //just for debug purpose
+  // it stores the position
+  //read for the first time
+//  Printf("@@@@@@@@@ fPos: %x fPosition: %d",fPos,fPosition);
+  if(fPos[fPosition]!=0) {
+//    Printf("Position already stored!!! Value %i at address %i",fPos[fPosition],fPosition); 
+    return;
+  }
+  fiPos++;
+  fPos[fPosition] = fiPos;
+//  if(stDeb)Printf("%i - Actual position %i",iPos,fPosition); 
+}
+//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++