]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDrawStream.h
added slewing correction by data
[u/mrichter/AliRoot.git] / TRD / AliTRDrawStream.h
index bd9d71d52ff55bdb866549ecd02808f41d7fe3ed..98d70d88afa4eef4e0e0e9fd476e6a7599e4a6b5 100644 (file)
 #ifndef ALITRDRAWSTREAM_H
 #define ALITRDRAWSTREAM_H
 
-#include "TObject.h"
-#include "TClonesArray.h"
 #include "TTree.h"
-
 #include "AliLog.h"
+#include "TClonesArray.h"
 
+class TObject;
 class TObjArray;
 class TString;
 class TBranch;
@@ -27,6 +26,7 @@ class AliTRDdigitsParam;
 class AliTRDarrayADC;
 class AliTRDSignalIndex;
 class AliTRDtrackletContainer;
+class AliESDTrdTrack;
 
 class AliTRDrawStream : public TObject
 {
@@ -48,16 +48,13 @@ class AliTRDrawStream : public TObject
   TClonesArray* GetMarkerArray() const { return fMarkers; }
 
   AliTRDdigitsManager* GetDigitsManager() const { return fDigitsManager; }
-  TTree *GetTrackletTree() const { return fTrackletTree; }
 
-  Bool_t ReadEvent(TTree *trackletTree = 0x0);
+  Bool_t ReadEvent();
 
   Bool_t NextDDL();
   Int_t NextChamber(AliTRDdigitsManager *digMgr);
   Int_t NextChamber(AliTRDdigitsManager *digMgr,
-                     UInt_t ** /* trackletContainer */, UShort_t ** /* errorContainer */) { AliError("Deprecated, use NextChamber(AliTRDdigitsManger*) instead!"); return NextChamber(digMgr); }
-
-  Bool_t ConnectTracklets(TTree *trklTree);
+                     UInt_t ** /* trackletContainer */, UShort_t ** /* errorContainer */) { AliError("Deprecated, use NextChamber(AliTRDdigitsManger*) instead!"); return NextChamber(digMgr); }
 
   void StoreErrorsInTree()   { fStoreError = &AliTRDrawStream::StoreErrorTree; }
   void StoreErrorsInArray()  { fStoreError = &AliTRDrawStream::StoreErrorArray; }
@@ -68,10 +65,11 @@ class AliTRDrawStream : public TObject
   enum ErrorCode_t {
     kUnknown = 0,
     kLinkMonitor,
-    kPtrgCntMismatch,
+    kEvCntMismatch,
     kNonTrdEq,
     kStackHeaderInvalid,
     kInvalidDetector,
+    kInvalidPadRow,
     kNoDigits,
     kHCmismatch,
     kHCcheckFailed,
@@ -84,14 +82,17 @@ class AliTRDrawStream : public TObject
     kAdcDataAbort,
     kAdcChannelsMiss,
     kMissMcmHeaders,
+    kMissTpData,
+    kCRCmismatch,
     kLastErrorCode
   };
 
   enum ErrorBehav_t {
     kTolerate = 0,
-    kAbort = 1,
-    kDiscardMCM = 2,
-    kDiscardHC = 4
+    kDiscardMCM = 1,
+    kDiscardHC = 2,
+    kDiscardDDL = 4,
+    kAbort = 8
   };
 
   enum MarkerCode_t {
@@ -155,24 +156,59 @@ class AliTRDrawStream : public TObject
     ClassDef(AliTRDrawStats, 1);
   };
 
-  AliTRDrawStats fStats;            // event statistics, clearing must be done by the user
-
   AliTRDrawStats* GetStats() { return &fStats; }
   Int_t GetEventSize(Int_t sector)  const { return fStats.fStatsSector[sector].fBytes; }
+  Int_t GetEventSize(Int_t sector, Int_t stack)  const {
+    Int_t size = 0;
+    for (Int_t iHC = 0; iHC < 12; iHC++) {
+      size += fStats.fStatsSector[sector].fStatsHC[12*stack + iHC].fBytes;
+    }
+    return size;
+  }
+  Int_t GetEventSize(Int_t sector, Int_t stack, Int_t hc) const {
+    return fStats.fStatsSector[sector].fStatsHC[12*stack + hc].fBytes;
+  }
   Int_t GetNTracklets(Int_t sector) const { return fStats.fStatsSector[sector].fNTracklets; }
   Int_t GetNMCMs(Int_t sector)      const { return fStats.fStatsSector[sector].fNMCMs; }
   Int_t GetNChannels(Int_t sector)  const { return fStats.fStatsSector[sector].fNChannels; }
 
