psc_service.hpp

/* psc_service.hpp */
/* $Id$ */
#ifndef __psc_service_hpp
#define __psc_service_hpp

#include <vector>
#include <string>

#ifdef GE_BUILD
#include <trace.hpp>
#include <curve_smooth_sb.hpp>
#else
#include <s2proc/trace.hpp>
#include <mth/curve_smooth_sb.hpp>
#endif

namespace psc
{
  static const float floatNumEpsilon = 1.e-8;
  
  struct CorrDirMask
  {
    static const unsigned int ACF    =  (1 << 0); 
    // 2D observation system
    static const unsigned int SHOT   =  (1 << 1); 
    static const unsigned int RECV   =  (1 << 2); 
    static const unsigned int CMP    =  (1 << 3); 
    static const unsigned int OFFSET =  (1 << 4); 
    static const unsigned int COMPLETE_2D = ACF | SHOT | RECV | CMP | OFFSET; 
    
    // 3D observation system
    static const unsigned int NSS     = (1 << 5);  
    static const unsigned int NSL     = (1 << 6);  
    static const unsigned int NSS_NSL = (1 << 7);  
    static const unsigned int NRS     = (1 << 8);  
    static const unsigned int NRL     = (1 << 9);  
    static const unsigned int NRS_NRL = (1 << 10); 
    static const unsigned int COMPLETE_3D = ACF | NSS | NSL | NSS_NSL | NRS | NRL | NRS_NRL;
    
    static std::string toString(const unsigned int dir_mask);
    
    static std::vector<unsigned int> toVec(const unsigned int dir_mask);
  };
  
  bool applyStatics(const std::map<int, Trace::TimeT> &shot_st,
            const std::map<int, Trace::TimeT> &recv_st,
            Trace &ccfs);

  bool removeLowFreq(const double sb_size, rbm_smooth::Curve &st_curve);
  
  bool selectCCF(const unsigned int dir, 
         const unsigned int stride, 
         const Trace &ccfs); 

  Trace::TimeT getTimeShiftOfMaxValue(const Trace &ccf, 
                      Trace::AmplT *pmax_val = NULL);
  
  Trace::TimeT getTimeShiftOfMinValue(const Trace &ccf, 
                      Trace::AmplT *pmin_val = NULL);
  
  bool isZeroTrace(const Trace &tr);


  double getPartOfZeros(const Trace &tr, Trace::AmplT eps = 0.);
};

#endif /* psc_service.hpp */