]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - RAW/AliRawReaderFile.cxx
Run number is stored as a file <run#> among the DDL files. AliRawReaderFile looks...
[u/mrichter/AliRoot.git] / RAW / AliRawReaderFile.cxx
index ad1201977c13f3def7c3dadda74887c8c07335e7..22f29c556b12f9a04ecd56b1975ab83a2a8b13cd 100644 (file)
@@ -32,7 +32,9 @@
 
 #include <cstdlib>
 #include "AliRawReaderFile.h"
+#include "AliDAQ.h"
 #include <TSystem.h>
+#include <TArrayC.h>
 
 
 ClassImp(AliRawReaderFile)
@@ -45,7 +47,14 @@ AliRawReaderFile::AliRawReaderFile(Int_t eventNumber) :
   fStream(NULL),
   fEquipmentId(-1),
   fBuffer(NULL),
-  fBufferSize(0)
+  fBufferSize(0),
+  fEquipmentSize(0),
+  fDDLIndex(NULL),
+  fDDLCurrent(-1),
+  fType(7),
+  fRunNb(0),
+  fDetectorPattern(0),
+  fTimestamp(0)
 {
 // create an object to read digits from the given event
 // in the current directory
@@ -53,6 +62,9 @@ AliRawReaderFile::AliRawReaderFile(Int_t eventNumber) :
   fDirectory = OpenDirectory();
   OpenNextFile();
   fHeader = new AliRawDataHeader;
+
+  fId[0] = fId[1] = 0;
+  fTriggerPattern[0] = fTriggerPattern[1] = 0;
 }
 
 AliRawReaderFile::AliRawReaderFile(const char* dirName, Int_t eventNumber) :
@@ -62,13 +74,23 @@ AliRawReaderFile::AliRawReaderFile(const char* dirName, Int_t eventNumber) :
   fStream(NULL),
   fEquipmentId(-1),
   fBuffer(NULL),
-  fBufferSize(0)
+  fBufferSize(0),
+  fEquipmentSize(0),
+  fDDLIndex(NULL),
+  fDDLCurrent(-1),
+  fType(7),
+  fRunNb(0),
+  fDetectorPattern(0),
+  fTimestamp(0)
 {
 // create an object to read digits from the given directory
 
   fDirectory = OpenDirectory();
   OpenNextFile();
   fHeader = new AliRawDataHeader;
+
+  fId[0] = fId[1] = 0;
+  fTriggerPattern[0] = fTriggerPattern[1] = 0;
 }
 
 AliRawReaderFile::~AliRawReaderFile()
@@ -84,8 +106,9 @@ AliRawReaderFile::~AliRawReaderFile()
 #endif
     delete fStream;
   }
-  delete fHeader;
+  if (fHeader) delete fHeader;
   if (fBuffer) delete[] fBuffer;
+  if (fDDLIndex) delete fDDLIndex; fDDLIndex=NULL;
 }
 
 void AliRawReaderFile::RequireHeader(Bool_t required)
@@ -93,8 +116,14 @@ void AliRawReaderFile::RequireHeader(Bool_t required)
   // Reading of raw data in case of missing
   // raw data header is not implemented for
   // this class
-  if (!required)
-    Fatal("AliRawReaderFile","Reading of raw data without raw data header is not implemented !");
+  if (!required) {
+    Warning("AliRawReaderFile","Reading of raw data without raw data header!");
+    if (fHeader) delete fHeader;
+    fHeader = NULL;
+  }
+  else {
+    if (!fHeader) fHeader = new AliRawDataHeader;
+  }
 
   AliRawReader::RequireHeader(required);
 }
@@ -123,11 +152,51 @@ void* AliRawReaderFile::OpenDirectory()
   return directory;
 }
 
