CaUvConverter.h

Go to the documentation of this file.

/* Copyright (C) 2023 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sergey Gorbunov [committer] */
 
#ifndef CaUvConverter_h
#define CaUvConverter_h 1
 
 
#include "CaDefs.h"
 
namespace cbm::ca
{
 
  class CaUvConverter {
 
   public:
    CaUvConverter(double phiU, double phiV) { SetFromUV(phiU, phiV); }
 
    CaUvConverter(double phiU, double dx2, double dxy, double dy2) { SetFromXYCovMatrix(phiU, dx2, dxy, dy2); }
 
    void SetFromUV(double phiU, double phiV);
 
    void SetFromXYCovMatrix(double phiU, double dx2, double dxy, double dy2);
 
    std::pair<double, double> ConvertXYtoUV(double x, double y) const
    {
      return std::make_pair(fcosU * x + fsinU * y, fcosV * x + fsinV * y);
    }
 
    std::pair<double, double> ConvertUVtoXY(double u, double v) const
    {
      return std::make_pair(fcosX * u + fsinX * v, fcosY * u + fsinY * v);
    }
 
    std::tuple<double, double, double> ConvertCovMatrixUVtoXY(double du2, double duv, double dv2) const
    {
      return std::make_tuple(fcosX * fcosX * du2 + 2. * fsinX * fcosX * duv + fsinX * fsinX * dv2,
                             fcosX * fcosY * du2 + (fcosX * fsinY + fsinX * fcosY) * duv + fsinX * fsinY * dv2,
                             fcosY * fcosY * du2 + 2. * fsinY * fcosY * duv + fsinY * fsinY * dv2);
    }
 
    std::tuple<double, double, double> ConvertCovMatrixXYtoUV(double dx2, double dxy, double dy2) const
    {
      return std::make_tuple(fcosU * fcosU * dx2 + 2. * fsinU * fcosU * dxy + fsinU * fsinU * dy2,
                             fcosU * fcosV * dx2 + (fcosU * fsinV + fsinU * fcosV) * dxy + fsinU * fsinV * dy2,
                             fcosV * fcosV * dx2 + 2. * fsinV * fcosV * dxy + fsinV * fsinV * dy2);
    }
 
 
   private:
    double fcosU{cbm::algo::ca::constants::Undef<double>};  
    double fsinU{cbm::algo::ca::constants::Undef<double>};
 
    double fcosV{cbm::algo::ca::constants::Undef<double>};  
    double fsinV{cbm::algo::ca::constants::Undef<double>};
 
    double fcosX{cbm::algo::ca::constants::Undef<double>};  
    double fsinX{cbm::algo::ca::constants::Undef<double>};
 
    double fcosY{cbm::algo::ca::constants::Undef<double>};  
    double fsinY{cbm::algo::ca::constants::Undef<double>};
  };
 
}  // namespace cbm::ca
 
 
#endif