]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - RAW/AliAltroRawStreamV3.cxx
fix finding of pad neighbours; remove methods to write them in OCDB
[u/mrichter/AliRoot.git] / RAW / AliAltroRawStreamV3.cxx
index 73f883075fc2bce26da6ec037c07ba318e3821bc..3baa34f0da542de0238ed23db5733df5f765f421 100644 (file)
@@ -30,6 +30,7 @@
 #include "AliRawReader.h"
 #include "AliLog.h"
 #include "AliAltroRawStream.h"
+#include "AliRawEventHeaderBase.h"
 
 ClassImp(AliAltroRawStreamV3)
 
@@ -42,12 +43,14 @@ AliAltroRawStreamV3::AliAltroRawStreamV3(AliRawReader* rawReader) :
   fHWAddress(-1),
   fRawReader(rawReader),
   fData(NULL),
+  fChannelStartPos(-1),
   fPosition(-1),
   fCount(-1),
   fStartTimeBin(-1),
   fBunchLength(-1),
   fBadChannel(kFALSE),
   fPayloadSize(-1),
+  fChannelPayloadSize(-1),
   fBunchDataPointer(NULL),
   fBunchDataIndex(-1),
   fRCUTrailerData(NULL),
@@ -61,7 +64,8 @@ AliAltroRawStreamV3::AliAltroRawStreamV3(AliRawReader* rawReader) :
   fAltroCFG1(0),
   fAltroCFG2(0),
   fOldStream(NULL),
