Seismogram.h

#ifndef _SEISMOGRAM_H_
#define _SEISMOGRAM_H_
#include "mspass/seismic/CoreSeismogram.h"
#include "mspass/utility/ErrorLogger.h"
#include "mspass/utility/ProcessingHistory.h"
 
namespace mspass::seismic {
class Seismogram : public mspass::seismic::CoreSeismogram,
                   public mspass::utility::ProcessingHistory {
public:
  Seismogram()
      : mspass::seismic::CoreSeismogram(),
        mspass::utility::ProcessingHistory() {};
  Seismogram(const size_t nsamples);
  Seismogram(const mspass::seismic::CoreSeismogram &d);
  Seismogram(const mspass::seismic::CoreSeismogram &d, const std::string alg);
  Seismogram(const mspass::seismic::BasicTimeSeries &bts,
             const mspass::utility::Metadata &md);
 
  Seismogram(const mspass::seismic::BasicTimeSeries &b,
             const mspass::utility::Metadata &m,
             const mspass::utility::ProcessingHistory &his, const bool card,
             const bool ortho, const mspass::utility::dmatrix &tm,
             const mspass::utility::dmatrix &uin);
  Seismogram(const Metadata &md,
             const std::string jobname = std::string("test"),
             const std::string jobid = std::string("UNDEFINED"),
             const std::string readername = std::string("load3C"),
             const std::string algid = std::string("0"));
 
  Seismogram(const Metadata &md, bool load_data)
      : mspass::seismic::CoreSeismogram(md, load_data),
        mspass::utility::ProcessingHistory() {};
  Seismogram(const Seismogram &parent)
      : mspass::seismic::CoreSeismogram(parent),
        mspass::utility::ProcessingHistory(parent) {};
  virtual ~Seismogram() {};
  Seismogram &operator=(const Seismogram &parent);
  Seismogram &operator+=(const Seismogram &d) {
    dynamic_cast<CoreSeismogram &>(*this) +=
        dynamic_cast<const CoreSeismogram &>(d);
    this->elog += d.elog;
    return (*this);
  };
  Seismogram &operator+=(const CoreSeismogram &d) {
    dynamic_cast<CoreSeismogram &>(*this) += d;
    return (*this);
  };
  const Seismogram operator+(const Seismogram &other) const {
    Seismogram result(*this);
    result += other;
    return result;
  };
  const CoreSeismogram operator+(const CoreSeismogram &other) const {
    CoreSeismogram result(dynamic_cast<const CoreSeismogram &>(*this));
    result += other;
    return result;
  };
  Seismogram &operator*=(const double scale) {
    dynamic_cast<CoreSeismogram &>(*this) *= scale;
    return *this;
  };
  Seismogram &operator-=(const Seismogram &d) {
    dynamic_cast<CoreSeismogram &>(*this) -=
        dynamic_cast<const CoreSeismogram &>(d);
    this->elog += d.elog;
    return (*this);
  };
  Seismogram &operator-=(const CoreSeismogram &d) {
    dynamic_cast<CoreSeismogram &>(*this) -= d;
    return (*this);
  };
  const Seismogram operator-(const Seismogram &other) const {
    Seismogram result(*this);
    result -= other;
    return result;
  };
  const CoreSeismogram operator-(const CoreSeismogram &other) const {
    CoreSeismogram result(dynamic_cast<const CoreSeismogram &>(*this));
    result -= other;
    return result;
  };
  void load_history(const mspass::utility::ProcessingHistory &h);
  size_t memory_use() const;
};
} // namespace mspass::seismic
#endif