/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
* *
* Author: The ALICE Off-line Project. *
* Contributors are mentioned in the code where appropriate. *
#include "AliITSRawStreamSPD.h"
#include "AliRawReader.h"
+#include "AliLog.h"
ClassImp(AliITSRawStreamSPD)
-
const Int_t AliITSRawStreamSPD::fgkDDLModuleMap[kDDLsNumber][kModulesPerDDL] = {
- { 0, 1, 4, 5, 80, 81, 84, 85, 88, 89, 92, 93},
- { 8, 9,12,13, 96, 97,100,101,104,105,108,109},
- {16,17,20,21,112,113,116,117,120,121,124,125},
- {24,25,28,29,128,129,132,133,136,137,140,141},
- {32,33,36,37,144,145,148,149,152,153,156,157},
- {40,41,44,45,160,161,164,165,168,169,172,173},
- {48,49,52,53,176,177,180,181,184,185,188,189},
- {56,57,60,61,192,193,196,197,200,201,204,205},
- {64,65,68,69,208,209,212,213,216,217,220,221},
- {72,73,76,77,224,225,228,229,232,233,236,237},
- { 2, 3, 6, 7, 82, 83, 86, 87, 90, 91, 94, 95},
- {10,11,14,15, 98, 99,102,103,106,107,110,111},
- {18,19,22,23,114,115,118,119,122,123,126,127},
- {26,27,30,31,130,131,134,135,138,139,142,143},
- {34,35,38,39,146,147,150,151,154,155,158,159},
- {42,43,46,47,162,163,166,167,170,171,174,175},
- {50,51,54,55,178,179,182,183,186,187,190,191},
- {58,59,62,63,194,195,198,199,202,203,206,207},
- {66,67,70,71,210,211,214,215,218,219,222,223},
- {74,75,78,79,226,227,230,231,234,235,238,239}
+ { 4, 5, 0, 1, 80, 81, 84, 85, 88, 89, 92, 93},
+ {12,13, 8, 9, 96, 97,100,101,104,105,108,109},
+ {20,21,16,17,112,113,116,117,120,121,124,125},
+ {28,29,24,25,128,129,132,133,136,137,140,141},
+ {36,37,32,33,144,145,148,149,152,153,156,157},
+ {44,45,40,41,160,161,164,165,168,169,172,173},
+ {52,53,48,49,176,177,180,181,184,185,188,189},
+ {60,61,56,57,192,193,196,197,200,201,204,205},
+ {68,69,64,65,208,209,212,213,216,217,220,221},
+ {76,77,72,73,224,225,228,229,232,233,236,237},
+ { 7, 6, 3, 2, 83, 82, 87, 86, 91, 90, 95, 94},
+ {15,14,11,10, 99, 98,103,102,107,106,111,110},
+ {23,22,19,18,115,114,119,118,123,122,127,126},
+ {31,30,27,26,131,130,135,134,139,138,143,142},
+ {39,38,35,34,147,146,151,150,155,154,159,158},
+ {47,46,43,42,163,162,167,166,171,170,175,174},
+ {55,54,51,50,179,178,183,182,187,186,191,190},
+ {63,62,59,58,195,194,199,198,203,202,207,206},
+ {71,70,67,66,211,210,215,214,219,218,223,222},
+ {79,78,75,74,227,226,231,230,235,234,239,238}
};
-
-
+//__________________________________________________________________________
AliITSRawStreamSPD::AliITSRawStreamSPD(AliRawReader* rawReader) :
AliITSRawStream(rawReader),
- fData(0),
- fDDLNumber(-1),
- fEventNumber(-1),
- fOffset(0),
- fHitCount(0),
- fDataChar1(0),
- fDataChar2(0),
- fDataChar3(0),
- fDataChar4(0),
- fFirstWord(kTRUE)
+ fEventCounter(-1),fChipAddr(0),fHalfStaveNr(0),fCol(0),fRow(0),fCalHeadLen(0),
+ fData(0),fOffset(0),fHitCount(0),
+ fDataChar1(0),fDataChar2(0),fDataChar3(0),fDataChar4(0),
+ fFirstWord(kTRUE),fPrevEventId(0xffffffff),
+ fEqPLBytesRead(0),fEqPLChipHeadersRead(0),fEqPLChipTrailersRead(0),
+ fHeaderOrTrailerReadLast(kFALSE),fExpectedHeaderTrailerCount(0),
+ fFillOutOfSynch(kFALSE),fDDLID(-1),fLastDDLID(-1),fAdvancedErrorLog(kFALSE),
+ fAdvLogger(NULL)
{
-// create an object to read ITS SPD raw digits
-
+ // create an object to read ITS SPD raw digits
+ fRawReader->Reset();
fRawReader->Select("ITSSPD");
+ // reset calib header words
+ for (UInt_t iword=0; iword<kCalHeadLenMax; iword++) {
+ fCalHeadWord[iword]=0xffffffff;
+ }
+ for (UInt_t eq=0; eq<20; eq++) {
+ fActiveEq[eq]=kFALSE;
+ for (UInt_t hs=0; hs<6; hs++) {
+ fActiveHS[eq][hs]=kFALSE;
+ for (UInt_t chip=0; chip<10; chip++) {
+ fActiveChip[eq][hs][chip]=kFALSE;
+ fEventCounterFull[eq][hs][chip] = -1;
+ }
+ }
+ }
+ NewEvent();
}
-
-Bool_t AliITSRawStreamSPD::ReadNextShort()
+//__________________________________________________________________________
+AliITSRawStreamSPD::AliITSRawStreamSPD(const AliITSRawStreamSPD& rstream) :
+ AliITSRawStream(rstream.fRawReader),
+ fEventCounter(-1),fChipAddr(0),fHalfStaveNr(0),fCol(0),fRow(0),fCalHeadLen(0),
+ fData(0),fOffset(0),fHitCount(0),
+ fDataChar1(0),fDataChar2(0),fDataChar3(0),fDataChar4(0),
+ fFirstWord(kTRUE),fPrevEventId(0xffffffff),
+ fEqPLBytesRead(0),fEqPLChipHeadersRead(0),fEqPLChipTrailersRead(0),
+ fHeaderOrTrailerReadLast(kFALSE),fExpectedHeaderTrailerCount(0),
+ fFillOutOfSynch(kFALSE),fDDLID(-1),fLastDDLID(-1),fAdvancedErrorLog(kFALSE),
+ fAdvLogger(NULL)
{
+ // copy constructor
+ AliError("Copy constructor should not be used.");
+}
+//__________________________________________________________________________
+AliITSRawStreamSPD& AliITSRawStreamSPD::operator=(const AliITSRawStreamSPD& rstream) {
+ // assignment operator
+ if (this!=&rstream) {}
+ AliError("Assignment opertator should not be used.");
+ return *this;
+}
+//__________________________________________________________________________
+Bool_t AliITSRawStreamSPD::ReadNextShort() {
+ // read next 16 bit word into fData
if (fFirstWord) {
- fFirstWord=kFALSE;
Bool_t b1 = fRawReader->ReadNextChar(fDataChar1);
if (!b1) return kFALSE;
Bool_t b2, b3, b4;
return kFALSE;
}
fData = fDataChar3+(fDataChar4<<8);
+ const UInt_t *intPtr = fRawReader->GetEventId();
+ if (intPtr!=0) {
+ if (*intPtr!=fPrevEventId) { // if new event...
+ NewEvent();
+ fPrevEventId=*intPtr;
+ }
+ }
+ fFirstWord=kFALSE;
}
else {
fFirstWord=kTRUE;
fData = fDataChar1+(fDataChar2<<8);
}
-
+ fEqPLBytesRead+=2;
+
return kTRUE;
}
+//__________________________________________________________________________
+Bool_t AliITSRawStreamSPD::ReadNextInt() {
+ // reads next 32 bit into fDataChar1..4
+ // (if first 16 bits read already, just completes the present word)
+ if (fFirstWord) {
+ if (ReadNextShort() && ReadNextShort()) {
+ return kTRUE;
+ }
+ }
+ else {
+ if (ReadNextShort()) {
+ return kTRUE;
+ }
+ }
+ return kFALSE;
+}
+//__________________________________________________________________________
+void AliITSRawStreamSPD::NewEvent() {
+ // call this to reset flags for a new event
+ for (UInt_t eqId=0; eqId<20; eqId++) {
+ fCalHeadRead[eqId]=kFALSE;
+ }
+ fEventCounter = -1;
+ fDDLID = -1;
+ fLastDDLID = -1;
+ for (UInt_t eq=0; eq<20; eq++) {
+ for (UInt_t hs=0; hs<6; hs++) {
+ for (UInt_t chip=0; chip<10; chip++) {
+ fFastOrSignal[eq][hs][chip] = kFALSE;
+ }
+ }
+ }
+}
+//__________________________________________________________________________
+Int_t AliITSRawStreamSPD::ReadCalibHeader() {
+ // needs to be called in the beginning of each equipment data
+ // read the extra calibration header
+ // returns the length of the header if it is present, -1 otherwise
-void AliITSRawStreamSPD::SkipCalibHeader()
-{
- // Checks if there is an extra calibration header
- // present in the raw data. Reads past this in that case.
-
- fRawReader->ReadHeader(); // need this to get access to the block attributes
- UChar_t attr = fRawReader->GetBlockAttributes();
- if ((attr & 0x40) == 0x40) { // is the header present?
- Bool_t b1, b2, b3, b4;
- b1 = fRawReader->ReadNextChar(fDataChar1);
- b2 = fRawReader->ReadNextChar(fDataChar2);
- b3 = fRawReader->ReadNextChar(fDataChar3);
- b4 = fRawReader->ReadNextChar(fDataChar4);
- if (b1 && b2 && b3 && b4) {
+ Int_t ddlID = fRawReader->GetDDLID();
+ if (ddlID==-1) { // we may need to read one word to get the blockAttr
+ if (!ReadNextShort()) return -1;
+ ddlID = fRawReader->GetDDLID();
+ if(ddlID<0){
+ AliError("fRawReader->GetDDLID() returns a negative value");
+ ddlID=0;
+ }
+ }
+ // reset flags and counters
+ fEqPLBytesRead = 2;
+ fEqPLChipHeadersRead = 0;
+ fEqPLChipTrailersRead = 0;
+ fHeaderOrTrailerReadLast = kFALSE;
+ fFillOutOfSynch = kFALSE;
+
+ // check router error bits:
+ UInt_t statusBits = fRawReader->GetStatusBits();
+ if ((statusBits >> 5) & 1) { // linkrx/detector fatal error bit
+ TString errMess = "LinkRx Error Bit Set";
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kLinkRxDetectorFatalErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kLinkRxDetectorFatalErr,ddlID,-1,-1,errMess.Data());
+ }
+ if ((statusBits >> 12) & 1) { // trigger sequence monitor error bit
+ TString errMess = "TSM Trigger Error Bit Set";
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kTSMtriggerErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kTSMtriggerErr,ddlID,-1,-1,errMess.Data());
+ }
+ if (fRawReader->TestBlockAttribute(7)) { // bunch crossing difference error bit
+ TString errMess = "High Multiplicity Event Flag Set";
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kHighMultiplicityFlag,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHighMultiplicityFlag,ddlID,-1,-1,errMess.Data());
+ }
+
+ // set eq active and the participating half-staves (for access from outside this class),
+ // set what number of chip headers/trailers to expect (for the moment not used)
+ fActiveEq[ddlID]=kTRUE;
+ fExpectedHeaderTrailerCount = 0;
+ for (UInt_t hs=0; hs<6; hs++) {
+ if (!fRawReader->TestBlockAttribute(hs)) {
+ fActiveHS[ddlID][hs]=kTRUE;
+ fExpectedHeaderTrailerCount+=10;
+ }
+ else {
+ fActiveHS[ddlID][hs]=kFALSE;
+ }
+ }
+
+ if (ddlID>=0 && ddlID<20) fCalHeadRead[ddlID]=kTRUE;
+ if (fRawReader->TestBlockAttribute(6)) { // is the calib header present?
+ if (ReadNextInt()) {
// length of cal header:
- UInt_t calLen = fDataChar1+(fDataChar2<<8)+(fDataChar3<<16)+(fDataChar4<<24);
- // read past the cal header:
- UInt_t tmpData;
- for (UInt_t iword=0; iword<calLen; iword++) {
- fRawReader->ReadNextInt(tmpData);
+ fCalHeadLen = fDataChar1+(fDataChar2<<8)+(fDataChar3<<16)+(fDataChar4<<24);
+ if (fCalHeadLen>kCalHeadLenMax) {
+ TString errMess = Form("Header length %d > max = %d",fCalHeadLen,kCalHeadLenMax);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kCalHeaderLengthErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kCalHeaderLengthErr,ddlID,-1,-1,errMess.Data());
+ return -1;
+ }
+ else {
+ for (UInt_t iword=0; iword<fCalHeadLen; iword++) {
+ if (ReadNextInt()) {
+ fCalHeadWord[iword] = fDataChar1+(fDataChar2<<8)+(fDataChar3<<16)+(fDataChar4<<24);
+ }
+ else {
+ TString errMess = "Header length problem";
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kCalHeaderLengthErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kCalHeaderLengthErr,ddlID,-1,-1,errMess.Data());
+ return -1;
+ }
+ }
+ return fCalHeadLen;
}
}
}
-}
-Bool_t AliITSRawStreamSPD::Next()
-{
+ return -1;
+}
+//__________________________________________________________________________
+Bool_t AliITSRawStreamSPD::Next() {
// read the next raw digit
// returns kFALSE if there is no digit left
fPrevModuleID = fModuleID;
- // while (fRawReader->ReadNextShort(fData)) {
+
while (ReadNextShort()) {
+
+ fLastDDLID = fDDLID;
+ fDDLID = fRawReader->GetDDLID();
+ if (fDDLID>=0 && fDDLID<20) {
+ if (!fCalHeadRead[fDDLID]) {
+ if (fLastDDLID!=-1) { // if not the first equipment for this event
+ fEqPLBytesRead -= 2;
+ // The next line is commented, since we do not have the information how many chips to expect
+ // CheckHeaderAndTrailerCount(fLastDDLID);
+ }
+ 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)
+ }
+ }
+ else {
+ TString errMess = Form("Error in DDL number (=%d) , setting it to 19",fDDLID);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kDDLNumberErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->AddMessage(errMess.Data());
+ fDDLID=19;
+ }
+
+ Short_t hs = ((fData & 0x3800)>>11) ; Short_t chip=(fData & 0x000F); // -> for printouts
+
if ((fData & 0xC000) == 0x4000) { // header
+ if (fHeaderOrTrailerReadLast) {
+ TString errMess = Form("Chip trailer missing - (eq,hs,chip) = (%d,%d,%d)",fDDLID,hs,chip);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kTrailerMissingErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kTrailerMissingErr,fLastDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ }
+ fHeaderOrTrailerReadLast = kTRUE;
+ fEqPLChipHeadersRead++;
fHitCount = 0;
- UShort_t eventNumber = (fData >> 4) & 0x007F;
- if (fEventNumber < 0) {
- fEventNumber = eventNumber;
+ UShort_t eventCounter = (fData >> 4) & 0x007F;
+ if (fEventCounter < 0) {
+ fEventCounter = eventCounter;
}
- else if (eventNumber != fEventNumber) {
- Warning("Next", "mismatching event numbers: %d != %d",
- eventNumber, fEventNumber);
- }
- UShort_t chipAddr = fData & 0x000F;
- if (chipAddr>9) {
- Error("Next", "overflow chip addr (= %d) , setting it to 9", chipAddr);
- chipAddr=9;
- }
- UShort_t halfStaveNr = (fData & 0x3800)>>11;
- if (halfStaveNr>5 || fRawReader->TestBlockAttribute(halfStaveNr)) {
- Error("Next", "half stave number error(= %d) , setting it to 5", halfStaveNr);
- halfStaveNr=5;
- }
- fDDLNumber = fRawReader->GetDDLID();
- if (fDDLNumber>19 || fDDLNumber<0) {
- Error("Next", "DDL number error (= %d) , setting it to 19", fDDLNumber);
- fDDLNumber=19;
- }
- fModuleID = fgkDDLModuleMap[fDDLNumber][halfStaveNr*2+chipAddr/5];
- fOffset = 32 * (chipAddr % 5);
- }
+ else if (eventCounter != fEventCounter) {
+ TString errMess;
+ if (fEqPLChipHeadersRead==1) {
+ errMess = Form("Mismatching event counters between this equipment and the previous: %d != %d - (eq,hs,chip) = (%d,%d,%d)",
+ eventCounter,fEventCounter,fDDLID,hs,chip);
+ }
+ else {
+ errMess = Form("Mismatching event counters between this chip header and the previous: %d != %d - (eq,hs,chip) = (%d,%d,%d)",
+ eventCounter,fEventCounter,fDDLID,hs,chip);
+ }
+ fEventCounter = eventCounter;
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kEventCounterErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kEventCounterErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ }
+ fChipAddr = fData & 0x000F;
+ if (fChipAddr>9) {
+ TString errMess = Form("Overflow chip address %d - set to 0 (eq,hs) = (%d,%d)",fChipAddr,fDDLID,hs);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kChipAddrErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kChipAddrErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ fChipAddr=0;
+ }
+ fHalfStaveNr = (fData & 0x3800)>>11;
+ if (fHalfStaveNr>5 || fRawReader->TestBlockAttribute(fHalfStaveNr)) {
+ TString errMess = Form("Half stave number error: %d - set to 0 - eq %d",fHalfStaveNr,fDDLID);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kHSNumberErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHSNumberErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ fHalfStaveNr=0;
+ }
+ fActiveChip[fDDLID][fHalfStaveNr][fChipAddr]=kTRUE;
+ fEventCounterFull[fDDLID][fHalfStaveNr][fChipAddr] = eventCounter;
+ // translate ("online") ddl, hs, chip nr to ("offline") module id :
+ fModuleID = GetOfflineModuleFromOnline(fDDLID,fHalfStaveNr,fChipAddr);
+ }
+
else if ((fData & 0xC000) == 0x0000) { // trailer
- UShort_t hitCount = fData & 0x1FFF;
- if (hitCount != fHitCount) Error("Next", "wrong number of hits: %d != %d", fHitCount, hitCount);
- }
+ if ( (fEqPLBytesRead+fFillOutOfSynch*2)%4 != 0 ) {
+ TString errMess = Form("Fill word is missing - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kFillMissingErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kFillMissingErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ if (fFillOutOfSynch) fFillOutOfSynch = kFALSE;
+ else fFillOutOfSynch = kTRUE;
+ }
+ if (!fHeaderOrTrailerReadLast) {
+ TString errMess = Form("Trailer without previous header - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kTrailerWithoutHeaderErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kTrailerWithoutHeaderErr,fLastDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ }
+ fHeaderOrTrailerReadLast = kFALSE;
+ fEqPLChipTrailersRead++;
+ UShort_t hitCount = fData & 0x0FFF;
+ if (hitCount != fHitCount){
+ TString errMess = Form("Number of hits %d, while %d expected - (eq,hs,chip) = (%d,%d,%d)",fHitCount,hitCount,fDDLID,fHalfStaveNr,fChipAddr);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kNumberHitsErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kNumberHitsErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ }
+ Bool_t errorBit = fData & 0x1000;
+ if (errorBit) {
+ TString errMess = Form("Trailer error bit set for chip (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kTrailerErrorBitErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kTrailerErrorBitErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ }
+ Bool_t fastorBit = fData & 0x2000;
+ fFastOrSignal[fDDLID][fHalfStaveNr][fChipAddr] = fastorBit;
+ }
+
else if ((fData & 0xC000) == 0x8000) { // pixel hit
+ if (!fHeaderOrTrailerReadLast) {
+ TString errMess = Form("Chip header missing - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kHeaderMissingErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHeaderMissingErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ fHeaderOrTrailerReadLast = kTRUE;
+ }
fHitCount++;
- fCoord1 = (fData & 0x001F) + fOffset;
- fCoord2 = (fData >> 5) & 0x00FF;
+ fCol = (fData & 0x001F);
+ fRow = (fData >> 5) & 0x00FF;
+
+ fCoord1 = GetOfflineColFromOnline(fDDLID,fHalfStaveNr,fChipAddr,fCol);
+ fCoord2 = GetOfflineRowFromOnline(fDDLID,fHalfStaveNr,fChipAddr,fRow);
+
return kTRUE;
- }
+ }
+
else { // fill word
- if ((fData & 0xC000) != 0xC000) Error("Next", "wrong fill word!");
+ if ((fData & 0xC000) != 0xC000) {
+ TString errMess = Form("Wrong fill word! - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kWrongFillWordErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kWrongFillWordErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ }
+ if ( (fEqPLBytesRead+fFillOutOfSynch*2)%4 != 2 ) {
+ TString errMess = Form("Fill word is unexpected - (eq,hs,chip) = (%d,%d,%d)",fDDLID,fHalfStaveNr,fChipAddr);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kFillUnexpectErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kFillUnexpectErr,fDDLID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ if (fFillOutOfSynch) fFillOutOfSynch = kFALSE;
+ else fFillOutOfSynch = kTRUE;
+ }
}
}
-
+ if (fDDLID>=0 && fDDLID<20) {
+ // The next line is commented, since we do not have the information how many chips to expect
+ // CheckHeaderAndTrailerCount(fDDLID);
+ }
return kFALSE;
}
+//__________________________________________________________________________
+void AliITSRawStreamSPD::CheckHeaderAndTrailerCount(Int_t ddlID) {
+ // Checks that the number of header and trailers found for the ddl are as expected
+ if (fEqPLChipHeadersRead != fExpectedHeaderTrailerCount) {
+ TString errMess = Form("Chip header count inconsistent %d != %d (expected) for ddl %d",
+ fEqPLChipHeadersRead,fExpectedHeaderTrailerCount,ddlID);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kHeaderCountErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHeaderCountErr,ddlID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ }
+ if (fEqPLChipTrailersRead != fExpectedHeaderTrailerCount) {
+ TString errMess = Form("Chip trailer count inconsistent %d != %d (expected) for ddl %d",
+ fEqPLChipTrailersRead,fExpectedHeaderTrailerCount,ddlID);
+ AliError(errMess.Data());
+ fRawReader->AddMajorErrorLog(kHeaderCountErr,errMess.Data());
+ if (fAdvancedErrorLog) fAdvLogger->ProcessError(kHeaderCountErr,ddlID,fEqPLBytesRead,fEqPLChipHeadersRead,errMess.Data());
+ }
+}
+//__________________________________________________________________________
+void AliITSRawStreamSPD::ActivateAdvancedErrorLog(Bool_t activate, AliITSRawStreamSPDErrorLog* advLogger) {
+ // Activate the advanced error logging.
+ // Logger object has to be created elsewhere and a pointer sent here.
+ fAdvancedErrorLog=activate;
+ if (activate && advLogger!=NULL) {
+ fAdvLogger = advLogger;
+ }
+ if (fAdvLogger==NULL) {
+ fAdvancedErrorLog=kFALSE;
+ }
+}
+//__________________________________________________________________________
+const Char_t* AliITSRawStreamSPD::GetErrorName(UInt_t errorCode) {
+ // Returns a string for each error code
+ if (errorCode==kTotal) return "All Errors";
+ else if (errorCode==kHeaderMissingErr) return "Header Missing";
+ else if (errorCode==kTrailerMissingErr) return "Trailer Missing";
+ else if (errorCode==kTrailerWithoutHeaderErr) return "Trailer Unexpected";
+ else if (errorCode==kHeaderCountErr) return "Header Count Wrong";
+ else if (errorCode==kTrailerCountErr) return "Trailer Count Wrong";
+ else if (errorCode==kFillUnexpectErr) return "Fill Unexpected";
+ else if (errorCode==kFillMissingErr) return "Fill Missing";
+ else if (errorCode==kWrongFillWordErr) return "Fill Word Wrong";
+ else if (errorCode==kNumberHitsErr) return "Hit Count Wrong";
+ else if (errorCode==kEventCounterErr) return "Event Counter Error";
+ else if (errorCode==kDDLNumberErr) return "DDL Number Error";
+ else if (errorCode==kHSNumberErr) return "HS Number Error";
+ else if (errorCode==kChipAddrErr) return "Chip Address Error";
+ else if (errorCode==kCalHeaderLengthErr) return "Calib Header Length Error";
+ else if (errorCode==kAdvEventCounterErr) return "Event Counter Error (Adv)";
+ else if (errorCode==kAdvEventCounterOrderErr) return "Event Counter Jump Error (Adv)";
+ else if (errorCode==kTrailerErrorBitErr) return "Trailer Error Bit Set";
+ else if (errorCode==kLinkRxDetectorFatalErr) return "LinkRx/Detector Fatal Error Bit Set";
+ else if (errorCode==kTSMtriggerErr) return "TSM Trigger Error Bit Set";
+ else if (errorCode==kHighMultiplicityFlag) return "High Multiplicity Event Flag Set";
+ else return "";
+}
+//__________________________________________________________________________
+Bool_t AliITSRawStreamSPD::GetFastOrSignal(UInt_t eq, UInt_t hs, UInt_t chip) {
+ // returns if there was a fastor signal from this chip
+ if (eq>=20 || hs>=6 || chip>=10) {
+ TString errMess = Form("eq,hs,chip = %d,%d,%d out of bounds. Return kFALSE.",eq,hs,chip);
+ AliError(errMess.Data());
+ return kFALSE;
+ }
+ return fFastOrSignal[eq][hs][chip];
+}
+//__________________________________________________________________________
+Bool_t AliITSRawStreamSPD::IsActiveEq(UInt_t eq) const {
+ // returns if the eq is active (seen in data)
+ if (eq>=20) {
+ TString errMess = Form("eq = %d out of bounds. Return kFALSE.",eq);
+ AliError(errMess.Data());
+ return kFALSE;
+ }
+ return fActiveEq[eq];
+}
+//__________________________________________________________________________
+Bool_t AliITSRawStreamSPD::IsActiveHS(UInt_t eq, UInt_t hs) const {
+ // returns if the hs is active (info from block attr)
+ if (eq>=20 || hs>=6) {
+ TString errMess = Form("eq,hs = %d,%d out of bounds. Return kFALSE.",eq,hs);
+ AliError(errMess.Data());
+ return kFALSE;
+ }
+ return fActiveHS[eq][hs];
+}
+//__________________________________________________________________________
+Bool_t AliITSRawStreamSPD::IsActiveChip(UInt_t eq, UInt_t hs, UInt_t chip) const {
+ // returns if the chip is active (seen in data)
+ if (eq>=20 || hs>=6 || chip>=10) {
+ TString errMess = Form("eq,hs,chip = %d,%d,%d out of bounds. Return kFALSE.",eq,hs,chip);
+ AliError(errMess.Data());
+ return kFALSE;
+ }
+ return fActiveChip[eq][hs][chip];
+}
+//__________________________________________________________________________
+Bool_t AliITSRawStreamSPD::GetHalfStavePresent(UInt_t hs) {
+ // Reads the half stave present status from the block attributes
+ // This is not really needed anymore (kept for now in case it is still used somewhere).
+ // The same information can be reached through the "IsActiveHS" method instead.
+ Int_t ddlID = fRawReader->GetDDLID();
+ if (ddlID==-1) {
+ AliWarning("DDL ID = -1. Cannot read block attributes. Return kFALSE.");
+ return kFALSE;
+ }
+ else {
+ if (hs>=6) {
+ AliWarning(Form("HS >= 6 requested (%d). Return kFALSE.",hs));
+ return kFALSE;
+ }
+ UChar_t attr = fRawReader->GetBlockAttributes();
+ if (((attr>>hs) & 0x01) == 0x01) { // bit set means not present
+ return kFALSE;
+ }
+ else {
+ return kTRUE;
+ }
+ }
+}
+//__________________________________________________________________________
+Bool_t AliITSRawStreamSPD::IsEventCounterFullConsistent() const {
+ // checks if the event counter values are consistent within the event
+ Short_t reference = -1;
+ for (UInt_t eq=0; eq<20; eq++) {
+ if (IsActiveEq(eq)) {
+ for (UInt_t hs=0; hs<6; hs++) {
+ if (IsActiveHS(eq,hs)) {
+ for (UInt_t chip=0; chip<10; chip++) {
+ if (fEventCounterFull[eq][hs][chip]!=-1) {
+ if (reference==-1) reference = fEventCounterFull[eq][hs][chip];
+ if (fEventCounterFull[eq][hs][chip] != reference) return kFALSE;
+ }
+ }
+ }
+ }
+ }
+ }
+ return kTRUE;
+}
+//__________________________________________________________________________
+Short_t AliITSRawStreamSPD::GetEventCounterFullEq(UInt_t eq) const {
+ // if the eq is active; returns the event counter value
+ if (eq>=20) {
+ TString errMess = Form("eq (%d) out of bounds",eq);
+ AliError(errMess.Data());
+ return -1;
+ }
+ if (IsActiveEq(eq)) {
+ for (UInt_t hs=0; hs<6; hs++) {
+ if (IsActiveHS(eq,hs)) {
+ for (UInt_t chip=0; chip<10; chip++) {
+ if (fEventCounterFull[eq][hs][chip]!=-1) {
+ return fEventCounterFull[eq][hs][chip];
+ }
+ }
+ }
+ }
+ }
+ return -1;
+}
+//__________________________________________________________________________
+Short_t AliITSRawStreamSPD::GetEventCounterFullHS(UInt_t eq, UInt_t hs) const {
+ // if the eq,hs is active; returns the event counter value
+ if (eq>=20 || hs>=6) {
+ TString errMess = Form("eq,hs (%d,%d) out of bounds",eq,hs);
+ AliError(errMess.Data());
+ return -1;
+ }
+ if (IsActiveEq(eq)) {
+ if (IsActiveHS(eq,hs)) {
+ for (UInt_t chip=0; chip<10; chip++) {
+ if (fEventCounterFull[eq][hs][chip]!=-1) {
+ return fEventCounterFull[eq][hs][chip];
+ }
+ }
+ }
+ }
+ return -1;
+}
+//__________________________________________________________________________
+Short_t AliITSRawStreamSPD::GetEventCounterFullChip(UInt_t eq, UInt_t hs, UInt_t chip) const {
+ // if the eq,hs,chip is active; returns the event counter value
+ if (eq>=20 || hs>=6 || chip>=10) {
+ TString errMess = Form("eq,hs,chip (%d,%d,%d) out of bounds",eq,hs,chip);
+ AliError(errMess.Data());
+ return -1;
+ }
+ if (IsActiveEq(eq)) {
+ if (IsActiveHS(eq,hs)) {
+ if (IsActiveChip(eq,hs,chip)) {
+ return fEventCounterFull[eq][hs][chip];
+ }
+ }
+ }
+ return -1;
+}
+//__________________________________________________________________________
+Int_t AliITSRawStreamSPD::GetHword(UInt_t index) {
+ if (index<kCalHeadLenMax) return fCalHeadWord[index];
+ else return 0;
+}
+Bool_t AliITSRawStreamSPD::GetHhalfStaveScanned(UInt_t hs) const {
+ if (hs<6) return (Bool_t)((fCalHeadWord[0]>>(6+hs)) & (0x00000001));
+ else return kFALSE;
+}
+Bool_t AliITSRawStreamSPD::GetHchipPresent(UInt_t hs, UInt_t chip) const {
+ if (hs<6 && chip<10) return ((( fCalHeadWord[hs/3+3]>>((hs%3)*10+chip)) & 0x00000001) == 1);
+ else return kFALSE;
+}
+UInt_t AliITSRawStreamSPD::GetHdacHigh(UInt_t hs) const {
+ if (hs<6) return (fCalHeadWord[hs/2+7]>>(24-16*(hs%2)) & 0x000000ff);
+ else return 0;
+}
+UInt_t AliITSRawStreamSPD::GetHdacLow(UInt_t hs) const {
+ if (hs<6) return (fCalHeadWord[hs/2+7]>>(16-16*(hs%2)) & 0x000000ff);
+ else return 0;
+}
+UInt_t AliITSRawStreamSPD::GetHTPAmp(UInt_t hs) const {
+ if (hs<6) return fCalHeadWord[hs+10];
+ else return 0;
+}
+Bool_t AliITSRawStreamSPD::GetHminTHchipPresent(UInt_t chip) const {
+ if (chip<10) return ((( fCalHeadWord[7]>>(16+chip)) & 0x00000001) == 1);
+ else return kFALSE;
+}
+UInt_t AliITSRawStreamSPD::GetFOHnumDacs() const {
+ return (fCalHeadLen-37)/2;
+}
+UInt_t AliITSRawStreamSPD::GetFOHdacIndex(UInt_t index) const {
+ if (index>=GetFOHnumDacs()) {
+ TString errMess = Form("Only %d DACs in this run, returning 0",GetFOHnumDacs());
+ AliError(errMess.Data());
+ return 0;
+ }
+ return fCalHeadWord[7+index*2];
+}
+UInt_t AliITSRawStreamSPD::GetFOHdacValue(UInt_t index) const {
+ if (index>=GetFOHnumDacs()) {
+ TString errMess = Form("Only %d DACs in this run, returning 0",GetFOHnumDacs());
+ AliError(errMess.Data());
+ return 0;
+ }
+ return fCalHeadWord[7+1+index*2];
+}
+UInt_t AliITSRawStreamSPD::GetFOHchipCount(UInt_t hs, UInt_t chip) const {
+ if (hs<6 && chip<10) {
+ if (chip%2==0) {
+ return ((fCalHeadWord[7 + GetFOHnumDacs()*2 + (hs*10 + chip)/2] >> 16) & 0x0000ffff);
+ }
+ else {
+ return (fCalHeadWord[7 + GetFOHnumDacs()*2 + (hs*10 + chip)/2] & 0x0000ffff);
+ }
+ }
+ else return 0;
+}
+//__________________________________________________________________________
+Int_t AliITSRawStreamSPD::GetModuleNumber(UInt_t iDDL, UInt_t iModule) {
+ if (iDDL<20 && iModule<12) return fgkDDLModuleMap[iDDL][iModule];
+ else return 240;
+}
+//__________________________________________________________________________
+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) {
+ // converts offline coordinates to online
+ eq = GetOnlineEqIdFromOffline(module);
+ hs = GetOnlineHSFromOffline(module);
+ chip = GetOnlineChipFromOffline(module,colM);
+ col = GetOnlineColFromOffline(module,colM);
+ row = GetOnlineRowFromOffline(module,rowM);
+ if (eq>=20 || hs>=6 || chip>=10 || col>=32 || row>=256) return kFALSE;
+ else return kTRUE;
+}
+//__________________________________________________________________________
+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) {
+ // converts online coordinates to offline
+ module = GetOfflineModuleFromOnline(eq,hs,chip);
+ colM = GetOfflineColFromOnline(eq,hs,chip,col);
+ rowM = GetOfflineRowFromOnline(eq,hs,chip,row);
+ if (module>=240 || colM>=160 || rowM>=256) return kFALSE;
+ else return kTRUE;
+}
+//__________________________________________________________________________
+UInt_t AliITSRawStreamSPD::GetOnlineEqIdFromOffline(UInt_t module) {
+ // offline->online (eq)
+ for (UInt_t eqId=0; eqId<20; eqId++) {
+ for (UInt_t iModule=0; iModule<12; iModule++) {
+ if (GetModuleNumber(eqId,iModule)==(Int_t)module) return eqId;
+ }
+ }
+ return 20; // error
+}
+//__________________________________________________________________________
+UInt_t AliITSRawStreamSPD::GetOnlineHSFromOffline(UInt_t module) {
+ // offline->online (hs)
+ for (UInt_t eqId=0; eqId<20; eqId++) {
+ for (UInt_t iModule=0; iModule<12; iModule++) {
+ if (GetModuleNumber(eqId,iModule)==(Int_t)module) return iModule/2;
+ }
+ }
+ return 6; // error
+}
+//__________________________________________________________________________
+UInt_t AliITSRawStreamSPD::GetOnlineChipFromOffline(UInt_t module, UInt_t colM) {
+ // offline->online (chip)
+ for (UInt_t eq=0; eq<20; eq++) {
+ for (UInt_t iModule=0; iModule<12; iModule++) {
+ if (GetModuleNumber(eq,iModule)==(Int_t)module) {
+ if (module<80) {
+ if (eq<10) { // side A
+ return (159-colM)/32 + 5*(iModule%2);
+ }
+ else { // side C
+ return colM/32 + 5*(iModule%2);
+ }
+ }
+ else if (module<240) {
+ if (eq<10) { // side A
+ return colM/32 + 5*(iModule%2);
+ }
+ else { // side C
+ return (159-colM)/32 + 5*(iModule%2);
+ }
+ }
+ }
+ }
+ }
+ return 10; // error
+}
+//__________________________________________________________________________
+UInt_t AliITSRawStreamSPD::GetOnlineColFromOffline(UInt_t module, UInt_t colM) {
+ // offline->online (col)
+ if (module<80) { // inner layer
+ return colM%32;
+ }
+ else if (module<240) { // outer layer
+ return colM%32;
+ }
+ return 32; // error
+}
+//__________________________________________________________________________
+UInt_t AliITSRawStreamSPD::GetOnlineRowFromOffline(UInt_t module, UInt_t rowM) {
+ // offline->online (row)
+ if (module<80) { // inner layer
+ return (255-rowM);
+ }
+ else if (module<240) { // outer layer
+ return (255-rowM);
+ }
+ return 256; // error
+}
+//__________________________________________________________________________
+UInt_t AliITSRawStreamSPD::GetOfflineModuleFromOnline(UInt_t eqId, UInt_t hs, UInt_t chip) {
+ // online->offline (module)
+ if (eqId<20 && hs<6 && chip<10) return fgkDDLModuleMap[eqId][hs*2+chip/5];
+ else return 240;
+}
+//__________________________________________________________________________
+UInt_t AliITSRawStreamSPD::GetOfflineChipKeyFromOnline(UInt_t eqId, UInt_t hs, UInt_t chip) {
+ // online->offline (chip key: 0-1199)
+ if (eqId<20 && hs<6 && chip<10) {
+ UInt_t module = GetOfflineModuleFromOnline(eqId,hs,chip);
+ UInt_t chipInModule = ( chip>4 ? chip-5 : chip );
+ if(eqId>9) chipInModule = 4 - chipInModule; // side C only
+ return (module*5 + chipInModule);
+ } else return 1200;
+}
+//__________________________________________________________________________
+UInt_t AliITSRawStreamSPD::GetOfflineColFromOnline(UInt_t eqId, UInt_t hs, UInt_t chip, UInt_t col) {
+ // online->offline (col)
+ if (eqId>=20 || hs>=6 || chip>=10 || col>=32) return 160; // error
+ UInt_t offset = 32 * (chip % 5);
+ if (hs<2) {
+ if (eqId<10) {
+ return 159 - (31-col + offset); // inner layer, side A
+ }
+ else {
+ return col + offset; // inner layer, side C
+ }
+ }
+ else {
+ if (eqId<10) {
+ return (col + offset); // outer layer, side A
+ }
+ else {
+ return 159 - (31-col + offset); // outer layer, side C
+ }
+ }
+}
+//__________________________________________________________________________
+UInt_t AliITSRawStreamSPD::GetOfflineRowFromOnline(UInt_t eqId, UInt_t hs, UInt_t chip, UInt_t row) {
+ // online->offline (row)
+ if (eqId>=20 || hs>=6 || chip>=10 || row>=256) return 256; // error
+ return 255-row;
+}
+