X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSRawStreamSPD.cxx;h=6c0a72d06e8b0009b9dfb7152976be902dc250c7;hb=05da1b4e2007170ee10c3a94cfef582320a3d258;hp=7f3149046e4a9f2dd5079af0e12454ec82e13a93;hpb=72600597e34e128f6b1e58cc3ebd220aeec02a5b;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSRawStreamSPD.cxx b/ITS/AliITSRawStreamSPD.cxx index 7f3149046e4..6c0a72d06e8 100644 --- a/ITS/AliITSRawStreamSPD.cxx +++ b/ITS/AliITSRawStreamSPD.cxx @@ -1,5 +1,5 @@ /************************************************************************** - * 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. * @@ -23,56 +23,89 @@ #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; iwordReadNextChar(fDataChar1); if (!b1) return kFALSE; Bool_t b2, b3, b4; @@ -83,91 +116,696 @@ Bool_t AliITSRawStreamSPD::ReadNextShort() 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; iwordReadNextInt(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; iwordAddMajorErrorLog(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>(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; +} +