-  fCheckAltroPayload(kTRUE)
+  fCheckAltroPayload(kTRUE),
+  fFormatVersion(0)
 {
   // Constructor
   // Create an object to read Altro raw digits in
@@ -86,12 +90,14 @@ AliAltroRawStreamV3::AliAltroRawStreamV3(const AliAltroRawStreamV3& stream) :
   fHWAddress(stream.fHWAddress),
   fRawReader(stream.fRawReader),
   fData(stream.fData),
+  fChannelStartPos(stream.fChannelStartPos),
   fPosition(stream.fPosition),
   fCount(stream.fCount),
   fStartTimeBin(stream.fStartTimeBin),
   fBunchLength(stream.fBunchLength),
   fBadChannel(stream.fBadChannel),
   fPayloadSize(stream.fPayloadSize),
+  fChannelPayloadSize(stream.fChannelPayloadSize),
   fBunchDataPointer(stream.fBunchDataPointer),
   fBunchDataIndex(stream.fBunchDataIndex),
   fRCUTrailerData(stream.fRCUTrailerData),
@@ -105,7 +111,8 @@ AliAltroRawStreamV3::AliAltroRawStreamV3(const AliAltroRawStreamV3& stream) :
   fAltroCFG1(stream.fAltroCFG1),
   fAltroCFG2(stream.fAltroCFG2),
   fOldStream(NULL),
-  fCheckAltroPayload(stream.fCheckAltroPayload)
+  fCheckAltroPayload(stream.fCheckAltroPayload),
+  fFormatVersion(0)
 {
   // Copy constructor
   // Copy the bunch data array
@@ -128,12 +135,14 @@ AliAltroRawStreamV3& AliAltroRawStreamV3::operator = (const AliAltroRawStreamV3&
   fHWAddress         = stream.fHWAddress;
   fRawReader         = stream.fRawReader;
   fData              = stream.fData;
+  fChannelStartPos   = stream.fChannelStartPos;
   fPosition          = stream.fPosition;
   fCount             = stream.fCount;
   fStartTimeBin      = stream.fStartTimeBin;
   fBunchLength       = stream.fBunchLength;
   fBadChannel        = stream.fBadChannel;
   fPayloadSize       = stream.fPayloadSize;
+  fChannelPayloadSize= stream.fChannelPayloadSize;
   fBunchDataPointer  = stream.fBunchDataPointer;
   fBunchDataIndex    = stream.fBunchDataIndex;
   fRCUTrailerData    = stream.fRCUTrailerData;
@@ -146,6 +155,7 @@ AliAltroRawStreamV3& AliAltroRawStreamV3::operator = (const AliAltroRawStreamV3&
   fActiveFECsB       = stream.fActiveFECsB;
   fAltroCFG1         = stream.fAltroCFG1;
   fAltroCFG2         = stream.fAltroCFG2;
+  fFormatVersion     = stream.fFormatVersion;
 
   for(Int_t i = 0; i < kMaxNTimeBins; i++) fBunchData[i] = stream.fBunchData[i];
 
@@ -167,10 +177,12 @@ void AliAltroRawStreamV3::Reset()
 // Reset of the raw-reader as well
 
   fDDLNumber = fRCUId = fHWAddress = -1;
+  fChannelStartPos = -1;
   fPosition = fCount = -1;
   fBunchLength = fStartTimeBin = -1;
   fBadChannel = kFALSE;
   fPayloadSize = -1;
+  fChannelPayloadSize = -1;
   fBunchDataPointer = NULL;
   fBunchDataIndex = -1;
 
@@ -190,7 +202,7 @@ Bool_t AliAltroRawStreamV3::NextDDL()
 // Read the next DDL payload (CDH + RCU trailer)
 // Updates the information which is coming from these
 // two sources
-
+  fFormatVersion = 0;
   fPosition = 0;
   // Get next DDL payload
   // return wtih false in case no more data payloads
@@ -200,15 +212,17 @@ Bool_t AliAltroRawStreamV3::NextDDL()
   } while (fRawReader->GetDataSize() == 0);
 
   fDDLNumber = fRawReader->GetDDLID();
+  fChannelPayloadSize = -1;
+  fChannelStartPos = -1;
 
-  UChar_t rcuVer = fRawReader->GetBlockAttributes();
+  fFormatVersion = (fRawReader->GetBlockAttributes() & 0xF);
 
-  if (rcuVer < 2) {
+  if (fFormatVersion < 2) {
     // old altro format data
     if (!fOldStream) {
       fOldStream = new AliAltroRawStream(fRawReader);
       AliInfo(Form("RCU firmware verion %d detected. Using AliAltroRawStream to decode the data.",
-                  rcuVer));
+                  fFormatVersion));
     }
     Bool_t status = fOldStream->NextDDL(fData);
     if (status) {
@@ -224,11 +238,13 @@ Bool_t AliAltroRawStreamV3::NextDDL()
       fActiveFECsB = fOldStream->GetActiveFECsB();
       fAltroCFG1 = fOldStream->GetAltroCFG1();
       fAltroCFG2 = fOldStream->GetAltroCFG2();
+      if (fRawReader->GetType() == AliRawEventHeaderBase::kStartOfData)
+       fPayloadSize = fOldStream->GetRCUPayloadSizeInSOD();
     }
     return status;
   }
 
-  return ReadRCUTrailer(rcuVer);
+  return ReadRCUTrailer(fFormatVersion);
 }
 
 //_____________________________________________________________________________
@@ -241,10 +257,15 @@ Bool_t AliAltroRawStreamV3::NextChannel()
   // RCU signals readout error in this channel
   if (fOldStream) {
     Bool_t status = fOldStream->NextChannel();
-    if (status) fHWAddress = fOldStream->GetHWAddress();
+    if (status) {
+      fHWAddress = fOldStream->GetHWAddress();
+      fChannelPayloadSize = fOldStream->GetChannelPayloadSize();
+    }
     return status;
   }
 
+  Int_t channelStartPos=fPosition;
+  fChannelStartPos = -1;
   fCount = -1;
   fBadChannel = kFALSE;
   fBunchDataIndex = 0;
@@ -262,6 +283,7 @@ Bool_t AliAltroRawStreamV3::NextChannel()
 
   // extract channel payload and hw address
   fCount = (word >> 16) & 0x3FF; 
+  fChannelPayloadSize = fCount;
   fHWAddress = word & 0xFFF;
 
   // Now unpack the altro data
@@ -275,9 +297,10 @@ Bool_t AliAltroRawStreamV3::NextChannel()
     word = Get32bitWord(fPosition++);
     if ((word >> 30) != 0) {
       // Unexpected end of altro channel payload
-      AliWarning(Form("Unexpected end of payload in altro channel payload! Address=0x%x, word=0x%x",
-                     fHWAddress,word));
+      AliWarning(Form("Unexpected end of payload in altro channel payload! DDL=%03d, Address=0x%x, word=0x%x",
+                     fDDLNumber,fHWAddress,word));
       fRawReader->AddMinorErrorLog(kAltroPayloadErr,Form("hw=0x%x",fHWAddress));
+      if (AliDebugLevel() > 0) HexDumpChannel();
       fCount = -1;
       fPosition--;
       return kFALSE;
@@ -287,6 +310,7 @@ Bool_t AliAltroRawStreamV3::NextChannel()
     fBunchData[isample++] = word & 0x3FF;
   }  
 
+  fChannelStartPos=channelStartPos;
   return kTRUE;
 }
 
