Fix for FMD DA
[u/mrichter/AliRoot.git] / ITS / ITSbase / AliITSRawStreamSPD.cxx
1 /**************************************************************************
2  * Copyright(c) 2007-2009, 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 /* $Id$ */
17
18 ///////////////////////////////////////////////////////////////////////////////
19 ///
20 /// This class provides access to ITS SPD digits in raw data.
21 ///
22 ///////////////////////////////////////////////////////////////////////////////
23
24 #include "AliITSRawStreamSPD.h"
25 #include "AliRawReader.h"
26 #include "AliLog.h"
27
28 ClassImp(AliITSRawStreamSPD)
29
30 const Int_t AliITSRawStreamSPD::fgkDDLModuleMap[kDDLsNumber][kModulesPerDDL] = {
31   { 4, 5, 0, 1, 80, 81, 84, 85, 88, 89, 92, 93},
32   {12,13, 8, 9, 96, 97,100,101,104,105,108,109},
33   {20,21,16,17,112,113,116,117,120,121,124,125},
34   {28,29,24,25,128,129,132,133,136,137,140,141},
35   {36,37,32,33,144,145,148,149,152,153,156,157},
36   {44,45,40,41,160,161,164,165,168,169,172,173},
37   {52,53,48,49,176,177,180,181,184,185,188,189},
38   {60,61,56,57,192,193,196,197,200,201,204,205},
39   {68,69,64,65,208,209,212,213,216,217,220,221},
40   {76,77,72,73,224,225,228,229,232,233,236,237},
41   { 7, 6, 3, 2, 83, 82, 87, 86, 91, 90, 95, 94},
42   {15,14,11,10, 99, 98,103,102,107,106,111,110},
43   {23,22,19,18,115,114,119,118,123,122,127,126},
44   {31,30,27,26,131,130,135,134,139,138,143,142},
45   {39,38,35,34,147,146,151,150,155,154,159,158},
46   {47,46,43,42,163,162,167,166,171,170,175,174},
47   {55,54,51,50,179,178,183,182,187,186,191,190},
48   {63,62,59,58,195,194,199,198,203,202,207,206},
49   {71,70,67,66,211,210,215,214,219,218,223,222},
50   {79,78,75,74,227,226,231,230,235,234,239,238}
51 };
52 //__________________________________________________________________________
53 AliITSRawStreamSPD::AliITSRawStreamSPD(AliRawReader* rawReader) :
54   AliITSRawStream(rawReader),
55   fEventCounter(-1),fChipAddr(0),fHalfStaveNr(0),fCol(0),fRow(0),fCalHeadLen(0),
56   fData(0),fOffset(0),fHitCount(0),
57   fDataChar1(0),fDataChar2(0),fDataChar3(0),fDataChar4(0),
58   fFirstWord(kTRUE),fPrevEventId(0xffffffff),
59   fEqPLBytesRead(0),fEqPLChipHeadersRead(0),fEqPLChipTrailersRead(0),
60   fHeaderOrTrailerReadLast(kFALSE),fExpectedHeaderTrailerCount(0),
61   fFillOutOfSynch(kFALSE),fDDLID(-1),fLastDDLID(-1),fAdvancedErrorLog(kFALSE),
62   fAdvLogger(NULL)
63 {
64   // create an object to read ITS SPD raw digits
65   fRawReader->Reset();
66   fRawReader->Select("ITSSPD");
67   // reset calib header words
68   for (UInt_t iword=0; iword<kCalHeadLenMax; iword++) {
69     fCalHeadWord[iword]=0xffffffff;
70   }
71   for (UInt_t eq=0; eq<20; eq++) {
72     fActiveEq[eq]=kFALSE;
73     for (UInt_t hs=0; hs<6; hs++) {
74       fActiveHS[eq][hs]=kFALSE;
75       for (UInt_t chip=0; chip<10; chip++) {
76         fActiveChip[eq][hs][chip]=kFALSE;
77         fEventCounterFull[eq][hs][chip] = -1;
78       }
79     }
80   }
81   NewEvent();
82 }
83 //__________________________________________________________________________
84 AliITSRawStreamSPD::AliITSRawStreamSPD(const AliITSRawStreamSPD& rstream) :
85   AliITSRawStream(rstream.fRawReader),
86   fEventCounter(-1),fChipAddr(0),fHalfStaveNr(0),fCol(0),fRow(0),fCalHeadLen(0),
87   fData(0),fOffset(0),fHitCount(0),
88   fDataChar1(0),fDataChar2(0),fDataChar3(0),fDataChar4(0),
89   fFirstWord(kTRUE),fPrevEventId(0xffffffff),
90   fEqPLBytesRead(0),fEqPLChipHeadersRead(0),fEqPLChipTrailersRead(0),
91   fHeaderOrTrailerReadLast(kFALSE),fExpectedHeaderTrailerCount(0),
92   fFillOutOfSynch(kFALSE),fDDLID(-1),fLastDDLID(-1),fAdvancedErrorLog(kFALSE),
93   fAdvLogger(NULL)
94 {
95   // copy constructor
96   AliError("Copy constructor should not be used.");
97 }
98 //__________________________________________________________________________
99 AliITSRawStreamSPD& AliITSRawStreamSPD::operator=(const AliITSRawStreamSPD& rstream) {
100   // assignment operator
101   if (this!=&rstream) {}
102   AliError("Assignment opertator should not be used.");
103   return *this;
104 }
105 //__________________________________________________________________________
106 Bool_t AliITSRawStreamSPD::ReadNextShort() {
107   // read next 16 bit word into fData
108   if (fFirstWord) {
109     Bool_t b1 = fRawReader->ReadNextChar(fDataChar1);
110     if (!b1) return kFALSE;
111     Bool_t  b2, b3, b4;
112     b2 = fRawReader->ReadNextChar(fDataChar2);
113     b3 = fRawReader->ReadNextChar(fDataChar3);
114     b4 = fRawReader->ReadNextChar(fDataChar4);
115     if (!(b2 && b3 && b4)) {
116       return kFALSE;
117     }
118     fData = fDataChar3+(fDataChar4<<8);
119     const UInt_t *intPtr = fRawReader->GetEventId();
120     if (intPtr!=0) {
121       if (*intPtr!=fPrevEventId) { // if new event...
122         NewEvent();
123         fPrevEventId=*intPtr;
124       }
125     }
126     fFirstWord=kFALSE;
127   }
128   else {
129     fFirstWord=kTRUE;
130     fData = fDataChar1+(fDataChar2<<8);
131   }
132   fEqPLBytesRead+=2;
133   
134   return kTRUE;
135 }
136 //__________________________________________________________________________
137 Bool_t AliITSRawStreamSPD::ReadNextInt() {
138   // reads next 32 bit into fDataChar1..4 
139   // (if first 16 bits read already, just completes the present word)
140   if (fFirstWord) {
141     if (ReadNextShort() && ReadNextShort()) {
142       return kTRUE;
143     }
144   }
145   else {
146     if (ReadNextShort()) {
147       return kTRUE;
148     }
149   }
150   return kFALSE;
151 }
152 //__________________________________________________________________________
153 void AliITSRawStreamSPD::NewEvent() {
154   // call this to reset flags for a new event
155   for (UInt_t eqId=0; eqId<20; eqId++) {
156     fCalHeadRead[eqId]=kFALSE;
157   }
158   fEventCounter = -1;
159   fDDLID = -1;
160   fLastDDLID = -1;
161   for (UInt_t eq=0; eq<20; eq++) {
162     for (UInt_t hs=0; hs<6; hs++) {
163       for (UInt_t chip=0; chip<10; chip++) {
164         fFastOrSignal[eq][hs][chip] = kFALSE;
165       }
166     }
167   }
168 }
169 //__________________________________________________________________________
170 Int_t AliITSRawStreamSPD::ReadCalibHeader() {
171   // needs to be called in the beginning of each equipment data
172   // read the extra calibration header
173   // returns the length of the header if it is present, -1 otherwise
174
175   Int_t ddlID = fRawReader->GetDDLID();
176   if (ddlID==-1) { // we may need to read one word to get the blockAttr
177     if (!ReadNextShort()) return -1;
178     ddlID = fRawReader->GetDDLID();
179     if(ddlID<0){
180       AliError("fRawReader->GetDDLID() returns a negative value");
181       ddlID=0;
182     }
183   }
184   // reset flags and counters
185   fEqPLBytesRead = 2;
186   fEqPLChipHeadersRead = 0;
187   fEqPLChipTrailersRead = 0;
188   fHeaderOrTrailerReadLast = kFALSE;
189   fFillOutOfSynch = kFALSE;
190
191   // check router error bits:
192   UInt_t statusBits = fRawReader->GetStatusBits();
193   if ((statusBits >> 5) & 1) { // linkrx/detector fatal error bit
194     TString errMess = "LinkRx Error Bit Set";
195     AliError(errMess.Data());
196     fRawReader->AddMajorErrorLog(kLinkRxDetectorFatalErr,errMess.Data());
197     if (fAdvancedErrorLog) fAdvLogger->ProcessError(kLinkRxDetectorFatalErr,ddlID,-1,-1,errMess.Data());
198   }
199   if ((statusBits >> 12) & 1) { // trigger sequence monitor error bit
200     TString errMess = "TSM Trigger Error Bit Set";
201     AliError(errMess.Data());
202     fRawReader->AddMajorErrorLog(kTSMtriggerErr,errMess.Data());
203     if (fAdvancedErrorLog) fAdvLogger->ProcessError(kTSMtriggerErr,ddlID,-1,-1,errMess.Data());
204   }
205   if (fRawReader->TestBlockAttribute(7)) { // bunch crossing difference error bit
206     TString errMess = "High Multiplicity Event Flag Set";
207     AliError(errMess.Data());
208     fRawReader->AddMajorErrorLog(kHighMultiplicityFlag,errMess.Data());
209     if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHighMultiplicityFlag,ddlID,-1,-1,errMess.Data());
210   }
211
212   // set eq active and the participating half-staves (for access from outside this class), 
213   // set what number of chip headers/trailers to expect (for the moment not used)
214   fActiveEq[ddlID]=kTRUE;
215   fExpectedHeaderTrailerCount = 0;
216   for (UInt_t hs=0; hs<6; hs++) {
217     if (!fRawReader->TestBlockAttribute(hs)) {
218       fActiveHS[ddlID][hs]=kTRUE;
219       fExpectedHeaderTrailerCount+=10;
220     }
221     else {
222       fActiveHS[ddlID][hs]=kFALSE;
223     }
224   }
225
226   if (ddlID>=0 && ddlID<20) fCalHeadRead[ddlID]=kTRUE;
227   if (fRawReader->TestBlockAttribute(6)) { // is the calib header present?
228     if (ReadNextInt()) {
229       // length of cal header:
230       fCalHeadLen = fDataChar1+(fDataChar2<<8)+(fDataChar3<<16)+(fDataChar4<<24);
231       if (fCalHeadLen>kCalHeadLenMax) {
232         TString errMess = Form("Header length %d > max = %d",fCalHeadLen,kCalHeadLenMax);
233         AliError(errMess.Data());
234         fRawReader->AddMajorErrorLog(kCalHeaderLengthErr,errMess.Data());
235         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kCalHeaderLengthErr,ddlID,-1,-1,errMess.Data());
236         return -1;
237       }
238       else {
239         for (UInt_t iword=0; iword<fCalHeadLen; iword++) {
240           if (ReadNextInt()) {
241             fCalHeadWord[iword] = fDataChar1+(fDataChar2<<8)+(fDataChar3<<16)+(fDataChar4<<24);
242           }
243           else {
244             TString errMess = "Header length problem";
245             AliError(errMess.Data());
246             fRawReader->AddMajorErrorLog(kCalHeaderLengthErr,errMess.Data());
247             if (fAdvancedErrorLog) fAdvLogger->ProcessError(kCalHeaderLengthErr,ddlID,-1,-1,errMess.Data());
248             return -1;
249           }
250         }
251         return fCalHeadLen;
252       }
253     }
254   }
255
256   return -1;
257 }
258 //__________________________________________________________________________
259 Bool_t AliITSRawStreamSPD::Next() {
260 // read the next raw digit
261 // returns kFALSE if there is no digit left
262
263   fPrevModuleID = fModuleID;
264
265   while (ReadNextShort()) {
266
267     fLastDDLID = fDDLID;
268     fDDLID = fRawReader->GetDDLID();
269     if (fDDLID>=0 && fDDLID<20) {
270       if (!fCalHeadRead[fDDLID]) {
271         if (fLastDDLID!=-1) {  // if not the first equipment for this event
272           fEqPLBytesRead -= 2;
273           // The next line is commented, since we do not have the information how many chips to expect
274           // CheckHeaderAndTrailerCount(fLastDDLID);
275         }
276         if (ReadCalibHeader()>=0) continue; // skip to next word if we found a calib header, otherwise parse this word as regular data (go on with this while loop)
277       }
278     }
279     else {
280       TString errMess = Form("Error in DDL number (=%d) , setting it to 19",fDDLID);
281       AliError(errMess.Data());
282       fRawReader->AddMajorErrorLog(kDDLNumberErr,errMess.Data()); 
283       if (fAdvancedErrorLog) fAdvLogger->AddMessage(errMess.Data());
284       fDDLID=19;
285     }
286     
287     Short_t hs = ((fData & 0x3800)>>11) ; Short_t chip=(fData & 0x000F); // -> for printouts
288     
289     if ((fData & 0xC000) == 0x4000) {         // header
290       if (fHeaderOrTrailerReadLast) {
291         TString errMess = Form("Chip trailer missing - (eq,hs,chip) = (%d,%d,%d)",fDDLID,hs,chip);
292         AliError(errMess.Data());
293         fRawReader->AddMajorErrorLog(kTrailerMissingErr,errMess.Data());
294         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kTrailerMissingErr,fLastDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
295       }
296       fHeaderOrTrailerReadLast = kTRUE;
297       fEqPLChipHeadersRead++;
298       fHitCount = 0;
299       UShort_t eventCounter = (fData >> 4) & 0x007F;
300       if (fEventCounter < 0) {
301         fEventCounter = eventCounter;
302       } 
303       else if (eventCounter != fEventCounter) {
304         TString errMess;
305         if (fEqPLChipHeadersRead==1) {
306           errMess = Form("Mismatching event counters between this equipment and the previous: %d != %d - (eq,hs,chip) = (%d,%d,%d)",
307                          eventCounter,fEventCounter,fDDLID,hs,chip);
308         }
309         else {
310           errMess = Form("Mismatching event counters between this chip header and the previous: %d != %d - (eq,hs,chip) = (%d,%d,%d)",
311                          eventCounter,fEventCounter,fDDLID,hs,chip);
312         }
313         fEventCounter = eventCounter;
314         AliError(errMess.Data());
315         fRawReader->AddMajorErrorLog(kEventCounterErr,errMess.Data());
316         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kEventCounterErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
317       }
318       fChipAddr = fData & 0x000F;
319       if (fChipAddr>9) {
320         TString errMess = Form("Overflow chip address %d - set to 0 (eq,hs) = (%d,%d)",fChipAddr,fDDLID,hs);
321         AliError(errMess.Data());
322         fRawReader->AddMajorErrorLog(kChipAddrErr,errMess.Data());
323         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kChipAddrErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
324         fChipAddr=0;
325       }
326       fHalfStaveNr = (fData & 0x3800)>>11;
327       if (fHalfStaveNr>5 || fRawReader->TestBlockAttribute(fHalfStaveNr)) {
328         TString errMess = Form("Half stave number error: %d - set to 0 - eq %d",fHalfStaveNr,fDDLID);
329         AliError(errMess.Data());
330         fRawReader->AddMajorErrorLog(kHSNumberErr,errMess.Data());
331         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHSNumberErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
332         fHalfStaveNr=0;
333       }
334       fActiveChip[fDDLID][fHalfStaveNr][fChipAddr]=kTRUE;
335       fEventCounterFull[fDDLID][fHalfStaveNr][fChipAddr] = eventCounter;
336       // translate  ("online") ddl, hs, chip nr  to  ("offline") module id :
337       fModuleID = GetOfflineModuleFromOnline(fDDLID,fHalfStaveNr,fChipAddr);
338     }
339
340     else if ((fData & 0xC000) == 0x0000) {    // trailer
341       if ( (fEqPLBytesRead+fFillOutOfSynch*2)%4 != 0 ) {
342         TString errMess = Form("Fill word is missing - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
343         AliError(errMess.Data());
344         fRawReader->AddMajorErrorLog(kFillMissingErr,errMess.Data());
345         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kFillMissingErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
346         if (fFillOutOfSynch) fFillOutOfSynch = kFALSE;
347         else fFillOutOfSynch = kTRUE;
348       }
349       if (!fHeaderOrTrailerReadLast) {
350         TString errMess = Form("Trailer without previous header - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
351         AliError(errMess.Data());
352         fRawReader->AddMajorErrorLog(kTrailerWithoutHeaderErr,errMess.Data());
353         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kTrailerWithoutHeaderErr,fLastDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
354       }
355       fHeaderOrTrailerReadLast = kFALSE;
356       fEqPLChipTrailersRead++;
357       UShort_t hitCount = fData & 0x0FFF;
358       if (hitCount != fHitCount){
359         TString errMess = Form("Number of hits %d, while %d expected - (eq,hs,chip) = (%d,%d,%d)",fHitCount,hitCount,fDDLID,fHalfStaveNr,fChipAddr);
360         AliError(errMess.Data());
361         fRawReader->AddMajorErrorLog(kNumberHitsErr,errMess.Data());
362         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kNumberHitsErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
363       }
364       Bool_t errorBit = fData & 0x1000;
365       if (errorBit) {
366         TString errMess = Form("Trailer error bit set for chip (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
367         AliError(errMess.Data());
368         fRawReader->AddMajorErrorLog(kTrailerErrorBitErr,errMess.Data());
369         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kTrailerErrorBitErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
370       }
371       Bool_t fastorBit = fData & 0x2000;
372       fFastOrSignal[fDDLID][fHalfStaveNr][fChipAddr] = fastorBit;
373     }
374
375     else if ((fData & 0xC000) == 0x8000) {    // pixel hit
376       if (!fHeaderOrTrailerReadLast) {
377         TString errMess = Form("Chip header missing - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
378         AliError(errMess.Data());
379         fRawReader->AddMajorErrorLog(kHeaderMissingErr,errMess.Data());
380         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHeaderMissingErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
381         fHeaderOrTrailerReadLast = kTRUE;
382       }
383       fHitCount++;
384       fCol = (fData & 0x001F);
385       fRow = (fData >> 5) & 0x00FF;
386
387       fCoord1 = GetOfflineColFromOnline(fDDLID,fHalfStaveNr,fChipAddr,fCol);
388       fCoord2 = GetOfflineRowFromOnline(fDDLID,fHalfStaveNr,fChipAddr,fRow);
389
390       return kTRUE;
391     }
392
393     else {                                    // fill word
394       if ((fData & 0xC000) != 0xC000) {
395         TString errMess = Form("Wrong fill word! - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
396         AliError(errMess.Data());
397         fRawReader->AddMajorErrorLog(kWrongFillWordErr,errMess.Data());
398         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kWrongFillWordErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
399       }
400       if ( (fEqPLBytesRead+fFillOutOfSynch*2)%4 != 2 ) {
401         TString errMess = Form("Fill word is unexpected - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
402         AliError(errMess.Data());
403         fRawReader->AddMajorErrorLog(kFillUnexpectErr,errMess.Data());
404         if (fAdvancedErrorLog) fAdvLogger->ProcessError(kFillUnexpectErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
405         if (fFillOutOfSynch) fFillOutOfSynch = kFALSE;
406         else fFillOutOfSynch = kTRUE;
407       }
408     }
409
410   }
411   if (fDDLID>=0 && fDDLID<20) {
412     // The next line is commented, since we do not have the information how many chips to expect
413     // CheckHeaderAndTrailerCount(fDDLID);
414   }
415   return kFALSE;
416 }
417 //__________________________________________________________________________
418 void AliITSRawStreamSPD::CheckHeaderAndTrailerCount(Int_t ddlID) {
419   // Checks that the number of header and trailers found for the ddl are as expected
420   if (fEqPLChipHeadersRead != fExpectedHeaderTrailerCount) {
421     TString errMess = Form("Chip header count inconsistent %d != %d (expected) for ddl %d",
422                            fEqPLChipHeadersRead,fExpectedHeaderTrailerCount,ddlID);
423     AliError(errMess.Data());
424     fRawReader->AddMajorErrorLog(kHeaderCountErr,errMess.Data());
425     if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHeaderCountErr,ddlID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
426   }
427   if (fEqPLChipTrailersRead != fExpectedHeaderTrailerCount) {
428     TString errMess = Form("Chip trailer count inconsistent %d != %d (expected) for ddl %d",
429                            fEqPLChipTrailersRead,fExpectedHeaderTrailerCount,ddlID);
430     AliError(errMess.Data());
431     fRawReader->AddMajorErrorLog(kHeaderCountErr,errMess.Data());
432     //   if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHeaderCountErr,ddlID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
433     if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHeaderCountErr,ddlID,fEqPLBytesRead,fEqPLChipTrailersRead,errMess.Data());
434   }
435 }
436 //__________________________________________________________________________
437 void AliITSRawStreamSPD::ActivateAdvancedErrorLog(Bool_t activate, AliITSRawStreamSPDErrorLog* advLogger) {
438   // Activate the advanced error logging.
439   // Logger object has to be created elsewhere and a pointer sent here.
440   fAdvancedErrorLog=activate;
441   if (activate && advLogger!=NULL) {
442     fAdvLogger = advLogger;
443   }
444   if (fAdvLogger==NULL) {
445     fAdvancedErrorLog=kFALSE;
446   }
447 }
448 //__________________________________________________________________________
449 const Char_t* AliITSRawStreamSPD::GetErrorName(UInt_t errorCode) {
450   // Returns a string for each error code
451   if      (errorCode==kTotal)                   return "All Errors";
452   else if (errorCode==kHeaderMissingErr)        return "Header Missing";
453   else if (errorCode==kTrailerMissingErr)       return "Trailer Missing";
454   else if (errorCode==kTrailerWithoutHeaderErr) return "Trailer Unexpected";
455   else if (errorCode==kHeaderCountErr)          return "Header Count Wrong";
456   else if (errorCode==kTrailerCountErr)         return "Trailer Count Wrong";
457   else if (errorCode==kFillUnexpectErr)         return "Fill Unexpected";
458   else if (errorCode==kFillMissingErr)          return "Fill Missing";
459   else if (errorCode==kWrongFillWordErr)        return "Fill Word Wrong";
460   else if (errorCode==kNumberHitsErr)           return "Hit Count Wrong";
461   else if (errorCode==kEventCounterErr)         return "Event Counter Error";
462   else if (errorCode==kDDLNumberErr)            return "DDL Number Error";
463   else if (errorCode==kHSNumberErr)             return "HS Number Error";
464   else if (errorCode==kChipAddrErr)             return "Chip Address Error";
465   else if (errorCode==kCalHeaderLengthErr)      return "Calib Header Length Error";
466   else if (errorCode==kAdvEventCounterErr)      return "Event Counter Error (Adv)";
467   else if (errorCode==kAdvEventCounterOrderErr) return "Event Counter Jump Error (Adv)";
468   else if (errorCode==kTrailerErrorBitErr)      return "Trailer Error Bit Set";
469   else if (errorCode==kLinkRxDetectorFatalErr)  return "LinkRx/Detector Fatal Error Bit Set";
470   else if (errorCode==kTSMtriggerErr)           return "TSM Trigger Error Bit Set";
471   else if (errorCode==kHighMultiplicityFlag)    return "High Multiplicity Event Flag Set";
472   else return "";
473 }
474 //__________________________________________________________________________
475 Bool_t AliITSRawStreamSPD::GetFastOrSignal(UInt_t eq, UInt_t hs, UInt_t chip) {
476   // returns if there was a fastor signal from this chip
477   if (eq>=20 || hs>=6 || chip>=10) {
478     TString errMess = Form("eq,hs,chip = %d,%d,%d out of bounds. Return kFALSE.",eq,hs,chip);
479     AliError(errMess.Data());
480     return kFALSE;
481   }
482   return fFastOrSignal[eq][hs][chip];
483 }
484 //__________________________________________________________________________
485 Bool_t AliITSRawStreamSPD::IsActiveEq(UInt_t eq) const {
486   // returns if the eq is active (seen in data)
487   if (eq>=20) {
488     TString errMess = Form("eq = %d out of bounds. Return kFALSE.",eq);
489     AliError(errMess.Data());
490     return kFALSE;
491   }
492   return fActiveEq[eq];
493 }
494 //__________________________________________________________________________
495 Bool_t AliITSRawStreamSPD::IsActiveHS(UInt_t eq, UInt_t hs) const {
496   // returns if the hs is active (info from block attr)
497   if (eq>=20 || hs>=6) {
498     TString errMess = Form("eq,hs = %d,%d out of bounds. Return kFALSE.",eq,hs);
499     AliError(errMess.Data());
500     return kFALSE;
501   }
502   return fActiveHS[eq][hs];
503 }
504 //__________________________________________________________________________
505 Bool_t AliITSRawStreamSPD::IsActiveChip(UInt_t eq, UInt_t hs, UInt_t chip) const {
506   // returns if the chip is active (seen in data)
507   if (eq>=20 || hs>=6 || chip>=10) {
508     TString errMess = Form("eq,hs,chip = %d,%d,%d out of bounds. Return kFALSE.",eq,hs,chip);
509     AliError(errMess.Data());
510     return kFALSE;
511   }
512   return fActiveChip[eq][hs][chip];
513 }
514 //__________________________________________________________________________
515 Bool_t AliITSRawStreamSPD::GetHalfStavePresent(UInt_t hs) {
516   // Reads the half stave present status from the block attributes
517   // This is not really needed anymore (kept for now in case it is still used somewhere).
518   // The same information can be reached through the "IsActiveHS" method instead.
519   Int_t ddlID = fRawReader->GetDDLID();
520   if (ddlID==-1) {
521     AliWarning("DDL ID = -1. Cannot read block attributes. Return kFALSE.");
522     return kFALSE;
523   }
524   else {
525     if (hs>=6) {
526       AliWarning(Form("HS >= 6 requested (%d). Return kFALSE.",hs));
527       return kFALSE;
528     }
529     UChar_t attr = fRawReader->GetBlockAttributes();
530     if (((attr>>hs) & 0x01) == 0x01) { // bit set means not present
531       return kFALSE;
532     }
533     else {
534       return kTRUE;
535     }
536   }
537 }
538 //__________________________________________________________________________
539 Bool_t AliITSRawStreamSPD::IsEventCounterFullConsistent() const {
540   // checks if the event counter values are consistent within the event
541   Short_t reference = -1;
542   for (UInt_t eq=0; eq<20; eq++) {
543     if (IsActiveEq(eq)) {
544       for (UInt_t hs=0; hs<6; hs++) {
545         if (IsActiveHS(eq,hs)) {
546           for (UInt_t chip=0; chip<10; chip++) {
547             if (fEventCounterFull[eq][hs][chip]!=-1) {
548               if (reference==-1) reference = fEventCounterFull[eq][hs][chip];
549               if (fEventCounterFull[eq][hs][chip] != reference) return kFALSE;
550             }
551           }
552         }
553       }
554     }
555   }
556   return kTRUE;
557 }
558 //__________________________________________________________________________
559 Short_t AliITSRawStreamSPD::GetEventCounterFullEq(UInt_t eq) const {
560   // if the eq is active; returns the event counter value
561   if (eq>=20) {
562     TString errMess = Form("eq (%d) out of bounds",eq);
563     AliError(errMess.Data());
564     return -1;
565   }
566   if (IsActiveEq(eq)) {
567     for (UInt_t hs=0; hs<6; hs++) {
568       if (IsActiveHS(eq,hs)) {
569         for (UInt_t chip=0; chip<10; chip++) {
570           if (fEventCounterFull[eq][hs][chip]!=-1) {
571             return fEventCounterFull[eq][hs][chip];
572           }
573         }
574       }
575     }
576   }
577   return -1;
578 }
579 //__________________________________________________________________________
580 Short_t AliITSRawStreamSPD::GetEventCounterFullHS(UInt_t eq, UInt_t hs) const {
581   // if the eq,hs is active; returns the event counter value
582   if (eq>=20 || hs>=6) {
583     TString errMess = Form("eq,hs (%d,%d) out of bounds",eq,hs);
584     AliError(errMess.Data());
585     return -1;
586   }
587   if (IsActiveEq(eq)) {
588     if (IsActiveHS(eq,hs)) {
589       for (UInt_t chip=0; chip<10; chip++) {
590         if (fEventCounterFull[eq][hs][chip]!=-1) {
591           return fEventCounterFull[eq][hs][chip];
592         }
593       }
594     }
595   }
596   return -1;
597 }
598 //__________________________________________________________________________
599 Short_t AliITSRawStreamSPD::GetEventCounterFullChip(UInt_t eq, UInt_t hs, UInt_t chip) const {
600   // if the eq,hs,chip is active; returns the event counter value
601   if (eq>=20 || hs>=6 || chip>=10) {
602     TString errMess = Form("eq,hs,chip (%d,%d,%d) out of bounds",eq,hs,chip);
603     AliError(errMess.Data());
604     return -1;
605   }
606   if (IsActiveEq(eq)) {
607     if (IsActiveHS(eq,hs)) {
608       if (IsActiveChip(eq,hs,chip)) {
609         return fEventCounterFull[eq][hs][chip];
610       }
611     }
612   }
613   return -1;
614 }
615 //__________________________________________________________________________
616 Int_t AliITSRawStreamSPD::GetHword(UInt_t index) {
617   if (index<kCalHeadLenMax) return fCalHeadWord[index];
618   else return 0;
619 }
620 Bool_t AliITSRawStreamSPD::GetHhalfStaveScanned(UInt_t hs) const {
621   if (hs<6) return (Bool_t)((fCalHeadWord[0]>>(6+hs)) & (0x00000001));
622   else return kFALSE;
623 }
624 Bool_t AliITSRawStreamSPD::GetHchipPresent(UInt_t hs, UInt_t chip) const {
625   if (hs<6 && chip<10) return ((( fCalHeadWord[hs/3+3]>>((hs%3)*10+chip)) & 0x00000001) == 1);
626   else return kFALSE;
627 }
628 UInt_t AliITSRawStreamSPD::GetHdacHigh(UInt_t hs) const {
629   if (hs<6) return (fCalHeadWord[hs/2+7]>>(24-16*(hs%2)) & 0x000000ff);
630   else return 0;
631 }
632 UInt_t AliITSRawStreamSPD::GetHdacLow(UInt_t hs) const {
633   if (hs<6) return (fCalHeadWord[hs/2+7]>>(16-16*(hs%2)) & 0x000000ff);
634   else return 0;
635 }
636 UInt_t AliITSRawStreamSPD::GetHTPAmp(UInt_t hs) const {
637   if (hs<6) return fCalHeadWord[hs+10];
638   else return 0;
639 }
640 Bool_t AliITSRawStreamSPD::GetHminTHchipPresent(UInt_t chip) const {
641   if (chip<10) return ((( fCalHeadWord[7]>>(16+chip)) & 0x00000001) == 1);
642   else return kFALSE;
643 }
644 UInt_t AliITSRawStreamSPD::GetFOHnumDacs() const {
645   return (fCalHeadLen-37)/2;
646 }
647 UInt_t AliITSRawStreamSPD::GetFOHdacIndex(UInt_t index) const {
648   if (index>=GetFOHnumDacs()) {
649     TString errMess = Form("Only %d DACs in this run, returning 0",GetFOHnumDacs());
650     AliError(errMess.Data());
651     return 0;
652   }
653   return fCalHeadWord[7+index*2];
654 }
655 UInt_t AliITSRawStreamSPD::GetFOHdacValue(UInt_t index) const {
656   if (index>=GetFOHnumDacs()) {
657     TString errMess = Form("Only %d DACs in this run, returning 0",GetFOHnumDacs());
658     AliError(errMess.Data());
659     return 0;
660   }
661   return fCalHeadWord[7+1+index*2];
662 }
663 UInt_t AliITSRawStreamSPD::GetFOHchipCount(UInt_t hs, UInt_t chip) const {
664   if (hs<6 && chip<10) {
665     if (chip%2==0) {
666       return ((fCalHeadWord[7 + GetFOHnumDacs()*2 + (hs*10 + chip)/2] >> 16) & 0x0000ffff);
667     }
668     else {
669       return (fCalHeadWord[7 + GetFOHnumDacs()*2 + (hs*10 + chip)/2] & 0x0000ffff);
670     }
671   }
672   else return 0;
673 }
674 //__________________________________________________________________________
675 Int_t AliITSRawStreamSPD::GetModuleNumber(UInt_t iDDL, UInt_t iModule) {
676   if (iDDL<20 && iModule<12) return fgkDDLModuleMap[iDDL][iModule];
677   else return 240;
678 }
679 //__________________________________________________________________________
680 Bool_t AliITSRawStreamSPD::OfflineToOnline(UInt_t module, UInt_t colM, UInt_t rowM, UInt_t& eq, UInt_t& hs, UInt_t& chip, UInt_t& col, UInt_t& row) {
681   // converts offline coordinates to online
682   eq = GetOnlineEqIdFromOffline(module);
683   hs = GetOnlineHSFromOffline(module);
684   chip = GetOnlineChipFromOffline(module,colM);
685   col = GetOnlineColFromOffline(module,colM);
686   row = GetOnlineRowFromOffline(module,rowM);
687   if (eq>=20 || hs>=6 || chip>=10 || col>=32 || row>=256) return kFALSE;
688   else return kTRUE;
689 }
690 //__________________________________________________________________________
691 Bool_t AliITSRawStreamSPD::OnlineToOffline(UInt_t eq, UInt_t hs, UInt_t chip, UInt_t col, UInt_t row, UInt_t& module, UInt_t& colM, UInt_t& rowM) {
692   // converts online coordinates to offline
693   module = GetOfflineModuleFromOnline(eq,hs,chip);
694   colM = GetOfflineColFromOnline(eq,hs,chip,col);
695   rowM = GetOfflineRowFromOnline(eq,hs,chip,row);
696   if (module>=240 || colM>=160 || rowM>=256) return kFALSE;
697   else return kTRUE;
698 }
699 //__________________________________________________________________________
700 UInt_t AliITSRawStreamSPD::GetOnlineEqIdFromOffline(UInt_t module) {
701   // offline->online (eq)
702   for (UInt_t eqId=0; eqId<20; eqId++) {
703     for (UInt_t iModule=0; iModule<12; iModule++) {
704       if (GetModuleNumber(eqId,iModule)==(Int_t)module) return eqId;
705     }
706   }
707   return 20; // error
708 }
709 //__________________________________________________________________________
710 UInt_t AliITSRawStreamSPD::GetOnlineHSFromOffline(UInt_t module) {
711   // offline->online (hs)
712   for (UInt_t eqId=0; eqId<20; eqId++) {
713     for (UInt_t iModule=0; iModule<12; iModule++) {
714       if (GetModuleNumber(eqId,iModule)==(Int_t)module) return iModule/2;
715     }
716   }
717   return 6; // error
718 }
719 //__________________________________________________________________________
720 UInt_t AliITSRawStreamSPD::GetOnlineChipFromOffline(UInt_t module, UInt_t colM) {
721   // offline->online (chip)
722   for (UInt_t eq=0; eq<20; eq++) {
723     for (UInt_t iModule=0; iModule<12; iModule++) {
724       if (GetModuleNumber(eq,iModule)==(Int_t)module) {
725         if (module<80) {
726           if (eq<10) { // side A
727             return (159-colM)/32 + 5*(iModule%2);
728           }
729           else { // side C
730             return colM/32 + 5*(iModule%2);
731           }
732         }
733         else if (module<240) {
734           if (eq<10) { // side A
735             return colM/32 + 5*(iModule%2);
736           }
737           else { // side C
738             return (159-colM)/32 + 5*(iModule%2);
739           }
740         }
741       }
742     }
743   }
744   return 10; // error
745 }
746 //__________________________________________________________________________
747 UInt_t AliITSRawStreamSPD::GetOnlineColFromOffline(UInt_t module, UInt_t colM) {
748   // offline->online (col)
749   if (module<80) { // inner layer
750     return colM%32;
751   }
752   else if (module<240) { // outer layer
753     return colM%32;
754   }
755   return 32; // error
756 }
757 //__________________________________________________________________________
758 UInt_t AliITSRawStreamSPD::GetOnlineRowFromOffline(UInt_t module, UInt_t rowM) {
759   // offline->online (row)
760   if (module<80) { // inner layer
761     return (255-rowM);
762   }
763   else if (module<240) { // outer layer
764     return (255-rowM);
765   }
766   return 256; // error
767 }
768 //__________________________________________________________________________
769 UInt_t AliITSRawStreamSPD::GetOfflineModuleFromOnline(UInt_t eqId, UInt_t hs, UInt_t chip) {
770   // online->offline (module)
771   if (eqId<20 && hs<6 && chip<10) return fgkDDLModuleMap[eqId][hs*2+chip/5];
772   else return 240;
773 }
774 //__________________________________________________________________________
775 UInt_t AliITSRawStreamSPD::GetOfflineChipKeyFromOnline(UInt_t eqId, UInt_t hs, UInt_t chip) {
776   // online->offline (chip key: 0-1199)
777   if (eqId<20 && hs<6 && chip<10) {
778     UInt_t module = GetOfflineModuleFromOnline(eqId,hs,chip);
779     UInt_t chipInModule = ( chip>4 ? chip-5 : chip ); 
780     if(eqId>9) chipInModule = 4 - chipInModule;  // side C only
781     return (module*5 + chipInModule);
782   } else return 1200;
783 }
784 //__________________________________________________________________________
785 UInt_t AliITSRawStreamSPD::GetOfflineColFromOnline(UInt_t eqId, UInt_t hs, UInt_t chip, UInt_t col) {
786   // online->offline (col)
787   if (eqId>=20 || hs>=6 || chip>=10 || col>=32) return 160; // error
788   UInt_t offset = 32 * (chip % 5);
789   if (hs<2) {
790     if (eqId<10) {
791       return 159 - (31-col + offset); // inner layer, side A
792     }
793     else {
794       return col + offset; // inner layer, side C
795     }
796   }
797   else {
798     if (eqId<10) {
799       return (col + offset); // outer layer, side A
800     }
801     else {
802       return 159 - (31-col + offset); // outer layer, side C
803     }
804   }
805 }
806 //__________________________________________________________________________
807 UInt_t AliITSRawStreamSPD::GetOfflineRowFromOnline(UInt_t eqId, UInt_t hs, UInt_t chip, UInt_t row) {
808   // online->offline (row)
809   if (eqId>=20 || hs>=6 || chip>=10 || row>=256) return 256; // error
810   return 255-row;
811 }
812