MTDataComp.h

Go to the documentation of this file.

#ifndef CMTDATACOMP_H
#define CMTDATACOMP_H

#include <string>
#include "types.h"

/** \addtogroup mttools MT data analysis, processing and inversion */
/* @{ */

/*! The class MTDataComp stores the values of the impedance tensor,
 * the error estimates along with the apparent resistivity and phase.
 * Units are \f$ km/s\f$ for Z, \f$ \Omega m \f$ for \f$ \rho_a \f$ and
 * degree for the phase. The corresponding frequencies are not stored
 * because normally it is used by CMTData or a derived class.
 */
 
class MTDataComp
{
 private:
        //! Whenever an impedance value changes, Update() recalculates apparent resistvity and phase. 
    void Update(trealdata frequency); 
    void Assign(const int nfreq);
    //!  name can be either Zxx, Zxy, Zyx, Zyy, Tx or Ty
    std::string name;
    //! The vector containing apparent resistivities
    trealdata rhoa; 
    //! The vector containing the phase
    trealdata phi;
    //! The corresponding errors in apparent resistivities
    trealdata drhoa;
    //! The corresponding errors in phase
    trealdata dphi;
    //! A complex vector containing the impedance values
    tcompdata Z;
    //! A real vector containing the error estimates
    trealdata dZ;
 public:
    std::string GetName()const{return name;}
    trealdata GetRhoa()const{return rhoa;} 
    trealdata GetPhi()const{return phi;}
    trealdata GetdRhoa()const{return drhoa;}
    trealdata GetdPhi()const{return dphi;}
    tcompdata GetZ()const{return Z;}
    trealdata GetdZ()const{return dZ;}
    //! The constructor is empty
    MTDataComp(const int nfreq = 0); 
    //! The destructor is empty,too
    virtual ~MTDataComp();
    //! CMTStation needs to access the Update method, so we declare it a friend
    friend class MTStation;
    friend class C1DMTSynthData;
    friend class JParser;
    friend class EDIParser;
    virtual MTDataComp& operator= (const MTDataComp& source);
    virtual void erase(const int index); 
};
/* @} */
#endif

---