UtcTprhData.h

Go to the documentation of this file.

#ifndef _anal_UtcTprhData_h_
#define _anal_UtcTprhData_h_

#include <iostream>


namespace anal {

  /*
    This is the part of the data that actually stays in memory.
    Move any items needed for energy corrections from class UtcTprhRawRow
    to this class UtcTprhData.
  */

  struct UtcTprhData {

    UtcTprhData() = default;

    UtcTprhData(const unsigned long int gps) : fGPSSecond(gps) { }

    void Interpolate(const UtcTprhData& g1, const UtcTprhData& g2,
                     const unsigned long int gps);

    bool operator<(const UtcTprhData& d) const { return fGPSSecond < d.fGPSSecond; }

    friend
    std::ostream&
    operator<<(std::ostream& os, const UtcTprhData& d)
    {
      return
        os << d.fGPSSecond << ' '
           << d.fTemperature << ' '
           << d.fPressure << ' '
           << d.fDensity << ' '
           << d.fDailyAvgDensity << ' '
           << d.fDensity2HoursBefore;
    }

    unsigned long int fGPSSecond = 0;
    float fTemperature = 0;  // C
    float fPressure = 0;  // hPa
    float fDensity = 0;  // kg / m^3
    float fDailyAvgDensity = 0;  // kg / m^3
    float fDensity2HoursBefore = 0;  // kg / m^3

  };


  //


  /*
    This is the format of the line in utctprh.dat.bz2
    file that can be obtained from Hernan Asorey's web page
    http://auger.uis.edu.co/data/private.html

    http://auger.uis.edu.co/data/files/v2r0/utctprh.dat.bz2
  */

  struct UtcTprhRawRow : public UtcTprhData {

    void Clear() { *this = UtcTprhRawRow(); }

    friend
    std::istream&
    operator>>(std::istream& is, UtcTprhRawRow& g)
    {
      is >> g.fUTCSecond >> std::ws
         >> g.fTemperature >> std::ws
         >> g.fPressure >> std::ws
         >> g.fDensity >> std::ws
         >> g.fDailyAvgDensity >> std::ws
         >> g.fN6T5Hexagons >> std::ws
         >> g.fN5T5Hexagons >> std::ws
         >> g.fQuality >> std::ws
         >> g.fBadPeriod >> std::ws
         >> g.fUndocumentedQuantity1 >> std::ws
         >> g.fDensity2HoursBefore >> std::ws
         >> g.fDailyAvgDensityCentered;
      const unsigned long int unixGpsDiff =
        315964800ul - 25ul; // leap seconds as in 2014
      if (is && g.fUTCSecond)
        g.fGPSSecond = g.fUTCSecond - unixGpsDiff;
      return is;
    }

    unsigned long int fUTCSecond = 0;
    float fN6T5Hexagons = 0;
    float fN5T5Hexagons = 0;
    unsigned short fQuality = 0; // 0 = orig,
                                 // 1 = interpolated,
                                 // 2 = other station,
                                 // 3 = moni files,
                                 // 4 = long gap
    unsigned short fBadPeriod = 0;
    //
    float fUndocumentedQuantity1 = 0;
    float fDailyAvgDensityCentered = 0;  // kg / m^3

  };

}


#endif