Macro to read PHOS raw data and print the PHOS trigger events

diff --git a/PHOS/macros/RawProcessing/ReadRawCaloTrigger.C b/PHOS/macros/RawProcessing/ReadRawCaloTrigger.C
new file mode 100644 (file)
index 0000000..a7174f7
--- /dev/null
+++ b/PHOS/macros/RawProcessing/ReadRawCaloTrigger.C
@@ -0,0 +1,125 @@
+#if !defined(__CINT__) || defined(__MAKECINT__)
+  #include <TStopwatch.h>
+  #include <TGrid.h>
+  #include <TString.h>
+  #include <TMath.h>
+
+  #include "AliRawReader.h"
+  #include "AliCaloRawStreamV3.h"
+  #include "AliLog.h"
+  #include "AliPHOSRawFitterv1.h"
+  #include "AliCentralTrigger.h"
+  #include "AliTriggerConfiguration.h"
+  #include "AliTriggerClass.h"
+  #include "AliCDBManager.h"
+#endif
+
+//-----------------------------------------------------------------------------
+static Bool_t              firstEvent = kTRUE;
+static Int_t               runNum;
+static UInt_t              period;
+static UInt_t              orbitID;
+static UInt_t              bcID;
+static AliRawReader       *reader;
+TString GetTriggerClass(ULong64_t);
+
+//-----------------------------------------------------------------------------
+void ReadRawCaloTrigger(const TString fileName = "./", const TString calo="PHOS")
+{
+
+  if (fileName.BeginsWith("alien://")) {
+    TGrid::Connect("alien://");
+  }
+
+  AliCDBManager *man = AliCDBManager::Instance();
+  man->SetDefaultStorage("raw://");
+
+  reader = AliRawReader::Create(fileName);
+  reader ->Reset();
+  AliCaloRawStreamV3 *stream = new AliCaloRawStreamV3(reader,calo);
+
+  Int_t iev = 0;
+
+  TStopwatch stopwatch;
+  stopwatch.Start();
+  Int_t    module,cellX,cellZ,caloFlag;
+
+  while (reader->NextEvent()) {
+    if (firstEvent) {
+      firstEvent = kFALSE;
+      runNum = reader->GetRunNumber();
+      man = AliCDBManager::Instance();
+      man ->SetRun(runNum);
+    }
+    ULong64_t triggerMask  = reader->GetClassMask();
+    TString trclasses = GetTriggerClass(triggerMask);
+    
+    period  = reader->GetPeriod();
+    orbitID = reader->GetOrbitID();
+    bcID    = reader->GetBCID();
+    iev++;
+    if (!trclasses.Contains("CPHI1-B-NOPF-ALLNOTRD")) continue;
+
+    AliInfoGeneral("",Form("Reading event %d of type %d, time %d, trig.class \"%s\"",
+                          iev,reader->GetType(), reader->GetTimestamp(), trclasses.Data()));
+
+    while (stream->NextDDL()) {
+      while (stream->NextChannel()) {
+       // Get geometry of the channel
+       module   = stream->GetModule();
+       cellX    = stream->GetCellX();
+       cellZ    = stream->GetCellZ();
+       caloFlag = stream->GetCaloFlag();
+
+       Int_t nBunches = 0;
+       while (stream->NextBunch()) {
+         const UShort_t *sig = stream->GetSignals();
+         Int_t sigLength = stream->GetBunchLength();
+         Int_t maxAmp  = 0;
+         for (Int_t i = 0; i < sigLength; i++) {
+           if (sig[i] >  maxAmp) {
+             maxAmp = sig[i];
+           }
+         }
+         printf("(m,x,z,c)=(%d,%2d,%2d,%d), bunch=%d, amplutide=%d\n",
+                module,cellX,cellZ,caloFlag,nBunches,maxAmp);
+         nBunches++;
+       }
+      }
+    }
+    stream->Reset();
+  }
+  stopwatch.Print();
+
+}
+//-----------------------------------------------------------------------------
+
+TString GetTriggerClass(ULong64_t triggerMask)
+{
+  // Convert a trigger mask to a trigger class
+  
+  AliCentralTrigger aCTP;
+  TString configstr("");
+  TString trclasses;
+  if (!aCTP.LoadConfiguration(configstr)) { // Load CTP config from OCDB
+    AliInfoGeneral("","No trigger configuration found in OCDB! The trigger configuration information will not be used!");
+    return trclasses;
+  }
+  aCTP.SetClassMask(triggerMask);
+  AliTriggerConfiguration *config = aCTP.GetConfiguration();
+  const TObjArray& classesArray = config->GetClasses();
+  Int_t nclasses = classesArray.GetEntriesFast();
+  for( Int_t iclass=0; iclass < nclasses; iclass++ ) {
+    AliTriggerClass* trclass = (AliTriggerClass*)classesArray.At(iclass);
+    if (trclass && trclass->GetMask()>0) {
+      Int_t trindex = TMath::Nint(TMath::Log2(trclass->GetMask()));
+      reader->LoadTriggerClass(trclass->GetName(),trindex);
+      if (triggerMask & (1ull << trindex)) {
+       trclasses += " ";
+       trclasses += trclass->GetName();
+       trclasses += " ";
+      }
+    }
+  }
+  return trclasses;
+}