PMD raw digits reading

diff --git a/PMD/AliPMDRawStream.cxx b/PMD/AliPMDRawStream.cxx
new file mode 100644 (file)
index 0000000..c6e1241
--- /dev/null
+++ b/PMD/AliPMDRawStream.cxx
@@ -0,0 +1,229 @@
+/**************************************************************************
+ * 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$ */
+
+///////////////////////////////////////////////////////////////////////////////
+///
+/// This class provides access to PMD digits in raw data.
+///
+/// It loops over all PMD digits in the raw data given by the AliRawReader.
+/// The Next method goes to the next digit. If there are no digits left
+/// it returns kFALSE.
+/// Several getters provide information about the current digit.
+///
+///////////////////////////////////////////////////////////////////////////////
+
+#include "AliPMDRawStream.h"
+#include "AliRawReader.h"
+
+ClassImp(AliPMDRawStream)
+
+
+//_____________________________________________________________________________
+AliPMDRawStream::AliPMDRawStream(AliRawReader* rawReader) :
+  fRawReader(rawReader),
+  fModule(-1),
+  fPrevModule(-1),
+  fMCM(-1),
+  fChannel(-1),
+  fRow(-1),
+  fColumn(-1),
+  fSignal(-1),
+  fDetector(-1),
+  fSMN(-1)
+{
+// create an object to read PMD raw digits
+
+  fRawReader->Select(12);
+}
+
+//_____________________________________________________________________________
+AliPMDRawStream::AliPMDRawStream(const AliPMDRawStream& stream) :
+  TObject(stream),
+  fRawReader(NULL),
+  fModule(-1),
+  fPrevModule(-1),
+  fMCM(-1),
+  fChannel(-1),
+  fRow(-1),
+  fColumn(-1),
+  fSignal(-1),
+  fDetector(-1),
+  fSMN(-1)
+{
+// copy constructor
+
+  Fatal("AliPMDRawStream", "copy constructor not implemented");
+}
+
+//_____________________________________________________________________________
+AliPMDRawStream& AliPMDRawStream::operator = (const AliPMDRawStream& 
+                                             /* stream */)
+{
+// assignment operator
+
+  Fatal("operator =", "assignment operator not implemented");
+  return *this;
+}
+
+//_____________________________________________________________________________
+AliPMDRawStream::~AliPMDRawStream()
+{
+// destructor
+
+}
+
+
+//_____________________________________________________________________________
+Bool_t AliPMDRawStream::Next()
+{
+// read the next raw digit
+// returns kFALSE if there is no digit left
+
+  fPrevModule = fModule;
+
+  UInt_t data;
+  if (!fRawReader->ReadNextInt(data)) return kFALSE;
+
+  fSignal  = data & 0x0FFF;
+  fChannel = (data >> 12) & 0x003F;
+  fMCM     = (data >> 18) & 0x07FF;
+
+  Int_t  iddl  = fRawReader->GetDDLID();
+  Int_t  ium;
+  GetRowCol(iddl, fMCM, fChannel, ium, fRow, fColumn);
+
+  if (iddl < 4)
+    {
+      fModule = iddl*6 + ium;
+      fDetector = 0;
+      fSMN = iddl*6 + ium;
+    }
+  else if (iddl == 4)
+    {
+      if (ium < 6)
+       {
+         fModule = 24 + ium;
+         fSMN    = ium;
+       }
+      else if (ium >= 6)
+       {
+         fModule = 30 + ium;
+         fSMN    = 6 + ium;
+       }
+      fDetector = 1;
+    }
+  else if (iddl == 5)
+    {
+
+      if (ium < 6)
+       {
+         fModule = 30 + ium;
+         fSMN    = 6 + ium;
+       }
+      else if (ium >= 6)
+       {
+         fModule = 36 + ium;
+         fSMN    = 12 + ium;
+       }
+      fDetector = 1;
+    }
+
+
+  return kTRUE;
+}
+
+
+//_____________________________________________________________________________
+void AliPMDRawStream::GetRowCol(Int_t ddlno, UInt_t mcmno, UInt_t chno,
+                               Int_t &um, Int_t &row, Int_t &col) const
+{
+// decode: ddlno, mcmno, chno -> um, row, col
+
+
+  Int_t remmcm = 0;
+  Int_t divmcm = 0;
+
+  static const UInt_t kCh[64] = { 21, 25, 29, 28, 17, 24, 20, 16,
+                                 12, 13, 8, 4, 0, 1, 9, 5,
+                                 10, 6, 2, 3, 14, 7, 11, 15,
+                                 19, 18, 23, 27, 31, 30, 22, 26,
+                                 53, 57, 61, 60, 49, 56, 52, 48,
+                                 44, 45, 40, 36, 32, 33, 41, 37,
+                                 42, 38, 34, 35, 46, 39, 43, 47,
+                                 51, 50, 55, 59, 63, 62, 54, 58 };
+
+  if (ddlno == 0 || ddlno == 1)
+    {
+      um  = mcmno/72;
+      Int_t mcmnonew = mcmno - 72*um;
+      Int_t rowcol  = kCh[chno];
+      Int_t irownew = rowcol/4;
+      Int_t icolnew = rowcol%4;
+      
+      remmcm  = mcmnonew%12;
+      divmcm  = mcmnonew/12;
+      
+      row = 16*divmcm + irownew;
+      col =  4*remmcm + icolnew;
+      // This obtatined row and col (0,0) are the top left corner.
+      // Needs transformation to get the Geant (0,0)
+
+    }
+  else   if (ddlno == 2 || ddlno == 3)
+    {
+      um  = mcmno/72;
+      Int_t mcmnonew = mcmno - 72*um;
+      Int_t rowcol  = kCh[chno];
+      Int_t irownew = rowcol/4;
+      Int_t icolnew = rowcol%4;
+      
+      remmcm  = mcmnonew%24;
+      divmcm  = mcmnonew/24;
+      
+      row = 16*divmcm + irownew;
+      col =  4*remmcm + icolnew;
+      // This obtatined row and col (0,0) are the top left corner.
+      // Needs transformation to get the Geant (0,0)
+
+    }
+  else if (ddlno == 4 || ddlno == 5)
+    {
+      um  = mcmno/72;
+      Int_t mcmnonew = mcmno - 72*um;
+      Int_t rowcol  = kCh[chno];
+      Int_t irownew = rowcol/4;
+      Int_t icolnew = rowcol%4;
+      
+      if (um < 6)
+       {
+         remmcm  = mcmnonew%12;
+         divmcm  = mcmnonew/12;
+       }
+      else if (um >= 6)
+       {
+         remmcm  = mcmnonew%24;
+         divmcm  = mcmnonew/24;
+       }
+
+      row = 16*divmcm + irownew;
+      col =  4*remmcm + icolnew;
+      // This obtatined row and col (0,0) are the top left corner.
+      // Needs transformation to get the Geant (0,0)
+
+    }
+
+}

