ovtaper.h

Go to the documentation of this file.

// include guard
#ifndef TSXX_OVTAPER_H_VERSION

#define TSXX_OVTAPER_H_VERSION \
  "TSXX_OVTAPER_H   V1.0   "

#include <iostream>
#include <list>
#include <tsxx/tsxx.h>
#include <tsxx/tapers.h>
#include <tsxx/dttimeseries.h>

namespace ts {

  namespace tapers {

    namespace ovtaper {

      struct Pick {
        Pick(): t(0), x(0) { }
        Pick(const double& offset): t(0), x(offset) { }
        double t; 
        double x; 
      }; // struct Pick

      /*----------------------------------------------------------------------*/

      inline
      bool operator<(const Pick& p1, const Pick& p2) { return (p1.x<p2.x); }
      inline
      bool operator>(const Pick& p1, const Pick& p2) { return (p1.x>p2.x); }

      /*----------------------------------------------------------------------*/

      class Picks {
        public:
          typedef std::list<Pick> Tlistofpick;
          Picks(const bool& debug): Mdebug(debug) { }
          void read(std::istream& is);
          Pick pick(const double& offset) const;
          double time(const double& offset) const;
        private:
          bool Mdebug;
          Tlistofpick Mpicks;
      }; // class Picks

    } // namespace ovtaper

    /*----------------------------------------------------------------------*/

    class OffsetVariableTaper {
      public:
        typedef double Tvalue;
        typedef ts::TDtimeseries Ttimeseries;
        typedef Ttimeseries::Tseries Tseries;
        typedef Ttimeseries::Theader Theader;
        OffsetVariableTaper(const bool& debug=false)
          : Mdebug(debug), Mvalid(false),
            Mt1(debug), Mt2(debug), Mt3(debug), Mt4(debug)
      { }

        ts::tapers::FourPoint taper(const double& offset,
                                    const double& T0,
                                    const double& T) const;

        void read(std::istream& is);

        void read(const std::string& filename);

        ovtaper::Picks t1() const { return Mt1; }
        ovtaper::Picks t2() const { return Mt2; }
        ovtaper::Picks t3() const { return Mt3; }
        ovtaper::Picks t4() const { return Mt4; }
      private:
        bool Mdebug;
        bool Mvalid;
        ovtaper::Picks Mt1, Mt2, Mt3, Mt4;
    }; // class OffsetVariableTaper

  } // namespace tapers

} // namespace ts

#endif // TSXX_OVTAPER_H_VERSION (includeguard)

/* ----- END OF ovtaper.h ----- */