RFVelCalc.h

Go to the documentation of this file.

#ifndef RFVELCALC_H_
#define RFVELCALC_H_
#include "RecCalc.h"
#include "SeismicDataComp.h"
#include <string>

/** \addtogroup seistools Seismic data analysis and modeling */
/* @{ */

//! This class implements the method to calculate absolute S-Wave velocities from Receiver function data
// as described by Sevnningsen and Jacobsen, GJI 2007
class RFVelCalc
{
private:
        RecCalc RFCalculator;
        ttsdata Velocities;
        ttsdata Periods;
        //! The core estimation routine, take radial and vertical receiver functions
        void AbsVelCalc(const double slowness,const SeismicDataComp &RadRec, const SeismicDataComp &VerRec,ttsdata &AppVel);
public:
        const ttsdata &GetVelocities(){return Velocities;}
        //! Calculate absolute velocities from the radial and vertical components of the seismogram
        void CalcRFVel(const double slowness, const SeismicDataComp &RadComp, const SeismicDataComp &VerComp,ttsdata &AppVel);
        //!Calculate absolute velocities from the radial receiver function and the vertical component of the seismogram
        void CalcRFVelFromRec(const double slowness, const SeismicDataComp &RRec, const SeismicDataComp &VerComp,ttsdata &AppVel);
        //! Write the velocity estiamtes and corresponding periods into an ascii file
        void WriteVelocities(const std::string filename);
        RFVelCalc(const double myomega, const double mysigma,const double myc, const RecCalc::trfmethod themethod = RecCalc::specdiv);
        RFVelCalc& operator=(const RFVelCalc& source);
        RFVelCalc(const RFVelCalc &Old);
        virtual ~RFVelCalc();
};
/* @} */
#endif /*RFVELCALC_H_*/

---