ZHAireSFile.h

Go to the documentation of this file.

#ifndef _io_ZHAireSFile_h_
#define _io_ZHAireSFile_h_

#include <utl/Math.h>
#include <utl/AugerUnits.h>
#include <io/ZHAireSIOException.h>
#include <io/VEventFile.h>
#include <fstream>
#include <vector>


namespace evt {
  class Event;
  class ShowerSimData;
}

namespace io {

#ifdef HAVE_AIRES
  class AiresShowerFile;
#endif


  class ZHAireSFile : public VEventFile {

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

    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 { return 1; }  // there is currently one event in a ZHAireS file

    double
    ZHAireSAzimuthToAuger(const double airesAzimuth)
    {
      //const double kMagneticFieldDeclination = 4.233*utl::deg;
      const double kMagneticFieldDeclination = fBDeclination / utl::rad;
      return utl::NormalizeAngleMinusPiPi(airesAzimuth + 90*utl::deg - kMagneticFieldDeclination);
    }

    static
    void
    Rotatez(const double theta, const double vi[], double vf[])
    {
      const double ct = cos(theta);
      const double st = sin(theta);
      vf[0] = vi[0]*ct - vi[1]*st;
      vf[1] = vi[0]*st + vi[1]*ct;
      vf[2] = vi[2];
    }

  private:
    ZHAireSFile(const ZHAireSFile&);
    ZHAireSFile& operator=(const ZHAireSFile&);

    unsigned int fCurrentPosition = 0;
    double fBDeclination = 0;
    std::vector<std::string> fBinFileList;
    std::ifstream* fZHAireSin = nullptr;
    std::ifstream* fsry = nullptr;
    std::ifstream* fdef = nullptr;
    //for "reas reader" path finder
    std::string fOrigDirectory;
    std::string fZHAireSDirectory;
    // from REASFile.h
    std::string fHostName;
    unsigned long fProcessId;
    //TEST LOG
    std::ofstream logfile;
    
#ifdef HAVE_AIRES
    AiresShowerFile* fAiresShowerFile = nullptr;
#endif

  };

}


#endif