+  ULong64_t GetTrkFlags(Int_t sector, Int_t stack) const { return (fCurrTrgFlags[sector] & (1 << (27 + stack))) ? fCurrTrkFlags[sector*fgkNstacks + stack] : 0; }
+  UInt_t GetTriggerFlags(Int_t sector) const { return fCurrTrgFlags[sector]; }
+  UInt_t GetLinkMonitorFlags(Int_t sector, Int_t stack) const {
+    UInt_t temp = 0;
+    for (Int_t iLink = 0; iLink < 12; iLink++)
+      temp |= (fCurrLinkMonitorFlags[(sector * fgkNstacks + stack) * fgkNlinks + iLink] & 0x3) << (iLink * 2);
+    temp |= ((GetMatchFlagsBP(sector)   >> stack) & 0x1) << 24;
+    temp |= ((GetMatchFlagsSRAM(sector) >> stack) & 0x1) << 25;
+    return temp;
+  }
+  UInt_t GetMatchFlagsSRAM(Int_t sector) const { return fCurrMatchFlagsSRAM[sector]; }
+  UInt_t GetMatchFlagsBP(Int_t sector)   const { return fCurrMatchFlagsPostBP[sector]; }
+
+#ifdef TRD_RAW_DEBUG
+  UInt_t GetBC(Int_t hc)               const { return fCurrBC[hc]; }
+  UInt_t GetEvCount(Int_t det)     const { return fCurrEvCount[det]; }
+  UInt_t GetL0Count(Int_t sector)  const { return fCurrL0Count[sector]; }
+  UInt_t GetL1aCount(Int_t sector) const { return fCurrL1aCount[sector]; }
+  UInt_t GetL1rCount(Int_t sector) const { return fCurrL1rCount[sector]; }
+  UInt_t GetL2aCount(Int_t sector) const { return fCurrL2aCount[sector]; }
+  UInt_t GetL2rCount(Int_t sector) const { return fCurrL2rCount[sector]; }
+  UInt_t GetChecksumStack(Int_t sector, Int_t stack) const { return fCurrChecksumStack[sector][stack]; }
+#endif
+
   // raw data dumping
   void SetDumpMCM(Int_t det, Int_t rob, Int_t mcm, Bool_t dump = kTRUE);
 
   Bool_t IsDumping() const { return (fNDumpMCMs > 0); }
   Bool_t DumpingMCM(Int_t det, Int_t rob, Int_t mcm) const;
 
-  TString DumpRaw(TString title, UInt_t *start, Int_t length, UInt_t endmarker = 0xffffffff);
+  TString DumpRaw(TString title, const UInt_t *start, Int_t length, UInt_t endmarker = 0xffffffff);
   TString DumpMcmHeader(TString title, UInt_t word);
   TString DumpAdcMask(TString title, UInt_t word);
 
+  static void SortTracklets(TClonesArray *trklArray, TList &sortedTracklets, Int_t *indices);
+  static void AssignTracklets(AliESDTrdTrack *trdTrack, Int_t *trackletIndex, Int_t refIndex[6]);
+
   // temporary: allow to change expected readout order
   static void SetMCMReadoutPos(Int_t mcm, Int_t pos) { if (mcm > -1 && mcm < 16) fgMcmOrder[mcm] = pos; }
   static void SetROBReadoutPos(Int_t robpair, Int_t pos) { if (robpair > -1 && robpair < 4) fgMcmOrder[robpair] = pos; }
@@ -193,6 +229,9 @@ class AliTRDrawStream : public TObject
   Int_t ReadZSData();
   Int_t ReadNonZSData();
 
+  UShort_t CalcLinkChecksum(UInt_t *data, Int_t size);
+
+  Int_t SeekNextStack();
   Int_t SeekNextLink();
 
   // MCM header decoding
@@ -215,11 +254,11 @@ class AliTRDrawStream : public TObject
   Int_t CouldBeADCmask(UInt_t adcmask) const { return ((0xf & adcmask) == 0xc && (0x3 & adcmask >> 30) == 0x1); }
 
   // error message generation
-  void EquipmentError(ErrorCode_t err = kUnknown, const char *const msg = "", ...);
-  void StackError    (ErrorCode_t err = kUnknown, const char *const msg = "", ...);
-  void LinkError     (ErrorCode_t err = kUnknown, const char *const msg = "", ...);
-  void ROBError      (ErrorCode_t err = kUnknown, const char *const msg = "", ...);
-  void MCMError      (ErrorCode_t err = kUnknown, const char *const msg = "", ...);
+  void EquipmentError(ErrorCode_t err = kUnknown, const char *const msg = " ", ...);
+  void StackError    (ErrorCode_t err = kUnknown, const char *const msg = " ", ...);
+  void LinkError     (ErrorCode_t err = kUnknown, const char *const msg = " ", ...);
+  void ROBError      (ErrorCode_t err = kUnknown, const char *const msg = " ", ...);
+  void MCMError      (ErrorCode_t err = kUnknown, const char *const msg = " ", ...);
   void StoreErrorTree() { fErrors->Fill(); }
   void StoreErrorArray() { new ((*fMarkers)[fMarkers->GetEntriesFast()]) AliTRDrawStreamError(fLastError); }
   void ForgetError() { return; }
@@ -239,6 +278,8 @@ class AliTRDrawStream : public TObject
   UInt_t fErrorFlags;                           // error flags used to steer subsequent reading
   char   fErrorBuffer[100];                     // buffer for error message
 
