SELFASFile.h

Go to the documentation of this file.

#ifndef _io_SELFASFile_h_
#define _io_SELFASFile_h_

#include <io/SELFASIOException.h>
#include <io/VEventFile.h>
#include <dirent.h>
#include <fstream>
#include <vector>
#include <utl/Vector.h>
#include <revt/Station.h>
#include <utl/VParticleProperties.h>
#include <cmath>
#include <iostream>
#include <iomanip>
#include <sstream>
#include <string>
#include <boost/tokenizer.hpp>
#include <boost/filesystem/path.hpp>


namespace evt {
  class Event;
  class ShowerSimData;
}

namespace io {

  class SELFASFile : public VEventFile {

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

    void Write(const evt::Event& theEvent) override;
    Status Read(evt::Event& theEvent) override;
    Status FindEvent(const unsigned int eventId) override;
    int GetNEvents() override { return 1; }
    void Close() override;
    void Open(const std::string& fileName, const Mode mode, utl::Branch* const b = nullptr) override;
    Status GotoPosition(const unsigned int pos) override;

  private: 
    template<class T> bool from_string(T& t, const std::string& s, std::ios_base &(*f)(std::ios_base&));
    std::vector<std::string> StringSplitter(const std::string& line, const std::string& sep);

    // logistics for reading the file
    unsigned int fCurrentPosition = 0;
    std::string FileName;
    std::string FileType;
    std::string name;
    std::string path_name;
    std::string version_string;
    std::string file_extension;
    std::ifstream finput_file;
    bool file_open = false;
    unsigned int events_remaining = 1;
    
    // This is all data that is read from the .selfas file.
    //utl::Vector3D core_position; // unused. LN.
    //utl::Vector3D refcoord; // unused. LN.
    std::string fprimary_particle; 
    double fprimary_energy;
    std::string fparticle_id;
    //bool fhas_magnetic_field = false; // unused. LN.
    double fzenith_angle;
    double fazimuth_angle;
    double fx_core_position;
    double fy_core_position;
    double fz_core_position;
    double ffirst_interaction_lenght;
    double ftime_step;
    double fsimulated_particles;
    bool fhas_charge_excess = false;
    bool fhas_multiple_scattering = false;
    double charge_excess;
    double de_dx;
    double multiple_scattering;
    double Cerenkov;
    bool fhas_dE_dx = false;
    int fantennas_number;
    double fx_core_relative_position;
    double fy_core_relative_position;
    double fz_core_relative_position;
    double fxefield;
    double fyefield;
    double fzefield;
    double fB0;
    double fBx;
    double fBy;
    double fBz;
    double fSiteGround_m;
    double fSiteGround_gcm2;
    double fxmax;
    double fxinflexion;
    bool fhas_Cerenkov;

  };

}


#endif