SdReconstruction/SdHorizontalReconstruction/Utilities.h

Go to the documentation of this file.

#ifndef _SdHorizontalReconstruction_Utilities_h_
#define _SdHorizontalReconstruction_Utilities_h_

#include <utl/CoordinateSystemPtr.h>
#include <utl/Point.h>
#include <utl/Vector.h>

namespace SdHorizontalReconstructionNS {

  inline
  utl::CoordinateSystemPtr
  GetShowerCoordinateSystem(const double theta, const double phi,
                            const utl::CoordinateSystemPtr& coreCS)
  {
    return utl::CoordinateSystem::RotationY(
      theta,
      utl::CoordinateSystem::RotationZ(
        phi,
        coreCS));
  }

  inline
  void
  GetRhoAndDelta(double& rho, double& delta,
                 const utl::Point& pos, const utl::Point& core, const utl::Point& origin)
  {
    using namespace utl;
    const Vector coreToPos = pos - core;
    const Vector coreToOrigin = origin - core;
    delta = coreToPos*coreToOrigin/coreToOrigin.GetMag();
    rho = std::sqrt(coreToPos.GetMag2() - delta*delta);
  }

  inline
  double
  GetRho(const utl::Point& pos, const utl::Point& core, const utl::Point& origin)
  {
    double rho, delta;
    GetRhoAndDelta(rho, delta, pos, core, origin);
    return rho;
  }

} // NS SdHorizontalReconstructionNS

#endif