1st version of EMCAL STU decoding (not yet with final format)

diff --git a/EMCAL/AliEMCALSTURawStream.cxx b/EMCAL/AliEMCALSTURawStream.cxx
new file mode 100644 (file)
index 0000000..2a7d278
--- /dev/null
+++ b/EMCAL/AliEMCALSTURawStream.cxx
@@ -0,0 +1,129 @@
+/**************************************************************************
+ * 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.                  *
+ **************************************************************************/
+
+///////////////////////////////////////////////////////////////////////////////
+///
+/// This class provides access to EMCALSTU DDL raw data.
+///
+///////////////////////////////////////////////////////////////////////////////
+
+#include "AliEMCALSTURawStream.h"
+#include "AliRawReader.h"
+#include "AliLog.h"
+
+ClassImp(AliEMCALSTURawStream)
+
+//_____________________________________________________________________________
+AliEMCALSTURawStream::AliEMCALSTURawStream(AliRawReader* rawReader) :
+  fRawReader(rawReader),
+  fNum2x2Words(0)
+{
+  // reset arrays
+  memset(fJetPatchWords, 0, sizeof(fJetPatchWords));
+
+  // Reset RawStream
+  //
+  fRawReader->Reset();
+
+  // select the raw data corresponding to
+  // the EMCALSTU detector id
+  fRawReader->Select("EMCAL", kEMCALSTUDDL);
+}
+
+//_____________________________________________________________________________
+AliEMCALSTURawStream::AliEMCALSTURawStream(const AliEMCALSTURawStream& stream) :
+  TObject(stream),
+  fRawReader(NULL),
+  fNum2x2Words(0)
+{
+  // Copy constructor
+  AliFatal("Copy constructor not implemented");
+}
+
+//_____________________________________________________________________________
+AliEMCALSTURawStream& AliEMCALSTURawStream::operator = (const AliEMCALSTURawStream& source
+                                                       /* stream */)
+{
+  // assignment operator; use copy ctor
+  if (&source == this) return *this;
+
+  new (this) AliEMCALSTURawStream(source);
+  return *this;
+}
+
+//_____________________________________________________________________________
+AliEMCALSTURawStream::~AliEMCALSTURawStream()
+{
+  // destructor
+}
+
+//_____________________________________________________________________________
+void AliEMCALSTURawStream::Reset()
+{
+  // reset raw stream params
+  if (fRawReader) fRawReader->Reset();
+}
+
+//_____________________________________________________________________________
+Bool_t AliEMCALSTURawStream::Next()
+{
+  // read the whole EMCALSTU raw data stream
+  // return kFALSE in case of error
+
+  UChar_t *data = NULL;
+
+  // EMCALSTU raw data should contain CDH so we don't need any mods a la
+  //  fRawReader->RequireHeader(kTRUE);
+
+  // let's pick up the data block, if we can
+  if (!fRawReader->ReadNextData(data)) {
+    return kFALSE;
+  }
+
+  // check that the amount of data we picked up is what we expected
+  // 20080721 (DS): the rest of this method needs to be verified as the exact data format 
+  // is defined
+  int nwordsExtra = fRawReader->GetDataSize() - (kNumJetPatchWords + kNumGammaJetPatchWords);
+  if ( nwordsExtra < 0 ) {
+    AliError(Form("Wrong EMCALSTU raw data size: %d", fRawReader->GetDataSize()));
+    return kFALSE;
+  }
+
+  // unpack the data into our local arrays
+  // jet-patch
+  int iword = 0;
+  int ioffset = 0;
+  for (iword = 0; iword<kNumJetPatchWords; iword++) {
+    fJetPatchWords[iword] = (data[iword] & 0xffffffff);
+  }
+  // gamma-jet patch
+  ioffset += kNumJetPatchWords;
+  for (iword = 0; iword<kNumGammaJetPatchWords; iword++) {
+    fGammaJetPatchWords[iword] = (data[iword+ioffset] & 0xffffffff);
+  }
+
+  // also 2x2sum's, if there appears to be enough/more data..
+  fNum2x2Words = 0;
+  if ( nwordsExtra > kNum2x2Words ) {
+    ioffset += kNumGammaJetPatchWords;
+    for (iword = 0; iword<kNum2x2Words; iword++) {
+      f2x2Words[iword] = (data[iword+ioffset] & 0xffffffff);
+    }
+    fNum2x2Words = kNum2x2Words;
+  }
+
+  return kTRUE;
+}
+

diff --git a/EMCAL/AliEMCALSTURawStream.h b/EMCAL/AliEMCALSTURawStream.h
new file mode 100644 (file)
index 0000000..4b401ad
--- /dev/null
+++ b/EMCAL/AliEMCALSTURawStream.h
@@ -0,0 +1,61 @@
+#ifndef ALIEMCALSTURAWSTREAM_H
+#define ALIEMCALSTURAWSTREAM_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+///////////////////////////////////////////////////////////////////////////////
+///
+/// This is a class for reading the EMCAL STU (trigger) DDL raw data
+///
+///  Two different formats foreseen (different by their length and Format Version):
+///
+///  - The regular one, containing jet patch indices (7x 32 bit words) and
+///  gamma jet patch indices (~96 combinations per region x 32 regions / 32
+///  bit per words = 96 words).
+///
+///  - The debug one containing the previous informations (103 words) plus
+///  all the time integrated 2x2 sums used for generating the accepted
+///  trigger (96 x 32 x 16 bit / 32 = 1536 words).
+///
+///////////////////////////////////////////////////////////////////////////////
+
+#include <TObject.h>
+
+class AliRawReader;
+
+class AliEMCALSTURawStream: public TObject {
+
+  public :
+
+    enum { kNumJetPatchWords = 7, // max number of jet patch indices words
+          kNumGammaJetPatchWords = 96, // max number of gamma-jet patch indices words
+          kNum2x2Words = 1536,     //  max number of 2x2 indices words
+          kEMCALSTUDDL = 0x12c0 }; // 4800, in EMCAL DDL range but separated from regular FEE
+
+    AliEMCALSTURawStream(AliRawReader* rawReader);
+    virtual ~AliEMCALSTURawStream();
+  
+    virtual void             Reset();
+    virtual Bool_t           Next();
+
+    UInt_t    GetJetPatchWords(int i) const { return fJetPatchWords[i]; }
+    UInt_t    GetGammaJetPatchWords(int i) const { return fGammaJetPatchWords[i]; }
+    UInt_t    Get2x2Words(int i) const { return f2x2Words[i]; }
+    UInt_t    GetNum2x2Words() const { return fNum2x2Words; }
+
+  protected:
+    AliEMCALSTURawStream(const AliEMCALSTURawStream& stream);
+    AliEMCALSTURawStream& operator = (const AliEMCALSTURawStream& stream);
+
+  private:
+
+    AliRawReader*    fRawReader;   // object for reading the raw data
+    UInt_t           fJetPatchWords[kNumJetPatchWords]; // jet patch indices
+    UInt_t           fGammaJetPatchWords[kNumGammaJetPatchWords]; // gamma jet patch indices
+    UInt_t           f2x2Words[kNum2x2Words]; // 2x2 sums
+    UInt_t           fNum2x2Words; // how many 2x2 sums did we actually read?
+
+    ClassDef(AliEMCALSTURawStream, 0)   // class for reading EMCAL STU DDL raw data
+};
+
+#endif