+  AliTRDrawStats fStats;            // event statistics, clearing must be done by the user
+
   UInt_t *fPayloadStart;                        // pointer to start of data payload
   UInt_t *fPayloadCurr;                         // pointer to current reading position in the payload
   Int_t   fPayloadSize;                         // size of the payload (in UInt_t words)
@@ -249,6 +290,7 @@ class AliTRDrawStream : public TObject
   static const Int_t fgkNtriggers;              // number of triggers in data stream
   static const UInt_t fgkDataEndmarker;         // data endmarker
   static const UInt_t fgkTrackletEndmarker;     // tracklet endmarker
+  static const UInt_t fgkStackEndmarker[];      // stack endmarker (used from version 0xd on)
   static       Int_t fgMcmOrder [];             // expected readout order of the MCMs
   static       Int_t fgRobOrder [];             // expected readout order of the ROBs
 
@@ -274,20 +316,32 @@ class AliTRDrawStream : public TObject
   UInt_t fCurrTrgHeaderAvail;                  // current trigger information availability
   UInt_t fCurrTrgHeaderReadout;                 // current readout mode for the trigger headers
   UInt_t fCurrTrkHeaderAvail;                  // current tracking information availability
+  UInt_t fCurrStackEndmarkerAvail;              // current stack endmarker availability
   UInt_t fCurrEvType;                          // current event type
   UInt_t fCurrTriggerEnable;                    // current trigger enable
   UInt_t fCurrTriggerFired;                     // current trigger fired
   UInt_t fCurrTrackEnable;                     // current value of track enable
   UInt_t fCurrTrackletEnable;                  // current value of tracklet enable
   UInt_t fCurrStackMask;                       // current mask of active stacks
+#ifdef TRD_RAW_DEBUG
+  UInt_t *fCurrL0Count;                        // number of received L0 triggers
+  UInt_t *fCurrL1aCount;                       // number of received L1a triggers
+  UInt_t *fCurrL1rCount;                       // number of received L1r triggers
+  UInt_t *fCurrL2aCount;                       // number of received L2a triggers
+  UInt_t *fCurrL2rCount;                       // number of received L2r triggers
+   Int_t fCurrL0offset[540];                   // current offset for L0 accepts from GTU and chambers
+  UInt_t fCurrEvCount[540];                    // current event count from the MCMs
+#endif
 
   // Tracking header
   UInt_t *fCurrTrkHeaderIndexWord;              // current tracking header index word
   UInt_t *fCurrTrkHeaderSize;                   // current tracking header index word
+  ULong64_t *fCurrTrkFlags;                     // current tracking done flags
 
   // Trigger header
   UInt_t *fCurrTrgHeaderIndexWord;              // current tracking header index word
   UInt_t *fCurrTrgHeaderSize;                   // current tracking header index word
+  UInt_t *fCurrTrgFlags;                       // current trigger flags of all sectors
 
   // Stack header
   UInt_t *fCurrStackIndexWord;                 // current stack index words
@@ -302,6 +356,12 @@ class AliTRDrawStream : public TObject
   UInt_t *fCurrLinkDataTypeFlags;              // current link data flags
   UInt_t *fCurrLinkDebugFlags;                 // current link debug flags
 
+  // CRC checks from trailer
+  Char_t fCurrMatchFlagsSRAM[18];
+  Char_t fCurrMatchFlagsPostBP[18];
+  UInt_t fCurrChecksumStack[18][5];
+  UInt_t fCurrChecksumSIU;
+
   // HC information
   Int_t fCurrSpecial;                          // current value of the special flag
   Int_t fCurrMajor;                            // current major version
@@ -314,21 +374,19 @@ class AliTRDrawStream : public TObject
   Int_t fCurrHC;                               // current HC
   Int_t fCurrCheck;                            // current check bits
   Int_t fCurrNtimebins;                                // current number of timebins
-  Int_t fCurrBC;                               // current BC
   Int_t fCurrPtrgCnt;                          // current pretrigger count
   Int_t fCurrPtrgPhase;                                // current pretrigger phase
+#ifdef TRD_RAW_DEBUG
+  Int_t fCurrBC[1080];                         // current BC
+#endif
 
   // settings for dumping
   Int_t fDumpMCM[100];                         // MCMs to dump
   Int_t fNDumpMCMs;                             // number of MCMs to dump
 
-  // tracklet information
-  TClonesArray *fTrackletArray;                        // pointer to array for tracklet storage
-
   // output data
   AliTRDarrayADC *fAdcArray;                   // pointer to ADC array
   AliTRDSignalIndex *fSignalIndex;             // pointer to the signal index
-  TTree *fTrackletTree;                        // pointer to the tree for tracklet storage
   TClonesArray *fTracklets;                    // pointer to array of tracklets
   TClonesArray *fTracks;                       // pointer to array of GTU tracks
   TClonesArray *fMarkers;                      // pointer to array of markers (data present, errors, ...)