DataGap.h

#ifndef _MSPASS_SEISMIC_DATAGAP_H_
#define _MSPASS_SEISMIC_DATAGAP_H_
#include "mspass/algorithms/TimeWindow.h"
#include <list>
#include <set>
 
namespace mspass::seismic {
class TimeWindowCmp {
public:
  bool operator()(const mspass::algorithms::TimeWindow ti1,
                  const mspass::algorithms::TimeWindow ti2) const {
    return (ti1.end < ti2.start);
  };
};
 
class DataGap {
public:
  DataGap() {};
  DataGap(const std::list<mspass::algorithms::TimeWindow> &twlist);
  DataGap(const DataGap &parent) : gaps(parent.gaps) {};
  virtual ~DataGap() {};
  bool is_gap(const double ttest); // query by time
  bool has_gap(const mspass::algorithms::TimeWindow twin);
  bool has_gap() { return (!gaps.empty()); };
  void add_gap(const mspass::algorithms::TimeWindow tw);
  std::list<mspass::algorithms::TimeWindow> get_gaps() const;
  void clear_gaps() {
    if (!gaps.empty())
      gaps.clear();
  };
  int number_gaps() const { return gaps.size(); };
  DataGap subset(const mspass::algorithms::TimeWindow tw) const;
  void translate_origin(double time_of_new_origin);
  DataGap &operator=(const DataGap &parent);
  DataGap &operator+=(const DataGap &other);
 
protected:
  std::set<mspass::algorithms::TimeWindow, mspass::algorithms::TimeWindowCmp>
      gaps;
};
} // namespace mspass::seismic
#endif // end guard