@@ -313,17 +337,19 @@ Bool_t AliAltroRawStreamV3::NextBunch()
   fBunchLength = fBunchData[fBunchDataIndex];
   if (fBunchLength <= 2) {
     // Invalid bunch size
-    AliWarning(Form("Too short bunch length (%d) in Address=0x%x !",
-                   fBunchLength,fHWAddress));
+     AliDebug(1,Form("Too short bunch length (%d) @ %d in Address=0x%x (DDL=%03d)!",
+                   fBunchLength,fBunchDataIndex,fHWAddress,fDDLNumber));
     fRawReader->AddMinorErrorLog(kAltroBunchHeadErr,Form("hw=0x%x",fHWAddress));
+    if (AliDebugLevel() > 0) HexDumpChannel();
     fCount = fBunchLength = -1;
     return kFALSE;
   }
   if ((fBunchDataIndex + fBunchLength) > fCount) {
     // Too long bunch detected
-    AliWarning(Form("Too long bunch detected in Address=0x%x ! Expected <= %d 10-bit words, found %d !",
-                   fHWAddress,fCount-fBunchDataIndex,fBunchLength));
+    AliDebug(1,Form("Too long bunch detected @ %d in Address=0x%x (DDL=%03d) ! Expected <= %d 10-bit words, found %d !", fBunchDataIndex,
+                   fHWAddress,fDDLNumber,fCount-fBunchDataIndex,fBunchLength));
     fRawReader->AddMinorErrorLog(kAltroBunchHeadErr,Form("hw=0x%x",fHWAddress));
+    if (AliDebugLevel() > 0) HexDumpChannel();
     fCount = fBunchLength = -1;
     return kFALSE;
   }
@@ -333,16 +359,18 @@ Bool_t AliAltroRawStreamV3::NextBunch()
   fStartTimeBin = fBunchData[fBunchDataIndex++];
   if (fCheckAltroPayload) {
     if ((fStartTimeBin-fBunchLength+1) < 0) {
-      AliWarning(Form("Invalid start time-bin in Address=0x%x ! (%d-%d+1) < 0",
-                     fHWAddress,fStartTimeBin,fBunchLength));
+      AliWarning(Form("Invalid start time-bin @ %d in Address=0x%x (DDL=%03d)! (%d-%d+1) < 0", fBunchDataIndex-1,
+                     fHWAddress,fDDLNumber,fStartTimeBin,fBunchLength));
       fRawReader->AddMinorErrorLog(kAltroPayloadErr,Form("hw=0x%x",fHWAddress));
+      if (AliDebugLevel() > 0) HexDumpChannel();
       fCount = fBunchLength = -1;
       return kFALSE;
     }
     if (fStartTimeBin >= prevTimeBin) {
-      AliWarning(Form("Invalid start time-bin in Address=0x%x ! (%d>=%d)",
-                     fHWAddress,fStartTimeBin,prevTimeBin));
+      AliWarning(Form("Invalid start time-bin @ %d in Address=0x%x (DDL=%03d)! (%d>=%d)", fBunchDataIndex-1,
+                     fHWAddress,fDDLNumber,fStartTimeBin,prevTimeBin));
       fRawReader->AddMinorErrorLog(kAltroPayloadErr,Form("hw=0x%x",fHWAddress));
+      if (AliDebugLevel() > 0) HexDumpChannel();
       fCount = fBunchLength = -1;
       return kFALSE;
     }
@@ -355,6 +383,16 @@ Bool_t AliAltroRawStreamV3::NextBunch()
   return kTRUE;
 }
 
