Removing centrality check from default:
[u/mrichter/AliRoot.git] / 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 ClassImp(AliHMPIDRawStream)
29
30 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
31 AliHMPIDRawStream::AliHMPIDRawStream(AliRawReader* rawReader) :
32   fNPads(0),
33   fCharge(0x0),
34   fPad(0x0),
35   fDDLNumber(-1),
36   fnDDLInStream(0x0),
37   fnDDLOutStream(0x0),
38   fLDCNumber( 0),
39   fTimeStamp( 0),
40   fRawReader(rawReader),
41   fData(0x0),
42   fNumOfErr(0x0),
43   fPosition(-1),
44   fWord(0),
45   fZeroSup(kTRUE),
46   fPos(0x0),
47   fiPos(0),
48   fTurbo(kFALSE),
49   fRawDataSize(0)
50 {
51   //
52   // Constructor
53   //
54   fNumOfErr = new Int_t*[kNDDL];                                 // Store the numner of errors for a given error type and a given DD
55   for(Int_t i=0;i<kNDDL;i++) {
56     fNumOfErr[i] = new Int_t [kSumErr];
57   }
58   
59   fnDDLInStream=new Int_t[kNDDL];
60   fnDDLOutStream=new Int_t[kNDDL];
61   for(Int_t iddl=0;iddl<kNDDL;iddl++) { fnDDLInStream[iddl]=-1;fnDDLOutStream[iddl]=-1;}
62  
63   for(Int_t iddl=0;iddl<kNDDL;iddl++) 
64     for(Int_t ierr=0; ierr < kSumErr; ierr++) fNumOfErr[iddl][ierr]=0;               //reset errors  
65   fRawReader->Reset();
66   fRawReader->Select("HMPID");
67   
68   
69   
70 }
71 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
72 AliHMPIDRawStream::AliHMPIDRawStream() :
73   fNPads(0),
74   fCharge(0x0),
75   fPad(0x0),
76   fDDLNumber(-1),
77   fnDDLInStream(0x0),
78   fnDDLOutStream(0x0),
79   fLDCNumber( 0),
80   fTimeStamp( 0),
81   fRawReader(0x0),
82   fData(0x0),
83   fNumOfErr(0x0),  
84   fPosition(-1),
85   fWord(0),
86   fZeroSup(kTRUE),
87   fPos(0x0),
88   fiPos(0),
89   fTurbo(kFALSE) ,
90   fRawDataSize(0)
91 {
92   //
93   // Constructor
94   //
95   fNumOfErr = new Int_t*[kNDDL];                                 // Store the numner of errors for a given error type and a given DD
96   for(Int_t i=0;i<kNDDL;i++) {
97     fNumOfErr[i] = new Int_t [kSumErr];
98   }
99   fnDDLInStream=new Int_t[kNDDL];
100   fnDDLOutStream=new Int_t[kNDDL];
101   for(Int_t iddl=0;iddl<kNDDL;iddl++) { fnDDLInStream[iddl]=-1;fnDDLOutStream[iddl]=-1;}
102
103   
104   for(Int_t iddl=0;iddl<kNDDL;iddl++) 
105     for(Int_t ierr=0; ierr < kSumErr; ierr++) fNumOfErr[iddl][ierr]=0;               //reset errors  
106
107   
108 }
109 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
110 AliHMPIDRawStream::~AliHMPIDRawStream()
111 {
112   //
113   // destructor
114   //
115   DelVars();
116   
117   
118   fDDLNumber=0;
119   fLDCNumber=0;
120   fTimeStamp=0;
121   fPosition=0;
122   fWord=0;
123   fZeroSup=0;
124   fTurbo=0;
125   fRawDataSize=0;
126   for(Int_t i=0;i<kNDDL;i++) delete [] fNumOfErr[i]; 
127   delete [] fNumOfErr; 
128
129   delete [] fnDDLInStream;
130   delete [] fnDDLOutStream;
131 }
132 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
133 void AliHMPIDRawStream::Reset()
134 {
135   // reset raw stream params
136   // Reinitalize the containers
137   fDDLNumber = -1;
138   fLDCNumber =  0;
139   fTimeStamp =  0;
140   fPosition = -1;
141   fData = NULL;
142   if (fRawReader) fRawReader->Reset();
143 }
144 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
145 Bool_t AliHMPIDRawStream::Turbo()
146 {
147   
148   Int_t row,dilogic;UInt_t pad;
149   Int_t cntGlob = fRawReader->GetDataSize()/4;
150   fPosition=0;
151   fNPads=0;
152 //  Int_t gw=0;
153   for(Int_t i=1;i<cntGlob;i++) {
154     if(!GetWord(1)) return kFALSE;
155     if (((fWord >> kbit27) & 1)) continue;
156     UInt_t statusControlRow = 0x32a8; 
157     UInt_t rowControlWord = fWord >> kbit0 & 0xfbff;
158     if(rowControlWord == statusControlRow) continue;
159
160     row = (fWord >> kbit22) & 0x1f;
161     dilogic = (fWord >> kbit18) & 0xf;                                              //dilogic info in raw word is between bits: 18...21
162     
163     pad = (fWord >> kbit12) & 0x3f;                                          //pad info in raw word is between bits: 12...17
164     if(!CheckPad(pad)) continue;
165     Int_t charge = fWord & 0xfff;
166     if(GetPad(fDDLNumber,row,dilogic,pad)<0) continue;
167     fPad[fNPads] = GetPad(fDDLNumber,row,dilogic,pad);
168     fCharge[fNPads] = charge; 
169     fNPads++;
170     if(charge==0) fNumOfErr[fDDLNumber][kPedQZero]++;
171   }//word loop
172   //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++);
173   return kTRUE;
174 }//Turbo()
175 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
176 Bool_t AliHMPIDRawStream::Next()
177 {
178   // read next DDL raw data from the HMPID raw data stream
179   // return kFALSE in case of error or no data left
180   AliDebug(1,"Start.");
181   do {
182     if (!fRawReader->ReadNextData(fData)) return kFALSE;
183   } while (fRawReader->GetDataSize() == 0);
184  
185   Int_t runNumber = fRawReader->GetRunNumber();
186   
187   Int_t ddlArray[] = {0,1,2,3,4,5,6,7,8,9,10,11,12,13}; 
188
189   // next line: fix for link chamber2 cables inversion for period LHC11d, LHC11e, LHC11f and LHC11h. 
190   
191   if(runNumber>=156620 && runNumber<=170593) {ddlArray[4] = 5; ddlArray[5] = 4;}  
192   else {ddlArray[4] = 4; ddlArray[5] = 5;}
193      
194   /*
195   Event type is selected as in $ALICE_ROOT/RAW/event.h  
196   #define START_OF_RUN                    ((eventTypeType) 1)
197   #define END_OF_RUN                      ((eventTypeType) 2)
198   #define START_OF_RUN_FILES              ((eventTypeType) 3)
199   #define END_OF_RUN_FILES                ((eventTypeType) 4)
200   #define START_OF_BURST                  ((eventTypeType) 5)
201   #define END_OF_BURST                    ((eventTypeType) 6)
202   #define PHYSICS_EVENT                   ((eventTypeType) 7) <<---------------  
203   #define CALIBRATION_EVENT               ((eventTypeType) 8)
204   #define EVENT_FORMAT_ERROR              ((eventTypeType) 9)
205   #define START_OF_DATA                   ((eventTypeType)10)
206   #define END_OF_DATA                     ((eventTypeType)11)
207   #define SYSTEM_SOFTWARE_TRIGGER_EVENT   ((eventTypeType)12)
208   #define DETECTOR_SOFTWARE_TRIGGER_EVENT ((eventTypeType)13)
209   #define EVENT_TYPE_MIN                  1
210   #define EVENT_TYPE_MAX                  13 
211   */
212
213   fPosition = 0;
214   Bool_t status=kFALSE;
215   fRawDataSize=0;        
216   fDDLNumber = -1;
217   if (fRawReader->GetDDLID()>=0)
218     fDDLNumber = ddlArray[fRawReader->GetDDLID()];
219   
220   if(fDDLNumber<0) {
221     AliWarning(Form("fDDLNumber not a acceptable value %i",fDDLNumber));
222     return kFALSE;
223   }
224     
225   if(fRawReader->GetType() == 7 || fRawReader->GetType() == 8 )  {           //New: Select Physics events, Old: Raw data size is not 0 and not 47148 (pedestal)
226     fnDDLInStream[fDDLNumber]=1; fnDDLOutStream[fDDLNumber]=0;
227     
228     fLDCNumber = fRawReader->GetLDCId();
229     fTimeStamp = fRawReader->GetTimestamp();
230     
231     AliDebug(1,Form("DDL %i started to be decoded!",fDDLNumber));
232     fRawDataSize=fRawReader->GetDataSize()/4;
233     DelVars();                                         //We have to delete the variables initialized in the InitVars before recall IntiVars!!!
234     InitVars(fRawDataSize);                            //To read the charge and pads we cannot delete before the status return
235     
236     if(fTurbo==kTRUE) status=Turbo();
237     else status = ReadHMPIDRawData();
238    
239     if(status) AliDebug(1,Form("Event DDL %i successfully decoded!.",fDDLNumber));
240     else AliDebug(1,Form("Event DDL %i ERROR in decoding!.",fDDLNumber));
241     //DumpData(fRawReader->GetDataSize());
242   
243    
244     }
245     if(status==kTRUE) {fnDDLOutStream[fDDLNumber]++; }//Printf("fnDDLOutStream[%d]=%d",fDDLNumber,fnDDLOutStream[fDDLNumber]); } //Count the number of events when the DDL was succesfully decoded
246    
247 //    return status;  // temporary solution...
248    return kTRUE;
249 }
250 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
251 void AliHMPIDRawStream::InitVars(Int_t n)
252 {
253   //
254   //
255   //
256   fNPads = 0;
257   fCharge = new Int_t[n]; 
258   fPad = new Int_t[n];
259   //for debug purpose
260   fPos = new Int_t[4*n+4];                     //reset debug
261   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
262   fiPos = 0;
263 }    
264 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
265 void AliHMPIDRawStream::DelVars()
266 {
267   //
268   //
269   // 
270   //Clean the initvars!!!!!!!!
271   fNPads = 0; 
272   if (fCharge)     { delete [] fCharge;    fCharge = 0x0; }
273   if (fPad)        { delete [] fPad;       fPad = 0x0;       }   
274   if (fPos)        { delete [] fPos;       fPos = 0x0;       }     
275  
276   fiPos=0;
277   
278      
279 }    
280 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
281 Bool_t AliHMPIDRawStream::ReadHMPIDRawData()
282 {
283   //Here the loop on the decoding the raw bank 
284   //for one ddl starts.
285   //It returns: kFALSE if any error occurs
286   //            kTRUE  if all OK
287   Int_t cntGlob = fRawReader->GetDataSize()/4;
288   if(cntGlob==0) {fNumOfErr[fDDLNumber][kRawDataSizeErr]++; return kFALSE; }
289   
290   Int_t cnt = cntGlob;
291   Int_t nwSeg;
292   Int_t cntSegment;
293
294   if(!GetWord(cnt)) return kFALSE;
295   cnt--;
296
297   
298   while (cnt>20) { //counter limit is changed from 0 to 20 to take into account (skipp) the 5 extra words in the equipment header
299     nwSeg = (fWord >> kbit8) & 0xfff;
300     if(!CheckSegment()) return kFALSE;
301     if(!ReadSegment(cntSegment)) return kFALSE;
302
303     if(nwSeg != cntSegment) {AliDebug(1,Form("Error in Segment counters: %i different wrt %i",nwSeg,cntSegment)); return kFALSE;}
304     if(!GetWord(cntSegment+1,kBwd)) return kFALSE;
305     cnt-=cntSegment+1;
306   }
307   
308   return kTRUE;
309   
310 }
311 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
312 Bool_t AliHMPIDRawStream::ReadSegment(Int_t &cntSegment)
313 {
314   //Read the segment
315   //It returns: kFALSE if any error occurs
316   //            kTRUE  if all OK
317   cntSegment = (fWord >> kbit8) & 0xfff;
318   Int_t cnt = cntSegment;
319   Int_t cntRow;
320   Int_t nwRow;
321
322   if(!GetWord(cnt,kBwd)) return kFALSE;
323   
324   while (cnt>0) {
325
326     cntRow  = (fWord >> kbit16) & 0xfff;
327     if(!CheckRowMarker()) return kFALSE;
328     if(!ReadRow(nwRow)) return kFALSE;
329
330     if(nwRow != cntRow) {AliDebug(1,Form("Error in Row counters: %i different wrt %i",nwRow,cntRow)); return kFALSE;}
331     if(!GetWord(cntRow+1)) return kFALSE;
332     cnt -= cntRow+1;
333     
334   }
335
336   cntSegment -= cnt;
337   
338   return kTRUE;
339     
340 }    
341 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
342 Bool_t AliHMPIDRawStream::ReadRow(Int_t &cntRow)
343 {
344   // Read the row
345   //It returns: kFALSE if any error occurs
346   //            kTRUE  if all OK
347
348   Int_t cnt;
349   Int_t cntDilogic;
350   Int_t nwDil;
351   
352   cntRow  = (fWord >> kbit16) & 0xfff;
353   cnt = cntRow;  
354   
355   if(!GetWord(cntRow)) return kFALSE;
356   
357   while (cnt>0) {
358     
359     if(!CheckEoE(nwDil)) return kFALSE;
360     if(!ReadDilogic(cntDilogic)) return kFALSE;
361     
362     if(nwDil != cntDilogic) {AliDebug(1,Form("Error in Dilogic counters: %i different wrt %i",nwDil,cntDilogic));return kFALSE;}
363     cnt -= cntDilogic;
364     if(!GetWord(1,kBwd)) return kFALSE; // go to next Dilogic bank...
365     cnt--;
366 //    Printf(" cnt %i cntDilogic %i ",cnt,cntDilogic);
367   }
368   
369   cntRow -= cnt;  
370   
371   return kTRUE;
372   
373 }
374 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
375 Bool_t AliHMPIDRawStream::ReadDilogic(Int_t &cntDilogic)
376 {
377   // Read the dilogic bank
378   //It returns: kFALSE if any error occurs
379   //            kTRUE  if all OK
380
381   cntDilogic = fWord & 0x7f;
382   
383   Int_t cnt = cntDilogic;
384   
385 //  Printf(" cnt DILOGIC %i at %i word %08X",cnt,fPosition,fWord);
386
387   for(Int_t iDil=0;iDil<cntDilogic;iDil++) {
388     UInt_t dilogic = 0, row = 0;
389     if(!GetWord(1,kBwd)) return kFALSE;
390 //check on row number      
391     cnt--;
392     row = (fWord >> kbit22) & 0x1f;
393     if(!CheckRow(row)) continue;
394 //check dilogic number     
395     dilogic = (fWord >> kbit18) & 0xf;                                              //dilogic info in raw word is between bits: 18...21
396     if(!CheckDilogic(dilogic)) continue;
397 //check pad number
398     UInt_t pad = (fWord >> kbit12) & 0x3f;                                          //pad info in raw word is between bits: 12...17
399     if(!CheckPad(pad)) continue;
400     Int_t charge = fWord & 0xfff;
401     if(GetPad(fDDLNumber,row,dilogic,pad)<0) continue;
402     fPad[fNPads] = GetPad(fDDLNumber,row,dilogic,pad);
403     fCharge[fNPads] = charge; 
404     fNPads++;
405     
406     if(charge==0) 
407     {
408       AliDebug(1,Form("If PEDESTAL run -> WARNING: ZERO charge is read from DDL: %d row: %d dil: %d pad: %d",fDDLNumber,row,dilogic,pad));
409       fNumOfErr[fDDLNumber][kPedQZero]++;
410     }
411     
412   }//iDil
413
414   cntDilogic -= cnt;  
415   return kTRUE;  
416 }
417 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
418 Bool_t AliHMPIDRawStream::CheckSegment()
419 {
420   // Check the segment marker
421   //It returns: kFALSE if any error occurs
422   //            kTRUE  if all OK
423
424   UInt_t markSegment = 0xAB0;
425   /*
426   if (iRow%8 == 0) {
427     UInt_t segWord = GetWord();          
428     if ((segWord >> 20) != statusSegWord) {
429       fRawReader->AddMajorErrorLog(kBadSegWordErr);
430       AliDebug(1,Form("Wrong segment word signature: %x, expected 0xab0!",(segWord >> 20)));
431       fNumOfErr[kBadSegWordErr]++;
432       return kTRUE;
433     }
434 */
435   UInt_t segMarker = (fWord >> kbit20) & 0xfff;
436   if (segMarker != markSegment ) {
437     //fRawReader->AddMajorErrorLog(kWrongSegErr,Form("Segment marker %0X wrong (expected %0X) at %i in word %0X!",segMarker,markSegment,fPosition,fWord));
438     AliDebug(1,Form("Segment marker %X wrong (expected %0X)! at %i in word %0X!",segMarker,markSegment,fPosition,fWord));
439     fNumOfErr[fDDLNumber][kWrongSegErr]++;
440     return kFALSE;
441   }
442   
443   UInt_t segAddress = fWord & 0xff;
444   if (segAddress<1 ||segAddress>3) {
445     //fRawReader->AddMajorErrorLog(kWrongSegErr,Form("Segment address %d not in the valid range [1-3] at %i in word %0X",segAddress,fPosition,fWord));
446     AliDebug(1,Form("Segment address %d not in the valid range [1-3]",segAddress));
447     fNumOfErr[fDDLNumber][kWrongSegErr]++;
448     return kFALSE;
449   }
450 //  Printf("Segment Marker found at %i! Number of segment is %i",fPosition,segAddress);
451   return kTRUE;
452 }
453 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
454 Bool_t AliHMPIDRawStream::CheckRow(UInt_t row)
455 {
456   //check on row number      
457   //It returns: kFALSE if any error occurs
458   //            kTRUE  if all OK
459 //  Printf("ROW %i word %0X",row,fWord);
460   if(row>=1 && row <=kNRows) return kTRUE;
461   
462   //fRawReader->AddMajorErrorLog(kWrongRowErr,Form("row %d",row));
463   AliDebug(1,Form("Wrong row index: %d, expected (1 -> %d) word %0X at %i...",row,kNRows,fWord,fPosition));
464   fNumOfErr[fDDLNumber][kWrongRowErr]++;
465   return kFALSE;
466 }
467 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
468 Bool_t AliHMPIDRawStream::CheckDilogic(UInt_t dilogic)
469 {
470 //check dilogic number     
471   //It returns: kFALSE if any error occurs
472   //            kTRUE  if all OK
473   if (dilogic>= 1 && dilogic <=kNDILOGICAdd) return kTRUE;
474
475   //fRawReader->AddMajorErrorLog(kWrongDilogicErr,Form("dil %d",dilogic));
476   AliDebug(1,Form("Wrong DILOGIC index: %d, expected (1 -> %d)!",dilogic,kNDILOGICAdd));
477   fNumOfErr[fDDLNumber][kWrongDilogicErr]++;
478   //dilogic = iDILOGIC;
479   return kFALSE;
480 }      
481 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
482 Bool_t AliHMPIDRawStream::CheckPad(UInt_t pad)
483 {
484 //check pad number     
485   //It returns: kFALSE if any error occurs
486   //            kTRUE  if all OK
487   if (pad < kNPadAdd) return kTRUE;
488   
489   //fRawReader->AddMajorErrorLog(kWrongPadErr,Form("pad %d",pad));
490   AliDebug(1,Form("Wrong pad index: %d, expected (0 -> %d)!",pad,kNPadAdd));
491   fNumOfErr[fDDLNumber][kWrongPadErr]++;
492   return kFALSE;
493 }    
494 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
495 Bool_t AliHMPIDRawStream::CheckEoE(Int_t &nDil)
496 {
497   //check the End of Event
498   //It returns: kFALSE if any error occurs
499   //            kTRUE  if all OK
500   if (!((fWord >> kbit27) & 0x1)) {                                                //check 27th bit in EoE. It must be 1!
501     //fRawReader->AddMajorErrorLog(kEoEFlagErr);
502     AliDebug(1,Form("Missing end-of-event flag! (%08X) at %i",fWord,fPosition));
503     fNumOfErr[fDDLNumber][kEoEFlagErr]++;
504     return kFALSE;
505   }
506   nDil = fWord & 0x7f;                                               //nDil=EoE word count
507   if(nDil < 0 || nDil > 48 ) { 
508
509     //fRawReader->AddMajorErrorLog(kEoESizeErr,Form("EoE size=%d",nDil));
510     AliDebug(1,Form("Wrong end-of-event word-count: %08X",fWord));
511     fNumOfErr[fDDLNumber][kEoESizeErr]++;
512     return kFALSE;
513   }
514 //  UInt_t da = (eOfEvent >> 18) & 0xf;
515 //  if (cntData!=0 && da != dilogic) {
516 //    fRawReader->AddMajorErrorLog(kEoEDILOGICErr,Form("eoe dil %d != %d",da,dilogic));
517 //    AliDebug(1,Form("Wrong DILOGIC address found in end-of-event: %d, expected %d!",da,dilogic));
518 //    fNumOfErr[kEoEDILOGICErr]++;
519 //    return kFALSE;  AliQAChecker::Instance()->Run(AliQAv1::kHMPID, task, obj) ;  
520
521 //  }
522 //  UInt_t ca = (eOfEvent >> 22) & 0x1f;
523 //  if (cntData!=0 &&  ca != row) {
524 //    fRawReader->AddMajorErrorLog(kEoERowErr,Form("eoe row %d != %d",ca,row));
525 //    AliDebug(1,Form("Wrong row index found in end-of-event: %d, expected %d!",ca,row));
526 //    fNumOfErr[kEoERowErr]++;
527 //    return kFALSE;
528 //  }
529   return kTRUE;
530 }
531 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
532 Bool_t AliHMPIDRawStream::CheckRowMarker()
533 {
534   //check the row marker
535   //It returns: kFALSE if any error occurs
536   //            kTRUE  if all OK
537   UInt_t nMAXwordsInRow = 0x1EA;
538   UInt_t statusControlRow = 0x32a8; // 0x36a8 for zero suppression
539 //First check on row marker    
540   UInt_t rowControlWord = fWord >> kbit0 & 0xfbff;
541   
542   if(rowControlWord != statusControlRow) {
543     //fRawReader->AddMajorErrorLog(kRowMarkerErr);
544     AliDebug(1,Form("Wrong row marker %x expected 0x32a8!",rowControlWord));
545               fNumOfErr[fDDLNumber][kRowMarkerErr]++;
546               return kFALSE;
547    }
548 //Second check on row marker    
549    UInt_t wordsInRow = fWord >> kbit16 & 0x0fff;                // Number of words after the row marker, bit 10 is skipped in this check
550   
551   if (wordsInRow > nMAXwordsInRow) {
552     //fRawReader->AddMajorErrorLog(kRowMarkerSizeErr);
553     AliDebug(1,Form(" FATAL: Number of words %x in a row exceeds the expected value: 0x1EA !",wordsInRow));
554     fNumOfErr[fDDLNumber][kRowMarkerSizeErr]++;
555     return kFALSE;
556   }
557   
558   return kTRUE;
559 }
560 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
561 Bool_t AliHMPIDRawStream::GetWord(Int_t n,EDirection dir)
562 {
563   // This method returns the n-th 32 bit word
564   // inside the raw data payload.
565   // The method is supposed to be endian (platform)
566   // independent.
567   
568   fWord = 0;
569   if (fPosition < 0) {
570     AliError("fPosition < 0 !!! Event skipped.");
571     fRawReader->AddMajorErrorLog(kRawDataSizeErr,"fPosition<0");
572     return kFALSE;
573   }
574
575   if(dir==kBwd) n = -n; 
576   fPosition+=4*n-4;
577
578   if(fPosition==-4) return kTRUE;
579   
580   if(fPosition<0 || fPosition > fRawReader->GetDataSize()) {
581     AliWarning(Form("fPosition out of boundaries %i",fPosition));
582     return kFALSE;
583   }
584     
585   StorePosition();
586   
587   fWord |= fData[fPosition++];
588   fWord |= fData[fPosition++] << 8;
589   fWord |= fData[fPosition++] << 16;
590   fWord |= fData[fPosition++] << 24;
591
592   return kTRUE;
593 }
594 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
595 void AliHMPIDRawStream::DumpData(Int_t nw)
596 {
597   //just a simple raw data dump
598   // in () is the position in bytes
599   //--
600    for(Int_t i=0;i<nw;i+=4) {
601      if(!(i%16)) printf(" \n %8i) ",i);
602      printf("%02X%02X%02X%02X [ %06i ] ",fData[i+3],fData[i+2],fData[i+1],fData[i+0],fPos[i]);
603    }
604    Printf(" \n -----end of dump ----------- ");
605 }
606 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
607 void AliHMPIDRawStream::StorePosition()
608 {
609   //just for debug purpose
610   // it stores the position
611   //read for the first time
612 //  Printf("@@@@@@@@@ fPos: %x fPosition: %d",fPos,fPosition);
613   if(fPos[fPosition]!=0) {
614 //    Printf("Position already stored!!! Value %i at address %i",fPos[fPosition],fPosition); 
615     return;
616   }
617   fiPos++;
618   fPos[fPosition] = fiPos;
619 //  if(stDeb)Printf("%i - Actual position %i",iPos,fPosition); 
620 }
621 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++