First version of the PHOS/EMCAL ALTRO decoder which derives from the new AliAltroRawS...

[u/mrichter/AliRoot.git] / RAW / readCaloV3.C

diff --git a/RAW/readCaloV3.C b/RAW/readCaloV3.C
new file mode 100644 (file)
index 0000000..a423465
--- /dev/null
+++ b/RAW/readCaloV3.C
@@ -0,0 +1,97 @@
+#if !defined(__CINT__) || defined(__MAKECINT__)
+  #include <TStopwatch.h>
+  #include <TStyle.h>
+  #include <TH1F.h>
+  #include <TH2F.h>
+  #include <TString.h>
+  #include <TCanvas.h>
+  #include "AliRawReader.h"
+  #include "AliCaloRawStreamV3.h"
+  #include "AliLog.h"
+#endif
+
+
+void readCaloV3(const char *fileName = "./", const TString calo="PHOS")
+{
+  TH1D *hAmplHG  = new TH1D("hAmplHG" ,"HG amplitude" ,1024,-0.5,1023.5);
+  TH1D *hAmplLG  = new TH1D("hAmplLG" ,"LG amplitude" ,1024,-0.5,1023.5);
+  TH1D *hAmplTRU = new TH1D("hAmplTRU","TRU amplitude",1024,-0.5,1023.5);
+  TH1F *hModule  = new TH1F("hModule" ,"Module number",10,-0.5,9.5);
+  TH2F *hXZHG    = new TH2F("hXZHG"   ,"XZ HG cells",64,-0.5,63.5,64,-0.5,63.5);
+  TH2F *hXZLG    = new TH2F("hXZLG"   ,"XZ LG cells",64,-0.5,63.5,64,-0.5,63.5);
+  TH2F *hHWaddr  = new TH2F("hHWaddr","DDL is vs HW addr",216,-0.5,215.5,4096,-0.5,4095.5);
+
+  //  AliLog::SetGlobalLogLevel(AliLog::kFatal);
+  //  AliLog::SetPrintRepetitions(kFALSE);
+
+  AliRawReader *reader = AliRawReader::Create(fileName);
+  reader->Reset();
+
+  TStopwatch timer;
+  timer.Start();
+
+  AliCaloRawStreamV3 *stream = new AliCaloRawStreamV3(reader,calo);
+
+  Int_t iev = 0;
+
+  while (reader->NextEvent()) {
+    AliInfoGeneral("",Form("Reading event %d\n",iev++));
+    while (stream->NextDDL()) {
+      while (stream->NextChannel()) {
+//     AliInfoGeneral("",Form("New channel: HW=%d, module=%d, row=%d, col=%d, caloflag=%d",
+//                         stream->GetHWAddress(),
+//                         stream->GetModule(),
+//                         stream->GetRow(),
+//                         stream->GetColumn(),
+//                         stream->GetCaloFlag()));
+       hHWaddr->Fill(stream->GetDDLNumber(),stream->GetHWAddress());
+       hModule->Fill(stream->GetModule());
+       if (stream->IsHighGain())
+         hXZHG->Fill(stream->GetRow(),stream->GetColumn());
+       if (stream->IsLowGain())
+         hXZLG->Fill(stream->GetRow(),stream->GetColumn());
+
+       while (stream->NextBunch()) {
+         const UShort_t *sig = stream->GetSignals();
+         Int_t startBin = stream->GetStartTimeBin();
+         for (Int_t i = 0; i < stream->GetBunchLength(); i++) {
+           if (stream->IsHighGain())
+             hAmplHG->Fill(sig[i]);
+           if (stream->IsLowGain())
+             hAmplLG->Fill(sig[i]);
+           if (stream->IsTRUData())
+             hAmplTRU->Fill(sig[i]);
+         }
+       }
+      }
+    }
+    stream->Reset();
+  }
+
+  timer.Stop();
+  timer.Print();
+
+  gStyle->SetPalette(1);
+  TCanvas *c1 = new TCanvas("c1","",0,0,900,600);
+  c1->Divide(3,1);
+  c1->cd(1);
+  hAmplHG->Draw();
+  c1->cd(2);
+  hAmplLG->Draw();
+  c1->cd(3);
+  hAmplTRU->Draw();
+
+  TCanvas *c2 = new TCanvas("c2","",0,0,900,600);
+  c2->Divide(2,1);
+  c2->cd(1);
+  hModule->Draw();
+  c2->cd(2);
+  hHWaddr->Draw("colz");
+
+  TCanvas *c3 = new TCanvas("c3","",0,0,900,600);
+  c3->Divide(2,1);
+  c3->cd(1);
+  hXZHG->Draw("colz");
+  c3->cd(2);
+  hXZLG->Draw("colz");
+}