]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDrawStream.h
Moving reconstruction settings to the reco params, new cosmics reco param OCDB entry...
[u/mrichter/AliRoot.git] / TRD / AliTRDrawStream.h
index 19d42f1d121980f122acc24b6875abce3b1232f7..14af8d7fe9e4c7dfe8dc5ad4b5898e11661abecf 100644 (file)
-#ifndef ALITRDRAWSTREAM_H\r\r
-#define ALITRDRAWSTREAM_H\r\r
-\r\r
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r\r
- * See cxx source for full Copyright notice                               */\r\r
-\r\r
-/* $Id: AliTRDrawStream.h 27696 2008-07-31 09:18:53Z cblume $ */\r\r
-\r\r
-///////////////////////////////////////////////////////////////////////////////\r\r
-//                                                                           //\r\r
-// This class provides access to TRD digits in raw data.                     //\r\r
-//                                                                           //\r\r
-///////////////////////////////////////////////////////////////////////////////\r\r
-\r\r
-#include "AliTRDrawStreamBase.h"\r\r
-\r\r
-class AliTRDgeometry;\r\r
-class AliRawReader;\r\r
-class AliTRDdigitsManager;\r\r
-class TTreeSRedirector;\r\r
-class AliTRDfeeParam;\r\r
-\r\r
-\r\r
-class AliTRDrawStream : public AliTRDrawStreamBase\r\r
-{ // class def begin\r\r
-\r\r
- public:\r\r
-    \r\r
-  //--------------------------------------------------------\r\r
-  // THE STRUCTURES\r\r
-  //--------------------------------------------------------\r\r
-\r\r
-  struct AliTRDrawADC\r\r
-  {//adc struct\r\r
-\r\r
-    UInt_t             *fPos;                    //! position of ADC 1st word in the buffer\r\r
-    Short_t             fADCnumber;              // number of the ADC 0 .. 20\r\r
-    Short_t             fExtendedCOL;            // column - row from MCM\r\r
-    Short_t             fCOL;                    // column - row from MCM\r\r
-    Int_t               fSignals[TRDMAXTBINS];   // signals for this adc\r\r
-    Bool_t              fIsShared;               // is pad chared between MCMs\r\r
-    Short_t             fCorrupted;              // is adc word corrupted\r\r
-\r\r
-    AliTRDrawADC()\r\r
-      : fPos(0)\r\r
-      , fADCnumber(0)\r\r
-      , fExtendedCOL(0)\r\r
-      , fCOL(0)\r\r
-      , fSignals()\r\r
-      , fIsShared(kFALSE)\r\r
-      , fCorrupted(0)\r\r
-    {\r\r
-      // default constructor\r\r
-    };\r\r
-\r\r
-    AliTRDrawADC(const AliTRDrawADC& p): \r\r
-        fPos(p.fPos)\r\r
-      , fADCnumber(p.fADCnumber)\r\r
-      , fExtendedCOL(p.fExtendedCOL)\r\r
-      , fCOL(p.fCOL)\r\r
-      , fSignals()\r\r
-      , fIsShared(p.fIsShared)\r\r
-      , fCorrupted(p.fCorrupted)\r\r
-    {\r\r
-      // copy constructor\r\r
-    };\r\r
-\r\r
-    AliTRDrawADC &operator=(const AliTRDrawADC &) \r\r
-    {\r\r
-      // assignment operator\r\r
-      // not implemented\r\r
-      return *this;\r\r
-    };\r\r
-\r\r
-  };\r\r
-  \r\r
-  //--------------------------------------------------------\r\r
-\r\r
-  struct AliTRDrawMCM\r\r
-  { // mcm struct\r\r
-    Int_t               fROB;                     // ROB number\r\r
-    Int_t               fMCM;                     // MCM number\r\r
-    Int_t               fROW;                     // row number filed during decoding\r\r
-      \r\r
-    UInt_t              fEvCounter;               // MCM event counter\r\r
-    UInt_t              fADCMask;                 // ADC mask\r\r
-    UInt_t              fADCMaskWord;             // word with ADC mask in\r\r
-    UInt_t              fADCchannel[TRDMAXADC]; // channels to be decoded accrording to ADC mask\r\r
-      \r\r
-    Int_t               fADCmax;                  // number of ADCs fired\r\r
-    Int_t               fADCcount;                // number of ADCs fired from double checking bit\r\r
-    Int_t               fMCMADCWords;             // mcm words to expect\r\r
-    Int_t               fSingleADCwords;          // n of words per ADC\r\r
-      \r\r
-    Int_t               fMCMhdCorrupted;          // is mcm header corrupted\r\r
-    Int_t               fADCmaskCorrupted;        // is mcm adc mask corrupted\r\r
-    Int_t               fCorrupted;               // is mcm data missing\r\r
-      \r\r
-    UInt_t             *fPos;                     //! position of mcm header in the buffer\r\r
-    UInt_t             *fAdcDataPos;              //! start of ADC data for this mcm\r\r
-\r\r
-    Int_t               fADCcounter;              // count the adcs decoded\r\r
-    AliTRDrawADC        fADCs[TRDMAXADC];       // 21 adcs\r\r
-      \r\r
-    AliTRDrawMCM()\r\r
-      : fROB(-1)\r\r
-      , fMCM(-1)\r\r
-      , fROW(-1)\r\r
-      , fEvCounter(0)\r\r
-      , fADCMask(0)\r\r
-      , fADCMaskWord(0)\r\r
-      , fADCchannel()      \r\r
-      , fADCmax(0)\r\r
-      , fADCcount(0)\r\r
-      , fMCMADCWords(0)      \r\r
-      , fSingleADCwords(0)\r\r
-      , fMCMhdCorrupted(0)      \r\r
-      , fADCmaskCorrupted(0)      \r\r
-      , fCorrupted(0)      \r\r
-      , fPos(0)\r\r
-      , fAdcDataPos(0)\r\r
-      , fADCcounter(0)\r\r
-      , fADCs()\r\r
-    {\r\r
-      // default constructor\r\r
-    };\r\r
-\r\r
-    AliTRDrawMCM(const AliTRDrawMCM & p):\r\r
-        fROB(p.fROB)\r\r
-      , fMCM(p.fMCM)\r\r
-      , fROW(p.fROW)\r\r
-      , fEvCounter(p.fEvCounter)\r\r
-      , fADCMask(p.fADCMask)\r\r
-      , fADCMaskWord(p.fADCMaskWord)\r\r
-      , fADCchannel()      \r\r
-      , fADCmax(p.fADCmax)\r\r
-      , fADCcount(p.fADCcount)\r\r
-      , fMCMADCWords(p.fMCMADCWords)      \r\r
-      , fSingleADCwords(p.fSingleADCwords)\r\r
-      , fMCMhdCorrupted(p.fMCMhdCorrupted)      \r\r
-      , fADCmaskCorrupted(p.fADCmaskCorrupted)      \r\r
-      , fCorrupted(p.fCorrupted)      \r\r
-      , fPos(p.fPos)\r\r
-      , fAdcDataPos(p.fAdcDataPos)\r\r
-      , fADCcounter(p.fADCcounter)\r\r
-      , fADCs()\r\r
-    {\r\r
-      // copy constructor\r\r
-    };\r\r
-\r\r
-    AliTRDrawMCM &operator=(const AliTRDrawMCM &)\r\r
-    {\r\r
-      // assignment operator\r\r
-      // not implemented\r\r
-      return *this;\r\r
-    };\r\r
-\r\r
-  };\r\r
-\r\r
-  //--------------------------------------------------------\r\r
-\r\r
-  struct AliTRDrawHC\r\r
-  { // hc struct\r\r
-      \r\r
-    //tacklet words of given HC\r\r
-    UInt_t              fTrackletWords[MAXTRACKLETSPERHC]; // array to keep tracklet words [mj]\r\r
-    Short_t             fTrackletError;                      // tracklet error \r\r
-    Short_t             fNTracklets;                    // number of tracklet  \r\r
-\r\r
-    // header word 0\r\r
-    Int_t               fSpecialRawV;       // Raw data version\r\r
-    Int_t               fRawVMajor;         // Raw data version\r\r
-    Int_t               fRawVMajorOpt;      // Raw data version\r\r
-    Int_t               fRawVMinor;         // Raw data version\r\r
-    Int_t               fNExtraWords;       // N extra HC header words\r\r
-    Int_t               fDCSboard;          // DCS board number\r\r
-    Int_t               fSM;                // Super Module number\r\r
-    Int_t               fStack;             // Stack number (some people might call it a chamber)\r\r
-    Int_t               fLayer;             // Layer number (some people might call it a plane)\r\r
-    Int_t               fSide;              // Side of HC\r\r
-      \r\r
-    // header word 1\r\r
-    Int_t               fTimeBins;          // N of t bins\r\r
-    UInt_t              fBunchCrossCounter; // Bunch crossing counter\r\r
-    UInt_t              fPreTriggerCounter; // Pre Trigger counter\r\r
-    UInt_t              fPreTriggerPhase;   // Pre Trigger phase\r\r
-      \r\r
-    // error \r\r
-    Int_t               fH0Corrupted;       // is hc header 0 corrupted \r\r
-    Int_t               fH1Corrupted;       // is hc header 1 corrupted\r\r
-    Int_t               fCorrupted;         // is hc data corrupted \r\r
-\r\r
-    UInt_t             *fPos[2];            //! position of the header words in buffer\r\r
-      \r\r
-    Int_t               fDET;               // filled while decoding\r\r
-    Int_t               fROC;               // filled while decoding\r\r
-    Int_t               fRowMax;            // filled while decoding\r\r
-    Int_t               fColMax;            // filled while decoding\r\r
-\r\r
-    // hc data\r\r
-    Int_t               fMCMmax;            // number of mcm found\r\r
-    AliTRDrawMCM        fMCMs[TRDMAXMCM]; // 4 ROBS 16 each \r\r
-\r\r
-    AliTRDrawHC()\r\r
-      : fTrackletWords() //[mj]\r\r
-      , fTrackletError(0)\r\r
-      , fNTracklets(0)\r\r
-      , fSpecialRawV(0)\r\r
-      , fRawVMajor(0)\r\r
-      , fRawVMajorOpt(0)\r\r
-      , fRawVMinor(0)\r\r
-      , fNExtraWords(0)\r\r
-      , fDCSboard(-1)\r\r
-      , fSM(-1)\r\r
-      , fStack(-1)\r\r
-      , fLayer(-1)\r\r
-      , fSide(-1)\r\r
-      , fTimeBins(0)\r\r
-      , fBunchCrossCounter(0)\r\r
-      , fPreTriggerCounter(0)\r\r
-      , fPreTriggerPhase(0)\r\r
-      , fH0Corrupted(0)\r\r
-      , fH1Corrupted(0)\r\r
-      , fCorrupted(0)\r\r
-      , fPos()\r\r
-      , fDET(-1)\r\r
-      , fROC(-1)\r\r
-      , fRowMax(-1)\r\r
-      , fColMax(-1)\r\r
-      , fMCMmax(0)\r\r
-      , fMCMs()\r\r
-    {\r\r
-      // default constructor hc info \r\r
-    };\r\r
-\r\r
-    AliTRDrawHC(const AliTRDrawHC & p):\r\r
-        fTrackletWords() //[mj]\r\r
-      , fTrackletError(p.fTrackletError)\r\r
-      , fNTracklets(p.fNTracklets)\r\r
-      , fSpecialRawV(p.fSpecialRawV)\r\r
-      , fRawVMajor(p.fRawVMajor)\r\r
-      , fRawVMajorOpt(p.fRawVMajorOpt)\r\r
-      , fRawVMinor(p.fRawVMinor)\r\r
-      , fNExtraWords(p.fNExtraWords)\r\r
-      , fDCSboard(p.fDCSboard)\r\r
-      , fSM(p.fSM)\r\r
-      , fStack(p.fStack)\r\r
-      , fLayer(p.fLayer)\r\r
-      , fSide(p.fSide)\r\r
-      , fTimeBins(p.fTimeBins)\r\r
-      , fBunchCrossCounter(p.fBunchCrossCounter)\r\r
-      , fPreTriggerCounter(p.fPreTriggerCounter)\r\r
-      , fPreTriggerPhase(p.fPreTriggerPhase)\r\r
-      , fH0Corrupted(p.fH0Corrupted)\r\r
-      , fH1Corrupted(p.fH1Corrupted)\r\r
-      , fCorrupted(p.fCorrupted)\r\r
-      , fPos()\r\r
-      , fDET(p.fDET)\r\r
-      , fROC(p.fROC)\r\r
-      , fRowMax(p.fRowMax)\r\r
-      , fColMax(p.fColMax)\r\r
-      , fMCMmax(p.fMCMmax)\r\r
-      , fMCMs()\r\r
-    {\r\r
-      // copy constructor\r\r
-    };\r\r
-\r\r
-    AliTRDrawHC &operator=(const AliTRDrawHC &)\r\r
-    {\r\r
-      // assignment operator\r\r
-      // not implemented\r\r
-      return *this;\r\r
-    };\r\r
-\r\r
-  };\r\r
-\r\r
-  //--------------------------------------------------------\r\r
-    \r\r
-  struct AliTRDrawStack\r\r
-  {\r\r
-    UInt_t           fHeaderSize;            // header size of the stack info\r\r
-    Bool_t           fLinksActive[12];       // data links active - 1 per half chamber\r\r
-    Short_t          fLinksDataType[12];     // 0 indicating real data for the front-end electronics \r\r
-    Short_t          fLinksMonitor[12];      // 0 indicating properly operating link \r\r
-    Short_t          fLinkMonitorError[12];  // record link monitor error\r\r
-    Int_t            fActiveLinks;           // number of active links\r\r
-    UInt_t          *fPos;                   //! position in the buffer\r\r
-                   \r\r
-    AliTRDrawHC      fHalfChambers[12];      // 12 half chambers in a stack\r\r
-      \r\r
-    AliTRDrawStack()\r\r
-      : fHeaderSize(0)\r\r
-      , fLinksActive()\r\r
-      , fLinksDataType()\r\r
-      , fLinksMonitor()\r\r
-      , fLinkMonitorError()\r\r
-      , fActiveLinks(0)\r\r
-      , fPos(0)\r\r
-      , fHalfChambers()\r\r
-    {\r\r
-      // default constructor\r\r
-    };      \r\r
-\r\r
-    AliTRDrawStack(const AliTRDrawStack & p):\r\r
-        fHeaderSize(p.fHeaderSize)\r\r
-      , fLinksActive()\r\r
-      , fLinksDataType()\r\r
-      , fLinksMonitor()\r\r
-      , fLinkMonitorError()\r\r
-      , fActiveLinks(p.fActiveLinks)\r\r
-      , fPos(p.fPos)\r\r
-      , fHalfChambers()\r\r
-    {\r\r
-      // copy constructor\r\r
-    };\r\r
-\r\r
-    AliTRDrawStack &operator=(const AliTRDrawStack &)\r\r
-    {\r\r
-      // assignment operator\r\r
-      // not implemented\r\r
-      return *this;\r\r
-    };\r\r
-\r\r
-  };\r\r
-\r\r
-  //--------------------------------------------------------\r\r
-\r\r
-  struct AliTRDrawSM\r\r
-  {\r\r
-    UInt_t            fHeaderSize;            // size of the header in words\r\r
-    Bool_t            fTrackletEnable;        // tracklet enable bit\r\r
-    Bool_t            fStackActive[5];        // map of active/expected stacks\r\r
-    Int_t             fActiveStacks;          // number of active stacks\r\r
-    Int_t             fCorrupted;             // is sm info corrupted\r\r
-    Int_t             fNexpectedHalfChambers; // number of half chambers to be read out in this sm\r\r
-    Bool_t            fClean;                 // true if everything went OK - false is some error occured\r\r
-    UInt_t           *fPos;                   // location of the sm info - should be the first word (after CDH if not DDL buffer)\r\r
-\r\r
-    AliTRDrawStack    fStacks[5];             // objects for 5 stacks\r\r
-\r\r
-    AliTRDrawSM()\r\r
-      : fHeaderSize(0)\r\r
-      , fTrackletEnable(0)\r\r
-      , fStackActive()\r\r
-      , fActiveStacks(0)\r\r
-      , fCorrupted(0)\r\r
-      , fNexpectedHalfChambers(0)\r\r
-      , fClean(kTRUE)\r\r
-      , fPos(0)\r\r
-      , fStacks()\r\r
-    {\r\r
-           // Default constructor\r\r
-    };      \r\r
-\r\r
-    AliTRDrawSM(const AliTRDrawSM & p):\r\r
-        fHeaderSize(p.fHeaderSize)\r\r
-      , fTrackletEnable(p.fTrackletEnable)\r\r
-      , fStackActive()\r\r
-      , fActiveStacks(p.fActiveStacks)\r\r
-      , fCorrupted(p.fCorrupted)\r\r
-      , fNexpectedHalfChambers(p.fNexpectedHalfChambers)\r\r
-      , fClean(p.fClean)\r\r
-      , fPos(p.fPos)\r\r
-      , fStacks()\r\r
-    {\r\r
-      // copy constructor\r\r
-    };\r\r
-\r\r
-    AliTRDrawSM &operator=(const AliTRDrawSM &)\r\r
-    {\r\r
-      // assignment operator\r\r
-      // not implemented\r\r
-      return *this;\r\r
-    };\r\r
-\r\r
-  };\r\r
-  \r\r
-  //--------------------------------------------------------\r\r
-     \r\r
-  AliTRDrawStream();\r\r
-  AliTRDrawStream(AliRawReader *rawReader);\r\r
-  virtual ~AliTRDrawStream();\r\r
-\r\r
-  //--------------------------------------------------------\r\r
-\r\r
-  virtual Bool_t       Next();                           // read the next data in the memory\r\r
-  //virtual Int_t      NextChamber(AliTRDdigitsManager *man); // read next chamber data in the momory\r\r
-  virtual Int_t        NextChamber(AliTRDdigitsManager *man, UInt_t **trackletContainer, UShort_t **/*errorCodeContainer*/=NULL); // read next chamber data in the memory\r\r
-  virtual Bool_t       Init();                           // initialize some internal variables\r\r
-\r\r
-  Int_t    NextBuffer(); // go and init next buffer if available - check the implementation file for return values\r\r
-\r\r
-  Bool_t   SetRawVersion(Int_t fraw); // set the raw version - used for backward compat.\r\r
-  \r\r
-  Bool_t   IsCurrentPadShared() const {return fADC->fIsShared;} // is current pad shared between mcms\r\r
-  void     SetSharedPadReadout(Bool_t fv) {fSharedPadsOn = fv;} // set the flag on if the reader should return the shared pads\r\r
-  \r\r
-  Bool_t   DecodeSM(void *buffer, UInt_t length); // decode a buffer\r\r
-  Int_t    DecodeSM();                            // used with raw reader\r\r
-  Int_t    DecodeSM(AliRawReader *reader);        // used with raw reader\r\r
-          \r\r
-  Bool_t   SetReader(AliRawReader *reader); // set the raw reader to use\r\r
-          \r\r
-  // info from Supermodule Index Word\r\r
-  Bool_t    IsTrackletEnableBitSet() const {return fSM.fTrackletEnable;} // get status of tracklet enable bit\r\r
-  Bool_t    IsStackActive(Int_t is) const {return fSM.fStackActive[is];} // get status of stack enable bit\r\r
-  Int_t     GetNofActiveStacks() const {return fSM.fActiveStacks;}       // get number of active stacks from stack mask\r\r
-  UInt_t   *GetGTUheaderWords() const {return fSM.fPos;}       // get number of active stacks from stack mask\r\r
-\r\r
-  // info from Stack Index Word\r\r
-  Int_t     GetNexpectedHalfChambers() const {return fSM.fNexpectedHalfChambers;}                    // get number of expected HC in a sm\r\r
-  Int_t     GetNofActiveLinksInStack(Int_t is) const {return fSM.fStacks[is].fActiveLinks;}          // get number of active links in a stack\r\r
-  Bool_t    IsLinkActiveInStack(Int_t is, Int_t il) const {return fSM.fStacks[is].fLinksActive[il];} // check whether the link is active\r\r
-\r\r
-  // info from Stack Header Word\r\r
-  Short_t    GetLinkMonitorError(Int_t is, Int_t il) const {return fSM.fStacks[is].fLinkMonitorError[il];} // get link monitor error\r\r
-\r\r
-  // info from Tracklet Data\r\r
-  Int_t     GetTrackletErrorCode(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fTrackletError;}\r\r
-  Int_t     GetNTracklets(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fNTracklets;} // get number of tracklets \r\r
-\r\r
-  // info from HC Header Word\r\r
-  Int_t     GetSM(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fSM;}\r\r
-  Int_t     GetLayer(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fLayer;}\r\r
-  Int_t     GetStack(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fStack;}\r\r
-  Int_t     GetSide(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fSide;}\r\r
-  Int_t     GetH0ErrorCode(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fH0Corrupted;}\r\r
-  Int_t     GetH1ErrorCode(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fH1Corrupted;}\r\r
-  Int_t     GetNumberOfTimeBins(Int_t is, Int_t il) const { return fSM.fStacks[is].fHalfChambers[il].fTimeBins;}\r\r
-  UInt_t   *GetTrackletWords(Int_t is, Int_t il) { return fSM.fStacks[is].fHalfChambers[il].fTrackletWords;}\r\r
-\r\r
-  // info from HC data\r\r
-  Int_t     GetHCErrorCode(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fCorrupted;}\r\r
-  Int_t     GetHCMCMmax(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fMCMmax;}\r\r
-\r\r
-  // from MCM Header Word\r\r
-  // rob and mcm ordering\r\r
-  // side 0(even link) - ROB: 0 2 4 6  MCM: 12 13 14 15 8 9 10 11 4 5 6 7 0 1 2 3  \r\r
-  // side 1( odd link) - ROB: 1 3 5 7  MCM: 12 13 14 15 8 9 10 11 4 5 6 7 0 1 2 3  \r\r
-  Int_t     GetMCM(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fMCM;}\r\r
-  Int_t     GetROB(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fROB;}\r\r
-  Int_t     GetMCMhdErrorCode(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fMCMhdCorrupted;}\r\r
-  Int_t     GetMCMADCMaskErrorCode(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCmaskCorrupted;}\r\r
-  Int_t     GetEventNumber(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fEvCounter;}\r\r
-  Int_t     GetADCcount(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCcount;}\r\r
-\r\r
-  // info from MCM data words\r\r
-  Int_t     GetMCMErrorCode(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fCorrupted;} // get MCM data error code\r\r
-  Int_t     GetADCErrorCode(Int_t stack, Int_t link, Int_t mcm, Int_t adc) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fCorrupted;} // get ADC error code\r\r
-  Int_t     GetADCnumber(Int_t stack, Int_t link, Int_t mcm, Int_t adc) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fADCnumber;} // get ADC error code\r\r
-\r\r
-  Int_t     GetRow(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fROW;}         // get current row number\r\r
-  Int_t     GetCol(Int_t stack, Int_t link, Int_t mcm, Int_t adc) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fCOL;}         // get current column number\r\r
-  Int_t     GetExtendedCol(Int_t stack, Int_t link, Int_t mcm, Int_t adc) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fExtendedCOL;}         // get current column number\r\r
-\r\r
-  // info from ADC data words\r\r
-  Int_t    *GetSignalDirect(Int_t stack, Int_t link, Int_t mcm, Int_t adc) {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fSignals;}\r\r
-\r\r
-\r\r
-  // from here, only works with returning ADC channel pointer using Next() \r\r
-  UInt_t   *GetTrackletWords() const { return fHC->fTrackletWords;}                 // return tracklet words pointer per hc [mj]\r\r
-  Int_t     GetTrackletErrorCode() const {return fHC ? fHC->fTrackletError : -1;}   // get tracklet error code\r\r
-  Int_t     GetNTracklets() const {return fHC ? fHC->fNTracklets : -1;}             // get number of tracklets \r\r
-\r\r
-  Int_t     GetSpecialRawVersion() const {return fHC ? fHC->fSpecialRawV : -1;}     // return special raw version\r\r
-  Int_t     GetMajorRawVersion() const {return fHC ? fHC->fRawVMajor : -1;}         // major raw version getter\r\r
-  Int_t     GetRawVersion() const {return fHC ? fHC->fRawVMajor : -1;}              // compatibility see funtion above\r\r
-  Int_t     GetMinorRawVersion() const {return fHC ? fHC->fRawVMinor : -1;}         // minor raw version\r\r
-\r\r
-  Int_t     GetSM() const {return fHC ? fHC->fSM : -1;}                              //  SM Position of CURRENT half chamber in full TRD\r\r
-  Int_t     GetLayer() const {return fHC ? fHC->fLayer : -1;}                        //  Layer Position of CURRENT half chamber in full TRD\r\r
-  Int_t     GetStack() const {return fHC ? fHC->fStack : -1;}                        //  Stack Position of CURRENT half chamber in full TRD\r\r
-  Int_t     GetSide() const {return fHC ? fHC->fSide : -1;}                          // get side\r\r
-  Int_t     GetDCS() const { return fHC ? fHC->fDCSboard : -1;}                      //  DCS board number read from data (HC header)\r\r
-  Int_t     GetROC() const { return fHC ? fHC->fROC : -1;}                           //  ROB Position of CURRENT half chamber in full TRD\r\r
-  Int_t     GetNumberOfTimeBins() const { return fHC ? fHC->fTimeBins : 0;}          // Get Ntime bins\r\r
-  UInt_t    GetBunchCrossCounter() const {return fHC ? fHC->fBunchCrossCounter : 0;} // get bunch cross counter\r\r
-  UInt_t    GetPreTriggerCounter() const {return fHC ? fHC->fPreTriggerCounter : 0;} // get pre trigger info\r\r
-  UInt_t    GetPreTriggerPhase() const {return fHC ? fHC->fPreTriggerPhase : 0;}     // get trigger phase\r\r
-\r\r
-  Int_t     GetRow() const {return fMCM ? fMCM->fROW : -1;}         // get current row number\r\r
-  Int_t     GetExtendedCol() const {return fADC ? fADC->fExtendedCOL : -1;}         // get current column number\r\r
-  Int_t     GetCol() const {return fADC ? fADC->fCOL : -1;}         // get current column number\r\r
-  Int_t     GetRowMax() const { return fHC ? fHC->fRowMax : -1;}    // Get maximum rows in the current HC\r\r
-  Int_t     GetColMax() const { return fHC ? fHC->fColMax : -1;}    // Get maximum cols in the current HC\r\r
-  // compatibility\r\r
-  Int_t     GetMaxRow() const { return fHC ? fHC->fRowMax : -1;}    // Get maximum rows in the current HC\r\r
-  Int_t     GetMaxCol() const { return fHC ? fHC->fColMax : -1;}    // Get maximum cols in the current HC\r\r
-\r\r
-  UInt_t    GetHCword0() const {return fHC ? *fHC->fPos[0] : 0;}    // get the HC word 0\r\r
-  UInt_t    GetHCword1() const {return fHC ? *fHC->fPos[1] : 0;}    // get the HC word 1\r\r
-           \r\r
-  Int_t     GetDET() const {return fHC ? fHC->fDET : -1;}           // get current det number\r\r
-  Int_t     GetDet() const {return fHC ? fHC->fDET : -1;}           // get current det number\r\r
-           \r\r
-  Int_t     GetROB() const {return fMCM ? fMCM->fROB : -1;}         // get current ROB number\r\r
-  Int_t     GetMCM() const {return fMCM ? fMCM->fMCM : -1;}         // get current MCM number\r\r
-  Int_t     GetEventNumber() const { return fMCM->fEvCounter;}      //  MCM Event number and position of current MCM on TRD chamber\r\r
-\r\r
-  Int_t     GetADC() const { return fADC ? fADC->fADCnumber : -1;}  //  MCM ADC channel and Time Bin of word 1\r\r
-  Int_t     GetTimeBin() const { return 0;}                         //  MCM ADC channel and Time Bin of word 1\r\r
-  Int_t    *GetSignals() const { return fADC ? fADC->fSignals : (Int_t *)fgEmptySignals;} // signals in the n-time bins from data word\r\r
-\r\r
-  Int_t     GetHCErrorCode() const {return fHC ? fHC->fCorrupted : -1;}    // get HC error code\r\r
-  Int_t     GetH0ErrorCode() const {return fHC ? fHC->fH0Corrupted : -1;}  // get HC header word0 error code\r\r
-  Int_t     GetH1ErrorCode() const {return fHC ? fHC->fH1Corrupted : -1;}  // get HC header word1 error code\r\r
-  Int_t     GetMCMErrorCode() const {return fMCM ? fMCM->fCorrupted : -1;} // get MCM data error code\r\r
-  Int_t     GetADCErrorCode() const {return fADC ? fADC->fCorrupted : -1;} // get ADC data error code\r\r
-  Int_t     GetMCMhdErrorCode() const {return fMCM ? fMCM->fMCMhdCorrupted: -1;} // get MCM header word error code\r\r
-  Int_t     GetMCMADCMaskErrorCode() const {return fMCM ? fMCM->fADCmaskCorrupted: -1;} // get MCM adc mask error code\r\r
-\r\r
-  UInt_t   *GetSMstreamPosition() const {return fSM.fPos;} // get position of the SM index word in the buffer\r\r
-\r\r
-  Bool_t    IsSMbufferClean() const {return fSM.fClean;}   // is data clean\r\r
-\r\r
-  //----------------------------------------------------------\r\r
\r\r
-  static void    SetNoDebug() {fgDebugFlag = kFALSE;} // allow debug info\r\r
-  static void    EnableMemoryReset() {fgEnableMemoryReset = kTRUE;} // allow memory reset\r\r
-  static void    SetNoErrorWarning() {fgWarnError = kFALSE;} // disable warning and error info\r\r
-  static void    SetForceCleanDataOnly() {fgCleanDataOnly = kTRUE;} // clean data only\r\r
-  static void    AllowCorruptedData() {fgCleanDataOnly = kFALSE;} // accept corrupted data\r\r
-\r\r
-  static void    SetExtraWordsFix() {fgExtraSkip = kTRUE;} // extra skip of 24 32-bit words \r\r
-  static void    SetSkipCDH() {fgSkipCDH = kTRUE;} // skip of 8 32-bit words \r\r
-  static void    SetDumpHead(Int_t iv) {fgDumpHead = iv;}\r\r
-  static void    DisableStackNumberChecker() {fgStackNumberChecker = kFALSE;}  // set false to cleanroom data \r\r
-  static void    DisableStackLinkNumberChecker() {fgStackLinkNumberChecker = kFALSE;}  \r\r
-  static void    DisableSkipData() {fgSkipData = kFALSE;} // keep reading next words even previous words were corrupted - debugging purpose  \r\r
-  static void    SetDumpingEnable() {fgDumpingEnable = kTRUE;} \r\r
-  static void    SetDumpingMCM(Int_t sm, Int_t stack, Int_t layer, Int_t rob, Int_t mcm) {fgDumpingSM = sm; fgDumpingStack = stack; fgDumpingLayer = layer; fgDumpingROB = rob; fgDumpingMCM = mcm;}\r\r
-\r\r
-  // this is a temporary solution!\r\r
-  // baseline should come with the HC header word 2 (count from 0!)\r\r
-  static void    SetSubtractBaseline(Int_t baseline) {fgCommonAdditive = baseline;}\r\r
-  Int_t          GetCommonAdditive() const {return fgCommonAdditive;} // return the common additive\r\r
-\r\r
-  static void    EnableDecodeConfigData() {fgEnableDecodeConfigData = kTRUE;} // allow configuration data decoding\r\r
-  static Bool_t  fgEnableDecodeConfigData; // flag to decode configuration data\r\r
-\r\r
-\r\r
-  //--------------------------------------------------------\r\r
-  // Decoding functions\r\r
-  //--------------------------------------------------------\r\r
-\r\r
-  void DecodeSMInfo(const UInt_t *word, struct AliTRDrawSM *sm) const ;\r\r
-  const char *DumpSMInfo(const struct AliTRDrawSM *sm);\r\r
-  void DecodeStackInfo(const UInt_t *word, struct AliTRDrawStack *st) const;\r\r
-  void DecodeStackHeader(const UInt_t *word, struct AliTRDrawStack *st, Int_t iword) const;\r\r
-  const char *DumpStackInfo(const struct AliTRDrawStack *st);\r\r
-  Bool_t DecodeHCwordH0(const UInt_t *word, struct AliTRDrawHC *hc) const;\r\r
-  Bool_t DecodeHCwordH1(const UInt_t *word, struct AliTRDrawHC *hc) const;\r\r
-  const char *DumpHCinfoH0(const struct AliTRDrawHC *hc);\r\r
-  const char *DumpHCinfoH1(const struct AliTRDrawHC *hc);\r\r
-  void DecodeMCMheader(const UInt_t *word, struct AliTRDrawMCM *mcm) const;\r\r
-  UInt_t GetMCMadcMask(const UInt_t *word, struct AliTRDrawMCM *mcm) const;\r\r
-  void DecodeMask(const UInt_t *word, struct AliTRDrawMCM *mcm) const;\r\r
-  void MCMADCwordsWithTbins(UInt_t fTbins, struct AliTRDrawMCM *mcm) const;\r\r
-  const char *DumpMCMinfo(const struct AliTRDrawMCM *mcm);\r\r
-  const char *DumpMCMadcMask(const struct AliTRDrawMCM *mcm);\r\r
-\r\r
-\r\r
- protected:\r\r
-\r\r
-  Bool_t InitBuffer(void *buffer, UInt_t length); // init the buffer - called by DecodeSM(void*, UInt_t)\r\r
-  Bool_t DumpWords(UInt_t *px, UInt_t iw, UInt_t marker = 0); // dump some words onto the screen\r\r
-\r\r
-  void   SwapOnEndian();         // swap if endian is BIG\r\r
-  Bool_t SkipWords(UInt_t iw);   // skip number of words\r\r
-  Bool_t DecodeGTUheader();      // decode data in GTU header\r\r
-  Bool_t DecodeTracklets();      // decode tracklets\r\r
-  Bool_t DecodeHC();             // decode data in HC\r\r
-\r\r
-  Bool_t DecodeADC();            // decode 10 ADC words\r\r
-  Bool_t DecodeADCExtended();    // decode TRAP data with ZS and more than 32 samples \r\r
-\r\r
-  Bool_t DecodeHCheader();       // decode HC  header\r\r
-  Bool_t SeekEndOfData();        // go to next end of raw data marker (actually 1 word after)\r\r
-  Bool_t SkipMCMdata(UInt_t iw);  // skip this mcm data due to mcm header corruption\r\r
-  Bool_t SeekNextMCMheader();    // go to next mcm header\r\r
-  Bool_t DecodeMCMheader();      // decode mcm header\r\r
-\r\r
-  Bool_t IsRowValid();       // check if row within the range\r\r
-  Bool_t IsHCheaderOK();     // check if current hc header data make sense\r\r
-  Bool_t IsMCMheaderOK();    // check if current mcm header data make sense\r\r
-    \r\r
-  void   ResetCounters();    // reset some counters\r\r
-  void   ResetIterators();   // needed for Next()\r\r
-  void   ResetPerSM();       // reset every SM \r\r
-  void   ResetPerStack();    // reset every Stack \r\r
-  void   ResetPerHC();       // reset every HC \r\r
-  void   ResetPerMCM();      // reset every MCM  \r\r
-  void   ResetPerADC();      // reset every ADC  \r\r
-  void   ResetMemory();      // reset all data members\r\r
-\r\r
-  AliTRDrawStream(const AliTRDrawStream& st);\r\r
-  AliTRDrawStream &operator=(const AliTRDrawStream &);\r\r
-\r\r
-  // ----------------- DATA MEMBERS START\r\r
-\r\r
-  struct AliTRDrawSM       fSM;    // one SM per buffer\r\r
-  struct AliTRDrawStack   *fStack; //! pointer to the current stack\r\r
-  struct AliTRDrawHC      *fHC;    //! current HC\r\r
-  struct AliTRDrawHC      *fLastHC; //! last HC\r\r
-  struct AliTRDrawMCM     *fMCM;   //! current MCM\r\r
-  struct AliTRDrawADC     *fADC;   //! current ADC\r\r
-  \r\r
-  UInt_t *fpPos;   // current position in the buffer\r\r
-  UInt_t *fpBegin; // begin - pointer to the buffer word 0\r\r
-  UInt_t *fpEnd;   // end of the buffer\r\r
-\r\r
-  UInt_t  fWordLength;  // length of the buffer in 32bit words\r\r
-\r\r
-  Bool_t  fIsGlobalDigitsParamSet; // set true if the global digit parameters are already set\r\r
-\r\r
-  Int_t   fStackNumber;     // current stack number\r\r
-  Int_t   fStackLinkNumber; // current link in the stack\r\r
-\r\r
-  Int_t   fhcMCMcounter;    // mcm counter inside single hc - used in Next()  \r\r
-  Int_t   fmcmADCcounter;   // adc counrer inside single adc - used in Next() \r\r
-\r\r
-  Int_t   fLinkTrackletCounter; // count the tracklets in the current HC\r\r
-  Int_t   fEndOfTrackletCount;  // count link by link (hc by hc) used for debug\r\r
-  Int_t   fNWordsCounter;       // counts words of given hc having link monitor error\r\r
-\r\r
-  UInt_t  fMaskADCword; // temp mask when decoding adcs\r\r
-  UInt_t  fTbinADC;     // temp adc \r\r
-  Int_t   fDecodedADCs; // counter of decoded adcs [mj] do we need?\r\r
-\r\r
-  UInt_t  fEventCounter;     // stores the valid/current MCM event counter\r\r
-  UInt_t  fLastEventCounter; // last known event counter of MCM\r\r
-\r\r
-  Bool_t  fSharedPadsOn; // do we want to output shared pads - default is off\r\r
-  Int_t   fMaxADCgeom;   // maximum ADC channels per mcm\r\r
-\r\r
-  Bool_t  fBufferRead; // set to be true if buffer is read for a given equipment\r\r
-\r\r
-  AliTRDgeometry *fGeometry;  //! TRD geometry\r\r
-  AliRawReader   *fRawReader; //! raw reader    \r\r
-\r\r
-  AliTRDfeeParam      *fTRDfeeParam; // pointer to the fee params\r\r
-\r\r
-  // STATIC \r\r
-\r\r
-  static Bool_t fgExtraSkip; // whether we should skip the leading 24 words\r\r
-  static Bool_t fgSkipCDH; // whether we should skip CDH (8 words)\r\r
-  static Bool_t fgWarnError; // no errors no warnings\r\r
-  static Bool_t fgCleanDataOnly; // release only clean events = no errors\r\r
-  static Bool_t fgDebugFlag; // allow debugging info\r\r
-  static Bool_t fgEnableMemoryReset; // allow memory reset\r\r
-  static Bool_t fgStackNumberChecker; // decide if we check stack number insanity - set false to cleanroom data\r\r
-  static Bool_t fgStackLinkNumberChecker; // decide if we check stack link number insanity - debuging purpose\r\r
-  static Bool_t fgSkipData; // decide if we skip corrupted data of given HC\r\r
-  static Int_t fgDumpHead; // number of words to dump (from the start of the buffer) on each Init\r\r
-  static Int_t  fgEmptySignals[30]; // empty signals in case of ADC pointer = NULL\r\r
-  static Short_t  fgMCMordering[16]; // mcm number odering for mcm header corruption check\r\r
-  static Short_t  fgROBordering[16]; // mcm number odering for mcm header corruption check\r\r
-  static Int_t  fgLastHC;  // HC number indicator to check HC numbering\r\r
-  static Int_t  fgLastROB; // ROB number indicator to check ROB numbering \r\r
-  static Int_t  fgLastIndex; // variable to check mcm number odering\r\r
-\r\r
-  static Bool_t fgDumpingEnable; // flag to dump MCM data words  \r\r
-\r\r
-  static Int_t  fgDumpingSM; // sm number to dump MCM data words\r\r
-  static Int_t  fgDumpingStack; // stack number to dump MCM data words\r\r
-  static Int_t  fgDumpingLayer; // layer number to dump MCM data words\r\r
-  static Int_t  fgDumpingROB; // rob number to dump MCM data words\r\r
-  static Int_t  fgDumpingMCM; // mcm number to dump MCM data words\r\r
-\r\r
-  // this is a temporary solution!\r\r
-  // baseline should come with the HC header word 2 (count from 0!)\r\r
-  static Int_t   fgCommonAdditive; // common additive  - should be decoded! from HC word2\r\r
-\r\r
-  // ----------------- DATA MEMBERS STOP\r\r
-\r\r
-  enum ETRDzRawStreamError \r\r
-    {\r\r
-      kDecodeStackInfo          = 1 //\r\r
-      , kMissingData             = 2 //\r\r
-      , kLinkDataMissing        = 3 //\r\r
-      , kHCdataMissing         = 4 //\r\r
-      , kTrackletOverflow       = 5 //\r\r
-      , kEOTrackeltsMissing     = 6 //\r\r
-      , kWrongPadrow             = 7 //\r\r
-      , kMCMheaderCorrupted     = 8 //\r\r
-      , kWrongMCMorROB         = 9 //\r\r
-      , kMCMeventMissmatch      = 10 //\r\r
-      , kMCMADCMaskMissing      = 11 //\r\r
-      , kHCHeaderCorrupt        = 12 //\r\r
-      , kHCHeaderWrongStack     = 13 //\r\r
-      , kHCHeaderWrongLayer     = 14 //\r\r
-      , kHCHeaderWrongSide      = 15 //\r\r
-      , kHCHeaderWrongSM        = 16 //\r\r
-      , kHCHeaderWrongDet       = 17 //\r\r
-      , kHCHeaderWrongROC       = 18 //\r\r
-      , kHCWordMissing         = 19 //\r\r
-      , kMCMdataMissing                = 20 //\r\r
-      , kMCMoverflow             = 21 //\r\r
-      , kADCdataMissing                = 22 //\r\r
-      , kADCmaskMissmatch       = 23 //\r\r
-      , kWrongPadcolumn                = 24 //\r\r
-    };                        \r\r
-\r\r
-  ClassDef(AliTRDrawStream, 0)\r\r
-}; //clas def end\r\r
-\r\r
-#endif\r\r
+#ifndef ALITRDRAWSTREAM_H\r\r\r
+#define ALITRDRAWSTREAM_H\r\r\r
+\r\r\r
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r\r\r
+ * See cxx source for full Copyright notice                               */\r\r\r
+\r\r\r
+/* $Id: AliTRDrawStream.h 27696 2008-07-31 09:18:53Z cblume $ */\r\r\r
+\r\r\r
+///////////////////////////////////////////////////////////////////////////////\r\r\r
+//                                                                           //\r\r\r
+// This class provides access to TRD digits in raw data.                     //\r\r\r
+//                                                                           //\r\r\r
+///////////////////////////////////////////////////////////////////////////////\r\r\r
+\r\r\r
+#include "AliTRDrawStreamBase.h"\r\r\r
+\r\r\r
+class AliTRDgeometry;\r\r\r
+class AliRawReader;\r\r\r
+class AliTRDdigitsManager;\r\r\r
+class TTreeSRedirector;\r\r\r
+class AliTRDfeeParam;\r\r\r
+\r\r\r
+\r\r\r
+class AliTRDrawStream : public AliTRDrawStreamBase\r\r\r
+{ // class def begin\r\r\r
+\r\r\r
+ public:\r\r\r
+    \r\r\r
+  //--------------------------------------------------------\r\r\r
+  // THE STRUCTURES\r\r\r
+  //--------------------------------------------------------\r\r\r
+\r\r\r
+  struct AliTRDrawADC\r\r\r
+  {//adc struct\r\r\r
+\r\r\r
+    UInt_t             *fPos;                    //! position of ADC 1st word in the buffer\r\r\r
+    Short_t             fADCnumber;              // number of the ADC 0 .. 20\r\r\r
+    Short_t             fExtendedCOL;            // column - row from MCM\r\r\r
+    Short_t             fCOL;                    // column - row from MCM\r\r\r
+    Int_t               fSignals[TRDMAXTBINS];   // signals for this adc\r\r\r
+    Bool_t              fIsShared;               // is pad chared between MCMs\r\r\r
+    Short_t             fCorrupted;              // is adc word corrupted\r\r\r
+\r\r\r
+    AliTRDrawADC()\r\r\r
+      : fPos(0)\r\r\r
+      , fADCnumber(0)\r\r\r
+      , fExtendedCOL(0)\r\r\r
+      , fCOL(0)\r\r\r
+      , fSignals()\r\r\r
+      , fIsShared(kFALSE)\r\r\r
+      , fCorrupted(0)\r\r\r
+    {\r\r\r
+      // default constructor\r\r\r
+    };\r\r\r
+\r\r\r
+    AliTRDrawADC(const AliTRDrawADC& p): \r\r\r
+        fPos(p.fPos)\r\r\r
+      , fADCnumber(p.fADCnumber)\r\r\r
+      , fExtendedCOL(p.fExtendedCOL)\r\r\r
+      , fCOL(p.fCOL)\r\r\r
+      , fSignals()\r\r\r
+      , fIsShared(p.fIsShared)\r\r\r
+      , fCorrupted(p.fCorrupted)\r\r\r
+    {\r\r\r
+      // copy constructor\r\r\r
+    };\r\r\r
+\r\r\r
+    AliTRDrawADC &operator=(const AliTRDrawADC &) \r\r\r
+    {\r\r\r
+      // assignment operator\r\r\r
+      // not implemented\r\r\r
+      return *this;\r\r\r
+    };\r\r\r
+\r\r\r
+  };\r\r\r
+  \r\r\r
+  //--------------------------------------------------------\r\r\r
+\r\r\r
+  struct AliTRDrawMCM\r\r\r
+  { // mcm struct\r\r\r
+    Int_t               fROB;                     // ROB number\r\r\r
+    Int_t               fMCM;                     // MCM number\r\r\r
+    Int_t               fROW;                     // row number filed during decoding\r\r\r
+      \r\r\r
+    UInt_t              fEvCounter;               // MCM event counter\r\r\r
+    UInt_t              fADCMask;                 // ADC mask\r\r\r
+    UInt_t              fADCMaskWord;             // word with ADC mask in\r\r\r
+    UInt_t              fADCchannel[TRDMAXADC]; // channels to be decoded accrording to ADC mask\r\r\r
+      \r\r\r
+    Int_t               fADCmax;                  // number of ADCs fired\r\r\r
+    Int_t               fADCcount;                // number of ADCs fired from double checking bit\r\r\r
+    Int_t               fMCMADCWords;             // mcm words to expect\r\r\r
+    Int_t               fSingleADCwords;          // n of words per ADC\r\r\r
+      \r\r\r
+    Int_t               fMCMhdCorrupted;          // is mcm header corrupted\r\r\r
+    Int_t               fADCmaskCorrupted;        // is mcm adc mask corrupted\r\r\r
+    Int_t               fCorrupted;               // is mcm data missing\r\r\r
+      \r\r\r
+    UInt_t             *fPos;                     //! position of mcm header in the buffer\r\r\r
+    UInt_t             *fAdcDataPos;              //! start of ADC data for this mcm\r\r\r
+\r\r\r
+    Int_t               fADCcounter;              // count the adcs decoded\r\r\r
+    AliTRDrawADC        fADCs[TRDMAXADC];       // 21 adcs\r\r\r
+      \r\r\r
+    AliTRDrawMCM()\r\r\r
+      : fROB(-1)\r\r\r
+      , fMCM(-1)\r\r\r
+      , fROW(-1)\r\r\r
+      , fEvCounter(0)\r\r\r
+      , fADCMask(0)\r\r\r
+      , fADCMaskWord(0)\r\r\r
+      , fADCchannel()      \r\r\r
+      , fADCmax(0)\r\r\r
+      , fADCcount(0)\r\r\r
+      , fMCMADCWords(0)      \r\r\r
+      , fSingleADCwords(0)\r\r\r
+      , fMCMhdCorrupted(0)      \r\r\r
+      , fADCmaskCorrupted(0)      \r\r\r
+      , fCorrupted(0)      \r\r\r
+      , fPos(0)\r\r\r
+      , fAdcDataPos(0)\r\r\r
+      , fADCcounter(0)\r\r\r
+      , fADCs()\r\r\r
+    {\r\r\r
+      // default constructor\r\r\r
+    };\r\r\r
+\r\r\r
+    AliTRDrawMCM(const AliTRDrawMCM & p):\r\r\r
+        fROB(p.fROB)\r\r\r
+      , fMCM(p.fMCM)\r\r\r
+      , fROW(p.fROW)\r\r\r
+      , fEvCounter(p.fEvCounter)\r\r\r
+      , fADCMask(p.fADCMask)\r\r\r
+      , fADCMaskWord(p.fADCMaskWord)\r\r\r
+      , fADCchannel()      \r\r\r
+      , fADCmax(p.fADCmax)\r\r\r
+      , fADCcount(p.fADCcount)\r\r\r
+      , fMCMADCWords(p.fMCMADCWords)      \r\r\r
+      , fSingleADCwords(p.fSingleADCwords)\r\r\r
+      , fMCMhdCorrupted(p.fMCMhdCorrupted)      \r\r\r
+      , fADCmaskCorrupted(p.fADCmaskCorrupted)      \r\r\r
+      , fCorrupted(p.fCorrupted)      \r\r\r
+      , fPos(p.fPos)\r\r\r
+      , fAdcDataPos(p.fAdcDataPos)\r\r\r
+      , fADCcounter(p.fADCcounter)\r\r\r
+      , fADCs()\r\r\r
+    {\r\r\r
+      // copy constructor\r\r\r
+    };\r\r\r
+\r\r\r
+    AliTRDrawMCM &operator=(const AliTRDrawMCM &)\r\r\r
+    {\r\r\r
+      // assignment operator\r\r\r
+      // not implemented\r\r\r
+      return *this;\r\r\r
+    };\r\r\r
+\r\r\r
+  };\r\r\r
+\r\r\r
+  //--------------------------------------------------------\r\r\r
+\r\r\r
+  struct AliTRDrawHC\r\r\r
+  { // hc struct\r\r\r
+      \r\r\r
+    //tacklet words of given HC\r\r\r
+    UInt_t              fTrackletWords[MAXTRACKLETSPERHC]; // array to keep tracklet words [mj]\r\r\r
+    Short_t             fTrackletError;                      // tracklet error \r\r\r
+    Short_t             fNTracklets;                    // number of tracklet  \r\r\r
+\r\r\r
+    // header word 0\r\r\r
+    Int_t               fSpecialRawV;       // Raw data version\r\r\r
+    Int_t               fRawVMajor;         // Raw data version\r\r\r
+    Int_t               fRawVMajorOpt;      // Raw data version\r\r\r
+    Int_t               fRawVMinor;         // Raw data version\r\r\r
+    Int_t               fNExtraWords;       // N extra HC header words\r\r\r
+    Int_t               fDCSboard;          // DCS board number\r\r\r
+    Int_t               fSM;                // Super Module number\r\r\r
+    Int_t               fStack;             // Stack number (some people might call it a chamber)\r\r\r
+    Int_t               fLayer;             // Layer number (some people might call it a plane)\r\r\r
+    Int_t               fSide;              // Side of HC\r\r\r
+      \r\r\r
+    // header word 1\r\r\r
+    Int_t               fTimeBins;          // N of t bins\r\r\r
+    UInt_t              fBunchCrossCounter; // Bunch crossing counter\r\r\r
+    UInt_t              fPreTriggerCounter; // Pre Trigger counter\r\r\r
+    UInt_t              fPreTriggerPhase;   // Pre Trigger phase\r\r\r
+      \r\r\r
+    // error \r\r\r
+    Int_t               fH0Corrupted;       // is hc header 0 corrupted \r\r\r
+    Int_t               fH1Corrupted;       // is hc header 1 corrupted\r\r\r
+    Int_t               fCorrupted;         // is hc data corrupted \r\r\r
+\r\r\r
+    UInt_t             *fPos[2];            //! position of the header words in buffer\r\r\r
+      \r\r\r
+    Int_t               fDET;               // filled while decoding\r\r\r
+    Int_t               fROC;               // filled while decoding\r\r\r
+    Int_t               fRowMax;            // filled while decoding\r\r\r
+    Int_t               fColMax;            // filled while decoding\r\r\r
+\r\r\r
+    // hc data\r\r\r
+    Int_t               fMCMmax;            // number of mcm found\r\r\r
+    AliTRDrawMCM        fMCMs[TRDMAXMCM]; // 4 ROBS 16 each \r\r\r
+\r\r\r
+    AliTRDrawHC()\r\r\r
+      : fTrackletWords() //[mj]\r\r\r
+      , fTrackletError(0)\r\r\r
+      , fNTracklets(0)\r\r\r
+      , fSpecialRawV(0)\r\r\r
+      , fRawVMajor(0)\r\r\r
+      , fRawVMajorOpt(0)\r\r\r
+      , fRawVMinor(0)\r\r\r
+      , fNExtraWords(0)\r\r\r
+      , fDCSboard(-1)\r\r\r
+      , fSM(-1)\r\r\r
+      , fStack(-1)\r\r\r
+      , fLayer(-1)\r\r\r
+      , fSide(-1)\r\r\r
+      , fTimeBins(0)\r\r\r
+      , fBunchCrossCounter(0)\r\r\r
+      , fPreTriggerCounter(0)\r\r\r
+      , fPreTriggerPhase(0)\r\r\r
+      , fH0Corrupted(0)\r\r\r
+      , fH1Corrupted(0)\r\r\r
+      , fCorrupted(0)\r\r\r
+      , fPos()\r\r\r
+      , fDET(-1)\r\r\r
+      , fROC(-1)\r\r\r
+      , fRowMax(-1)\r\r\r
+      , fColMax(-1)\r\r\r
+      , fMCMmax(0)\r\r\r
+      , fMCMs()\r\r\r
+    {\r\r\r
+      // default constructor hc info \r\r\r
+    };\r\r\r
+\r\r\r
+    AliTRDrawHC(const AliTRDrawHC & p):\r\r\r
+        fTrackletWords() //[mj]\r\r\r
+      , fTrackletError(p.fTrackletError)\r\r\r
+      , fNTracklets(p.fNTracklets)\r\r\r
+      , fSpecialRawV(p.fSpecialRawV)\r\r\r
+      , fRawVMajor(p.fRawVMajor)\r\r\r
+      , fRawVMajorOpt(p.fRawVMajorOpt)\r\r\r
+      , fRawVMinor(p.fRawVMinor)\r\r\r
+      , fNExtraWords(p.fNExtraWords)\r\r\r
+      , fDCSboard(p.fDCSboard)\r\r\r
+      , fSM(p.fSM)\r\r\r
+      , fStack(p.fStack)\r\r\r
+      , fLayer(p.fLayer)\r\r\r
+      , fSide(p.fSide)\r\r\r
+      , fTimeBins(p.fTimeBins)\r\r\r
+      , fBunchCrossCounter(p.fBunchCrossCounter)\r\r\r
+      , fPreTriggerCounter(p.fPreTriggerCounter)\r\r\r
+      , fPreTriggerPhase(p.fPreTriggerPhase)\r\r\r
+      , fH0Corrupted(p.fH0Corrupted)\r\r\r
+      , fH1Corrupted(p.fH1Corrupted)\r\r\r
+      , fCorrupted(p.fCorrupted)\r\r\r
+      , fPos()\r\r\r
+      , fDET(p.fDET)\r\r\r
+      , fROC(p.fROC)\r\r\r
+      , fRowMax(p.fRowMax)\r\r\r
+      , fColMax(p.fColMax)\r\r\r
+      , fMCMmax(p.fMCMmax)\r\r\r
+      , fMCMs()\r\r\r
+    {\r\r\r
+      // copy constructor\r\r\r
+    };\r\r\r
+\r\r\r
+    AliTRDrawHC &operator=(const AliTRDrawHC &)\r\r\r
+    {\r\r\r
+      // assignment operator\r\r\r
+      // not implemented\r\r\r
+      return *this;\r\r\r
+    };\r\r\r
+\r\r\r
+  };\r\r\r
+\r\r\r
+  //--------------------------------------------------------\r\r\r
+    \r\r\r
+  struct AliTRDrawStack\r\r\r
+  {\r\r\r
+    UInt_t           fHeaderSize;            // header size of the stack info\r\r\r
+    Bool_t           fLinksActive[12];       // data links active - 1 per half chamber\r\r\r
+    Short_t          fLinksDataType[12];     // 0 indicating real data for the front-end electronics \r\r\r
+    Short_t          fLinksMonitor[12];      // 0 indicating properly operating link \r\r\r
+    Short_t          fLinkMonitorError[12];  // record link monitor error\r\r\r
+    Int_t            fActiveLinks;           // number of active links\r\r\r
+    UInt_t          *fPos;                   //! position in the buffer\r\r\r
+                   \r\r\r
+    AliTRDrawHC      fHalfChambers[12];      // 12 half chambers in a stack\r\r\r
+      \r\r\r
+    AliTRDrawStack()\r\r\r
+      : fHeaderSize(0)\r\r\r
+      , fLinksActive()\r\r\r
+      , fLinksDataType()\r\r\r
+      , fLinksMonitor()\r\r\r
+      , fLinkMonitorError()\r\r\r
+      , fActiveLinks(0)\r\r\r
+      , fPos(0)\r\r\r
+      , fHalfChambers()\r\r\r
+    {\r\r\r
+      // default constructor\r\r\r
+    };      \r\r\r
+\r\r\r
+    AliTRDrawStack(const AliTRDrawStack & p):\r\r\r
+        fHeaderSize(p.fHeaderSize)\r\r\r
+      , fLinksActive()\r\r\r
+      , fLinksDataType()\r\r\r
+      , fLinksMonitor()\r\r\r
+      , fLinkMonitorError()\r\r\r
+      , fActiveLinks(p.fActiveLinks)\r\r\r
+      , fPos(p.fPos)\r\r\r
+      , fHalfChambers()\r\r\r
+    {\r\r\r
+      // copy constructor\r\r\r
+    };\r\r\r
+\r\r\r
+    AliTRDrawStack &operator=(const AliTRDrawStack &)\r\r\r
+    {\r\r\r
+      // assignment operator\r\r\r
+      // not implemented\r\r\r
+      return *this;\r\r\r
+    };\r\r\r
+\r\r\r
+  };\r\r\r
+\r\r\r
+  //--------------------------------------------------------\r\r\r
+\r\r\r
+  struct AliTRDrawSM\r\r\r
+  {\r\r\r
+    UInt_t            fHeaderSize;            // size of the header in words\r\r\r
+    Bool_t            fTrackletEnable;        // tracklet enable bit\r\r\r
+    Bool_t            fStackActive[5];        // map of active/expected stacks\r\r\r
+    Int_t             fActiveStacks;          // number of active stacks\r\r\r
+    Int_t             fCorrupted;             // is sm info corrupted\r\r\r
+    Int_t             fNexpectedHalfChambers; // number of half chambers to be read out in this sm\r\r\r
+    Bool_t            fClean;                 // true if everything went OK - false is some error occured\r\r\r
+    UInt_t           *fPos;                   // location of the sm info - should be the first word (after CDH if not DDL buffer)\r\r\r
+\r\r\r
+    AliTRDrawStack    fStacks[5];             // objects for 5 stacks\r\r\r
+\r\r\r
+    AliTRDrawSM()\r\r\r
+      : fHeaderSize(0)\r\r\r
+      , fTrackletEnable(0)\r\r\r
+      , fStackActive()\r\r\r
+      , fActiveStacks(0)\r\r\r
+      , fCorrupted(0)\r\r\r
+      , fNexpectedHalfChambers(0)\r\r\r
+      , fClean(kTRUE)\r\r\r
+      , fPos(0)\r\r\r
+      , fStacks()\r\r\r
+    {\r\r\r
+           // Default constructor\r\r\r
+    };      \r\r\r
+\r\r\r
+    AliTRDrawSM(const AliTRDrawSM & p):\r\r\r
+        fHeaderSize(p.fHeaderSize)\r\r\r
+      , fTrackletEnable(p.fTrackletEnable)\r\r\r
+      , fStackActive()\r\r\r
+      , fActiveStacks(p.fActiveStacks)\r\r\r
+      , fCorrupted(p.fCorrupted)\r\r\r
+      , fNexpectedHalfChambers(p.fNexpectedHalfChambers)\r\r\r
+      , fClean(p.fClean)\r\r\r
+      , fPos(p.fPos)\r\r\r
+      , fStacks()\r\r\r
+    {\r\r\r
+      // copy constructor\r\r\r
+    };\r\r\r
+\r\r\r
+    AliTRDrawSM &operator=(const AliTRDrawSM &)\r\r\r
+    {\r\r\r
+      // assignment operator\r\r\r
+      // not implemented\r\r\r
+      return *this;\r\r\r
+    };\r\r\r
+\r\r\r
+  };\r\r\r
+  \r\r\r
+  //--------------------------------------------------------\r\r\r
+     \r\r\r
+  AliTRDrawStream();\r\r\r
+  AliTRDrawStream(AliRawReader *rawReader);\r\r\r
+  virtual ~AliTRDrawStream();\r\r\r
+\r\r\r
+  //--------------------------------------------------------\r\r\r
+\r\r\r
+  virtual Bool_t       Next();                           // read the next data in the memory\r\r\r
+  //virtual Int_t      NextChamber(AliTRDdigitsManager *man); // read next chamber data in the momory\r\r\r
+  virtual Int_t        NextChamber(AliTRDdigitsManager *man, UInt_t **trackletContainer, UShort_t **/*errorCodeContainer*/=NULL); // read next chamber data in the memory\r\r\r
+  virtual Bool_t       Init();                           // initialize some internal variables\r\r\r
+\r\r\r
+  Int_t    NextBuffer(); // go and init next buffer if available - check the implementation file for return values\r\r\r
+\r\r\r
+  Bool_t   SetRawVersion(Int_t fraw); // set the raw version - used for backward compat.\r\r\r
+  \r\r\r
+  Bool_t   IsCurrentPadShared() const {return fADC->fIsShared;} // is current pad shared between mcms\r\r\r
+  void     SetSharedPadReadout(Bool_t fv) {fSharedPadsOn = fv;} // set the flag on if the reader should return the shared pads\r\r\r
+  \r\r\r
+  Bool_t   DecodeSM(void *buffer, UInt_t length); // decode a buffer\r\r\r
+  Int_t    DecodeSM();                            // used with raw reader\r\r\r
+  Int_t    DecodeSM(AliRawReader *reader);        // used with raw reader\r\r\r
+          \r\r\r
+  Bool_t   SetReader(AliRawReader *reader); // set the raw reader to use\r\r\r
+          \r\r\r
+  // info from Supermodule Index Word\r\r\r
+  Bool_t    IsTrackletEnableBitSet() const {return fSM.fTrackletEnable;} // get status of tracklet enable bit\r\r\r
+  Bool_t    IsStackActive(Int_t is) const {return fSM.fStackActive[is];} // get status of stack enable bit\r\r\r
+  Int_t     GetNofActiveStacks() const {return fSM.fActiveStacks;}       // get number of active stacks from stack mask\r\r\r
+  UInt_t   *GetGTUheaderWords() const {return fSM.fPos;}       // get number of active stacks from stack mask\r\r\r
+\r\r\r
+  // info from Stack Index Word\r\r\r
+  Int_t     GetNexpectedHalfChambers() const {return fSM.fNexpectedHalfChambers;}                    // get number of expected HC in a sm\r\r\r
+  Int_t     GetNofActiveLinksInStack(Int_t is) const {return fSM.fStacks[is].fActiveLinks;}          // get number of active links in a stack\r\r\r
+  Bool_t    IsLinkActiveInStack(Int_t is, Int_t il) const {return fSM.fStacks[is].fLinksActive[il];} // check whether the link is active\r\r\r
+\r\r\r
+  // info from Stack Header Word\r\r\r
+  Short_t    GetLinkMonitorError(Int_t is, Int_t il) const {return fSM.fStacks[is].fLinkMonitorError[il];} // get link monitor error\r\r\r
+\r\r\r
+  // info from Tracklet Data\r\r\r
+  Int_t     GetTrackletErrorCode(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fTrackletError;}\r\r\r
+  Int_t     GetNTracklets(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fNTracklets;} // get number of tracklets \r\r\r
+\r\r\r
+  // info from HC Header Word\r\r\r
+  Int_t     GetSM(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fSM;}\r\r\r
+  Int_t     GetLayer(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fLayer;}\r\r\r
+  Int_t     GetStack(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fStack;}\r\r\r
+  Int_t     GetSide(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fSide;}\r\r\r
+  Int_t     GetH0ErrorCode(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fH0Corrupted;}\r\r\r
+  Int_t     GetH1ErrorCode(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fH1Corrupted;}\r\r\r
+  Int_t     GetNumberOfTimeBins(Int_t is, Int_t il) const { return fSM.fStacks[is].fHalfChambers[il].fTimeBins;}\r\r\r
+  UInt_t   *GetTrackletWords(Int_t is, Int_t il) { return fSM.fStacks[is].fHalfChambers[il].fTrackletWords;}\r\r\r
+\r\r\r
+  // info from HC data\r\r\r
+  Int_t     GetHCErrorCode(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fCorrupted;}\r\r\r
+  Int_t     GetHCMCMmax(Int_t is, Int_t il) const {return fSM.fStacks[is].fHalfChambers[il].fMCMmax;}\r\r\r
+\r\r\r
+  // from MCM Header Word\r\r\r
+  // rob and mcm ordering\r\r\r
+  // side 0(even link) - ROB: 0 2 4 6  MCM: 12 13 14 15 8 9 10 11 4 5 6 7 0 1 2 3  \r\r\r
+  // side 1( odd link) - ROB: 1 3 5 7  MCM: 12 13 14 15 8 9 10 11 4 5 6 7 0 1 2 3  \r\r\r
+  Int_t     GetMCM(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fMCM;}\r\r\r
+  Int_t     GetROB(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fROB;}\r\r\r
+  Int_t     GetMCMhdErrorCode(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fMCMhdCorrupted;}\r\r\r
+  Int_t     GetMCMADCMaskErrorCode(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCmaskCorrupted;}\r\r\r
+  Int_t     GetEventNumber(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fEvCounter;}\r\r\r
+  Int_t     GetADCcount(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCcount;}\r\r\r
+\r\r\r
+  // info from MCM data words\r\r\r
+  Int_t     GetMCMErrorCode(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fCorrupted;} // get MCM data error code\r\r\r
+  Int_t     GetADCErrorCode(Int_t stack, Int_t link, Int_t mcm, Int_t adc) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fCorrupted;} // get ADC error code\r\r\r
+  Int_t     GetADCnumber(Int_t stack, Int_t link, Int_t mcm, Int_t adc) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fADCnumber;} // get ADC error code\r\r\r
+\r\r\r
+  Int_t     GetRow(Int_t stack, Int_t link, Int_t mcm) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fROW;}         // get current row number\r\r\r
+  Int_t     GetCol(Int_t stack, Int_t link, Int_t mcm, Int_t adc) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fCOL;}         // get current column number\r\r\r
+  Int_t     GetExtendedCol(Int_t stack, Int_t link, Int_t mcm, Int_t adc) const {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fExtendedCOL;}         // get current column number\r\r\r
+\r\r\r
+  // info from ADC data words\r\r\r
+  Int_t    *GetSignalDirect(Int_t stack, Int_t link, Int_t mcm, Int_t adc) {return fSM.fStacks[stack].fHalfChambers[link].fMCMs[mcm].fADCs[adc].fSignals;}\r\r\r
+\r\r\r
+\r\r\r
+  // from here, only works with returning ADC channel pointer using Next() \r\r\r
+  UInt_t   *GetTrackletWords() const { return fHC->fTrackletWords;}                 // return tracklet words pointer per hc [mj]\r\r\r
+  Int_t     GetTrackletErrorCode() const {return fHC ? fHC->fTrackletError : -1;}   // get tracklet error code\r\r\r
+  Int_t     GetNTracklets() const {return fHC ? fHC->fNTracklets : -1;}             // get number of tracklets \r\r\r
+\r\r\r
+  Int_t     GetSpecialRawVersion() const {return fHC ? fHC->fSpecialRawV : -1;}     // return special raw version\r\r\r
+  Int_t     GetMajorRawVersion() const {return fHC ? fHC->fRawVMajor : -1;}         // major raw version getter\r\r\r
+  Int_t     GetRawVersion() const {return fHC ? fHC->fRawVMajor : -1;}              // compatibility see funtion above\r\r\r
+  Int_t     GetMinorRawVersion() const {return fHC ? fHC->fRawVMinor : -1;}         // minor raw version\r\r\r
+\r\r\r
+  Int_t     GetSM() const {return fHC ? fHC->fSM : -1;}                              //  SM Position of CURRENT half chamber in full TRD\r\r\r
+  Int_t     GetLayer() const {return fHC ? fHC->fLayer : -1;}                        //  Layer Position of CURRENT half chamber in full TRD\r\r\r
+  Int_t     GetStack() const {return fHC ? fHC->fStack : -1;}                        //  Stack Position of CURRENT half chamber in full TRD\r\r\r
+  Int_t     GetSide() const {return fHC ? fHC->fSide : -1;}                          // get side\r\r\r
+  Int_t     GetDCS() const { return fHC ? fHC->fDCSboard : -1;}                      //  DCS board number read from data (HC header)\r\r\r
+  Int_t     GetROC() const { return fHC ? fHC->fROC : -1;}                           //  ROB Position of CURRENT half chamber in full TRD\r\r\r
+  Int_t     GetNumberOfTimeBins() const { return fHC ? fHC->fTimeBins : 0;}          // Get Ntime bins\r\r\r
+  UInt_t    GetBunchCrossCounter() const {return fHC ? fHC->fBunchCrossCounter : 0;} // get bunch cross counter\r\r\r
+  UInt_t    GetPreTriggerCounter() const {return fHC ? fHC->fPreTriggerCounter : 0;} // get pre trigger info\r\r\r
+  UInt_t    GetPreTriggerPhase() const {return fHC ? fHC->fPreTriggerPhase : 0;}     // get trigger phase\r\r\r
+\r\r\r
+  Int_t     GetRow() const {return fMCM ? fMCM->fROW : -1;}         // get current row number\r\r\r
+  Int_t     GetExtendedCol() const {return fADC ? fADC->fExtendedCOL : -1;}         // get current column number\r\r\r
+  Int_t     GetCol() const {return fADC ? fADC->fCOL : -1;}         // get current column number\r\r\r
+  Int_t     GetRowMax() const { return fHC ? fHC->fRowMax : -1;}    // Get maximum rows in the current HC\r\r\r
+  Int_t     GetColMax() const { return fHC ? fHC->fColMax : -1;}    // Get maximum cols in the current HC\r\r\r
+  // compatibility\r\r\r
+  Int_t     GetMaxRow() const { return fHC ? fHC->fRowMax : -1;}    // Get maximum rows in the current HC\r\r\r
+  Int_t     GetMaxCol() const { return fHC ? fHC->fColMax : -1;}    // Get maximum cols in the current HC\r\r\r
+\r\r\r
+  UInt_t    GetHCword0() const {return fHC ? *fHC->fPos[0] : 0;}    // get the HC word 0\r\r\r
+  UInt_t    GetHCword1() const {return fHC ? *fHC->fPos[1] : 0;}    // get the HC word 1\r\r\r
+           \r\r\r
+  Int_t     GetDET() const {return fHC ? fHC->fDET : -1;}           // get current det number\r\r\r
+  Int_t     GetDet() const {return fHC ? fHC->fDET : -1;}           // get current det number\r\r\r
+           \r\r\r
+  Int_t     GetROB() const {return fMCM ? fMCM->fROB : -1;}         // get current ROB number\r\r\r
+  Int_t     GetMCM() const {return fMCM ? fMCM->fMCM : -1;}         // get current MCM number\r\r\r
+  Int_t     GetEventNumber() const { return fMCM->fEvCounter;}      //  MCM Event number and position of current MCM on TRD chamber\r\r\r
+\r\r\r
+  Int_t     GetADC() const { return fADC ? fADC->fADCnumber : -1;}  //  MCM ADC channel and Time Bin of word 1\r\r\r
+  Int_t     GetTimeBin() const { return 0;}                         //  MCM ADC channel and Time Bin of word 1\r\r\r
+  Int_t    *GetSignals() const { return fADC ? fADC->fSignals : (Int_t *)fgEmptySignals;} // signals in the n-time bins from data word\r\r\r
+\r\r\r
+  Int_t     GetHCErrorCode() const {return fHC ? fHC->fCorrupted : -1;}    // get HC error code\r\r\r
+  Int_t     GetH0ErrorCode() const {return fHC ? fHC->fH0Corrupted : -1;}  // get HC header word0 error code\r\r\r
+  Int_t     GetH1ErrorCode() const {return fHC ? fHC->fH1Corrupted : -1;}  // get HC header word1 error code\r\r\r
+  Int_t     GetMCMErrorCode() const {return fMCM ? fMCM->fCorrupted : -1;} // get MCM data error code\r\r\r
+  Int_t     GetADCErrorCode() const {return fADC ? fADC->fCorrupted : -1;} // get ADC data error code\r\r\r
+  Int_t     GetMCMhdErrorCode() const {return fMCM ? fMCM->fMCMhdCorrupted: -1;} // get MCM header word error code\r\r\r
+  Int_t     GetMCMADCMaskErrorCode() const {return fMCM ? fMCM->fADCmaskCorrupted: -1;} // get MCM adc mask error code\r\r\r
+\r\r\r
+  UInt_t   *GetSMstreamPosition() const {return fSM.fPos;} // get position of the SM index word in the buffer\r\r\r
+\r\r\r
+  Bool_t    IsSMbufferClean() const {return fSM.fClean;}   // is data clean\r\r\r
+\r\r\r
+  //----------------------------------------------------------\r\r\r
\r\r\r
+  static void    SetNoDebug() {fgDebugFlag = kFALSE;} // allow debug info\r\r\r
+  static void    EnableMemoryReset() {fgEnableMemoryReset = kTRUE;} // allow memory reset\r\r\r
+  static void    SetNoErrorWarning() {fgWarnError = kFALSE;} // disable warning and error info\r\r\r
+  static void    SetForceCleanDataOnly() {fgCleanDataOnly = kTRUE;} // clean data only\r\r\r
+  static void    AllowCorruptedData() {fgCleanDataOnly = kFALSE;} // accept corrupted data\r\r\r
+\r\r\r
+  static void    SetExtraWordsFix() {fgExtraSkip = kTRUE;} // extra skip of 24 32-bit words \r\r\r
+  static void    SetSkipCDH() {fgSkipCDH = kTRUE;} // skip of 8 32-bit words \r\r\r
+  static void    SetDumpHead(Int_t iv) {fgDumpHead = iv;}\r\r\r
+  static void    DisableStackNumberChecker() {fgStackNumberChecker = kFALSE;}  // set false to cleanroom data \r\r\r
+  static void    DisableStackLinkNumberChecker() {fgStackLinkNumberChecker = kFALSE;}  \r\r\r
+  static void    DisableSkipData() {fgSkipData = kFALSE;} // keep reading next words even previous words were corrupted - debugging purpose  \r\r\r
+  static void    SetDumpingEnable() {fgDumpingEnable = kTRUE;} \r\r\r
+  static void    SetDumpingMCM(Int_t sm, Int_t stack, Int_t layer, Int_t rob, Int_t mcm) {fgDumpingSM = sm; fgDumpingStack = stack; fgDumpingLayer = layer; fgDumpingROB = rob; fgDumpingMCM = mcm;}\r\r\r
+\r\r\r
+  // this is a temporary solution!\r\r\r
+  // baseline should come with the HC header word 2 (count from 0!)\r\r\r
+\r\r\r
+  static void    EnableDecodeConfigData() {fgEnableDecodeConfigData = kTRUE;} // allow configuration data decoding\r\r\r
+  static Bool_t  fgEnableDecodeConfigData; // flag to decode configuration data\r\r\r
+\r\r\r
+\r\r\r
+  //--------------------------------------------------------\r\r\r
+  // Decoding functions\r\r\r
+  //--------------------------------------------------------\r\r\r
+\r\r\r
+  void DecodeSMInfo(const UInt_t *word, struct AliTRDrawSM *sm) const ;\r\r\r
+  const char *DumpSMInfo(const struct AliTRDrawSM *sm);\r\r\r
+  void DecodeStackInfo(const UInt_t *word, struct AliTRDrawStack *st) const;\r\r\r
+  void DecodeStackHeader(const UInt_t *word, struct AliTRDrawStack *st, Int_t iword) const;\r\r\r
+  const char *DumpStackInfo(const struct AliTRDrawStack *st);\r\r\r
+  Bool_t DecodeHCwordH0(const UInt_t *word, struct AliTRDrawHC *hc) const;\r\r\r
+  Bool_t DecodeHCwordH1(const UInt_t *word, struct AliTRDrawHC *hc) const;\r\r\r
+  const char *DumpHCinfoH0(const struct AliTRDrawHC *hc);\r\r\r
+  const char *DumpHCinfoH1(const struct AliTRDrawHC *hc);\r\r\r
+  void DecodeMCMheader(const UInt_t *word, struct AliTRDrawMCM *mcm) const;\r\r\r
+  UInt_t GetMCMadcMask(const UInt_t *word, struct AliTRDrawMCM *mcm) const;\r\r\r
+  void DecodeMask(const UInt_t *word, struct AliTRDrawMCM *mcm) const;\r\r\r
+  void MCMADCwordsWithTbins(UInt_t fTbins, struct AliTRDrawMCM *mcm) const;\r\r\r
+  const char *DumpMCMinfo(const struct AliTRDrawMCM *mcm);\r\r\r
+  const char *DumpMCMadcMask(const struct AliTRDrawMCM *mcm);\r\r\r
+\r\r\r
+\r\r\r
+ protected:\r\r\r
+\r\r\r
+  Bool_t InitBuffer(void *buffer, UInt_t length); // init the buffer - called by DecodeSM(void*, UInt_t)\r\r\r
+  Bool_t DumpWords(UInt_t *px, UInt_t iw, UInt_t marker = 0); // dump some words onto the screen\r\r\r
+\r\r\r
+  void   SwapOnEndian();         // swap if endian is BIG\r\r\r
+  Bool_t SkipWords(UInt_t iw);   // skip number of words\r\r\r
+  Bool_t DecodeGTUheader();      // decode data in GTU header\r\r\r
+  Bool_t DecodeTracklets();      // decode tracklets\r\r\r
+  Bool_t DecodeHC();             // decode data in HC\r\r\r
+\r\r\r
+  Bool_t DecodeADC();            // decode 10 ADC words\r\r\r
+  Bool_t DecodeADCExtended();    // decode TRAP data with ZS and more than 32 samples \r\r\r
+\r\r\r
+  Bool_t DecodeHCheader();       // decode HC  header\r\r\r
+  Bool_t SeekEndOfData();        // go to next end of raw data marker (actually 1 word after)\r\r\r
+  Bool_t SkipMCMdata(UInt_t iw);  // skip this mcm data due to mcm header corruption\r\r\r
+  Bool_t SeekNextMCMheader();    // go to next mcm header\r\r\r
+  Bool_t DecodeMCMheader();      // decode mcm header\r\r\r
+\r\r\r
+  Bool_t IsRowValid();       // check if row within the range\r\r\r
+  Bool_t IsHCheaderOK();     // check if current hc header data make sense\r\r\r
+  Bool_t IsMCMheaderOK();    // check if current mcm header data make sense\r\r\r
+    \r\r\r
+  void   ResetCounters();    // reset some counters\r\r\r
+  void   ResetIterators();   // needed for Next()\r\r\r
+  void   ResetPerSM();       // reset every SM \r\r\r
+  void   ResetPerStack();    // reset every Stack \r\r\r
+  void   ResetPerHC();       // reset every HC \r\r\r
+  void   ResetPerMCM();      // reset every MCM  \r\r\r
+  void   ResetPerADC();      // reset every ADC  \r\r\r
+  void   ResetMemory();      // reset all data members\r\r\r
+\r\r\r
+  AliTRDrawStream(const AliTRDrawStream& st);\r\r\r
+  AliTRDrawStream &operator=(const AliTRDrawStream &);\r\r\r
+\r\r\r
+  // ----------------- DATA MEMBERS START\r\r\r
+\r\r\r
+  struct AliTRDrawSM       fSM;    // one SM per buffer\r\r\r
+  struct AliTRDrawStack   *fStack; //! pointer to the current stack\r\r\r
+  struct AliTRDrawHC      *fHC;    //! current HC\r\r\r
+  struct AliTRDrawHC      *fLastHC; //! last HC\r\r\r
+  struct AliTRDrawMCM     *fMCM;   //! current MCM\r\r\r
+  struct AliTRDrawADC     *fADC;   //! current ADC\r\r\r
+  \r\r\r
+  UInt_t *fpPos;   // current position in the buffer\r\r\r
+  UInt_t *fpBegin; // begin - pointer to the buffer word 0\r\r\r
+  UInt_t *fpEnd;   // end of the buffer\r\r\r
+\r\r\r
+  UInt_t  fWordLength;  // length of the buffer in 32bit words\r\r\r
+\r\r\r
+  Bool_t  fIsGlobalDigitsParamSet; // set true if the global digit parameters are already set\r\r\r
+\r\r\r
+  Int_t   fStackNumber;     // current stack number\r\r\r
+  Int_t   fStackLinkNumber; // current link in the stack\r\r\r
+\r\r\r
+  Int_t   fhcMCMcounter;    // mcm counter inside single hc - used in Next()  \r\r\r
+  Int_t   fmcmADCcounter;   // adc counrer inside single adc - used in Next() \r\r\r
+\r\r\r
+  Int_t   fLinkTrackletCounter; // count the tracklets in the current HC\r\r\r
+  Int_t   fEndOfTrackletCount;  // count link by link (hc by hc) used for debug\r\r\r
+  Int_t   fNWordsCounter;       // counts words of given hc having link monitor error\r\r\r
+\r\r\r
+  UInt_t  fMaskADCword; // temp mask when decoding adcs\r\r\r
+  UInt_t  fTbinADC;     // temp adc \r\r\r
+  Int_t   fDecodedADCs; // counter of decoded adcs [mj] do we need?\r\r\r
+\r\r\r
+  UInt_t  fEventCounter;     // stores the valid/current MCM event counter\r\r\r
+  UInt_t  fLastEventCounter; // last known event counter of MCM\r\r\r
+\r\r\r
+  Bool_t  fSharedPadsOn; // do we want to output shared pads - default is off\r\r\r
+  Int_t   fMaxADCgeom;   // maximum ADC channels per mcm\r\r\r
+\r\r\r
+  Bool_t  fBufferRead; // set to be true if buffer is read for a given equipment\r\r\r
+\r\r\r
+  AliTRDgeometry *fGeometry;  //! TRD geometry\r\r\r
+  AliRawReader   *fRawReader; //! raw reader    \r\r\r
+\r\r\r
+  AliTRDfeeParam      *fTRDfeeParam; // pointer to the fee params\r\r\r
+\r\r\r
+  // STATIC \r\r\r
+\r\r\r
+  static Bool_t fgExtraSkip; // whether we should skip the leading 24 words\r\r\r
+  static Bool_t fgSkipCDH; // whether we should skip CDH (8 words)\r\r\r
+  static Bool_t fgWarnError; // no errors no warnings\r\r\r
+  static Bool_t fgCleanDataOnly; // release only clean events = no errors\r\r\r
+  static Bool_t fgDebugFlag; // allow debugging info\r\r\r
+  static Bool_t fgEnableMemoryReset; // allow memory reset\r\r\r
+  static Bool_t fgStackNumberChecker; // decide if we check stack number insanity - set false to cleanroom data\r\r\r
+  static Bool_t fgStackLinkNumberChecker; // decide if we check stack link number insanity - debuging purpose\r\r\r
+  static Bool_t fgSkipData; // decide if we skip corrupted data of given HC\r\r\r
+  static Int_t fgDumpHead; // number of words to dump (from the start of the buffer) on each Init\r\r\r
+  static Int_t  fgEmptySignals[30]; // empty signals in case of ADC pointer = NULL\r\r\r
+  static Short_t  fgMCMordering[16]; // mcm number odering for mcm header corruption check\r\r\r
+  static Short_t  fgROBordering[16]; // mcm number odering for mcm header corruption check\r\r\r
+  static Int_t  fgLastHC;  // HC number indicator to check HC numbering\r\r\r
+  static Int_t  fgLastROB; // ROB number indicator to check ROB numbering \r\r\r
+  static Int_t  fgLastIndex; // variable to check mcm number odering\r\r\r
+\r\r\r
+  static Bool_t fgDumpingEnable; // flag to dump MCM data words  \r\r\r
+\r\r\r
+  static Int_t  fgDumpingSM; // sm number to dump MCM data words\r\r\r
+  static Int_t  fgDumpingStack; // stack number to dump MCM data words\r\r\r
+  static Int_t  fgDumpingLayer; // layer number to dump MCM data words\r\r\r
+  static Int_t  fgDumpingROB; // rob number to dump MCM data words\r\r\r
+  static Int_t  fgDumpingMCM; // mcm number to dump MCM data words\r\r\r
+\r\r\r
+  // ----------------- DATA MEMBERS STOP\r\r\r
+\r\r\r
+  enum ETRDzRawStreamError \r\r\r
+    {\r\r\r
+      kDecodeStackInfo          = 1 //\r\r\r
+      , kMissingData             = 2 //\r\r\r
+      , kLinkDataMissing        = 3 //\r\r\r
+      , kHCdataMissing         = 4 //\r\r\r
+      , kTrackletOverflow       = 5 //\r\r\r
+      , kEOTrackeltsMissing     = 6 //\r\r\r
+      , kWrongPadrow             = 7 //\r\r\r
+      , kMCMheaderCorrupted     = 8 //\r\r\r
+      , kWrongMCMorROB         = 9 //\r\r\r
+      , kMCMeventMissmatch      = 10 //\r\r\r
+      , kMCMADCMaskMissing      = 11 //\r\r\r
+      , kHCHeaderCorrupt        = 12 //\r\r\r
+      , kHCHeaderWrongStack     = 13 //\r\r\r
+      , kHCHeaderWrongLayer     = 14 //\r\r\r
+      , kHCHeaderWrongSide      = 15 //\r\r\r
+      , kHCHeaderWrongSM        = 16 //\r\r\r
+      , kHCHeaderWrongDet       = 17 //\r\r\r
+      , kHCHeaderWrongROC       = 18 //\r\r\r
+      , kHCWordMissing         = 19 //\r\r\r
+      , kMCMdataMissing                = 20 //\r\r\r
+      , kMCMoverflow             = 21 //\r\r\r
+      , kADCdataMissing                = 22 //\r\r\r
+      , kADCmaskMissmatch       = 23 //\r\r\r
+      , kWrongPadcolumn                = 24 //\r\r\r
+    };                        \r\r\r
+\r\r\r
+  ClassDef(AliTRDrawStream, 0)\r\r\r
+}; //clas def end\r\r\r
+\r\r\r
+#endif\r\r\r