+Bool_t AliRawReaderFile::CreateFileIndex()
+{
+// scan the files of the directory and create index of all DDL files
+// returns kFALSE if no DDL files available
+  Bool_t result=kFALSE;
+  fDDLCurrent=-1;
+  if (fDDLIndex) return fDDLIndex->GetSize()>0;
+  if (!fDirectory) return kFALSE;
+  fDDLIndex=new TArrayC(0);
+  if (!fDDLIndex) return kFALSE;
+  TString entry;
+  while (entry = gSystem->GetDirEntry(fDirectory)) {
+    const char* filename=entry.Data();
+    if (!filename || entry.IsNull()) break;
+    if (entry.BeginsWith("run")) {
+      entry.ReplaceAll("run","");
+      fRunNb = entry.Atoi();
+      continue;
+    }
+    if (!entry.EndsWith(".ddl")) continue;
+    result=kTRUE;
+    entry.Remove(0, entry.Last('_')+1);
+    entry.Remove(entry.Length()-4);
+    Int_t equipmentId = atoi(entry.Data());
+    Int_t ddlIndex = -1;
+    fDetectorPattern |= (1 << AliDAQ::DetectorIDFromDdlID(equipmentId,ddlIndex));
+    if (fDDLIndex->GetSize()<=equipmentId) {
+      fDDLIndex->Set(equipmentId+1);
+    }
+    char* array=(char*)fDDLIndex->GetArray();
+    array[equipmentId]=1;
+  }
+
+  return result;
+}
+
 Bool_t AliRawReaderFile::OpenNextFile()
 {
 // open the next file
 // returns kFALSE if the current file is the last one
 
+  if (!fDDLIndex && !CreateFileIndex()) return kFALSE;
+  if (fSelectMinEquipmentId>=0 && fSelectMinEquipmentId>fEquipmentId)
+    fDDLCurrent=fSelectMinEquipmentId-1;
+
   if (fStream) {
 #if defined(__HP_aCC) || defined(__DECCXX)
     if (fStream->rdbuf()->is_open()) fStream->close();
@@ -137,26 +206,34 @@ Bool_t AliRawReaderFile::OpenNextFile()
     delete fStream;
     fStream = NULL;
     fEquipmentId = -1;
+    fEquipmentSize = 0;
   }
 
   if (!fDirectory) return kFALSE;
-  TString entry;
-  while (entry = gSystem->GetDirEntry(fDirectory)) {
-    if (entry.IsNull()) return kFALSE;
-    if (!entry.EndsWith(".ddl")) continue;
+  while (++fDDLCurrent<(fDDLIndex->GetSize()) && 
+        (fDDLCurrent<=fSelectMaxEquipmentId || fSelectMaxEquipmentId<0)) {
+    if (fDDLIndex->At(fDDLCurrent)==0) continue;
+    Int_t dummy=0;
+    TString entry;
+    entry.Form("%s_%d.ddl", AliDAQ::DetectorNameFromDdlID(fDDLCurrent, dummy), fDDLCurrent);
     char* fileName = gSystem->ConcatFileName(GetDirName(), entry);
+    if (!fileName) continue;
+    // read the timestamp
+    FileStat_t buf;
+    if (gSystem->GetPathInfo(fileName,buf) == 0) {
+      fTimestamp = buf.fMtime;
+    }
 #ifndef __DECCXX 
     fStream = new fstream(fileName, ios::binary|ios::in);
 #else
     fStream = new fstream(fileName, ios::in);
 #endif
+    delete [] fileName;
     break;
   }
 
   if (!fStream) return kFALSE;
-  entry.Remove(0, entry.Last('_')+1);
-  entry.Remove(entry.Length()-4);
-  fEquipmentId = atoi(entry.Data());
+  fEquipmentId = fDDLCurrent;
 #if defined(__HP_aCC) || defined(__DECCXX)
   return (fStream->rdbuf()->is_open());
 #else
@@ -170,13 +247,19 @@ Bool_t AliRawReaderFile::ReadHeader()
 // read a data header at the current stream position
 // returns kFALSE if the mini header could not be read
 
-  if (!fStream) return kFALSE;
+  if (!fStream && !OpenNextFile()) return kFALSE;
   do {
     if (fCount > 0) fStream->seekg(Int_t(fStream->tellg()) + fCount);
-    while (!fStream->read((char*) fHeader, sizeof(AliRawDataHeader))) {
-      if (!OpenNextFile()) return kFALSE;
+    if (fHeader) {
+      while (!fStream->read((char*) fHeader, sizeof(AliRawDataHeader))) {
+       if (!OpenNextFile()) return kFALSE;
+      }
+    }
+    else {
+      if (fStream->eof())
+       if (!OpenNextFile()) return kFALSE;
     }
-    if (fHeader->fSize != 0xFFFFFFFF) {
+    if (fHeader && fHeader->fSize != 0xFFFFFFFF) {
       fCount = fHeader->fSize - sizeof(AliRawDataHeader);
     } else {
       UInt_t currentPos = fStream->tellg();
@@ -184,6 +267,8 @@ Bool_t AliRawReaderFile::ReadHeader()
       fCount = UInt_t(fStream->tellg()) - currentPos;
       fStream->seekg(currentPos);
     }
+    fEquipmentSize = fCount;
+    if (fHeader) fEquipmentSize += sizeof(AliRawDataHeader);
   } while (!IsSelected());
   return kTRUE;
 }
@@ -245,7 +330,20 @@ Bool_t AliRawReaderFile::Reset()
   if (fDirectory) gSystem->FreeDirectory(fDirectory);
   fDirectory = directory;
 
+  // Matthias 05.06.2008
+  // do not open the next file. That might collide with a subsequent
+  // SelectEquipment call as the 'file pointer' is already set.
+  // This is important for the indexing of the DDL files.
+  // ---------------------------------------------------------
+  // All ReadNext functions first require the fCount member to be
+  // non zero or call ReadHeader. That allows to postpone the call
+  // to OpenNextFile to the next invocation of ReadHeader.
+  // ReadHeader has been mofified according to that.
+  /*
   OpenNextFile();
+  */
+  fEquipmentId=-1;
+  fDDLCurrent=-1;
   fCount = 0;
   return kTRUE;
 }
@@ -254,6 +352,9 @@ Bool_t AliRawReaderFile::NextEvent()
 {
 // go to the next event directory
 
+  if (fDDLIndex) delete fDDLIndex;
+  fDDLIndex=NULL;
+  fDetectorPattern = 0;
   if (fEventIndex < -1) return kFALSE;
 
   do {
@@ -267,6 +368,19 @@ Bool_t AliRawReaderFile::NextEvent()
     Reset();
   } while (!IsEventSelected());
 
+  // Read the header of the first payload
+  // in order to fill the 'fake' event header
+  if (ReadHeader() && fHeader) {
+    fId[0] = ((fHeader->GetEventID2() >> 20) & 0xf);
+    fId[1] = (fHeader->GetEventID1() & 0xfff) | ((fHeader->GetEventID2() & 0xfffff) << 12);
+    fTriggerPattern[0] = (fHeader->GetTriggerClasses() & 0xffffffff);
+    fTriggerPattern[1] = ((fHeader->GetTriggerClasses() >> 32) & 0x3ffff);
+  }
+  else {
+    Warning("AliRawReaderFile","Can not read CDH header! The event header fields will be empty!");
+  }
+  Reset();
+
   fEventNumber++;
 
   return kTRUE;