X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSRawStreamSPD.cxx;h=6c0a72d06e8b0009b9dfb7152976be902dc250c7;hb=70f90d01a5c015f7a32fdcb24bfa048ff1658fef;hp=bcf52869d0902cd6a595626dfa6d50b71780fa2a;hpb=1fd93b67b9abf6d9c0c4339333a65f6d38a85d41;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSRawStreamSPD.cxx b/ITS/AliITSRawStreamSPD.cxx index bcf52869d09..6c0a72d06e8 100644 --- a/ITS/AliITSRawStreamSPD.cxx +++ b/ITS/AliITSRawStreamSPD.cxx @@ -27,50 +27,83 @@ ClassImp(AliITSRawStreamSPD) - - // this map has to change, waiting for the new geometry 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), - fEventNumber(-1),fChipAddr(0),fHalfStaveNr(0),fCol(0),fRow(0), + 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) + 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 + fRawReader->Reset(); fRawReader->Select("ITSSPD"); // reset calib header words for (UInt_t iword=0; iwordReadNextChar(fDataChar1); @@ -83,9 +116,12 @@ Bool_t AliITSRawStreamSPD::ReadNextShort() return kFALSE; } fData = fDataChar3+(fDataChar4<<8); - if ((*fRawReader->GetEventId())!=fPrevEventId) { // if new event... - NewEvent(); - fPrevEventId=(*fRawReader->GetEventId()); + const UInt_t *intPtr = fRawReader->GetEventId(); + if (intPtr!=0) { + if (*intPtr!=fPrevEventId) { // if new event... + NewEvent(); + fPrevEventId=*intPtr; + } } fFirstWord=kFALSE; } @@ -93,12 +129,12 @@ Bool_t AliITSRawStreamSPD::ReadNextShort() fFirstWord=kTRUE; fData = fDataChar1+(fDataChar2<<8); } - + fEqPLBytesRead+=2; + return kTRUE; } - -Bool_t AliITSRawStreamSPD::ReadNextInt() -{ +//__________________________________________________________________________ +Bool_t AliITSRawStreamSPD::ReadNextInt() { // reads next 32 bit into fDataChar1..4 // (if first 16 bits read already, just completes the present word) if (fFirstWord) { @@ -113,150 +149,379 @@ Bool_t AliITSRawStreamSPD::ReadNextInt() } return kFALSE; } - -void AliITSRawStreamSPD::NewEvent() -{ +//__________________________________________________________________________ +void AliITSRawStreamSPD::NewEvent() { // call this to reset flags for a new event for (UInt_t eqId=0; eqId<20; eqId++) { fCalHeadRead[eqId]=kFALSE; } - fEventNumber=-1; + 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; + } + } + } } - -Bool_t AliITSRawStreamSPD::ReadCalibHeader() -{ +//__________________________________________________________________________ +Int_t AliITSRawStreamSPD::ReadCalibHeader() { + // needs to be called in the beginning of each equipment data // read the extra calibration header - // returns kTRUE if the header is present and has length > 0 + // returns the length of the header if it is present, -1 otherwise Int_t ddlID = fRawReader->GetDDLID(); if (ddlID==-1) { // we may need to read one word to get the blockAttr - if (!ReadNextShort()) return kFALSE; + 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()); } - UChar_t attr = fRawReader->GetBlockAttributes(); + 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 ((attr & 0x40) == 0x40) { // is the header present? + 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); - if (calLen>kCalHeadLenMax) { - fRawReader->AddMajorErrorLog(kCalHeaderLengthErr,Form("Header length %d > max = %d",calLen,kCalHeadLenMax)); - AliWarning(Form("Header length problem. %d > %d (max)",calLen,kCalHeadLenMax)); - return kFALSE; + 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 if (calLen>0) { - for (UInt_t iword=0; iwordAddMajorErrorLog(kCalHeaderLengthErr,"header length problem"); - AliWarning("header length problem"); - return kFALSE; + 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 kTRUE; + return fCalHeadLen; } } } - return kFALSE; + return -1; } - -Bool_t AliITSRawStreamSPD::Next() -{ +//__________________________________________________________________________ +Bool_t AliITSRawStreamSPD::Next() { // read the next raw digit // returns kFALSE if there is no digit left - Int_t ddlID=-1; fPrevModuleID = fModuleID; while (ReadNextShort()) { - ddlID = fRawReader->GetDDLID(); - if (ddlID>=0 && ddlID<20) { - if (!fCalHeadRead[ddlID]) { - ReadCalibHeader(); + 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 { - fRawReader->AddMajorErrorLog(kDDLNumberErr,Form("Wrong DDL number %d",ddlID)); - AliWarning(Form("DDL number error (= %d) , setting it to 19", ddlID)); - ddlID=19; + 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) { - fRawReader->AddMajorErrorLog(kEventNumberErr,Form("Reading event number %d instead of %d",eventNumber,fEventNumber)); - AliWarning(Form("mismatching event numbers: %d != %d",eventNumber, fEventNumber)); + 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) { - fRawReader->AddMajorErrorLog(kChipAddrErr,Form("Overflow chip address %d - set to 9",fChipAddr)); - AliWarning(Form("overflow chip addr (= %d) , setting it to 9", fChipAddr)); - 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)) { - fRawReader->AddMajorErrorLog(kStaveNumberErr,Form("Half stave number error %d - set to 5",fHalfStaveNr)); - AliWarning(Form("half stave number error(=%d) , setting it to 5", fHalfStaveNr)); - fHalfStaveNr=5; + 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 = fgkDDLModuleMap[ddlID][fHalfStaveNr*2+fChipAddr/5]; - fOffset = 32 * (fChipAddr % 5); - } + fModuleID = GetOfflineModuleFromOnline(fDDLID,fHalfStaveNr,fChipAddr); + } + else if ((fData & 0xC000) == 0x0000) { // trailer - UShort_t hitCount = fData & 0x1FFF; + 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){ - fRawReader->AddMajorErrorLog(kNumbHitsErr,Form("Number of hits %d instead of %d",hitCount,fHitCount)); - AliWarning(Form("wrong number of hits: %d != %d", fHitCount, hitCount)); + 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++; fCol = (fData & 0x001F); fRow = (fData >> 5) & 0x00FF; - // translate ("online") chipcol, chiprow to ("offline") col (coord1), row (coord2): - // This will change, waiting for new geometry!!! - fCoord1 = fCol; - // if (fModuleID < 80 && ddlID < 10) fCoord1=31-fCoord1; - // else if (fModuleID >=80 && ddlID >=10) fCoord1=31-fCoord1; - fCoord1 += fOffset; - // if (ddlID>=10) fCoord1=159-fCoord1; - fCoord2 = fRow; - // if (fModuleID<80) fCoord2=255-fCoord2; + fCoord1 = GetOfflineColFromOnline(fDDLID,fHalfStaveNr,fChipAddr,fCol); + fCoord2 = GetOfflineRowFromOnline(fDDLID,fHalfStaveNr,fChipAddr,fRow); return kTRUE; - } + } + else { // fill word if ((fData & 0xC000) != 0xC000) { - fRawReader->AddMajorErrorLog(kWrongWordErr,"Wrong fill word"); - AliWarning("wrong fill word!"); + 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) { - fRawReader->AddMinorErrorLog(kHalfStaveStatusErr,"DDL ID = -1. Cannot read block attributes."); AliWarning("DDL ID = -1. Cannot read block attributes. Return kFALSE."); return kFALSE; } else { if (hs>=6) { - fRawReader->AddMinorErrorLog(kHalfStaveStatusErr,Form( "HS >= 6 requested (%d). Return kFALSE.",hs)); AliWarning(Form("HS >= 6 requested (%d). Return kFALSE.",hs)); return kFALSE; } @@ -269,7 +534,88 @@ Bool_t AliITSRawStreamSPD::GetHalfStavePresent(UInt_t hs) { } } } - +//__________________________________________________________________________ +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>(6+hs)) & (0x00000001)); else return kFALSE; @@ -294,3 +640,172 @@ 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; +} +