Added PHOS/macros/Trigger/raw and contained code.

[u/mrichter/AliRoot.git] / PHOS / macros / Trigger / raw / AliTRUPedestalAnalysis.cxx

diff --git a/PHOS/macros/Trigger/raw/AliTRUPedestalAnalysis.cxx b/PHOS/macros/Trigger/raw/AliTRUPedestalAnalysis.cxx
new file mode 100644 (file)
index 0000000..485a241
--- /dev/null
+++ b/PHOS/macros/Trigger/raw/AliTRUPedestalAnalysis.cxx
@@ -0,0 +1,77 @@
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+// Author: Henrik Qvigstad <henrik.qvigstad@cern.ch>
+/* $Id$ */
+
+#include "AliCaloRawStreamV3.h"
+#include "AliTRUPedestalOutput.h"
+#include "AliPHOSTriggerRawReader.h"
+#include "AliPHOSTRURawReader.h"
+#include "TH1I.h"
+
+#include "AliTRUPedestalAnalysis.h"
+#include <TRandom.h>
+
+AliTRUPedestalAnalysis::AliTRUPedestalAnalysis()
+: fOutput(new AliTRUPedestalOutput),
+  fTriggerReader(new AliPHOSTriggerRawReader)
+{
+}
+
+AliTRUPedestalAnalysis::~AliTRUPedestalAnalysis()
+{
+  delete fOutput;
+}
+
+void AliTRUPedestalAnalysis::ProcessEvent ( AliCaloRawStreamV3* phosStream )
+{
+  fTriggerReader->Reset();
+
+  phosStream->Reset();
+  while (phosStream->NextDDL()) {
+    while (phosStream->NextChannel()) {
+      if (phosStream->IsTRUData()) {
+       fTriggerReader->ReadFromStream(phosStream);
+      }// IsTRUData
+    }// NextChannel
+  }//NextDDL
+
+  bool data_in_event = false;
+  for (unsigned int mod=0; mod<kNMods; ++mod) {
+    for (unsigned int row=0; row<kNTRURows; ++row) {
+      for (unsigned int branch=0; branch<kNBranches; ++branch) {
+       // if there is tru data in this branch, for this event.
+       if( fTriggerReader->GetTRU(mod, row, branch)->IsActive() ) {
+         data_in_event = true;
+         for (unsigned int xrow=0; xrow<kN2x2XPrTRURow; ++xrow) {
+           for (unsigned int zcol=0; zcol<kN2x2ZPrBranch; ++zcol) {
+             // if there is tru data for this timebin, for this branch.
+             for (unsigned int timeBin=0; timeBin < kNTRUTimeBins; ++timeBin) {
+               if( fTriggerReader->GetTRU(mod, row, branch)->IsActive(timeBin) ){
+                 double signal = fTriggerReader->GetTRU(mod, row, branch)->GetTriggerSignal(xrow, zcol, timeBin);
+                 fOutput->GetTRUSignals(mod, row, branch, xrow, zcol)->Fill(signal);
+               } // if timebin is active
+             } // timeBin
+           } // zcol
+         } // xrow
+       } // if branch is active
+      } // branch
+    } // row
+  } // mod
+
+  if( data_in_event )
+    fOutput->EventAdded();
+}
+
+
+UInt_t AliTRUPedestalAnalysis::Get2x2Max ( AliPHOSTriggerRawReader* reader, int mod, int row, int branch, int x, int z )
+{
+  UInt_t max = 0;
+  for(UInt_t timeBin = 0; timeBin < kNTRUTimeBins; ++timeBin) {
+    const UInt_t signal = reader->GetTRU(mod, row, branch)->GetTriggerSignal(x, z, timeBin);
+    if( max < signal )
+      max = signal;
+  }
+  return max;
+}