REvent/StationRecData.cc

Go to the documentation of this file.

#include <revt/StationRecData.h>

#include <utl/ErrorLogger.h>
#include <utl/AnalyticCoordinateTransformator.h>
#include <utl/Math.h>

using namespace revt;
using namespace utl;


double
StationRecData::GetSignalArrivalAzimuth()
  const
{
  if (HasParameter(eSignalArrivalDirectionX) && HasParameter(eSignalArrivalDirectionY)) {
    return
      AnalyticCoordinateTransformator::GetAzimuthFromCartesianCoordinates(
        GetParameter(eSignalArrivalDirectionX), GetParameter(eSignalArrivalDirectionY)
      );
  } else {
    WARNING("StationRecData::GetSignalArrivalAzimuth(): Parameter eSignalArrivalDirection was not set");
    return std::numeric_limits<double>::quiet_NaN();
  }
}


double
StationRecData::GetSignalArrivalAzimuthError()
  const
{
  if (HasParameter(eSignalArrivalDirectionX) &&
      HasParameter(eSignalArrivalDirectionY) &&
      HasParameter(eSignalArrivalDirectionZ) &&
      HasParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionX) &&
      HasParameterCovariance(eSignalArrivalDirectionY, eSignalArrivalDirectionY) &&
      HasParameterCovariance(eSignalArrivalDirectionZ, eSignalArrivalDirectionZ) &&
      HasParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionY) &&
      HasParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionZ) &&
      HasParameterCovariance(eSignalArrivalDirectionY, eSignalArrivalDirectionZ)) {
    const double covariance22 =
      AnalyticCoordinateTransformator::GetSphericalCovarianceFromCartesian(
        2,
        2,
        GetParameter(eSignalArrivalDirectionX),
        GetParameter(eSignalArrivalDirectionY),
        GetParameter(eSignalArrivalDirectionZ),
        GetParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionX),
        GetParameterCovariance(eSignalArrivalDirectionY, eSignalArrivalDirectionY),
        GetParameterCovariance(eSignalArrivalDirectionZ, eSignalArrivalDirectionZ),
        GetParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionY),
        GetParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionZ),
        GetParameterCovariance(eSignalArrivalDirectionY, eSignalArrivalDirectionZ)
      );
    return sqrt(covariance22);
  } else {
    WARNING("StationRecData::GetSignalArrivalAzimuthError(): Parameter or Covariances"
            " of eSignalArrivalDirection was not set, returning NaN");
    return std::numeric_limits<double>::quiet_NaN();
  }
}


double
StationRecData::GetSignalArrivalZenith()
  const
{
  if (HasParameter(eSignalArrivalDirectionX) &&
      HasParameter(eSignalArrivalDirectionY) &&
      HasParameter(eSignalArrivalDirectionZ)) {
    return
      AnalyticCoordinateTransformator::GetZenithFromCartesianCoordinates(
        GetParameter(eSignalArrivalDirectionX), GetParameter(eSignalArrivalDirectionY), GetParameter(eSignalArrivalDirectionZ)
      );
  } else {
    WARNING("StationRecData::GetSignalArrivalZenith(): Parameter eSignalArrivalDirection was not set");
    return std::numeric_limits<double>::quiet_NaN();
  }
}


double
StationRecData::GetSignalArrivalZenithError()
  const
{
  if (HasParameter(eSignalArrivalDirectionX) &&
      HasParameter(eSignalArrivalDirectionY) &&
      HasParameter(eSignalArrivalDirectionZ) &&
      HasParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionX) &&
      HasParameterCovariance(eSignalArrivalDirectionY, eSignalArrivalDirectionY) &&
      HasParameterCovariance(eSignalArrivalDirectionZ, eSignalArrivalDirectionZ) &&
      HasParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionY) &&
      HasParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionZ) &&
      HasParameterCovariance(eSignalArrivalDirectionY, eSignalArrivalDirectionZ)) {
    const double covariance11 =
      AnalyticCoordinateTransformator::GetSphericalCovarianceFromCartesian(
        1,
        1,
        GetParameter(eSignalArrivalDirectionX),
        GetParameter(eSignalArrivalDirectionY),
        GetParameter(eSignalArrivalDirectionZ),
        GetParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionX),
        GetParameterCovariance(eSignalArrivalDirectionY, eSignalArrivalDirectionY),
        GetParameterCovariance(eSignalArrivalDirectionZ, eSignalArrivalDirectionZ),
        GetParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionY),
        GetParameterCovariance(eSignalArrivalDirectionX, eSignalArrivalDirectionZ),
        GetParameterCovariance(eSignalArrivalDirectionY, eSignalArrivalDirectionZ)
      );
    return sqrt(covariance11);
  } else {
    WARNING("StationRecData::GetSignalArrivalZenithError(): Parameter or Covariances"
            " of eSignalArrivalDirection was not set, returning NaN");
    return std::numeric_limits<double>::quiet_NaN();
  }
}