REASH5File.h

Go to the documentation of this file.

#ifndef _io_REASH5File_h_
#define _io_REASH5File_h_

#include <io/REASH5IOException.h>
#include <io/VEventFile.h>
#include <fstream>
#include <vector>
#include <string>
#include "H5Cpp.h"
#include <limits>


namespace evt {
  class Event;
  class ShowerSimData;
}

namespace utl {
  class Branch;
}

namespace io {

  class CorsikaShowerFile;
  //class REASH5Header;
  //class REASH5Shower;


  class REASH5File : public VEventFile {

  public:
    REASH5File() { }
    REASH5File(const std::string& fileName, const Mode mode = eRead,
               utl::Branch* const b = nullptr);

    virtual ~REASH5File();

    void Open(const std::string& fileName, const Mode mode = eRead,
              utl::Branch* const b = nullptr) override;

    void Close() override;

    Status Read(evt::Event& event) override;

    void Write(const evt::Event& event) override;

    Status FindEvent(const unsigned int eventId) override;

    Status GotoPosition(const unsigned int position) override;

    int GetNEvents() override;

    //REASH5FileParticleIterator* GetREASH5FileParticleIterator() const;

  private:
    unsigned int fCurrentPosition = 0;
    std::vector<std::string> fBinFileList;
    H5::H5File fReasH5;
    std::string fOrigDirectory;
    std::string fReasDirectory;

    CorsikaShowerFile* fCorsikaShowerFile = nullptr;

    template<typename T>
    static
    T
    ReadAttribute(H5::Group group, const std::string& attributeName)
    {
      if (group.attrExists(attributeName)) {
        T tmp = T();
        H5::Attribute attr = group.openAttribute(attributeName);
        H5::DataType type = attr.getDataType();
        attr.read(type, &tmp);
        return tmp;
      }
      throw utl::NonExistentComponentException(attributeName + "does not exist");
    }

    static herr_t AddObserver(hid_t loc_id, const char* const name, void* const operator_data);

  };

}


#endif