diff --git a/PMD/AliPMDRawStream.h b/PMD/AliPMDRawStream.h
new file mode 100644 (file)
index 0000000..7ddde3f
--- /dev/null
+++ b/PMD/AliPMDRawStream.h
@@ -0,0 +1,62 @@
+#ifndef ALIPMDRAWSTREAM_H
+#define ALIPMDRAWSTREAM_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/* $Id$ */
+
+///////////////////////////////////////////////////////////////////////////////
+///
+/// This class provides access to PMD digits in raw data.
+///
+///////////////////////////////////////////////////////////////////////////////
+
+#include <TObject.h>
+
+class AliRawReader;
+
+
+class AliPMDRawStream: public TObject {
+  public :
+    AliPMDRawStream(AliRawReader* rawReader);
+    virtual ~AliPMDRawStream();
+
+    virtual Bool_t   Next();
+
+    Int_t            GetModule() const {return fModule;};
+    Int_t            GetPrevModule() const {return fPrevModule;};
+    Bool_t           IsNewModule() const {return fModule != fPrevModule;};
+    Int_t            GetMCM() const {return fMCM;};
+    Int_t            GetChannel() const {return fChannel;};
+    Int_t            GetRow() const {return fRow;};
+    Int_t            GetColumn() const {return fColumn;};
+    Int_t            GetSignal() const {return fSignal;};
+    Int_t            GetDetector() const {return fDetector;};
+    Int_t            GetSMN() const {return fSMN;};
+
+
+    enum {kDDLOffset = 0xC00};      // offset for DDL numbers
+
+  private :
+    AliPMDRawStream(const AliPMDRawStream& stream);
+    AliPMDRawStream& operator = (const AliPMDRawStream& stream);
+
+    void             GetRowCol(Int_t ddlno, UInt_t mcmno, UInt_t chno,
+                              Int_t &um, Int_t &row, Int_t &col) const;
+
+    AliRawReader*    fRawReader;    // object for reading the raw data
+
+    Int_t            fModule;       // index of current module
+    Int_t            fPrevModule;   // index of previous module
+    Int_t            fMCM;          // index of current MCM
+    Int_t            fChannel;      // index of current channel
+    Int_t            fRow;          // index of current row
+    Int_t            fColumn;       // index of current column
+    Int_t            fSignal;       // signal in ADC counts
+    Int_t            fDetector;     // PRE = 0, CPV = 1
+    Int_t            fSMN;          // serial module number (0-23)
+
+    ClassDef(AliPMDRawStream, 0)    // class for reading PMD raw digits
+};
+
+#endif