Class to produce PHOS digits from raw data stream.

diff --git a/PHOS/AliPHOSRawDigiProducer.cxx b/PHOS/AliPHOSRawDigiProducer.cxx
new file mode 100644 (file)
index 0000000..a3e60e7
--- /dev/null
+++ b/PHOS/AliPHOSRawDigiProducer.cxx
@@ -0,0 +1,74 @@
+#include "TClonesArray.h"
+
+#include "AliPHOSRawDigiProducer.h"
+#include "AliPHOSRawDecoder.h"
+#include "AliPHOSGeometry.h"
+#include "AliPHOSDigit.h"
+
+ClassImp(AliPHOSRawDigiProducer)
+
+void AliPHOSRawDigiProducer::MakeDigits(TClonesArray *digits, AliPHOSRawDecoder* decoder) 
+{
+  digits->Clear();
+ 
+  AliPHOSGeometry* geo = AliPHOSGeometry::GetInstance();
+  if(!geo) geo = AliPHOSGeometry::GetInstance("IHEP");
+
+  Int_t    iDigit   = 0 ;
+  Double_t time     = 0. ;
+  Int_t    iOldDigit;
+  Bool_t   seen,lowGainFlag;
+  Int_t relId[4], absId =0;
+
+  while (decoder->NextDigit()) {
+
+    lowGainFlag = decoder->IsLowGain();
+    time = decoder->GetTime();
+
+    relId[0] = decoder->GetModule();
+    relId[1] = 0;
+    relId[2] = decoder->GetRow();
+    relId[3] = decoder->GetColumn();
+    geo->RelToAbsNumbering(relId, absId);
+
+    // Add low gain digit only
+    //if the high gain digit does not exist in the digits array
+
+    seen = kFALSE;
+
+    if(lowGainFlag) {
+      for (iOldDigit=iDigit-1; iOldDigit>=0; iOldDigit--) {
+        if ((dynamic_cast<AliPHOSDigit*>(digits->At(iOldDigit)))->GetId() == absId) {
+          seen = kTRUE;
+          break;
+        }
+      }
+      if (!seen) {
+        new((*digits)[iDigit]) AliPHOSDigit(-1,absId,(Float_t)decoder->GetEnergy(),time);
+        iDigit++;
+      }
+    }
+
+    // Add high gain digit only if it is not saturated;
+    // replace low gain digit by a high gain one
+    else {
+      if (decoder->GetEnergy() >= 1023) continue;
+      for (iOldDigit=iDigit-1; iOldDigit>=0; iOldDigit--) {
+        if ((dynamic_cast<AliPHOSDigit*>(digits->At(iOldDigit)))->GetId() == absId) {
+          digits->RemoveAt(iOldDigit);
+          new((*digits)[iOldDigit]) AliPHOSDigit(-1,absId,(Float_t)decoder->GetEnergy(),time);
+          seen = kTRUE;
+          break;
+        }
+      }
+      if (!seen) {
+        new((*digits)[iDigit]) AliPHOSDigit(-1,absId,(Float_t)decoder->GetEnergy(),time);
+        iDigit++;
+      }
+    }
+
+  }
+
+  digits->Compress();
+  digits->Sort();
+}

diff --git a/PHOS/AliPHOSRawDigiProducer.h b/PHOS/AliPHOSRawDigiProducer.h
new file mode 100644 (file)
index 0000000..eda2a0f
--- /dev/null
+++ b/PHOS/AliPHOSRawDigiProducer.h
@@ -0,0 +1,20 @@
+#ifndef ALIPHOSRAWDIGIPRODUCER_H
+#define ALIPHOSRAWDIGIPRODUCER_H
+
+class AliPHOSRawDecoder;
+
+class AliPHOSRawDigiProducer {
+
+public:
+
+  AliPHOSRawDigiProducer() {}
+  AliPHOSRawDigiProducer(const AliPHOSRawDigiProducer& );
+  AliPHOSRawDigiProducer& operator = (const AliPHOSRawDigiProducer& );
+  virtual ~AliPHOSRawDigiProducer() {}
+
+  virtual void MakeDigits(TClonesArray *digits, AliPHOSRawDecoder* decoder);
+
+  ClassDef(AliPHOSRawDigiProducer,1)
+};
+
+#endif

diff --git a/PHOS/PHOSbaseLinkDef.h b/PHOS/PHOSbaseLinkDef.h
index 8ce8d693ca7f058f8018f4fe1866bb8e4158a058..808b82454fe53ea7fbe1c6b2afc8fca8cdd5fc16 100644 (file)
--- a/PHOS/PHOSbaseLinkDef.h
+++ b/PHOS/PHOSbaseLinkDef.h
@@ -28,5 +28,6 @@
 #pragma link C++ class AliPHOSEmcCalibData+;
 #pragma link C++ class AliPHOSPulseGenerator+;
 #pragma link C++ class AliPHOSRawDecoder+;
+#pragma link C++ class AliPHOSRawDigiProducer+;
 
 #endif

diff --git a/PHOS/libPHOSbase.pkg b/PHOS/libPHOSbase.pkg
index 3bd306cb4531bc6042ebb3a5fe5f6d82ec6f7c7f..4eda387ccc97574800b102698290cb5dde1c90d5 100644 (file)
--- a/PHOS/libPHOSbase.pkg
+++ b/PHOS/libPHOSbase.pkg
@@ -14,7 +14,8 @@ SRCS          =  AliPHOS.cxx \
                  AliPHOSCpvCalibData.cxx \
                  AliPHOSEmcCalibData.cxx \
                 AliPHOSPulseGenerator.cxx \
-                AliPHOSRawDecoder.cxx
+                AliPHOSRawDecoder.cxx \
+                AliPHOSRawDigiProducer.cxx
 
 HDRS:= $(SRCS:.cxx=.h)