-/**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * *
- * Author: The ALICE Off-line Project. *
- * Contributors are mentioned in the code where appropriate. *
- * *
- * Permission to use, copy, modify and distribute this software and its *
- * documentation strictly for non-commercial purposes is hereby granted *
- * without fee, provided that the above copyright notice appears in all *
- * copies and that both the copyright notice and this permission notice *
- * appear in the supporting documentation. The authors make no claims *
- * about the suitability of this software for any purpose. It is *
- * provided "as is" without express or implied warranty. *
- **************************************************************************/
+#/**************************************************************************
+* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+* *
+* Author: The ALICE Off-line Project. *
+* Contributors are mentioned in the code where appropriate. *
+* *
+* Permission to use, copy, modify and distribute this software and its *
+* documentation strictly for non-commercial purposes is hereby granted *
+* without fee, provided that the above copyright notice appears in all *
+* copies and that both the copyright notice and this permission notice *
+* appear in the supporting documentation. The authors make no claims *
+* about the suitability of this software for any purpose. It is *
+* provided "as is" without express or implied warranty. *
+**************************************************************************/
/* $Id$ */
#include "AliTRDgeometry.h"
#include "AliTRDfeeParam.h"
#include "AliTRDdigitsManager.h"
+#include "AliTRDdataArrayS.h"
#include "AliTRDdataArrayI.h"
+#include "AliTRDdataArrayDigits.h"
#include "AliTRDSignalIndex.h"
+#include "AliTRDrecoParam.h"
+#include "AliTRDcalibDB.h"
+#include "Cal/AliTRDCalPadStatus.h"
+
+//[mj tracklet writing] #include "AliTRDrawTracklet.h"
#include "AliLog.h"
#include "AliRawReader.h"
// Stack word masks
#define STACK_HEADER_SIZE(w) GET_VALUE_AT(w,0xffff,16)
#define STACK_LINK_WORD(w) ((w) & 0xfff)
+#define LINK0_DATA_TYPE_FLAG(w) (GET_VALUE_AT(w,0x3,4) == (0x0) ? 0 : 1) // 0 if physics data
+#define LINK1_DATA_TYPE_FLAG(w) (GET_VALUE_AT(w,0x3,20) == (0x0) ? 0 : 1) // 0 if physics data
+#define LINK0_MONITOR_FLAG(w) (GET_VALUE_AT(w,0xf,0) == (0x0) ? 0 : 1) // 0 if OK
+#define LINK1_MONITOR_FLAG(w) (GET_VALUE_AT(w,0xf,16) == (0x0) ? 0 : 1) // 0 if OK
// HC word masks
//#define HC_HEADER_MASK_ERR(w) ( ((w) & (0x80000003)) == (0x80000001) ? 0 : 1) // 0 if OK!!!
// HC word 0
#define HC_SPECIAL_RAW_VERSION(w) IS_BIT_SET(w,31)
#define HC_MAJOR_RAW_VERSION(w) GET_VALUE_AT(w,0x7f,24)
+#define HC_MAJOR_RAW_VERSION_OPT(w) GET_VALUE_AT(w,0x7,24)
#define HC_MINOR_RAW_VERSION(w) GET_VALUE_AT(w,0x7f,17)
#define HC_EXTRA_WORDS(w) GET_VALUE_AT(w,0x7,14)
#define HC_DCS_BOARD(w) GET_VALUE_AT(w,0xfff<<20,20)
#define HC_PRETRIGGER_COUNTER(w) GET_VALUE_AT(w,0xf,6)
#define HC_PRETRIGGER_PHASE(w) GET_VALUE_AT(w,0xf,6)
-//--------------------------------------------------------
-
// MCM word and ADC mask
#define MCM_HEADER_MASK_ERR(w) ( ((w) & (0xf)) == (0xc) ? 0 : 1) // 0 if OK!!!
#define MCM_ADCMASK_MASK_ERR(w) ( ((w) & (0xf)) == (0xc) ? 0 : 1) // 0 if OK!!!
#define MCM_ROB_NUMBER(w) GET_VALUE_AT(w,0x7,28)
#define MCM_EVENT_COUNTER(w) GET_VALUE_AT(w,0x00fffff,4)
#define MCM_ADCMASK_VAL(w) GET_VALUE_AT(w,0x1fffff,4)
+#define MCM_ADCMASK_NADC(w) GET_VALUE_AT(w,0x1f,25)
-//this should give 0x1fffff as a mask
-#define MCM_DUMMY_ADCMASK_VAL 0x001fffffc
-//--------------------------------------------------------
-
-#define MAX_TRACKLETS_PERHC 256 // max number of tracklets per HC - large number for now
+#define MCM_DUMMY_ADCMASK_VAL 0x015fffffc // updated
+#define ADCDATA_VAL1 0x2 // updated
+#define ADCDATA_VAL2 0x3 // updated
//--------------------------------------------------------
#define ADC_WORD_MASK(w) ((w) & 0x3)
//--------------------------------------------------------
+
+
ClassImp(AliTRDrawStreamTB)
Bool_t AliTRDrawStreamTB::fgExtraSkip = kFALSE;
Bool_t AliTRDrawStreamTB::fgSkipCDH = kFALSE;
Bool_t AliTRDrawStreamTB::fgWarnError = kTRUE;
-Bool_t AliTRDrawStreamTB::fgCleanDataOnly = kTRUE;
+Bool_t AliTRDrawStreamTB::fgCleanDataOnly = kFALSE;
Bool_t AliTRDrawStreamTB::fgDebugFlag = kTRUE;
-Bool_t AliTRDrawStreamTB::fgDebugStreamFlag = kFALSE;
-TTreeSRedirector *AliTRDrawStreamTB::fgDebugStreamer = 0;
-UInt_t AliTRDrawStreamTB::fgStreamEventCounter = 0;
+Bool_t AliTRDrawStreamTB::fgEnableMemoryReset = kTRUE;
+Bool_t AliTRDrawStreamTB::fgStackNumberChecker = kTRUE;
+Bool_t AliTRDrawStreamTB::fgStackLinkNumberChecker = kTRUE;
+Bool_t AliTRDrawStreamTB::fgSkipData = kTRUE;
UInt_t AliTRDrawStreamTB::fgDumpHead = 0;
-
+Int_t AliTRDrawStreamTB::fgCommonAdditive = 0;
Int_t AliTRDrawStreamTB::fgEmptySignals[] =
{
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1
};
+Short_t AliTRDrawStreamTB::fgMCMordering[] =
+ {
+ 12, 13, 14, 15, 8, 9, 10, 11, 4, 5, 6, 7, 0, 1, 2, 3
+ };
+Short_t AliTRDrawStreamTB::fgROBordering[] =
+ {
+ 0, 1, 2, 3
+ };
+Int_t AliTRDrawStreamTB::fgLastHC = -1;
+Int_t AliTRDrawStreamTB::fgLastROB = -1;
+Int_t AliTRDrawStreamTB::fgLastIndex = -1;
+Bool_t AliTRDrawStreamTB::fDumpingEnable = kFALSE;
+Int_t AliTRDrawStreamTB::fDumpingSM = -1;
+Int_t AliTRDrawStreamTB::fDumpingStack = -1;
+Int_t AliTRDrawStreamTB::fDumpingLayer = -1;
+Int_t AliTRDrawStreamTB::fDumpingROB = -1;
+Int_t AliTRDrawStreamTB::fDumpingMCM = -1;
+
AliTRDrawStreamTB::AliTRDrawStreamTB()
- : TObject()
+ : AliTRDrawStreamBase()
, fSM()
, fStack(0)
, fHC(0)
, fmcmADCcounter(0)
, fLinkTrackletCounter(-1)
, fEndOfTrackletCount(-1)
+ , fNWordsCounter(-1)
, fMaskADCword(0)
, fTbinADC(0)
, fDecodedADCs(-1)
, fEventCounter(0)
, fLastEventCounter(0)
- , fCommonAdditive(0)
, fSharedPadsOn(kFALSE)
, fMaxADCgeom(0)
+ , fBufferRead(0)
, fGeometry(0)
, fRawReader(0)
, fTRDfeeParam(0)
- , fDebugStreamOwned(kFALSE)
{
//
// default constructor
//--------------------------------------------------------
AliTRDrawStreamTB::AliTRDrawStreamTB(AliRawReader *rawReader)
- : TObject()
+ : AliTRDrawStreamBase(rawReader)
, fSM()
, fStack(0)
, fHC(0)
, fmcmADCcounter(0)
, fLinkTrackletCounter(-1)
, fEndOfTrackletCount(-1)
+ , fNWordsCounter(-1)
, fMaskADCword(0)
, fTbinADC(0)
, fDecodedADCs(-1)
, fEventCounter(0)
, fLastEventCounter(0)
- , fCommonAdditive(0)
, fSharedPadsOn(kFALSE)
, fMaxADCgeom(0)
+ , fBufferRead(0)
, fGeometry(0)
, fRawReader(rawReader)
, fTRDfeeParam(0)
- , fDebugStreamOwned(kFALSE)
{
//
// default constructor
if (fRawReader)
{
if (Init() == kFALSE)
- {
- AliWarning("Unable to Init. Try setting up the reader with SetReader or buffer with Init(void *, UInt_t )");
- }
+ {
+ AliWarning("Unable to Init. Try setting up the reader with SetReader or buffer with Init(void *, UInt_t )");
+ }
}
else
{
//------------------------------------------------------------
AliTRDrawStreamTB::AliTRDrawStreamTB(const AliTRDrawStreamTB& /*st*/)
- : TObject()
+ : AliTRDrawStreamBase()
, fSM()
, fStack(0)
, fHC(0)
, fmcmADCcounter(0)
, fLinkTrackletCounter(-1)
, fEndOfTrackletCount(-1)
+ , fNWordsCounter(-1)
, fMaskADCword(0)
, fTbinADC(0)
, fDecodedADCs(-1)
, fEventCounter(0)
, fLastEventCounter(0)
- , fCommonAdditive(0)
, fSharedPadsOn(kFALSE)
, fMaxADCgeom(0)
+ , fBufferRead(0)
, fGeometry(0)
, fRawReader(0)
, fTRDfeeParam(0)
- , fDebugStreamOwned(kFALSE)
{
//
// Copy constructor
}
//------------------------------------------------------------
-void AliTRDrawStreamTB::SetRawVersion(Int_t fraw)
+Bool_t AliTRDrawStreamTB::SetRawVersion(Int_t fraw)
{
//
// function provided for backward compatibility
//
AliWarning("Raw data version is read from raw data stream! No point of setting it in here.");
fraw = 0; // avoid warnings
+ return kFALSE;
}
//------------------------------------------------------------
// destructor
//
delete fGeometry;
-
- if (fDebugStreamOwned == kTRUE)
- {
- if (fgDebugStreamer)
- {
- delete fgDebugStreamer;
- fgDebugStreamer = 0;
- }
- }
}
//------------------------------------------------------------
return *this;
}
-//------------------------------------------------------------
-
-void
-AliTRDrawStreamTB::EnableDebug(TTreeSRedirector *debugStream)
-{
- //
- // replace the current debug streamer or create a new one if owned
- //
-
- if (debugStream == 0)
- {
- if (fDebugStreamOwned == kTRUE)
- {
- if (fgDebugStreamer)
- delete fgDebugStreamer;
- }
-
- fDebugStreamOwned = kTRUE;
- fgDebugStreamer = new TTreeSRedirector("TRDrawDataDebug.root");
- fgStreamEventCounter = 0;
-
- }
- else
- {
- if (fDebugStreamOwned == kTRUE)
- {
- if (fgDebugStreamer)
- delete fgDebugStreamer;
- }
-
- fgDebugStreamer = debugStream;
- fDebugStreamOwned = kFALSE;
- }
-}
-
//___________________________________________________________
void
AliTRDrawStreamTB::SwapOnEndian()
//
// Check the endian and swap if needed
//
-
int itemp = 1;
char* ptemp = (char*) &itemp;
if (ptemp[0] != 1)
{
- // need to swap...
if (fgDebugFlag) AliDebug(8, "Swapping.");
fpPos = fpBegin;
UInt_t iutmp = 0;
while (fpPos < fpEnd)
- {
- fpPos += 1;
- iutmp = (((*fpPos & 0x000000ffU) << 24) | ((*fpPos & 0x0000ff00U) << 8) |
- ((*fpPos & 0x00ff0000U) >> 8) | ((*fpPos & 0xff000000U) >> 24));
- // here we override the value in the buffer!
- *fpPos = iutmp;
- }
+ {
+ fpPos += 1;
+ iutmp = (((*fpPos & 0x000000ffU) << 24) | ((*fpPos & 0x0000ff00U) << 8) |
+ ((*fpPos & 0x00ff0000U) >> 8) | ((*fpPos & 0xff000000U) >> 24));
+ // here we override the value in the buffer!
+ *fpPos = iutmp;
+ }
fpPos = fpBegin;
}
}
-//------------------------------------------------------------
-void
-AliTRDrawStreamTB::DumpErrorCount()
-{
- //
- // print the error statistics into the stdout
- //
-
- ;
- // AliInfo(Form("Error counts: HC0 %d HC1 %d MCM %d ADCmask %d ADC %d",
- // fHC->fH0ErrorCounter, fHC->fH1ErrorCounter, fMCM->fErrorCounter, fMCM->fMaskErrorCounter, adc.fErrorCounter));
-}
-
//------------------------------------------------------------
Bool_t
AliTRDrawStreamTB::DumpWords(UInt_t *px, UInt_t iw, UInt_t marker)
{
- //
- // skip few words
- // note: can be made faster
- //
-
+
TString tsreturn = Form("\n[ Dump Sequence at 0x%08x ] : ", px);
for (UInt_t i = 0; i < iw; i++)
{
if (px + iw >= fpEnd)
- return kFALSE;
+ return kFALSE;
if (i % 8 == 0)
- tsreturn += "\n ";
+ tsreturn += "\n ";
if (marker != 0 && marker == px[i])
- tsreturn += Form(" *>0x%08x<* ", px[i]);
+ tsreturn += Form(" *>0x%08x<* ", px[i]);
else
- tsreturn += Form("0x%08x ", px[i]);
+ tsreturn += Form("0x%08x ", px[i]);
}
tsreturn += "\n";
AliTRDrawStreamTB::SkipWords(UInt_t iw)
{
//
- // skip few words
- // note: can be made faster
+ // Skip words corresponding to iw
//
-
if ( fpPos + iw < fpEnd )
{
fpPos += iw;
else
{
if (fgWarnError) AliWarning(Form("Skip %d words failed. %d available", iw, fpEnd - fpPos - 1));
- if (fgDebugStreamer)
- {
- TTreeSRedirector &cstream = *fgDebugStreamer;
- cstream << "SkipWords"
- << "Event=" << fgStreamEventCounter
- << ".nwords=" << iw
-// << ".read=" << (Char_t*)(fpPos - fpBegin)
-// << ".togo=" << (Char_t*)(fpEnd - fpPos)
-// << ".length=" << (Char_t*)(fpEnd - fpBegin)
- << "\n";
- }
return kFALSE;
}
Bool_t
AliTRDrawStreamTB::SetReader(AliRawReader *reader)
{
- //
- //
- //
if (reader != 0)
{
fRawReader = reader;
if (fRawReader)
- {
- return Init();
- }
+ {
+ return Init();
+ }
else
- {
- AliWarning("Unable to setup reader.");
- return kFALSE;
- }
+ {
+ AliWarning("Unable to setup reader.");
+ return kFALSE;
+ }
}
else
{
{
//
// return -1 if no more buffers available
- // return 0 DecodeSM failed (clean data required for example) but still maybe more data to come
- // return 1 DecodeSM OK
+ // return 0 DecodeSM failed (clean data required for example) but still maybe more data to come
+ // return 1 DecodeSM OK
//
if (fRawReader != 0)
UInt_t length = 0;
Bool_t kBufferSet = fRawReader->ReadNextData(buffer);
if (kBufferSet == kTRUE)
- {
- if (fgDebugFlag) AliDebug(9, "Buffer is set.");
- length = fRawReader->GetDataSize();
- if (fgExtraSkip == kTRUE)
- {
- buffer += EXTRA_LEAD_WORDS * WORD_SIZE;
- length -= EXTRA_LEAD_WORDS * WORD_SIZE;
- }
-
- if (fgSkipCDH == kTRUE)
- {
- buffer += CDH_WORDS * WORD_SIZE;
- length -= CDH_WORDS * WORD_SIZE;
- }
-
- if (length > 0)
- {
- if (fgDebugFlag) AliDebug(9, Form("Buffer length : %d", length));
- //return InitBuffer((void*)buffer, length);
- if (DecodeSM((void*)buffer, length) == kTRUE)
- return 1;
- else
- return 0;
- }
- }
+ {
+ if (fgDebugFlag) AliDebug(9, "Buffer is set.");
+ length = fRawReader->GetDataSize();
+ if (fgExtraSkip == kTRUE)
+ {
+ buffer += EXTRA_LEAD_WORDS * WORD_SIZE;
+ length -= EXTRA_LEAD_WORDS * WORD_SIZE;
+ }
+
+ if (fgSkipCDH == kTRUE)
+ {
+ buffer += CDH_WORDS * WORD_SIZE;
+ length -= CDH_WORDS * WORD_SIZE;
+ }
+
+ if (length > 0)
+ {
+ if (fgDebugFlag) AliDebug(9, Form("Buffer length : %d", length));
+ if (fgEnableMemoryReset) ResetMemory(); //[mj]
+ if (DecodeSM((void*)buffer, length) == kTRUE)
+ return 1;
+ else
+ return 0;
+ }
+ }
else
- {
- return -1;
- }
+ {
+ return -1;
+ }
}
return -1;
//
// reset some global counters
//
+ fBufferRead = kFALSE; // important to read buffer
fStackNumber = 0;
fStackLinkNumber = 0;
-
fDecodedADCs = 0;
fSM.fActiveStacks = 0;
fSM.fNexpectedHalfChambers = 0;
fLinkTrackletCounter = 0;
-
fLastEventCounter = 0;
fEventCounter = 0;
}
AliTRDrawStreamTB::ResetIterators()
{
//
- // reset the data iterators used in the Next()
+ // reset data which should be reset every sm
//
- fStackNumber = 0;
- fStackLinkNumber = 0;
+ fStackNumber = 0; // reset for Next() function
+ fStackLinkNumber = 0; // reset for Next() function
fhcMCMcounter = 0;
fmcmADCcounter = 0;
}
+//------------------------------------------------------------
+void
+AliTRDrawStreamTB::ResetPerSM()
+{
+ //
+ // reset every SM
+ //
+
+ fSM.fHeaderSize = 0;
+ fSM.fTrackletEnable = kFALSE;
+ fSM.fCorrupted = 0;
+ fSM.fNexpectedHalfChambers = 0;
+ fSM.fNexpectedHalfChambers = 0;
+ fSM.fClean = kTRUE;
+ fSM.fPos = NULL;
+ for (Int_t i=0; i<5; i++){
+ fSM.fStackActive[i] = kFALSE;
+ }
+}
+
+//------------------------------------------------------------
+void
+AliTRDrawStreamTB::ResetPerStack()
+{
+ //
+ // reset every Stack
+ //
+
+ fStack->fHeaderSize = 0;
+ fStack->fActiveLinks = 0;
+ fStack->fPos = NULL;
+ for (Int_t i=0; i<12; i++){
+ fStack->fLinksActive[i] = kFALSE;
+ fStack->fLinksDataType[i] = 0;
+ fStack->fLinksMonitor[i] = 0;
+ fStack->fLinkMonitorError[i] = 0;
+ }
+}
+
+//------------------------------------------------------------
+void
+AliTRDrawStreamTB::ResetPerHC()
+{
+ //
+ // reset every HC
+ //
+ fEventCounter = 0;
+ fHC->fTrackletError = 0;
+ fHC->fNTracklets = 0;
+ fHC->fSpecialRawV = 0;
+ fHC->fRawVMajor = 0;
+ fHC->fRawVMajorOpt = 0;
+ fHC->fRawVMinor = 0;
+ fHC->fNExtraWords = 0;
+ fHC->fDCSboard = 0;
+ fHC->fSM = 0;
+ fHC->fStack = 0;
+ fHC->fLayer = 0;
+ fHC->fSide = 0;
+ fHC->fTimeBins = 0;
+ fHC->fBunchCrossCounter = 0;
+ fHC->fPreTriggerCounter = 0;
+ fHC->fPreTriggerPhase = 0;
+ fHC->fDET = 0;
+ fHC->fROC = 0;
+ fHC->fRowMax = 0;
+ fHC->fColMax = 0;
+ fHC->fMCMmax = 0;
+
+ fHC->fH0Corrupted = 0;
+ fHC->fH1Corrupted = 0;
+ fHC->fCorrupted = 0;
+}
+
+//------------------------------------------------------------
+void
+AliTRDrawStreamTB::ResetPerMCM()
+{
+ //
+ // reset every MCM
+ //
+ fMCM->fROB = 0;
+ fMCM->fMCM = 0;
+ fMCM->fROW = 0;
+ fMCM->fEvCounter = 0;
+ fMCM->fADCMask = 0;
+ fMCM->fADCMaskWord = 0;
+ fMCM->fADCmax = 0;
+ fMCM->fADCcount = 0;
+ fMCM->fMCMADCWords = 0;
+ fMCM->fSingleADCwords = 0;
+ fMCM->fMCMhdCorrupted = 0;
+ fMCM->fADCmaskCorrupted = 0;
+ fMCM->fCorrupted = 0;
+ fMCM->fPos = NULL;
+ fMCM->fAdcDataPos = NULL;
+ fMCM->fADCcounter = 0;
+
+ for (Int_t i=0; i<21; i++){
+ fMCM->fADCchannel[i] = 0;
+ }
+}
+
+//------------------------------------------------------------
+void
+AliTRDrawStreamTB::ResetPerADC()
+{
+ //
+ // reset every ADC
+ //
+ fADC->fPos = NULL;
+ fADC->fADCnumber = 0;
+ fADC->fCOL = 0;
+ fADC->fIsShared = kTRUE;
+ fADC->fCorrupted = 0;
+
+ for (Int_t i=0; i<30; i++){
+ fADC->fSignals[i] = 0;
+ }
+}
+
+//------------------------------------------------------------
+void
+AliTRDrawStreamTB::ResetMemory()
+{
+ //
+ // initialize all the data members to prevent read data
+ // from previous buffer
+ //
+ ResetPerSM();
+ for (Int_t istack=0; istack<5; istack++){
+ fStack = &fSM.fStacks[istack];
+ ResetPerStack();
+ for (Int_t ilink=0; ilink<12; ilink++){
+ fHC = &fStack->fHalfChambers[ilink];
+ ResetPerHC();
+ for (Int_t imcm=0; imcm<12; imcm++){
+ fMCM = &fHC->fMCMs[imcm];
+ ResetPerMCM();
+ for (Int_t iadc=0; iadc<12; iadc++){
+ fADC = &fMCM->fADCs[iadc];
+ ResetPerADC();
+ }
+ }
+ }
+ }
+}
+
//------------------------------------------------------------
Bool_t
// returns with true on next adc read
// returns false on errors and end of buffer
//
+if (fBufferRead)
+ {
- while (fStackNumber < 5 && fSM.fActiveStacks > 0)
- {
- if(fSM.fStackActive[fStackNumber] == kTRUE)
- {
- fStack = &fSM.fStacks[fStackNumber];
- while (fStackLinkNumber < 12)
- {
- if (fStack->fLinksActive[fStackLinkNumber] == kTRUE)
- {
- fHC = &fStack->fHalfChambers[fStackLinkNumber];
- if (!fHC)
- {
- AliError(Form("Super Strange. HC missing at stack %d link %d", fStackNumber, fStackLinkNumber));
- return kFALSE;
- }
- //if (fgDebugFlag) AliDebug(2, DumpHCinfoH0(fHC));
- if (fHC->fCorrupted == 0)
- {
- while (fhcMCMcounter < fHC->fMCMmax)
- {
- fMCM = &fHC->fMCMs[fhcMCMcounter];
- //if (fgDebugFlag) AliDebug(2, DumpMCMinfo(fMCM));
- if (!fMCM)
- {
- AliError(Form("Super Strange. HC missing at stack %d link %d atMCMslot %d",
- fStackNumber, fStackLinkNumber, fhcMCMcounter));
- return kFALSE;
- }
- while(fmcmADCcounter < fMCM->fADCmax)
- {
- //printf("%d %d %d %d\n", fStackNumber, fStackLinkNumber, fhcMCMcounter, fmcmADCcounter);
- fADC = &fMCM->fADCs[fmcmADCcounter];
- if (!fADC)
- {
- AliError(Form("Super Strange. ADC missing at stack %d link %d MCMslot %d ADCslot %d",
- fStackNumber, fStackLinkNumber, fhcMCMcounter, fmcmADCcounter));
- return kFALSE;
- }
- fmcmADCcounter++;
- //printf("ADC count %d \n", fmcmADCcounter);
- if (fSharedPadsOn)
- {
- return kTRUE;
- }
- else
- {
- if (fADC->fIsShared == kFALSE)
- return kTRUE;
- }
- } //while adc in MCM
- //printf("MCM count %d \n", fhcMCMcounter);
- fhcMCMcounter++;
- // next MCM should go through all active ADCs
- fmcmADCcounter = 0;
- } // while mcm
- } // if HC OK
- }// if link active
- fStackLinkNumber++;
- // next stack link (HC) should go through all active MCMs
- fhcMCMcounter = 0;
- }// while links
- }// if stack active
- fStackNumber++;
- // next stack should go through all links - start from 0
- fStackLinkNumber = 0;
- }
-
- // in case rawreader manages the mem buffers
- // lets go for the next buffer
+ while (fStackNumber < 5 && fSM.fActiveStacks > 0)
+ {
+ if(fSM.fStackActive[fStackNumber] == kTRUE)
+ {
+ fStack = &fSM.fStacks[fStackNumber];
+ while (fStackLinkNumber < 12)
+ {
+ if (fStack->fLinksActive[fStackLinkNumber] == kTRUE && fStack->fLinksMonitor[fStackLinkNumber] == 0)
+ {
+ fHC = &fStack->fHalfChambers[fStackLinkNumber];
+ if (!fHC)
+ {
+ AliError(Form("HC missing at stack %d link %d", fStackNumber, fStackLinkNumber));
+ return kFALSE;
+ }
+ if (fHC->fCorrupted == 0 && (fHC->fH0Corrupted == 0 && fHC->fH1Corrupted == 0)) // if HC data corrupted(in any case), we don't read data at all from this HC
+ {
+ while (fhcMCMcounter < fHC->fMCMmax)
+ {
+ fMCM = &fHC->fMCMs[fhcMCMcounter];
+ if (!fMCM)
+ {
+ AliError(Form("HC missing at stack %d link %d atMCMslot %d",
+ fStackNumber, fStackLinkNumber, fhcMCMcounter));
+ return kFALSE;
+ }
+ while(fmcmADCcounter < fMCM->fADCmax)
+ {
+ fADC = &fMCM->fADCs[fmcmADCcounter];
+ if (!fADC)
+ {
+ AliError(Form("ADC missing at stack %d link %d MCMslot %d ADCslot %d",
+ fStackNumber, fStackLinkNumber, fhcMCMcounter, fmcmADCcounter));
+ return kFALSE;
+ }
+ fmcmADCcounter++;
+ if (fSharedPadsOn)
+ {
+ return kTRUE;
+ }
+ else
+ {
+ if (fADC->fIsShared == kFALSE)
+ return kTRUE;
+ }
+ } //while adc in MCM
+ fhcMCMcounter++;
+ // next MCM should go through all active ADCs
+ fmcmADCcounter = 0;
+ } // while mcm
+ } // if HC OK
+ }// if link active
+ fStackLinkNumber++;
+ // next stack link (HC) should go through all active MCMs
+ fhcMCMcounter = 0;
+ }// while links
+ }// if stack active
+ fStackNumber++;
+ // next stack should go through all links - start from 0
+ fStackLinkNumber = 0;
+ }
+ } // fBufferRead
+
+ // in case rawreader manages the mem buffers, go for the next buffer
if (fRawReader)
{
Int_t nextBuff = NextBuffer();
while (nextBuff != -1)
- {
- if (nextBuff > 0)
- return Next();
- nextBuff = NextBuffer();
- }
+ {
+ if (nextBuff > 0)
+ {
+ fBufferRead = kTRUE;
+ return Next();
+ }
+ nextBuff = NextBuffer();
+ }
}
return kFALSE;
}
//------------------------------------------------------------
-
Int_t
-AliTRDrawStreamTB::NextChamber(AliTRDdigitsManager *digitsManager)
+AliTRDrawStreamTB::NextChamber(AliTRDdigitsManager *digitsManager, UInt_t **trackletContainer)
{
//
// Fills single chamber digit array
// Return value is the detector number
//
- AliTRDdataArrayI *digits = 0;
+ AliTRDcalibDB *cal = AliTRDcalibDB::Instance();
+ AliTRDdataArrayDigits *digits = 0;
AliTRDdataArrayI *track0 = 0;
AliTRDdataArrayI *track1 = 0;
AliTRDdataArrayI *track2 = 0;
// Loop through the digits
Int_t lastdet = -1;
- Int_t det = -1;
- // Int_t returnDet = -1;
+ Int_t det = -1;
+ Int_t lastside = -1;
+ Int_t side = -1;
Int_t it = 0;
+ Int_t ntracklets = 0;
+
+ if (trackletContainer){
+ for (Int_t i = 0; i < 2; i++)
+ for (Int_t j = 0; j < MAX_TRACKLETS_PERHC; j++)
+ trackletContainer[i][j] = 0;
+ }
+
while (Next())
{
det = GetDet();
-
+ side = GetSide();
+
+ if (trackletContainer)
+ {
+ if ((det + side*AliTRDgeometry::kNdet) != (lastdet + lastside*AliTRDgeometry::kNdet))
+ {
+ if (det != lastdet)
+ {
+ if (lastdet != -1)
+ {
+ return lastdet;
+ }
+ }
+ ntracklets = GetNTracklets();
+ if(ntracklets > 0) memcpy(trackletContainer[side], GetTrackletWords(), sizeof(UInt_t) * ntracklets); //copy tracklet words to trackletContainer array
+ lastside = side;
+ }
+ }
+
if (det != lastdet)
- {
- // If new detector found
- if (lastdet == -1)
- {
- lastdet = det;
- }
- else
- {
- return lastdet;
- }
-
- // Add a container for the digits of this detector
- digits = digitsManager->GetDigits(det);
- track0 = digitsManager->GetDictionary(det,0);
- track1 = digitsManager->GetDictionary(det,1);
- track2 = digitsManager->GetDictionary(det,2);
-
- Int_t rowMax = GetRowMax();
- Int_t colMax = GetColMax();
- Int_t ntbins = GetNumberOfTimeBins();
-
- // Allocate memory space for the digits buffer
- if (digits->GetNtime() == 0)
- {
- digits->Allocate(rowMax, colMax, ntbins);
- track0->Allocate(rowMax, colMax, ntbins);
- track1->Allocate(rowMax, colMax, ntbins);
- track2->Allocate(rowMax, colMax, ntbins);
- }
-
- indexes = digitsManager->GetIndexes(det);
- indexes->SetSM(GetSM());
- indexes->SetStack(GetStack());
- indexes->SetLayer(GetLayer());
- indexes->SetDetNumber(det);
- if (indexes->IsAllocated() == kFALSE)
- indexes->Allocate(rowMax, colMax, ntbins);
- }
-
+ {
+ // If new detector found
+ if (lastdet == -1)
+ {
+ lastdet = det;
+ }
+ else
+ {
+ return lastdet;
+ }
+
+ if (det < 0 || det >= AliTRDgeometry::kNdet)
+ {
+ if (fSM.fClean == kTRUE)
+ {
+ AliError(Form("Strange Det Number %d BUT event buffer seems to be clean.", det));
+ }
+ else
+ {
+ AliError(Form("Strange Det Number %d. Event buffer marked NOT clean!", det));
+ }
+ continue;
+ }
+
+ // Add a container for the digits of this detector
+ digits = (AliTRDdataArrayDigits *) digitsManager->GetDigits(det);
+
+ if (digitsManager->UsesDictionaries())
+ {
+ track0 = (AliTRDdataArrayI *) digitsManager->GetDictionary(det,0);
+ track1 = (AliTRDdataArrayI *) digitsManager->GetDictionary(det,1);
+ track2 = (AliTRDdataArrayI *) digitsManager->GetDictionary(det,2);
+ }
+
+ if (!digits)
+ {
+ if (fSM.fClean == kTRUE)
+ {
+ AliError(Form("Unable to get digits for det %d BUT event buffer seems to be clean.", det));
+ }
+ else
+ {
+ AliError(Form("Unable to get digits for det %d. Event buffer is NOT clean!", det));
+ }
+ return -1;
+ }
+
+ Int_t rowMax = GetRowMax();
+ Int_t colMax = GetColMax();
+ Int_t ntbins = GetNumberOfTimeBins();
+
+ // Allocate memory space for the digits buffer
+ if (digits->GetNtime() == 0)
+ {
+ digits->Allocate(rowMax, colMax, ntbins);
+ if (digitsManager->UsesDictionaries())
+ {
+ track0->Allocate(rowMax, colMax, ntbins);
+ track1->Allocate(rowMax, colMax, ntbins);
+ track2->Allocate(rowMax, colMax, ntbins);
+ }
+ }
+
+ indexes = digitsManager->GetIndexes(det);
+ indexes->SetSM(GetSM());
+ indexes->SetStack(GetStack());
+ indexes->SetLayer(GetLayer());
+ indexes->SetDetNumber(det);
+ if (indexes->IsAllocated() == kFALSE)
+ indexes->Allocate(rowMax, colMax, ntbins);
+ }
+
+ Char_t padStatus = cal->GetPadStatus(det, GetCol(), GetRow());
+
// ntimebins data are ready to read
for (it = 0; it < GetNumberOfTimeBins(); it++)
- {
- if (GetSignals()[it] > 0)
- {
- digits->SetDataUnchecked(GetRow(), GetCol(), it, GetSignals()[it]);
-
- indexes->AddIndexTBin(GetRow(), GetCol(), it);
- track0->SetDataUnchecked(GetRow(), GetCol(), it, 0);
- track1->SetDataUnchecked(GetRow(), GetCol(), it, 0);
- track2->SetDataUnchecked(GetRow(), GetCol(), it, 0);
- }
- } // tbins
+ {
+ if (GetSignals()[it] > 0)
+ {
+ digits->SetDataUnchecked(GetRow(), GetCol(), it, GetSignals()[it]);
+ if(padStatus)
+ digits->SetPadStatus(GetRow(), GetCol(), it, padStatus);
+
+ indexes->AddIndexTBin(GetRow(), GetCol(), it);
+ if (digitsManager->UsesDictionaries())
+ {
+ track0->SetDataUnchecked(GetRow(), GetCol(), it, 0);
+ track1->SetDataUnchecked(GetRow(), GetCol(), it, 0);
+ track2->SetDataUnchecked(GetRow(), GetCol(), it, 0);
+ }
+ }
+ } // tbins
}// while Next()
- // what happens if the last HC is turned off?
return det;
- //return -1;
}
//------------------------------------------------------------
AliTRDrawStreamTB::Init()
{
//
- // general init
+ // Initialize geometry and fee parameters
//
TDirectory *saveDir = gDirectory;
fMaxADCgeom = (Int_t)fGeometry->ADCmax();
- fCommonAdditive = 10;
-
- ResetCounters();
-
- // in case rawreader manages the mem buffers
- // lets go for the next buffer
- if (fRawReader)
- {
- Int_t nextBuff = NextBuffer();
- while (nextBuff != -1)
- {
- if (nextBuff > 0)
- return kTRUE;
- nextBuff = NextBuffer();
- }
- }
-
- if (fgDebugStreamFlag == kTRUE)
- {
- if (!fgDebugStreamer)
- {
- fgDebugStreamer = new TTreeSRedirector("TRDrawDataDebug.root");
- fgStreamEventCounter = 0;
-
- if (fgDebugStreamer)
- {
- AliInfo(Form("Debug Streamer Initialized! Remember to delete! %s::DeleteDebugStream()", this->IsA()->GetName()));
- //fDebugStreamOwned = kTRUE;
- }
- else
- {
- AliError("Unable to init debug stream");
- }
- }
- }
+ ResetCounters(); // fBufferRead is set to kFALSE - important
saveDir->cd();
return kTRUE;
}
-void AliTRDrawStreamTB::DeleteDebugStream()
-{
- //
- // static helper function
- //
-
- if (fgDebugStreamer)
- {
- delete fgDebugStreamer;
- fgDebugStreamer = 0;
- }
-
-}
-
//------------------------------------------------------------
Bool_t
AliTRDrawStreamTB::InitBuffer(void *buffer, UInt_t length)
{
//
- // init the class before reading from the buffer
- // and read SM heading info:
- // -- super module header
- // -- stack headers + link status words
+ // set initial information about the buffer
//
- // not in the pre scan mode
- fgStreamEventCounter++;
-
+ if (fgDebugFlag) AliDebug(5, Form("Equipment ID: %d",fRawReader->GetEquipmentId()));
+ if (fRawReader->GetEquipmentId()<1024 || fRawReader->GetEquipmentId()>1041) //tmp protection
+ return kFALSE;
+
ResetCounters();
fpBegin = (UInt_t *)buffer;
if (WORD_SIZE == 0)
{
- AliFatal("Strange! Size of UInt_t == 0");
+ AliFatal("Strange word size. size of UInt_t == 0");
return kFALSE;
}
if (fpBegin == 0 || length <= 0)
{
- AliError(Form("This will not work! Pointer to the buffer is 0x%08x of size %d", fpBegin, length));
+ AliError(Form("Buffer size or pointer is strange. pointer to the buffer is 0x%08x of size %d", fpBegin, length));
return kFALSE;
}
if (fgDumpHead > 0)
{
- AliInfo("------------------------------------------------");
+ AliInfo(Form("---------- Dumping %u words from the beginnig of the buffer ----------",fgDumpHead));
if (DumpWords(fpBegin, fgDumpHead) == kFALSE)
- {
- AliError("Dump failed. Not enough data.");
- }
- AliInfo("------------------------------------------------");
+ {
+ AliError("Dump failed. Not enough data.");
+ }
+ AliInfo(Form("---------- Dumping ended ----------------------------------------------"));
}
- //decode SM
+ return kTRUE;
+}
+
+//------------------------------------------------------------
+Bool_t
+AliTRDrawStreamTB::DecodeGTUheader()
+{
+ // Decode Supermodule Index Word
DecodeSMInfo(fpPos, &fSM);
if (fgDebugFlag) AliDebug(5, DumpSMInfo(&fSM));
fpPos++;
if (fpPos < fpEnd)
{
+ // fSM.fHeaderSize represent additional Supermodule header size which contains additional information regarding hardware design.
+ // For the moment, we skip decoding these words
if (SkipWords(fSM.fHeaderSize) == kTRUE)
- {
- //decode stacks info
- //for (Int_t istack = 0; istack < fSM.fActiveStacks; istack++)
- for (Int_t istack = 0; istack < 5; istack++)
- {
- if (fSM.fStackActive[istack] == kFALSE)
- continue;
-
- fStack = &fSM.fStacks[istack];
- DecodeStackInfo(fpPos, fStack);
- fpPos++;
-
- fSM.fNexpectedHalfChambers += fStack->fActiveLinks;
-
- if (fgDebugFlag) AliDebug(5, DumpStackInfo(fStack));
-
- if (SkipWords(fStack->fHeaderSize) == kFALSE)
- {
- if (fRawReader) fRawReader->AddMajorErrorLog(kDecodeStackInfo, "Stack head words missing");
- return kFALSE;
- }
- }
- }
+ {
+ for (Int_t istack = 0; istack < 5; istack++)
+ {
+ if (fSM.fStackActive[istack] == kFALSE)
+ continue;
+
+ fStack = &fSM.fStacks[istack];
+
+ // Decode Stack Index Word of given stack
+ DecodeStackInfo(fpPos, fStack);
+ fpPos++;
+
+ fSM.fNexpectedHalfChambers += fStack->fActiveLinks;
+
+ if (fgDebugFlag) AliDebug(5, DumpStackInfo(fStack));
+
+ if (SkipWords(fStack->fHeaderSize-6) == kFALSE) // 6 is the 6 stack header words for 12 links
+ {
+ if (fRawReader) fRawReader->AddMajorErrorLog(kDecodeStackInfo, "Stack header words missing");
+ return kFALSE;
+ }
+ for (Int_t iword=0; iword<6; iword++) // decode 6 stack header words
+ {
+ // Decode Stack Header Word of given stack
+ DecodeStackHeader(fpPos, fStack, iword);
+ fpPos++;
+ }
+ }
+ }
else
- {
- return kFALSE;
- }
+ {
+ return kFALSE;
+ }
}
else
{
- if (fgWarnError) AliWarning("No Stack info present. Strange.");
+ if (fgWarnError) AliWarning("No additional sm headers and stack index words present.");
if (fRawReader) fRawReader->AddMajorErrorLog(kDecodeStackInfo, "Stack info missing");
return kFALSE;
}
- if (fgDebugFlag) AliDebug(5, Form("Expected half chambers : %d", fSM.fNexpectedHalfChambers));
-
- //fpPos++;
if (fpPos < fpEnd)
{
- if (fgDebugFlag) AliDebug(5, "Init OK.");
+ if (fgDebugFlag) AliDebug(5, "GTU headers are OK.");
}
else
{
- if (fgWarnError) AliWarning("No data just after init. Strange.");
- if (fRawReader) fRawReader->AddMajorErrorLog(kMissingData, "Missing data");
+ if (fgWarnError) AliWarning("No data just after GTU headers.");
+ if (fRawReader) fRawReader->AddMajorErrorLog(kMissingData, "Missing sm data");
return kFALSE;
}
+ if (fgDebugFlag) AliDebug(5, Form("Expected half chambers from GTU header: %d", fSM.fNexpectedHalfChambers));
+
return kTRUE;
}
AliTRDrawStreamTB::DecodeSM(void *buffer, UInt_t length)
{
//
- // decode all sm at once
- // still here for devel and debug
+ // decode one sm data in buffer
//
- //memset(&fSM, 0, sizeof(fSM));
-
- ResetIterators();
+ ResetIterators();
fSM.fClean = kTRUE;
- if (InitBuffer(buffer, length) == kTRUE)
+ if (InitBuffer(buffer, length) == kFALSE)
{
- // no we are already set!
- // fpPos++;
- }
- else
- {
- if (fgWarnError) AliError("--- INIT failed. ---------------");
+ if (fgWarnError) AliError("InitBuffer failed.");
fSM.fClean = kFALSE;
return kFALSE;
}
- //decode data here
- //fSM.fActiveStacks
+ if (DecodeGTUheader()== kFALSE)
+ return kFALSE;
+
for (Int_t istack = 0; istack < 5; istack++)
{
- fStackNumber = istack;
+ fStackNumber = istack;
if (fSM.fStackActive[istack] == kFALSE)
- continue;
+ continue;
fStack = &fSM.fStacks[istack];
- //fStack[istack].fActiveLinks // max is 12
+ fgLastHC = -1; // to check rob number odering
for (Int_t ilink = 0; ilink < 12; ilink++)
- {
- fStackLinkNumber = ilink;
- if (fStack->fLinksActive[ilink] == kFALSE)
- continue;
-
- if (fpPos >= fpEnd)
- {
- if (fRawReader) fRawReader->AddMajorErrorLog(kLinkDataMissing, "Link data missing");
- fSM.fClean = kFALSE;
- break;
- }
-
- fHC = &fStack->fHalfChambers[ilink];
-
- if (fSM.fTrackletEnable == kTRUE)
- {
- if (DecodeTracklets() == kFALSE)
- {
- // if (fgWarnError) AliError(Form("--- Decode Tracklets failed. --------------- Link %d of %d",
- // ilink + 1, fStack->fActiveLinks));
- if (fgDebugStreamer)
- {
- TTreeSRedirector &cstream = *fgDebugStreamer;
- cstream << "TrackletDecodeError"
- << "Event=" << fgStreamEventCounter
- << ".Stack=" << fStackNumber
- << ".Link=" << fStackLinkNumber
- << ".EndOfTrackletCount=" << fEndOfTrackletCount
- << "\n";
- }
-
- fSM.fClean = kFALSE;
- SeekEndOfData();
-
- if (fgWarnError)
- {
- AliError(Form("Failed stack %d link %d", fStackNumber, fStackLinkNumber));
- AliError(Form("Debug Event Counter : %d", fgStreamEventCounter));
- }
- continue;
- //return kFALSE;
- }
- }
-
- if (fpPos >= fpEnd)
- {
- if (fRawReader) fRawReader->AddMajorErrorLog(kHCdataMissing, "HC data missing");
- fSM.fClean = kFALSE;
- break;
- }
-
- if (DecodeHC() == kFALSE)
- {
-// if (fgWarnError) AliError(Form("--- Decode HC failed. --------------- Link %d of %d",
-// ilink + 1, fStack->fActiveLinks));
- fSM.fClean = kFALSE;
- fHC->fCorrupted += 100;
- SeekEndOfData();
-
- if (fgWarnError)
- {
- AliError(Form("Failed HC : %s", DumpHCinfoH0(fHC)));
- AliError(Form("Failed HC : %s", DumpHCinfoH1(fHC)));
- AliError(Form("Debug Event Counter : %d", fgStreamEventCounter));
- }
- // let us assume that we have the HC data although link mask says different
- // recovery
- if (fStackLinkNumber == -1)
- {
- if (fgWarnError) AliWarning("Trying to recover to the right Link Mask.");
- ilink -= 1;
- }
-
- continue;
- //return kFALSE;
- }
- else
- {
- // just finish off with the end of data markers
- SeekEndOfData();
- }
-
- //if (fgDebugFlag) AliDebug(5, Form("++++ Done with link %d of %d", ilink + 1, fStack->fActiveLinks));
-
- } // ilink
+ {
+ fStackLinkNumber = ilink;
+ if (fStack->fLinksActive[ilink] == kFALSE)
+ continue;
+
+ // check GTU link monitor
+ if (!(fStack->fLinksDataType[ilink] == 0 && fStack->fLinksMonitor[ilink] == 0))
+ {
+ fStack->fLinkMonitorError[ilink] = 1;
+ SeekEndOfData(); // skip this HC data if GTU link monitor report error
+ fStack->fLinkMonitorError[ilink] += fNWordsCounter; // counts words of given hc having link monitor error
+ continue;
+ }
+
+ if (fpPos >= fpEnd)
+ {
+ if (fRawReader) fRawReader->AddMajorErrorLog(kLinkDataMissing, "Link data missing");
+ if (fgWarnError) AliError("Link data missing.");
+ fSM.fClean = kFALSE;
+ break;
+ }
+
+ fHC = &fStack->fHalfChambers[ilink];
+ ResetPerHC();
+
+ if (fSM.fTrackletEnable == kTRUE)
+ {
+ if (DecodeTracklets() == kFALSE)
+ {
+
+ fSM.fClean = kFALSE;
+ SeekEndOfData();
+
+ if (fgWarnError)
+ {
+ AliError(Form("Tracklet decoding failed stack %d link %d", fStackNumber, fStackLinkNumber));
+ }
+ continue;
+ }
+ }
+
+ if (fpPos >= fpEnd)
+ {
+ if (fRawReader) fRawReader->AddMajorErrorLog(kHCdataMissing, "HC data missing");
+ if (fgWarnError) AliError("HC data missing.");
+ fSM.fClean = kFALSE;
+ break;
+ }
+
+ fgLastROB = -1; // to check mcm number odering
+ fgLastIndex = -1 ; // to check mcm number odering
+ if (DecodeHC() == kFALSE)
+ {
+ fSM.fClean = kFALSE;
+ if (fHC->fCorrupted < 16) SeekEndOfData(); // In case that we meet END_OF_TRACKLET_MARKERNEW
+ // during ADC data decoding or MCM header decoding
+ // we don't seek ENDOFRAWDATAMARKER
+
+ if (fgWarnError)
+ {
+ AliError(Form("Failed HC : %s", DumpHCinfoH0(fHC)));
+ AliError(Form("Failed HC : %s", DumpHCinfoH1(fHC)));
+ }
+
+ continue;
+ }
+ else
+ {
+ SeekEndOfData(); // make sure that finish off with the end of data markers
+ }
+
+ } // ilink
} // istack
- // for next()
- ResetIterators();
+ ResetIterators(); // need to do it again for Next() function
if (fSM.fClean == kTRUE)
return kTRUE;
if (fgCleanDataOnly && (fSM.fClean == kFALSE))
{
- if (fgWarnError) AliWarning("Buffer with errors. Returning FALSE.");
- fSM.fActiveStacks = 0; // Next will not give data
+ if (fgWarnError)
+ {
+ AliWarning("Buffer with errors. Returning FALSE.");
+ AliWarning(Form("--- Failed SM : %s ---", DumpSMInfo(&fSM)));
+ }
+ fSM.fActiveStacks = 0; // Next() will not give data
return kFALSE;
}
{
//
// decode SM data in case AliRawReader is in use
-
+ //
if (fRawReader)
{
Int_t nextBuff = NextBuffer();
while (nextBuff != -1)
- {
- if (nextBuff > 0)
- return nextBuff;
- nextBuff = NextBuffer();
- }
+ {
+ if (nextBuff > 0)
+ return nextBuff;
+ nextBuff = NextBuffer();
+ }
return -1;
}
else
// go to end of data marker
//
Int_t fEndOfDataCount = 0;
+ fNWordsCounter = 0;
while ( *fpPos != ENDOFRAWDATAMARKER && fpPos < fpEnd )
{
fpPos++;
+ fNWordsCounter++;
}
-
while (*fpPos == ENDOFRAWDATAMARKER && fpPos < fpEnd )
{
fEndOfDataCount++;
fpPos++;
}
-// if (fEndOfDataCount == 0)
-// {
-// if (fgWarnError) AliWarning("End of buffer reached first. No end of data marker?");
-// return kFALSE;
-// }
-
return kTRUE;
}
+//------------------------------------------------------------
+Bool_t
+AliTRDrawStreamTB::SkipMCMdata(UInt_t iw)
+{
+
+ if (fgDebugFlag) AliDebug(11,Form("Skip %d words due to MCM header corruption.",iw));
+ UInt_t iwcounter = 0;
+ while ( *fpPos != ENDOFRAWDATAMARKER && iwcounter < iw)
+ {
+ if ( *fpPos == END_OF_TRACKLET_MARKERNEW)
+ {
+ if (fgDebugFlag) AliDebug(11,"Met END_OF_TRACKLET_MARKERNEW");
+ fMCM->fCorrupted += 16;
+ fHC->fCorrupted += 16;
+ return kFALSE;
+ }
+ fpPos++;
+ iwcounter++;
+ }
+
+ if (iwcounter == iw)
+ {
+ fpPos++;
+ return kTRUE;
+ }
+
+ if (fgDebugFlag) AliDebug(11,"Met ENDOFRAWDATAMARKER");
+ return kFALSE;
+}
+
//------------------------------------------------------------
Bool_t
AliTRDrawStreamTB::SeekNextMCMheader()
// go to mcm marker
//
- // move back to the mcm pos
- fpPos = fMCM->fPos + 1;
+ fpPos++;
while ( *fpPos != ENDOFRAWDATAMARKER && fpPos < fpEnd )
{
- //if (CheckMCMMask(fpPos) == 0)
- if (MCM_HEADER_MASK_ERR(*fpPos) == 0)
- {
- if (fgWarnError) AliWarning(Form("^^^ Found : Pos 0x%08x : Val 0x%08x", fpPos, *fpPos));
- return kTRUE;
- }
+ if (MCM_HEADER_MASK_ERR(*fpPos) == 0 && MCM_HEADER_MASK_ERR(*(fpPos+1)) == 0)
+ {
+ if (fgDebugFlag) AliDebug(11,Form("^^^ Found : Pos 0x%08x : Val 0x%08x", fpPos, *fpPos));
+ return kTRUE;
+ }
+ if ( *fpPos == END_OF_TRACKLET_MARKERNEW)
+ {
+ fMCM->fCorrupted += 16;
+ fHC->fCorrupted += 16;
+ return kFALSE;
+ }
fpPos++;
}
fLinkTrackletCounter = 0;
fEndOfTrackletCount = 0;
+ fHC->fNTracklets = 0;
+
+ for (Int_t i = 0; i < MAX_TRACKLETS_PERHC; i++)
+ fHC->fTrackletWords[i] = 0;
if (fgDebugFlag) AliDebug(10, Form("Decode tracklets at 0x%08x : 0x%08x", fpPos, *fpPos));
- //cout << Form("Decode tracklets..") << endl;;
- //printf("[I] Tracklet... ? \n");
while ( *fpPos != END_OF_TRACKLET_MARKEROLD && *fpPos != END_OF_TRACKLET_MARKERNEW && fpPos < fpEnd )
{
-
if (fgDebugFlag) AliDebug(10, Form("Tracklet found at 0x%08x : 0x%08x", fpPos, *fpPos));
- //decode tracklets here...
fLinkTrackletCounter++;
if (fLinkTrackletCounter > MAX_TRACKLETS_PERHC)
- {
- if (fgWarnError) AliWarning(Form("Max number of tracklets exceeded %d > %d. Something is wrong with the input data!",
- fLinkTrackletCounter, MAX_TRACKLETS_PERHC));
-
- if (fRawReader) fRawReader->AddMajorErrorLog(kTrackletOverflow,"Too many tracklets");
- return kFALSE;
- }
-
- if (fgDebugStreamer)
- {
- // jan!
- unsigned pid;
- unsigned row;
- signed defl;
- signed ypos;
-
- double deflm;
- double yposm;
-
- unsigned long val;
-
- val = *fpPos;
-
- pid = val >> 24;
- row = (val >> 20) & 0xF;
- defl = (val >> 13) & 0x7F;
- defl = defl << 25 >> 25;
- ypos = val & 0x1FFF;
- ypos = ypos << 19 >> 19;
-
- deflm = defl * 140.0e-6; /* m */
- yposm = ypos * 160.0e-6; /* m */
-
- TTreeSRedirector &cstream = *fgDebugStreamer;
- cstream << "MCMtracklets"
- << "Event=" << fgStreamEventCounter
- << ".stack=" << fStackNumber
- << ".link=" << fStackLinkNumber
- << ".pid=" << pid
- << ".row=" << row
- << ".defl=" << defl
- << ".ypos=" << ypos
- << ".deflm=" << deflm
- << ".yposm=" << yposm
- << "\n";
- }
+ {
+ if (fgDebugFlag) AliDebug(11,Form("Max number of tracklets exceeded %d > %d.",
+ fLinkTrackletCounter, MAX_TRACKLETS_PERHC));
+ if (fRawReader) fRawReader->AddMajorErrorLog(kTrackletOverflow,"Too many tracklets");
+ fHC->fTrackletError = 1;
+ return kFALSE;
+ }
+
+ fHC->fTrackletWords[fLinkTrackletCounter-1] = UInt_t(*fpPos); //store tracklet words into array
+ fHC->fNTracklets = fLinkTrackletCounter;
fpPos++;
}
while ( ( *fpPos == END_OF_TRACKLET_MARKEROLD || *fpPos == END_OF_TRACKLET_MARKERNEW ) && fpPos < fpEnd )
{
- //seek non end of tracklets
- //printf("[I] End of tracklets. \n");
if (fgDebugFlag) AliDebug(10, Form("EoTracklets found at 0x%08x : 0x%08x", fpPos, *fpPos));
fEndOfTrackletCount++;
if ( fEndOfTrackletCount < 2 )
{
- //if (fgWarnError) AliWarning(Form("End of tracklet marker words missing %d", 2 - fEndOfTrackletCount));
+ if (fgDebugFlag) AliDebug(11,"End of tracklets word missing");
if (fRawReader) fRawReader->AddMajorErrorLog(kEOTrackeltsMissing, "End of tracklets word missing");
+ fHC->fTrackletError += 2;
return kFALSE;
}
Bool_t
AliTRDrawStreamTB::IsRowValid()
{
- //SLOW GEOM
if ( (fHC->fStack == 2 && fMCM->fROW >= fGeometry->RowmaxC0()) ||
- (fHC->fStack != 2 && fMCM->fROW >= fGeometry->RowmaxC1()) || fMCM->fROW < 0 )
+ (fHC->fStack != 2 && fMCM->fROW >= fGeometry->RowmaxC1()) || fMCM->fROW < 0 )
{
- if (fgWarnError) AliWarning(Form("SM%d L%dS%d: Wrong Padrow (%d) fROB=%d, fSIDE=%d, fMCM=%02d"
- , fHC->fSM, fHC->fLayer, fHC->fStack, fMCM->fROW, fMCM->fROB, fHC->fSide, fMCM->fMCM ));
- if (fRawReader) fRawReader->AddMajorErrorLog(kWrongPadrow, "Wrong Row");
+ if (fgDebugFlag) AliDebug(11,Form("SM%d L%dS%d: Wrong Padrow (%d) fROB=%d, fSIDE=%d, fMCM=%02d"
+ , fHC->fSM, fHC->fLayer, fHC->fStack, fMCM->fROW, fMCM->fROB, fHC->fSide, fMCM->fMCM ));
return kFALSE;
}
return kTRUE;
// check the mcm header
//
- if ( fMCM->fCorrupted > 10 )
+ if (fgLastROB != fMCM->fROB)
{
- if (fRawReader) fRawReader->AddMajorErrorLog(kMCMheaderCorrupted,"ADC mask Corrupted");
- if (fgWarnError) AliWarning(Form("Wrong ADC Mask word 0x%08x %s. Error : %d", *fMCM->fPos, DumpMCMadcMask(fMCM), fMCM->fCorrupted));
- return kFALSE;
+ fgLastIndex = 0;
+ if (fgLastROB== -1) fgLastROB = fMCM->fROB;
+ }
+ else
+ {
+ Int_t matchingcounter = 0;
+ for (Int_t i=fgLastIndex+1; i<16; i++)
+ {
+ if ( fMCM->fMCM == fgMCMordering[i] )
+ {
+ fgLastIndex = i;
+ matchingcounter++;
+ break;
+ }
+ }
+ if (matchingcounter == 0)
+ {
+ fMCM->fMCMhdCorrupted += 2;
+ AliDebug(11,Form("MCM number from last MCM is larger: MCM # from last MCM %d, MCM # from current MCM %d \n",(fMCM-1)->fMCM, fMCM->fMCM));
+ }
}
- if ( fMCM->fCorrupted > 0 )
+ if ( fgLastHC == fHC->fLayer*2 + fHC->fSide )
{
- if (fRawReader) fRawReader->AddMajorErrorLog(kMCMheaderCorrupted,"Corrupted");
- if (fgWarnError) AliWarning(Form("Wrong MCM word 0x%08x %s. Error : %d", *fMCM->fPos, DumpMCMinfo(fMCM), fMCM->fCorrupted));
- return kFALSE;
+ if ( fMCM->fROB < (fMCM-1)->fROB )
+ {
+ fMCM->fMCMhdCorrupted += 2;
+ AliDebug(11,Form("ROB number from last MCM is larger: ROB # from last MCM %d, ROB # from current MCM %d \n",(fMCM-1)->fROB, fMCM->fROB));
+ }
+ else fgLastROB = fMCM->fROB;
}
- if ( fMCM->fMCM < 0 || fMCM->fMCM > 15 || fMCM->fROB < 0 || fMCM->fROB > 7 )
+ fgLastHC = fHC->fLayer*2 + fHC->fSide;
+
+ /*
+ // this check will come back later again when we have "patched MCM map"
+ int expectedROB = -1;
+ if(!fHC->fSide) expectedROB = int(fHC->fMCMmax/16)*2;
+ else expectedROB = int(fHC->fMCMmax/16)*2 + 1;
+ int expectedMCM = 4*(3-int((fHC->fMCMmax%16)/4)) + fHC->fMCMmax%4;
+
+ if ( expectedROB != fMCM->fROB || expectedMCM != fMCM->fMCM)
{
- if (fRawReader) fRawReader->AddMajorErrorLog(kWrongMCMorROB,"Wrong ROB or MCM");
- if (fgWarnError) AliWarning(Form("Wrong fMCM or fROB. %s Skip this data.", DumpMCMinfo(fMCM)));
- return kFALSE;
+ fMCM->fMCMhdCorrupted += 2;
+ AliDebug(11,Form("ROB expected %d ROB read %d, MCM expected %d MCM read %d\n",expectedROB, fMCM->fROB, expectedMCM, fMCM->fMCM));
}
+ */
+ // below two conditions are redundant
+ /*
+ if ( fMCM->fMCM < 0 || fMCM->fMCM > 15 || fMCM->fROB < 0 || fMCM->fROB > 7 )
+ {
+ fMCM->fMCMhdCorrupted += 8; // need to assign new number
+ if (fgDebugFlag) AliDebug(11,Form("ROB or MCM number is out of range. %s\n", DumpMCMinfo(fMCM)));
+ }
if (IsRowValid() == kFALSE)
- return kFALSE;
+ {
+ fMCM->fMCMhdCorrupted += 16; // need to assign new number
+ }
+ */
+
+ if (fEventCounter == 0)
+ {
+ fEventCounter = fMCM->fEvCounter;
+ }
+
+ if (fEventCounter != fMCM->fEvCounter)
+ {
+ fMCM->fMCMhdCorrupted += 4; if (fgDebugFlag) AliDebug(11,Form("Event number(%d) of current MCM is different from that(%d) of reference MCM %s.\n", fMCM->fEvCounter, fEventCounter, DumpMCMinfo(fMCM)));
+ }
+
+ if (fEventCounter < fLastEventCounter)
+ {
+ fMCM->fMCMhdCorrupted += 8; if (fgDebugFlag) AliDebug(11,Form("Event from the past? Current %d Last %d %s.\n", fEventCounter, fLastEventCounter, DumpMCMinfo(fMCM)));
+ }
+
+ if ( fMCM->fADCmaskCorrupted > 0 )
+ return kFALSE;
+
+ if ( fMCM->fMCMhdCorrupted > 0 )
+ return kFALSE;
return kTRUE;
}
+
//------------------------------------------------------------
Bool_t
AliTRDrawStreamTB::DecodeMCMheader()
// decode the mcm header
//
- DecodeMCMheader(fpPos, fMCM);
+ DecodeMCMheader(fpPos, fMCM);
+
+ if (fDumpingEnable)
+ {
+ if (fMCM->fMCM == fDumpingMCM)
+ {
+ if (fMCM->fROB == fDumpingROB && fHC->fLayer == fDumpingLayer)
+ {
+ if (fHC->fSM == fDumpingSM && fHC->fStack == fDumpingStack)
+ {
+ if (fgDebugFlag) {
+ AliDebug(5,DumpHCinfoH0(fHC));
+ AliDebug(5,DumpMCMinfo(fMCM));
+ }
+ DumpWords(fpPos, 212);
+ }
+ }
+ }
+ }
+
+ if (fHC->fCorrupted >= 16)
+ {
+ fpPos--;
+ return kFALSE;
+ }
- fMCM->fROW = fTRDfeeParam->GetPadRowFromMCM(fMCM->fROB, fMCM->fMCM);
+ fMCM->fROW = fTRDfeeParam->GetPadRowFromMCM(fMCM->fROB, fMCM->fMCM);
- if (fHC->fRawVMajor > 2)
+ if ((fHC->fRawVMajor > 2 && fHC->fRawVMajor <5) || (fHC->fRawVMajor > 31 && fHC->fRawVMajor < 64)) //cover old and new version definition of ZS data
{
fpPos++;
if ( fpPos < fpEnd )
- {
- DecodeMask(fpPos, fMCM);
- MCMADCwordsWithTbins(fHC->fTimeBins, fMCM);
- fMCM->fAdcDataPos = fpPos + 1;
- }
+ {
+ DecodeMask(fpPos, fMCM);
+ if (fHC->fCorrupted >= 16)
+ {
+ fpPos--;
+ return kFALSE;
+ }
+ MCMADCwordsWithTbins(fHC->fTimeBins, fMCM);
+ fMCM->fAdcDataPos = fpPos + 1;
+ }
else
- {
- if (fgWarnError) AliWarning("Expected ADC mask word. Fail due to buffer END.");
- if (fRawReader) fRawReader->AddMajorErrorLog(kMCMADCMaskMissing,"Missing");
- return kFALSE;
- }
+ {
+ if (fgDebugFlag) AliDebug(11,"Expected ADC mask word. Fail due to buffer END.");
+ if (fRawReader) fRawReader->AddMajorErrorLog(kMCMADCMaskMissing,"Missing");
+ fHC->fCorrupted += 32;
+ return kFALSE;
+ }
}
else
{
UInt_t dummyMask = MCM_DUMMY_ADCMASK_VAL;
- DecodeMask(&dummyMask, fMCM);
+ DecodeMask(&dummyMask, fMCM);
MCMADCwordsWithTbins(fHC->fTimeBins, fMCM);
fMCM->fAdcDataPos = fpPos + 1;
}
}
if (IsMCMheaderOK() == kFALSE)
- {
- return kFALSE;
- }
-
- if (fEventCounter == 0)
- {
- fEventCounter = fMCM->fEvCounter;
- }
-
- if (fEventCounter < fLastEventCounter)
- {
- if (fgWarnError) AliWarning(Form("Weird. Event from the past? Current %d Last %d %s", fEventCounter, fLastEventCounter, DumpMCMinfo(fMCM)));
- if (fRawReader) fRawReader->AddMajorErrorLog(kMCMeventMissmatch, "Wrong MCM event counter ? Past-Future");
- return kFALSE;
- }
-
- if (fEventCounter != fMCM->fEvCounter)
- {
- if (fgWarnError) AliWarning(Form("Weird. Event missmatch? %d %s", fEventCounter, DumpMCMinfo(fMCM)));
- if (fRawReader) fRawReader->AddMajorErrorLog(kMCMeventMissmatch, "Wrong MCM event counter ?");
return kFALSE;
- }
-
- //fLastEventCounter = fEventCounter;
- //AliInfo(Form("Current MCM %s", DumpMCMinfo(fMCM)));
return kTRUE;
- //return IsMCMheaderOK();
}
//------------------------------------------------------------
AliTRDrawStreamTB::IsHCheaderOK()
{
//
- // check the half chamber header
+ // check insanity of half chamber header
//
- if (fHC->fCorrupted > 0)
- {
- if (fgWarnError) AliWarning(Form("Wrong HC Header word. Word 0x%08x Error : %d", *fHC->fPos, fHC->fCorrupted));
- if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderCorrupt, "Corrupted");
-
- return kFALSE;
- }
-
if (fHC->fStack < 0 || fHC->fStack > 4)
{
- if (fgWarnError) AliWarning(Form("Wrong Stack %d", fHC->fStack));
- if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderWrongStack, "Wrong Stack");
+ if (fgDebugFlag) AliDebug(11,Form("Wrong Stack %d", fHC->fStack));
return kFALSE;
}
- if (fHC->fStack != fStackNumber)
+ if (fHC->fLayer < 0 || fHC->fLayer >= AliTRDgeometry::kNlayer)
{
- if (fgWarnError) AliWarning(Form("Missmatch: Stack in HC header %d HW-stack %d",
- fHC->fStack, fStackNumber));
- if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderWrongStack, "Stack-HWstack");
- // Try this for recovery in DecodeSM(void*,UInt_t) after DecodeHC failed
- // buffer will still will be marked as NOT clean
- fStackNumber = -1;
+ if (fgDebugFlag) AliDebug(11,Form("Wrong layer %d", fHC->fLayer));
return kFALSE;
}
- if (fHC->fLayer < 0 || fHC->fLayer >= AliTRDgeometry::kNplan)
+ if (fHC->fSide < 0 || fHC->fSide > 1)
{
- if (fgWarnError) AliWarning(Form("Wrong plane(layer) %d", fHC->fLayer));
- if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderWrongLayer, "Wrong Plane");
+ if (fgDebugFlag) AliDebug(11,Form("Wrong Side %d", fHC->fSide));
return kFALSE;
}
- if ((fHC->fLayer * 2 != fStackLinkNumber) && (fHC->fLayer * 2 != fStackLinkNumber - 1))
+ if (fgStackNumberChecker)
{
- if (fgWarnError) AliWarning(Form("Missmatch: plane(layer) in HCheader %d HW-Link %d | %s",
- fHC->fLayer, fStackLinkNumber, DumpStackInfo(fStack)));
- if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderWrongLayer, "Plane-Link missmatch");
- // Try this for recovery in DecodeSM(void*,UInt_t) after DecodeHC failed
- // buffer will still will be marked as NOT clean
- fStackLinkNumber = -1;
- return kFALSE;
+ if (fHC->fStack != fStackNumber)
+ {
+ if (fgDebugFlag) AliDebug(11,Form("Missmatch: Stack number between HC header %d and GTU link mask %d",
+ fHC->fStack, fStackNumber));
+ fStackNumber = -1;
+ return kFALSE;
+ }
}
- if (fHC->fSide < 0 || fHC->fSide > 1)
+ if (fgStackLinkNumberChecker)
{
- if (fgWarnError) AliWarning(Form("Wrong Side %d", fHC->fSide));
- if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderWrongSide, "Wrong Side");
- return kFALSE;
+ //if (fHC->fLayer * 2 + fHC->fSide != fStackLinkNumber)
+ // let it make flexible to consider known fiber swapping
+ if ((fHC->fLayer * 2 != fStackLinkNumber) && (fHC->fLayer * 2 != fStackLinkNumber - 1))
+ {
+ if (fgDebugFlag) AliDebug(11,Form("Missmatch: Layer number between HC header %d and GTU link mask %d | %s",
+ fHC->fLayer, fStackLinkNumber, DumpStackInfo(fStack)));
+ fStackLinkNumber = -1;
+ return kFALSE;
+ }
}
-
- // SLOW GEOM
+
+ // SLOW GEOM : consistancy check with geometry
fHC->fDET = fGeometry->GetDetector(fHC->fLayer, fHC->fStack, fHC->fSM);
if (fHC->fDET < 0 || fHC->fDET >= AliTRDgeometry::kNdet)
{
- if (fgWarnError) AliWarning(Form("Wrong detector %d", fHC->fDET));
+ if (fgDebugFlag) AliDebug(11,Form("Wrong detector %d", fHC->fDET));
if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderWrongDet, "Wrong Det");
return kFALSE;
}
- // SLOW GEOM
- // this check fails - raw data inconsistent with geometry?
if (fHC->fSM != fGeometry->GetSector(fHC->fDET)
- || fHC->fSM <0 || fHC->fSM >= AliTRDgeometry::kNsect)
+ || fHC->fSM <0 || fHC->fSM >= AliTRDgeometry::kNsector)
{
- if (fgWarnError) AliWarning(Form("Wrong SM(sector) %d (Geometry says: %d) Stack=%d Layer=%d Det=%d",
- fHC->fSM, fGeometry->GetSector(fHC->fDET),
- fHC->fStack, fHC->fLayer, fHC->fDET));
+ if (fgDebugFlag) AliDebug(11,Form("Wrong SM(sector) %d (Geometry says: %d) Stack=%d Layer=%d Det=%d",
+ fHC->fSM, fGeometry->GetSector(fHC->fDET),
+ fHC->fStack, fHC->fLayer, fHC->fDET));
if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderWrongSM, "Wrong SM");
return kFALSE;
}
- // SLOW GEOM
- // CPU EXPENSIVE!!!
fHC->fROC = fGeometry->GetDetectorSec(fHC->fLayer, fHC->fStack);
if (fHC->fROC < 0)
{
if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderWrongROC, "Wrong ROC Col Max");
return kFALSE;
}
- // extra - not needed
- // if (fHC->fSM <0 || fHC->fSM >= AliTRDgeometry::kNsect)
- // {
- // if (fgWarnError) AliWarning(Form("Wrong SM(sector) %d (Geometry says: %d) Stack=%d Layer=%d",
- // fHC->fSM, fGeometry->GetDetectorSec(fHC->fLayer, fHC->fStack),
- // fHC->fStack, fHC->fLayer));
- // return kFALSE;
- // }
return kTRUE;
}
// decode the half chamber header
//
- DecodeHCwordH0(fpPos, fHC);
+ if (DecodeHCwordH0(fpPos, fHC) == kFALSE)
+ return kFALSE;
if (fHC->fNExtraWords > 0)
{
fpPos++;
if (fpPos < fpEnd)
- {
- DecodeHCwordH1(fpPos, fHC);
- }
+ {
+ if (DecodeHCwordH1(fpPos, fHC) == kFALSE)
+ return kFALSE;
+ }
else
- {
- if (fgWarnError) AliWarning("Expected HC header word 1. Fail due to buffer END.");
- if (fRawReader) fRawReader->AddMajorErrorLog(kHCWordMissing,"Next HC word 1 (count from 0) missing");
- return kFALSE;
- }
+ {
+ if (fgDebugFlag) AliDebug(11,"Expected HC header word 1. Fail due to buffer END.");
+ if (fRawReader) fRawReader->AddMajorErrorLog(kHCWordMissing,"Next HC word 1 (count from 0) missing");
+ return kFALSE;
+ }
}
if (fgDebugFlag) AliDebug(5, DumpHCinfoH0(fHC));
fHC->fDET = -1;
if (IsHCheaderOK() == kFALSE)
{
+ fHC->fH0Corrupted += 2;
+ if (fgDebugFlag) AliDebug(11,Form("H0 Header Insane. Word 0x%08x", *fHC->fPos));
return kFALSE;
}
AliTRDrawStreamTB::DecodeHC()
{
//
- // decode the hc data
- // function for devel & debug
+ // decode hc header and data
//
if (DecodeHCheader() == kFALSE)
{
- if (fgDebugStreamer)
- {
- TTreeSRedirector &cstream = *fgDebugStreamer;
- cstream << "HCDecodeError"
- << "Event=" << fgStreamEventCounter
- << ".stack=" << fStackNumber
- << ".link=" << fStackLinkNumber
- << ".hcCorrupted" << fHC->fCorrupted
- << ".hcRVmajor=" << fHC->fRawVMajor
- << ".hcDCSboard=" << fHC->fDCSboard
- << ".hcSM=" << fHC->fSM
- << ".hcStack=" << fHC->fStack
- << ".hcLayer=" << fHC->fLayer
- << ".hcSide=" << fHC->fSide
- << ".hcTbins=" << fHC->fTimeBins
- << ".hcBunchCross=" << fHC->fBunchCrossCounter
- << ".hcPreTrig=" << fHC->fPreTriggerCounter
- << ".hcPreTrigPh=" << fHC->fPreTriggerPhase
- << ".hcDet=" << fHC->fDET
- << ".hcROC=" << fHC->fROC
- << ".hc=RowMax" << fHC->fRowMax
- << ".hc=ColMax" << fHC->fColMax
- << "\n";
- }
-
- if (fgWarnError) AliWarning("Header decode failed.");
+ if (fgWarnError) AliWarning(Form("HC Header decode failed. H0 Error: %d H1 Error: %d",fHC->fH0Corrupted,fHC->fH1Corrupted));
+ if (fRawReader) fRawReader->AddMajorErrorLog(kHCHeaderCorrupt, "HC header corrupted");
return kFALSE;
}
else
{
fpPos++;
if (fpPos >= fpEnd)
- {
- if (fgWarnError) AliWarning("No MCM data? Not enough data in the buffer.");
- if (fRawReader) fRawReader->AddMajorErrorLog(kMCMdataMissing, "MCM data missing");
- return kFALSE;
- }
+ {
+ fHC->fCorrupted += 1;
+ if (fgDebugFlag) AliDebug(11,"No MCM data? Not enough data in the buffer.");
+ if (fRawReader) fRawReader->AddMajorErrorLog(kMCMdataMissing, "MCM data missing");
+ return kFALSE;
+ }
}
fHC->fMCMmax = 0;
while (*fpPos != ENDOFRAWDATAMARKER && fpPos < fpEnd)
{
if (fHC->fMCMmax > TRD_MAX_MCM)
- {
- if (fgWarnError) AliWarning("More mcm data than expected!");
- if (fRawReader) fRawReader->AddMajorErrorLog(kMCMoverflow, "Too many mcms found!");
- return kFALSE;
- }
+ {
+ fHC->fCorrupted += 2;
+ if (fgDebugFlag) AliDebug(11,"More mcm data than expected!");
+ if (fRawReader) fRawReader->AddMajorErrorLog(kMCMoverflow, "Too many mcms found!");
+ return kFALSE;
+ }
fMCM = &fHC->fMCMs[fHC->fMCMmax];
if (DecodeMCMheader() == kFALSE)
- {
- if (fgDebugStreamer)
- {
- TTreeSRedirector &cstream = *fgDebugStreamer;
- cstream << "MCMDecodeError"
- << "Event=" << fgStreamEventCounter
- << ".stack=" << fStackNumber
- << ".link=" << fStackLinkNumber
- << ".hcDCSboard=" << fHC->fDCSboard
- << ".hcSM=" << fHC->fSM
- << ".hcStack=" << fHC->fStack
- << ".hcLayer=" << fHC->fLayer
- << ".hcSide=" << fHC->fSide
- << ".hcBunchCross=" << fHC->fBunchCrossCounter
- << ".hcPreTrig=" << fHC->fPreTriggerCounter
- << ".hcPreTrigPh=" << fHC->fPreTriggerPhase
- << ".hcDet=" << fHC->fDET
-
- << ".mcmCorrupted=" << fMCM->fCorrupted
- << ".mcmROB=" << fMCM->fROB
- << ".mcmMCM=" << fMCM->fMCM
- << ".mcmADCMaskWord=" << fMCM->fADCMaskWord
- << ".mcmEventCounter=" << fMCM->fEvCounter
-
- << ".fEventCounter=" << fEventCounter
- << ".fLastEventCounter=" << fLastEventCounter
- << "\n";
- }
-
- return kFALSE;
- }
+ {
+ if (fHC->fCorrupted < 4) fHC->fCorrupted += 4; // benchmark hc data corruption as 4
+
+ if (fgSkipData == kTRUE || fHC->fCorrupted >= 16)
+ return kFALSE; // stop HC data reading
+
+ fHC->fMCMmax++; // increase mcm counter to match with expected rob/mcm number
+
+ // in case we decide to keep reading data, skip this mcm data and find next mcm header
+ if (fMCM->fADCmaskCorrupted < 2)
+ {
+ if (SkipMCMdata(fMCM->fADCcount*fMCM->fSingleADCwords) == kFALSE)
+ return kFALSE;
+ continue;
+ }
+ else
+ {
+ if (SeekNextMCMheader() == kFALSE)
+ return kFALSE;
+ continue;
+ }
+ }
fHC->fMCMmax++;
if (fMCM->fADCmax > 0)
- {
- fpPos++;
- if (fpPos >= fpEnd)
- {
- if (fgDebugFlag) AliDebug(9, Form("Buffer short of data. ADC data expected."));
- if (fRawReader) fRawReader->AddMajorErrorLog(kADCdataMissing, "ADC data missing");
- }
-
- // for (Int_t iadc = 0; iadc < TRD_MAX_ADC; iadc++)
- for (Int_t iadc = 0; iadc < fMCM->fADCmax; iadc++)
- {
- fADC = &fMCM->fADCs[iadc];
- fADC->fADCnumber = fMCM->fADCchannel[iadc];
-
- if (fgDebugFlag) AliDebug(9, Form("This is ADC %d of %d. ADC number is %d.",
- iadc+1, fMCM->fADCmax, fMCM->fADCchannel[iadc]));
-
- if (fpPos + fMCM->fSingleADCwords >= fpEnd)
- {
- if (fgDebugStreamer)
- {
- TTreeSRedirector &cstream = *fgDebugStreamer;
- cstream << "ADCDecodeError"
- << "Event=" << fgStreamEventCounter
- << ".hcSM=" << fHC->fSM
- << ".hcStack=" << fHC->fStack
- << ".hcLayer=" << fHC->fLayer
- << ".mcmROB=" << fMCM->fROB
- << ".mcmMCM=" << fMCM->fMCM
- << ".mcmADCMaskWord=" << fMCM->fADCMaskWord
- << "\n";
- }
-
- if (fgWarnError) AliWarning(Form("This is ADC %d of %d. ADC number is %d.", iadc+1, fMCM->fADCmax, fMCM->fADCchannel[iadc]));
- if (fgWarnError) AliWarning("--> ADC (10 words) expected. Not enough data in the buffer.");
- if (fgWarnError) AliWarning(Form("--> ADC mask : %s", DumpMCMadcMask(fMCM)));
-
- if (fRawReader) fRawReader->AddMajorErrorLog(kADCdataMissing, "ADC data missing - less than expected");
- return kFALSE;
- }
-
- //DECODE the data here
- if (DecodeADC() == kFALSE)
- {
- // check if we are out of the det when the pad is shared
- if (fADC->fIsShared && fADC->fCorrupted == 111)
- {
- fADC->fCOL = -1;
- fpPos = fADC->fPos + fMCM->fSingleADCwords;
- }
- else
- {
- if (fgWarnError) AliWarning(Form("ADC decode failed."));
- //fpPos = fMCM->fPpos + fMCMADCWords;
- fpPos = fADC->fPos + fMCM->fSingleADCwords;
- return kFALSE;
- }
- }
- //decode the ADC words here
- }
- } //mcm data present
+ {
+ fpPos++;
+ if (fpPos >= fpEnd)
+ {
+ fMCM->fCorrupted += 1;
+ if (fHC->fCorrupted < 4) fHC->fCorrupted += 4; // benchmark hc data corruption as 4
+ if (fgDebugFlag) AliDebug(9, Form("Buffer short of data. ADC data expected."));
+ return kFALSE;
+ }
+
+ for (Int_t iadc = 0; iadc < fMCM->fADCmax; iadc++)
+ {
+ fADC = &fMCM->fADCs[iadc];
+ fADC->fADCnumber = fMCM->fADCchannel[iadc];
+
+ if (fgDebugFlag) AliDebug(9, Form("This is ADC %d of %d. ADC number is %d.",
+ iadc+1, fMCM->fADCmax, fMCM->fADCchannel[iadc]));
+
+ if (fpPos + fMCM->fSingleADCwords >= fpEnd)
+ {
+
+ fMCM->fCorrupted += 2;
+ if (fHC->fCorrupted < 4) fHC->fCorrupted += 4; // benchmark hc data corruption as 4
+ if (fgDebugFlag) AliDebug(11,"ADC (10 words) expected. Not enough data in the buffer.");
+ return kFALSE;
+ }
+
+ if (fHC->fRawVMajor < 64) // normal(real) ADC data
+ {
+ if (DecodeADC() == kFALSE)
+ {
+ if (fMCM->fCorrupted < 4) fMCM->fCorrupted += 4; // benchmark mcm data corruption as 4
+ if (fHC->fCorrupted < 4) fHC->fCorrupted += 4; // benchmark hc data corruption as 4
+ if (fADC->fIsShared && fADC->fCorrupted == 16) // check if we are out of the det when the pad is shared
+ {
+ fADC->fCOL = -1;
+ fpPos = fADC->fPos + fMCM->fSingleADCwords;
+ }
+ else
+ {
+ if (fgDebugFlag) AliDebug(11,Form("ADC decode failed."));
+ if (fgSkipData == kTRUE || fHC->fCorrupted >= 16)
+ return kFALSE; // stop HC data reading
+ }
+ }
+ }
+ else // test pattern data
+ {
+ if (fgWarnError) AliError("These are test pattern data. You need other reader"); // will be served in other class
+ }
+ }
+ }
else
- {
- fpPos++;
- }
+ {
+ fpPos++;
+ }
}//while eof data
if (fpPos >= fpEnd)
{
- if (fgWarnError) AliWarning("We are at the end of buffer. There should be one more word left.");
- //SeekEndOfData();
+ if (fgDebugFlag) AliDebug(11,"We are at the end of buffer. There should be one more word left.");
return kFALSE;
}
return kTRUE;
}
-
//------------------------------------------------------------
+
Bool_t
AliTRDrawStreamTB::DecodeADC()
{
//
fADC->fCorrupted = 0;
- fMaskADCword = ADC_WORD_MASK(*fpPos);
- fADC->fPos = fpPos;
+ if(fADC->fADCnumber%2==1) fMaskADCword = ADC_WORD_MASK(ADCDATA_VAL1);
+ if(fADC->fADCnumber%2==0) fMaskADCword = ADC_WORD_MASK(ADCDATA_VAL2);
+ fADC->fPos = fpPos;
fTbinADC = 0;
for (Int_t i = 0; i < TRD_MAX_TBINS; i++)
for (Int_t iw = 0; iw < fMCM->fSingleADCwords; iw++)
{
+ if (HC_HEADER_MASK_ERR(*fpPos) == 0 || *fpPos == END_OF_TRACKLET_MARKERNEW)
+ {
+ if (fgWarnError) AliError(Form("There should be ADC data. We meet HC header or END_OF_TRACKLET_MARKER 0x%08x",*fpPos));
+ fADC->fCorrupted += 16;
+ fHC->fCorrupted += 16;
+ fpPos--;
+
+ return kFALSE;
+ }
if (fMaskADCword != ADC_WORD_MASK(*fpPos))
- {
- //this is corrupted data
- fADC->fCorrupted += 1;
- if (fgWarnError) AliWarning(Form("Mask Change in ADC data 0x%08x 0x%08x", *(fpPos-1), fMaskADCword, *fpPos, ADC_WORD_MASK(*fpPos)));
- if (fRawReader) fRawReader->AddMajorErrorLog(kADCmaskMissmatch, "Mask change inside single channel");
-
- break;
- }
-
- fADC->fSignals[fTbinADC + 0] = (*fpPos & 0x00000ffc) >> 2;
- fADC->fSignals[fTbinADC + 1] = (*fpPos & 0x003ff000) >> 12;
- fADC->fSignals[fTbinADC + 2] = (*fpPos & 0xffc00000) >> 22;
+ {
+ fADC->fCorrupted += 1;
+ if (fgDebugFlag) AliDebug(11,Form("Wrong ADC data mask! ADC channel number: %02d [Expected mask: 0x%08x Current mask: 0x%08x] MCM= %s Error : %d",
+ fADC->fADCnumber, fMaskADCword, ADC_WORD_MASK(*fpPos),DumpMCMinfo(fMCM),fADC->fCorrupted));
+ fpPos++;
+ continue;
+ }
+
+ // here we subtract the baseline ( == common additive)
+ fADC->fSignals[fTbinADC + 0] = ((*fpPos & 0x00000ffc) >> 2) - fgCommonAdditive;
+ fADC->fSignals[fTbinADC + 1] = ((*fpPos & 0x003ff000) >> 12) - fgCommonAdditive;
+ fADC->fSignals[fTbinADC + 2] = ((*fpPos & 0xffc00000) >> 22) - fgCommonAdditive;
fTbinADC += 3;
fpPos++;
if ( fADC->fADCnumber >= fMaxADCgeom - 1)
{
- // let us guess the Column
- // take the one before last ADC and shift by one column
- // later we check if we are inside the limits of the chamber
fADC->fCOL = AliTRDfeeParam::Instance()->GetPadColFromADC(fMCM->fROB, fMCM->fMCM, fADC->fADCnumber - 1);
fADC->fCOL--;
}
if (fADC->fCOL >= fHC->fColMax || fADC->fCOL < 0)
{
if (fADC->fIsShared == kFALSE)
- {
- fADC->fCorrupted += 100;
-
- if (fgWarnError) AliWarning(Form("Wrong column! ADCnumber %d MaxIs %d Col %d MaxIs %d MCM= %s",
- fADC->fADCnumber, fMaxADCgeom, fADC->fCOL, fHC->fColMax, DumpMCMinfo(fMCM)));
- if (fRawReader) fRawReader->AddMajorErrorLog(kWrongPadcolumn, "Wrong column");
- }
- else
- {
- // flag it - we are out of the det when the pad is shared
- if (fgDebugFlag) AliDebug(10, Form("Column out of the detector! ADCnumber %d MaxIs %d Col %d MaxIs %d MCM= %s",
- fADC->fADCnumber, fMaxADCgeom, fADC->fCOL, fHC->fColMax, DumpMCMinfo(fMCM)));
- fADC->fCorrupted += 111;
- }
+ {
+ fADC->fCorrupted += 32;
+ if (fgDebugFlag) AliDebug(11,Form("Wrong column! ADCnumber %d MaxIs %d Col %d MaxIs %d MCM= %s",
+ fADC->fADCnumber, fMaxADCgeom, fADC->fCOL, fHC->fColMax, DumpMCMinfo(fMCM)));
+ }
+ //else
+ //{
+ // we are out of the det when the pad is shared
+ //if (fgDebugFlag) AliDebug(11, Form("Column out of the detector! ADCnumber %d MaxIs %d Col %d MaxIs %d MCM= %s",
+ // fADC->fADCnumber, fMaxADCgeom, fADC->fCOL, fHC->fColMax, DumpMCMinfo(fMCM)));
+ //fADC->fCorrupted += 32;
+ //}
}
if (fADC->fCorrupted > 0)
}
//--------------------------------------------------------
+
+
void AliTRDrawStreamTB::DecodeSMInfo(const UInt_t *word, struct AliTRDrawSM *sm) const
{
//
- // Decode the data *word into SM info structure
+ // Decode Supermodule Index Word
+ // The Supermodule Index Word is a 32-Bit word wit following structure
+ // ssssssss ssssssss vvvv rrrr r d t mmmm
+ // s: Size of the Supermodule Header, v: Supermodule Header Version, r: Reserved for future use
+ // d: Track Data Enabled Bit, t: Tracklet Data Enabled Bit, m: Stack Mask
//
-
- sm->fPos = (UInt_t*)word;
+ sm->fPos = (UInt_t*)word;
- // do it once here
UInt_t vword = *word;
sm->fHeaderSize = SM_HEADER_SIZE(vword);
sm->fActiveStacks = 0;
for (Int_t i = 0; i < 5; i++)
{
- //if ((stackMask >> i) & 0x1 > 0)
if (IS_BIT_SET(stackMask,i) > 0)
- {
- sm->fStackActive[i] = kTRUE;
- sm->fActiveStacks++;
- }
+ {
+ sm->fStackActive[i] = kTRUE;
+ sm->fActiveStacks++;
+ }
else
- {
- sm->fStackActive[i] = kFALSE;
- }
+ {
+ sm->fStackActive[i] = kFALSE;
+ }
}
}
//
// Get SM structure into a const char
//
-
return Form("[ SM Info 0x%08x] : Hsize %d TrackletEnable %d Stacks %d %d %d %d %d",
- *sm->fPos,
- sm->fHeaderSize, sm->fTrackletEnable,
- sm->fStackActive[0], sm->fStackActive[1], sm->fStackActive[2],
- sm->fStackActive[3], sm->fStackActive[4]);
+ *sm->fPos,
+ sm->fHeaderSize, sm->fTrackletEnable,
+ sm->fStackActive[0], sm->fStackActive[1], sm->fStackActive[2],
+ sm->fStackActive[3], sm->fStackActive[4]);
}
//--------------------------------------------------------
void AliTRDrawStreamTB::DecodeStackInfo(const UInt_t *word, struct AliTRDrawStack *st) const
{
//
- // Decode stack info - active links
+ // Decode Stack #i Index Word
+ // The Stack #i Index Word is a 32-Bit word wit following structure
+ // ssssssss ssssssss vvvv mmmm mmmmmmmm
+ // s: Size of the Stack #i Header, v: Supermodule Header Version, m: Link Mask
//
-
st->fPos = (UInt_t*)word;
- // do it once here
UInt_t vword = *word;
st->fHeaderSize = STACK_HEADER_SIZE(vword);
st->fActiveLinks = 0;
for (Int_t i = 0; i < 12; i++)
{
- //if (linkMask & (0x1 << i) > 0)
if (IS_BIT_SET(linkMask,i) > 0)
- {
- st->fLinksActive[i] = kTRUE;
- st->fTrackletDecode[i] = kTRUE;
- st->fHCDecode[i] = kTRUE;
- st->fActiveLinks++;
- }
+ {
+ st->fLinksActive[i] = kTRUE;
+ st->fActiveLinks++;
+ }
else
- {
- st->fTrackletDecode[i] = kFALSE;
- st->fLinksActive[i] = kFALSE;
- st->fHCDecode[i] = kFALSE;
- }
+ {
+ st->fLinksActive[i] = kFALSE;
+ }
}
}
+//--------------------------------------------------------
+void AliTRDrawStreamTB::DecodeStackHeader(const UInt_t *word, struct AliTRDrawStack *st, Int_t iword) const
+{
+ st->fPos = (UInt_t*)word;
+
+ UInt_t vword = *word;
+ st->fLinksDataType[2*iword] = LINK0_DATA_TYPE_FLAG(vword);
+ st->fLinksMonitor[2*iword] = LINK0_MONITOR_FLAG(vword);
+ st->fLinksDataType[2*iword+1] = LINK1_DATA_TYPE_FLAG(vword);
+ st->fLinksMonitor[2*iword+1] = LINK1_MONITOR_FLAG(vword);
+}
+
//--------------------------------------------------------
const char *AliTRDrawStreamTB::DumpStackInfo(const struct AliTRDrawStack *st)
{
//
// format the string with the stack info
//
+
return Form("[ Stack Info 0x%08x ] : Hsize %d Links Active %d %d %d %d %d %d %d %d %d %d %d %d",
- *st->fPos,
- st->fHeaderSize,
- st->fLinksActive[0], st->fLinksActive[1], st->fLinksActive[2], st->fLinksActive[3],
- st->fLinksActive[4], st->fLinksActive[5], st->fLinksActive[6], st->fLinksActive[7],
- st->fLinksActive[8], st->fLinksActive[9], st->fLinksActive[10], st->fLinksActive[11]);
+ *st->fPos,
+ st->fHeaderSize,
+ st->fLinksActive[0], st->fLinksActive[1], st->fLinksActive[2], st->fLinksActive[3],
+ st->fLinksActive[4], st->fLinksActive[5], st->fLinksActive[6], st->fLinksActive[7],
+ st->fLinksActive[8], st->fLinksActive[9], st->fLinksActive[10], st->fLinksActive[11]);
+
}
//--------------------------------------------------------
-void AliTRDrawStreamTB::DecodeHCwordH0(const UInt_t *word, struct AliTRDrawHC *hc) const
+Bool_t AliTRDrawStreamTB::DecodeHCwordH0(const UInt_t *word, struct AliTRDrawHC *hc) const
{
//
// decode the hc header word 0
//
-
- // do it once here
UInt_t vword = *word;
- hc->fCorrupted = HC_HEADER_MASK_ERR(vword);
- if (hc->fCorrupted > 0)
+ hc->fH0Corrupted = HC_HEADER_MASK_ERR(vword);
+ if (hc->fH0Corrupted > 0)
{
- hc->fH0ErrorCounter++;
+ if (fgDebugFlag) AliDebug(11,Form("H0 Header Mask Error. Word 0x%08x", *fHC->fPos));
+ return kFALSE;
}
hc->fSpecialRawV = HC_SPECIAL_RAW_VERSION(vword);
hc->fRawVMajor = HC_MAJOR_RAW_VERSION(vword);
+ hc->fRawVMajorOpt = HC_MAJOR_RAW_VERSION_OPT(vword);
hc->fRawVMinor = HC_MINOR_RAW_VERSION(vword);
hc->fNExtraWords = HC_EXTRA_WORDS(vword);
hc->fDCSboard = HC_DCS_BOARD(vword);
hc->fPos[0] = (UInt_t*)word;
+ return kTRUE;
}
//--------------------------------------------------------
-void AliTRDrawStreamTB::DecodeHCwordH1(const UInt_t *word, struct AliTRDrawHC *hc) const
+Bool_t AliTRDrawStreamTB::DecodeHCwordH1(const UInt_t *word, struct AliTRDrawHC *hc) const
{
//
- //
+ // decode the hc header word 1
//
- // do it once here
UInt_t vword = *word;
- hc->fCorrupted += 10 * HC_HEADER_MASK_ERR(vword);
- if (hc->fCorrupted > 10)
- {
- hc->fH1ErrorCounter++;
+ hc->fH1Corrupted = HC_HEADER_MASK_ERR(vword);
+ if (hc->fH1Corrupted > 0)
+ {
+ if (fgDebugFlag) AliDebug(11,Form("H1 Header Mask Error. Word 0x%08x", *fHC->fPos));
+ return kFALSE;
}
hc->fTimeBins = HC_NTIMEBINS(vword);
hc->fPreTriggerPhase = HC_PRETRIGGER_PHASE(vword);
hc->fPos[1] = (UInt_t*)word;
+
+ return kTRUE;
}
//--------------------------------------------------------
const char *AliTRDrawStreamTB::DumpHCinfoH0(const struct AliTRDrawHC *hc)
{
//
- //
+ // dump the hc header word 0
//
if (!hc)
return Form("Unable to dump. Null received as parameter!?!");
else
return Form("[ HC[0] at 0x%08x ] : 0x%08x Info is : RawV %d SM %d Stack %d Layer %d Side %d DCSboard %d",
- hc->fPos[0], *(hc->fPos[0]), hc->fRawVMajor, hc->fSM, hc->fStack, hc->fLayer, hc->fSide, hc->fDCSboard);
+ hc->fPos[0], *(hc->fPos[0]), hc->fRawVMajor, hc->fSM, hc->fStack, hc->fLayer, hc->fSide, hc->fDCSboard);
}
//--------------------------------------------------------
const char *AliTRDrawStreamTB::DumpHCinfoH1(const struct AliTRDrawHC *hc)
{
//
- //
+ // dump the hc header word 1
//
if (!hc)
return Form("Unable to dump. Null received as parameter!?!");
else
return Form("[ HC[1] at 0x%08x ] : 0x%08x Info is : TBins %d BCcount %d PreTrigCount %d PreTrigPhase %d",
- hc->fPos[1], *(hc->fPos[1]), hc->fTimeBins, hc->fBunchCrossCounter, hc->fPreTriggerCounter, hc->fPreTriggerPhase);
+ hc->fPos[1], *(hc->fPos[1]), hc->fTimeBins, hc->fBunchCrossCounter, hc->fPreTriggerCounter, hc->fPreTriggerPhase);
}
//--------------------------------------------------------
//
// decode the mcm header
//
-
- // do it once here
UInt_t vword = *word;
- mcm->fCorrupted = MCM_HEADER_MASK_ERR(vword);
- //printf("0x%08x %d\n", word, MCM_HEADER_MASK_ERR(word));
- if (mcm->fCorrupted)
- mcm->fErrorCounter++;
+ if (vword == END_OF_TRACKLET_MARKERNEW)
+ {
+ if (fgWarnError) AliError(Form("There should be MCM header. We meet END_OF_TRACKLET_MARKER 0x%08x",vword));
+ mcm->fMCMhdCorrupted += 16;
+ fHC->fCorrupted += 16; //to finish data reading of this HC
+ }
+
+ mcm->fMCMhdCorrupted = MCM_HEADER_MASK_ERR(vword); //if MCM header mask has error
+ if (fgDebugFlag && mcm->fMCMhdCorrupted != 0) AliDebug(11,Form("Wrong MCM header mask 0x%08x.\n", *fpPos));
+
mcm->fROB = MCM_ROB_NUMBER(vword);
mcm->fMCM = MCM_MCM_NUMBER(vword);
mcm->fEvCounter = MCM_EVENT_COUNTER(vword);
-
mcm->fPos = (UInt_t*)word;
}
//
// get the adc mask
//
-
- // do it once here
UInt_t vword = *word;
- mcm->fADCindex = 0;
- mcm->fADCmax = 0;
- mcm->fADCMask = 0;
+ mcm->fADCmax = 0;
+ mcm->fADCMask = 0;
+ mcm->fADCcount = 0;
mcm->fADCMaskWord = vword;
- //memset(mcm->fADCchannel, 0, sizeof(UInt_t) * 30);
- //if ((word & 0x0000000F) == 0x0000000C)
+
+ if (vword == END_OF_TRACKLET_MARKERNEW)
+ {
+ if (fgWarnError) AliError(Form("There should be MCMadcMask. We meet END_OF_TRACKLET_MARKER 0x%08x",vword));
+ mcm->fADCmaskCorrupted += 16;
+ fHC->fCorrupted += 16; //to finish data reading of this HC
+ }
+
if ( MCM_ADCMASK_MASK_ERR(vword) == 0 )
{
- //mcm->fADCMask = (word >> 4) & skADCmaskBits;
- mcm->fADCMask = MCM_ADCMASK_VAL(vword);
+ mcm->fADCMask = MCM_ADCMASK_VAL(vword);
+ mcm->fADCcount = MCM_ADCMASK_NADC(~vword);
}
else
{
mcm->fADCMask = 0xffffffff;
- mcm->fCorrupted += 10;
- mcm->fMaskErrorCounter++;
+ mcm->fADCmaskCorrupted = 1; // mcm adc mask error
+ if (fgDebugFlag) AliDebug(11,Form("Wrong ADC Mask word 0x%08x.\n", *fpPos));
}
return mcm->fADCMask;
//
// decode the adc mask - adcs to be read out
//
-
mcm->fMCMADCWords = 0;
mcm->fSingleADCwords = 0;
mcm->fADCmax = 0;
mcm->fADCMask = GetMCMadcMask(word, mcm);
+
if (mcm->fADCMask > 0)
{
for (Int_t i = 0; i < TRD_MAX_ADC; i++)
- {
- mcm->fADCchannel[mcm->fADCmax] = 0;
- if( IS_BIT_SET(mcm->fADCMask,i) )
- {
- mcm->fADCchannel[mcm->fADCmax] = i;
- mcm->fADCmax++;
- }
- }
+ {
+ mcm->fADCchannel[mcm->fADCmax] = 0;
+ if( IS_BIT_SET(mcm->fADCMask,i) )
+ {
+ mcm->fADCchannel[mcm->fADCmax] = i;
+ mcm->fADCmax++;
+ }
+ }
+ }
+ if (mcm->fADCcount != mcm->fADCmax && fHC->fRawVMajor >= 32) // backward compatibility
+ {
+ mcm->fADCmaskCorrupted += 2;
+ if (fgDebugFlag) AliDebug(11,Form("ADC counts from ADCMask are different %d %d : ADCMask word 0x%08x\n", mcm->fADCcount, mcm->fADCmax, *fMCM->fPos));
}
}
}
//--------------------------------------------------------
-//TString DumpMCMadcMask(const struct AliTRDrawMCMInfo *mcm)
const char *AliTRDrawStreamTB::DumpMCMadcMask(const struct AliTRDrawMCM *mcm)
{
//
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff