]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - MUON/AliMUONPayloadTracker.cxx
add a flag to reconstruct from clusters instead of digits
[u/mrichter/AliRoot.git] / MUON / AliMUONPayloadTracker.cxx
index 23dca5b67d3e8fef13d3ee5cc15f10baedf47fcc..4733c66c994428929960cec0021bc1c22ac87828 100644 (file)
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
+// $Id$
 
-///////////////////////////////////////////////////////////////////////////////
-///
-/// Class Payload
-///
+//-----------------------------------------------------------------------------
+/// \class AliMUONPayloadTracker
 /// Decodes rawdata from buffer and stores in TClonesArray.
-/// 
 /// First version implement for Tracker
 ///
-///////////////////////////////////////////////////////////////////////////////
+/// \author Christian Finck
+//-----------------------------------------------------------------------------
 
-#include "AliMUONPayloadTracker.h"
+#include <TObjString.h>
 
-#include "AliRawReader.h"
-#include "AliRawDataHeader.h"
-#include "AliLog.h"
+#include "AliMUONPayloadTracker.h"
 
 #include "AliMUONDspHeader.h"
 #include "AliMUONBlockHeader.h"
 #include "AliMUONBusStruct.h"
 #include "AliMUONDDLTracker.h"
 
+#include "AliMUONLogger.h"
+#include "AliLog.h"
+
+/// \cond CLASSIMP
 ClassImp(AliMUONPayloadTracker)
+/// \endcond
 
 AliMUONPayloadTracker::AliMUONPayloadTracker()
   : TObject(),
+    fBusPatchId(0),
+    fDspId(0),
+    fBlkId(0),
+    fMaxDDL(20),
     fMaxBlock(2),
     fMaxDsp(5),
-    fMaxBus(5)
+    fMaxBus(5),
+    fDDLTracker(new AliMUONDDLTracker()),
+    fBusStruct(new AliMUONBusStruct()),
+    fBlockHeader(new AliMUONBlockHeader()),
+    fDspHeader(new AliMUONDspHeader()),
+    fLog(new AliMUONLogger(1000)),
+    fParityErrors(0),
+    fGlitchErrors(0),
+    fPaddingErrors(0),
+    fWarnings(kTRUE)
 {
-  //
-  // create an object to read MUON raw digits
-  // Default ctor for monitoring purposes
-  //
-  fBusPatchManager = new AliMpBusPatch();
-  fBusPatchManager->ReadBusPatchFile();
-
-  fDDLTracker      = new AliMUONDDLTracker();
-  fBusStruct       = new AliMUONBusStruct();
-  fBlockHeader     = new AliMUONBlockHeader();
-  fDspHeader       = new AliMUONDspHeader();
-}
+  ///
+  /// create an object to decode MUON payload
+  ///
 
-//_________________________________________________________________
-AliMUONPayloadTracker::AliMUONPayloadTracker(const AliMUONPayloadTracker& stream) :
-  TObject(stream)
-{ 
-  //
-  // copy ctor
-  //
-  AliFatal("copy constructor not implemented");
 }
 
-//______________________________________________________________________
-AliMUONPayloadTracker& AliMUONPayloadTracker::operator = (const AliMUONPayloadTracker& 
-                                             /* stream */)
-{
-  // 
-  // assignment operator
-  //
-  AliFatal("assignment operator not implemented");
-  return *this;
-}
-
-
 //___________________________________
 AliMUONPayloadTracker::~AliMUONPayloadTracker()
 {
-  //
-  // clean up
-  //
-  delete fBusPatchManager;
+  ///
+  /// clean up
+  ///
   delete fDDLTracker;
   delete fBusStruct;
   delete fBlockHeader;
   delete fDspHeader;
+  delete fLog;
 }
 
 //______________________________________________________