+//_____________________________________________________________________________
+const UChar_t *AliAltroRawStreamV3::GetChannelPayload() const
+{
+  //returns raw channel data, length 4+(fChannelPayloadSize+2)/3*4
+  if (fChannelStartPos<0 || fChannelPayloadSize<0) return NULL;
+  Int_t channelSize=1+(fChannelPayloadSize+2)/3; // nof 32bit words
+  if (fPosition<fChannelStartPos+channelSize) return NULL;
+  return fData+(fChannelStartPos<<2);
+}
+
 //_____________________________________________________________________________
 UInt_t AliAltroRawStreamV3::Get32bitWord(Int_t index) const
 {
@@ -574,25 +612,31 @@ void AliAltroRawStreamV3::PrintRCUTrailer() const
 {
   // Prints the contents of
   // the RCU trailer data
-  printf("RCU trailer:\n===========\n");
-  printf("FECERRA: 0x%x\nFECERRB: 0x%x\n",fFECERRA,fFECERRB);
-  printf("ERRREG2: 0x%x\n",fERRREG2);
+  printf("RCU trailer (Format version %d):\n"
+        "==================================================\n",  GetFormatVersion());
+  printf("FECERRA:                                   0x%x\n", fFECERRA);
+  printf("FECERRB:                                   0x%x\n", fFECERRB);
+  printf("ERRREG2:                                   0x%x\n", fERRREG2);
   printf("#channels skipped due to address mismatch: %d\n",GetNChAddrMismatch());
   printf("#channels skipped due to bad block length: %d\n",GetNChLengthMismatch());
-  printf("Active FECs (branch A): 0x%x\nActive FECs (branch B): 0x%x\n",fActiveFECsA,fActiveFECsB);
-  printf("Baseline corr: 0x%x\n",GetBaselineCorr());
-  printf("Number of presamples: %d\nNumber of postsamples: %d\n",GetNPresamples(),GetNPostsamples());
-  printf("Second baseline corr: %d\n",GetSecondBaselineCorr());
-  printf("GlitchFilter: %d\n",GetGlitchFilter());
-  printf("Number of non-ZS postsamples: %d\nNumber of non-ZS presamples: %d\n",GetNNonZSPostsamples(),GetNNonZSPresamples());
-  printf("Number of ALTRO buffers: %d\n",GetNAltroBuffers());
-  printf("Number of pretrigger samples: %d\n",GetNPretriggerSamples());
-  printf("Number of samples per channel: %d\n",GetNSamplesPerCh());
-  printf("Sparse readout: %d\n",GetSparseRO());
-  printf("Sampling time: %e s\n",GetTSample());
-  printf("L1 Phase: %e s\n",GetL1Phase());
-  printf("AltroCFG1: 0x%x\nAltroCFG2: 0x%x\n",GetAltroCFG1(),GetAltroCFG2());
-  printf("===========\n");
+  printf("Active FECs (branch A):                    0x%x\n", fActiveFECsA);
+  printf("Active FECs (branch B):                    0x%x\n", fActiveFECsB);
+  printf("Baseline corr:                             0x%x\n",GetBaselineCorr());
+  printf("Number of presamples:                      %d\n", GetNPresamples());
+  printf("Number of postsamples:                     %d\n",GetNPostsamples());
+  printf("Second baseline corr:                      %d\n",GetSecondBaselineCorr());
+  printf("GlitchFilter:                              %d\n",GetGlitchFilter());
+  printf("Number of non-ZS postsamples:              %d\n",GetNNonZSPostsamples());
+  printf("Number of non-ZS presamples:               %d\n",GetNNonZSPresamples());
+  printf("Number of ALTRO buffers:                   %d\n",GetNAltroBuffers());
+  printf("Number of pretrigger samples:              %d\n",GetNPretriggerSamples());
+  printf("Number of samples per channel:             %d\n",GetNSamplesPerCh());
+  printf("Sparse readout:                            %d\n",GetSparseRO());
+  printf("Sampling time:                             %e s\n",GetTSample());
+  printf("L1 Phase:                                  %e s\n",GetL1Phase());
+  printf("AltroCFG1:                                 0x%x\n",GetAltroCFG1());
+  printf("AltroCFG2:                                 0x%x\n",GetAltroCFG2());
+  printf("==================================================\n");
 }
 
 //_____________________________________________________________________________
@@ -660,7 +704,7 @@ Double_t AliAltroRawStreamV3::GetL1Phase() const
 
   Double_t tSample = GetTSample();
   if (phase >= tSample) {
-    AliWarning(Form("Invalid L1 trigger phase (%f >= %d) !",
+    AliWarning(Form("Invalid L1 trigger phase (%f >= %f) !",
                    phase,tSample));
     phase = 0.;
   }
@@ -677,3 +721,61 @@ void AliAltroRawStreamV3::AddMappingErrorLog(const char *message)
 
   if (fRawReader) fRawReader->AddMinorErrorLog(kBadAltroMapping,message);
 }
+
+//_____________________________________________________________________________
+UChar_t *AliAltroRawStreamV3::GetRCUPayloadInSOD() const
+{
+  // Get a pointer to the data in case
+  // of SOD events
+  if (fRawReader) {
+    if (fRawReader->GetType() == AliRawEventHeaderBase::kStartOfData) {
+      return fData;
+    }
+  }
+  return NULL;
+}
+
+//_____________________________________________________________________________
+Int_t AliAltroRawStreamV3::GetRCUPayloadSizeInSOD() const
+{
+  // Get the size of the RCU data in case
+  // of SOD events
+  if (fRawReader) {
+    if (fRawReader->GetType() == AliRawEventHeaderBase::kStartOfData) {
+      return fPayloadSize;
+    }
+  }
+  return -1;
+}
+
+//_____________________________________________________________________________
+
+void AliAltroRawStreamV3::HexDumpChannel() const
+{
+  // Print of the Hex Data of the current channel
+  // to decipher read-out warnings and errors
+  if (fCount>0 && fPosition>0) {
+    printf("Hex-Dump of DDL: %3d, RCU ID: %d, HWADDR: 0x%03x\n",
+           fDDLNumber,fRCUId,fHWAddress);
+    printf("32-bit     - 2bit 10bit 10bit 10bit - ");
+    printf("Index 10bit  10bit  10bit\n");
+    printf("**********   **** ***** ***** ***** - ");
+    printf("***** ****** ****** ******\n");
+    Int_t nwords = (fCount+2)/3+1;
+    for (Int_t iword = 0; iword < nwords; iword++) {
+      UInt_t word32 = Get32bitWord(fPosition-nwords+iword);
+      UInt_t marker  = word32 >> 30 & 0x3;
+      UInt_t word101 = word32 >> 20 & 0x3FF;
+      UInt_t word102 = word32 >> 10 & 0x3FF;
+      UInt_t word103 = word32 >> 00 & 0x3FF; // nice octal number
+      printf("0x%08x - 0b%1d%1d 0x%03x 0x%03x 0x%03x",
+            word32,marker>>1,marker&0x1,word101,word102,word103);
+      if (iword == 0) { printf(" - Channel Header\n"); continue; }
+      Int_t base = 3*(iword-1);
+      printf(" - %5d 0x%03x%c 0x%03x%c 0x%03x%c\n", base,
+            fBunchData[base],   (fBunchDataIndex == base   ? '*' : ' '),
+            fBunchData[base+1], (fBunchDataIndex == base+1 ? '*' : ' '),
+            fBunchData[base+2], (fBunchDataIndex == base+2 ? '*' : ' '));
+    }
+  }
+}