FileEventStream.cc

Go to the documentation of this file.

#include "FileEventStream.hh"

FileEventStream::FileEventStream(char * filename, int nTubes)
  :NumberOfTubes(nTubes)
{
  myFileStream = fopen(filename,"r");
  if(myFileStream==NULL)
    {
      fprintf (stderr,"Could not open input file %s, exiting.", filename);
      exit(1);
    }
  fileEmpty = false;
  fscanf(myFileStream,"%d",&upcomingEventID); //Read in the first event number, to make the event processing work.
}

FileEventStream::~FileEventStream()
{
  fclose(myFileStream);
}

Event * FileEventStream::getNextEvent()
{
  Event * nextEvent = new Event();
  if(fileEmpty)
    return NULL;
  nextEvent->eventno = upcomingEventID;
  int currentHitPtr=0;
  do
    {
      if(nextEvent->eventno!=upcomingEventID) //New event, abort!
        {
          return nextEvent;
        }
      
      nextEvent->hit.push_back(new EventHit(NumberOfTubes));

      nextEvent->hit.back()->eventno = upcomingEventID;

      fscanf(myFileStream, "%d %d", &nextEvent->hit.back()->Z, &nextEvent->hit.back()->A);
      
      fscanf(myFileStream, "%le", &nextEvent->hit.back()->time);
      
      fscanf(myFileStream,"%e %e %e %e", &nextEvent->hit.back()->UFSP, &nextEvent->hit.back()->LFSP, &nextEvent->hit.back()->BSP, &nextEvent->hit.back()->Ge);
      
      for(int j = 0; j<NumberOfTubes; j++)
        {
          fscanf(myFileStream,"%e",&nextEvent->hit.back()->DET[j]);
        }
      ++currentHitPtr;
      if(feof(myFileStream) || fscanf(myFileStream, "%d",&upcomingEventID)<1 )
        {
          fileEmpty = true;
          return nextEvent;
        }
    } while(upcomingEventID==nextEvent->eventno);

   return nextEvent;
}

bool FileEventStream::hasMoreEvents()
{
  return !fileEmpty;
}