ShapingWavelet.h

#ifndef __SHAPING_WAVELET_H__
#define __SHAPING_WAVELET_H__
#include "mspass/algorithms/deconvolution/ComplexArray.h"
#include "mspass/seismic/CoreTimeSeries.h"
#include "mspass/utility/Metadata.h"
#include <boost/archive/text_iarchive.hpp>
#include <boost/archive/text_oarchive.hpp>
#include <boost/serialization/vector.hpp>
#include <utility>
namespace mspass::algorithms::deconvolution {
class ShapingWavelet {
public:
  ShapingWavelet() : nfft(0), dt(-1), df(-1) {};
  ShapingWavelet(const mspass::utility::Metadata &md, int npts = 0);
  ShapingWavelet(mspass::seismic::CoreTimeSeries d, int nfft = 0);
  ShapingWavelet(const double fpeak, const double dtin, const int n);
  ShapingWavelet(const int npolelo, const double f3dblo, const int npolehi,
                 const double f3dbhi, const double dtin, const int n);
  ShapingWavelet(const ShapingWavelet &parent);
  ShapingWavelet &operator=(const ShapingWavelet &parent);
  void swap(ShapingWavelet &other) noexcept {
    std::swap(nfft, other.nfft);
    w.swap(other.w);
    std::swap(dt, other.dt);
    std::swap(df, other.df);
    wavelet_name.swap(other.wavelet_name);
  }
  ComplexArray *wavelet() { return &w; };
  mspass::seismic::CoreTimeSeries impulse_response();
  double freq_bin_size() { return df; };
  double sample_interval() { return dt; };
  std::string type() { return wavelet_name; };
  int size() const { return w.size(); };
 
private:
  int nfft;
  ComplexArray w;
  double dt, df;
  std::string wavelet_name;
  friend boost::serialization::access;
  template <class Archive>
  void serialize(Archive &ar, const unsigned int version) {
    ar & nfft;
    ar & dt;
    ar & df;
    ar & wavelet_name;
    ar & w;
  }
};
} // namespace mspass::algorithms::deconvolution
#endif