-Bool_t AliMUONPayloadTracker::Decode(UInt_t* buffer, Int_t ddl)
+Bool_t AliMUONPayloadTracker::Decode(UInt_t* buffer, Int_t totalDDLSize)
 {
-  // reading tracker DDL
-  // store buspatch info into Array
-  // store only non-empty structures (buspatch info with datalength !=0)
+
+  /// Each DDL is made with 2 Blocks each of which consists of 5 DSP's at most 
+  /// and each of DSP has at most 5 buspatches.
+  /// The different structures, Block (CRT), DSP (FRT) and Buspatch,
+  /// are identified by a key word 0xFC0000FC, 0xF000000F and 0xB000000B respectively.
+  /// (fBusPatchManager no more needed !)
 
 
   //Read Header Size of DDL,Block,DSP and BusPatch
@@ -104,143 +93,231 @@ Bool_t AliMUONPayloadTracker::Decode(UInt_t* buffer, Int_t ddl)
   static Int_t kDspHeaderSize      = fDspHeader->GetHeaderLength();
   static Int_t kBusPatchHeaderSize = fBusStruct->GetHeaderLength();
 
-  //  Int_t totalDDLSize;
+  // size structures
   Int_t totalBlockSize;
   Int_t totalDspSize;
   Int_t totalBusPatchSize;
   Int_t dataSize; 
-
-
-  Int_t iBusPerDSP[5]; // number of bus patches per DSP
-  Int_t iDspMax;       // number max of DSP per block
-
-  // minimum data size (only header's)
-  //  Int_t blankDDLSize;
-  Int_t blankBlockSize;
-  Int_t blankDspSize; 
-
-  AliDebug(3, Form("DDL Number %d\n", ddl ));
-
-  // getting DSP info
-  fBusPatchManager->GetDspInfo(ddl/2, iDspMax, iBusPerDSP);
-
-  // Each DDL is made with 2 Blocks each of which consists of 5 DSP's at most 
-  // and each of DSP has at most 5 buspatches.
-  // This information is used to calculate the size of headers (Block and DSP) 
-  // which has empty data.
-
-  //  blankDDLSize   = 2*kBlockHeaderSize + 2*iDspMax*kDspHeaderSize;
-  blankBlockSize = kBlockHeaderSize + iDspMax*kDspHeaderSize;
-  // totalDDLSize   = sizeof(buffer)/4;
-
-  for (Int_t i = 0; i < iDspMax; i++) {
-    //  blankDDLSize   += 2*iBusPerDSP[i]*kBusPatchHeaderSize;
-    blankBlockSize +=   iBusPerDSP[i]*kBusPatchHeaderSize;
+  Int_t bufSize;
+
+  // indexes
+  Int_t indexBlk;
+  Int_t indexDsp;
+  Int_t indexBusPatch;
+  Int_t index = 0;
+  Int_t iBlock = 0;
+
+  // CROCUS CRT
+  while (buffer[index] == fBlockHeader->GetDefaultDataKey()) {
+    
+    if (iBlock > fMaxBlock) break;
+    
+    // copy within padding words
+    memcpy(fBlockHeader->GetHeader(),&buffer[index], (kBlockHeaderSize)*4);
+    
+    totalBlockSize = fBlockHeader->GetTotalLength();
+    
+    indexBlk = index;
+    index += kBlockHeaderSize;
+    
+    // copy in TClonesArray
+    fDDLTracker->AddBlkHeader(*fBlockHeader);
+    
+    // Crocus FRT
+    Int_t iDsp = 0;
+    while (buffer[index] == fDspHeader->GetDefaultDataKey()) {
+      
+      if (iDsp > fMaxDsp) break; // if ever...
+      
+      memcpy(fDspHeader->GetHeader(),&buffer[index], kDspHeaderSize*4);
+      
+      totalDspSize = fDspHeader->GetTotalLength();
+      
+      if (fDspHeader->GetErrorWord()) {
+        //     fDspHeader->Print("");
+        if ( fDspHeader->GetErrorWord() == (0x000000B1 |  fBlockHeader->GetDspId())
+            ||  fDspHeader->GetErrorWord() == (0x00000091 |  fBlockHeader->GetDspId()) ){
+          // an event with a glitch in the readout  has been detected
+          // it means that somewhere a 1 byte word has been randomly inserted
+          // all the readout sequence is shifted  untill the next event 
+          
+          Char_t* msg = Form("Glitch error detected in DSP %d, skipping event ", fBlockHeader->GetDspId());
+          
+          if (fWarnings) AliWarning(msg); 
+          AddErrorMessage(msg);
+          fGlitchErrors++;
+          return kFALSE ; 
+        }      
+      }
+      
+      indexDsp = index;
+      index += kDspHeaderSize;
+      
+      // copy in TClonesArray
+      fDDLTracker->AddDspHeader(*fDspHeader, iBlock);
+      
+      // buspatch structure
+      Int_t iBusPatch = 0;
+      while (buffer[index] == fBusStruct->GetDefaultDataKey()) {
+        
+        if (iBusPatch > fMaxBus) break; // if ever
+        
+        //copy buffer into header structure
+        memcpy(fBusStruct->GetHeader(), &buffer[index], kBusPatchHeaderSize*4);
+        
+        totalBusPatchSize = fBusStruct->GetTotalLength();
+        indexBusPatch     = index;
+        
+        //Check Buspatch header, not empty events
+        if(totalBusPatchSize > kBusPatchHeaderSize) {    
+          
+          index   += kBusPatchHeaderSize;
+          dataSize = fBusStruct->GetLength();
+          bufSize  = fBusStruct->GetBufSize();
+          
+          if(dataSize > 0) { // check data present
+            if (dataSize > bufSize) // check buffer size
+              fBusStruct->SetAlloc(dataSize);
+            
+            //copy buffer into data structure
+            memcpy(fBusStruct->GetData(), &buffer[index], dataSize*4);
+            fBusStruct->SetBlockId(iBlock); // could be usefull in future applications ?
+            fBusStruct->SetDspId(iDsp);
+            
+            // check parity
+            if(!CheckDataParity()) {
+              fParityErrors++;
+              return kFALSE;
+            }
+            
+            // copy in TClonesArray
+            fDDLTracker->AddBusPatch(*fBusStruct, iBlock, iDsp);
+            
+          } // dataSize test
+          
+        } // testing buspatch
+        
+        index = indexBusPatch + totalBusPatchSize;
+        if (index >= totalDDLSize) {// check the end of DDL
+          index = totalDDLSize - 1; // point to the last element of buffer
+          break;
+        }
+        iBusPatch++;
+      }  // buspatch loop
+      
+      // skipping additionnal word if padding
+      if (fDspHeader->GetPaddingWord() == 1) {
+        if (buffer[index++] != fDspHeader->GetDefaultPaddingWord()) {
+          
+          Char_t *msg = Form("Padding word error for iBlock %d, iDsp %d, iBus %d\n", 
+                             iBlock, iDsp, iBusPatch);
+          
+          if (fWarnings) AliWarning(msg);
+          AddErrorMessage(msg);
+          fPaddingErrors++;
+        }
+      }
+      
+      index = indexDsp + totalDspSize;
+      if (index >= totalDDLSize) {
+        index = totalDDLSize - 1;
+        break;
+      }
+      iDsp++;
+    }  // dsp loop
+    
+    index = indexBlk + totalBlockSize;
+    if (index >= totalDDLSize) {
+      index = totalDDLSize - 1;
+      break;
+    }
+    
+    iBlock++;
+  }  // block loop
+  
+  if (buffer[index++] != fBlockHeader->GetDdlDataKey() || 
+      buffer[index++] != fBlockHeader->GetDdlDataKey()) {
+    
+    Char_t *msg = Form("Bad end of DDL data key\n");
+    
+    if (fWarnings) AliWarning(msg);
+    AddErrorMessage(msg);
   }
-
-  // Compare the DDL header with an empty DDL header size to read the file
-  //  if(totalDDLSize > blankDDLSize) {  //should not happen in real life    
-    // indexes
-    Int_t indexDsp;
-    Int_t indexBusPatch;
-    Int_t index = 0;
-
-    for(Int_t iBlock = 0; iBlock < 2 ;iBlock++){  // loop over 2 blocks
-
-      // copy within padding words
-      memcpy(fBlockHeader->GetHeader(),&buffer[index], (kBlockHeaderSize+1)*4);
-
-      totalBlockSize = fBlockHeader->GetTotalLength();
-
-      fDDLTracker->AddBlkHeader(*fBlockHeader);
-
-      if (fBlockHeader->GetPadding() == 0xDEAD) // skipping padding word
-       index++;
-
-      if(totalBlockSize > blankBlockSize) {        // compare block header
-       index += kBlockHeaderSize;
-
-       for(Int_t iDsp = 0; iDsp < iDspMax ;iDsp++){   //DSP loop
-
-         if (iDsp > fMaxDsp) break;
-         
-         memcpy(fDspHeader->GetHeader(),&buffer[index], kDspHeaderSize*4);
-
-         totalDspSize = fDspHeader->GetTotalLength();
-         indexDsp = index;
-
-         blankDspSize =  kDspHeaderSize + iBusPerDSP[iDsp]*kBusPatchHeaderSize; // no data just header
-
-         fDDLTracker->AddDspHeader(*fDspHeader, iBlock);
-
-         if(totalDspSize > blankDspSize) {       // Compare DSP Header
-           index += kDspHeaderSize;
-               
-           for(Int_t iBusPatch = 0; iBusPatch < iBusPerDSP[iDsp]; iBusPatch++) {  
-
-             if (iBusPatch > fMaxBus) break; 
-
-             //copy buffer into header structure
-             memcpy(fBusStruct->GetBusPatchHeader(), &buffer[index], kBusPatchHeaderSize*4);
-
-             totalBusPatchSize = fBusStruct->GetTotalLength();
-             indexBusPatch     = index;
-               
-             //Check Buspatch header, not empty events
-             if(totalBusPatchSize > kBusPatchHeaderSize) {    
-
-               index   += kBusPatchHeaderSize;
-               dataSize = fBusStruct->GetLength();
-               Int_t bufSize = fBusStruct->GetBufSize();
-
-               if(dataSize > 0) { // check data present
-                 if (dataSize > bufSize) fBusStruct->SetAlloc(dataSize);
-
-                 //copy buffer into data structure
-                 memcpy(fBusStruct->GetData(), &buffer[index], dataSize*4);
-                 fBusStruct->SetBlockId(iBlock); // could be usefull in future applications ?
-                 fBusStruct->SetDspId(iDsp);
-                 fDDLTracker->AddBusPatch(*fBusStruct, iBlock, iDsp);
-
-               } // dataSize test
-
-             } // testing buspatch
-
-             index = indexBusPatch + totalBusPatchSize;
-
-           }  //buspatch loop
-               
-         }  // dsp test
-
-         index = indexDsp + totalDspSize;
-             
-       }  // dsp loop
-
-      }   //block test
-
-      index = totalBlockSize;
-
-    }  //block loop
-
-    //  } //loop checking the header size of DDL
-
+  
+  
   return kTRUE;
 }
 
 //______________________________________________________
 void AliMUONPayloadTracker::ResetDDL()
 {
-  // reseting TClonesArray
-  // after each DDL
-  //
+  /// reseting TClonesArray
+  /// after each DDL
+  ///
   fDDLTracker->GetBlkHeaderArray()->Delete();
+  fGlitchErrors  = 0;
+  fPaddingErrors = 0;
+  fParityErrors  = 0;
+
 }
 
 //______________________________________________________
 void AliMUONPayloadTracker::SetMaxBlock(Int_t blk) 
 {
-  // set regional card number
+  /// set regional card number
   if (blk > 2) blk = 2;
   fMaxBlock = blk;
 }
+
+//______________________________________________________
+Bool_t AliMUONPayloadTracker::CheckDataParity()
+{
+  /// parity check
+  /// taken from MuTrkBusPatch.cxx (sotfware test for CROCUS)
+  /// A. Baldisseri
+
+  Int_t  parity;
+  UInt_t data;
+  
+  Int_t dataSize = fBusStruct->GetLength();
+  for (int idata = 0; idata < dataSize; idata++) {
+
+    data  = fBusStruct->GetData(idata);
+
+    // Compute the parity for each data word
+    parity  = data & 0x1;
+    for (Int_t i = 1; i <= 30; i++) 
+      parity ^= ((data >> i) & 0x1);
+
+    // Check
+    if (parity != fBusStruct->GetParity(idata)) {
+
+      Char_t* msg = Form("Parity error in word %d for manuId %d and channel %d in buspatch %d\n", 
+                          idata, fBusStruct->GetManuId(idata), fBusStruct->GetChannelId(idata),
+                          fBusStruct->GetBusPatchId());
+      
+      if (fWarnings) AliWarning(msg);
+      AddErrorMessage(msg);
+      fParityErrors++;
+      
+      return kFALSE;
+                    
+    }
+  }
+  return kTRUE;
+}
+
+//______________________________________________________
+void AliMUONPayloadTracker::AddErrorMessage(const Char_t* msg)
+{
+/// adding message to logger
+    TString tmp(msg);
+  
+    Int_t pos = tmp.First("\n");
+    tmp[pos] = 0;
+    
+    fLog->Log(tmp.Data());
+}
+
+