Helix.hpp

/*
 * (c) Copyright 2020 CORSIKA Project, corsika-project@lists.kit.edu
 *
 * This software is distributed under the terms of the GNU General Public
 * Licence version 3 (GPL Version 3). See file LICENSE for a full version of
 * the license.
 */

#pragma once

#include <corsika/framework/core/PhysicalUnits.hpp>
#include <corsika/framework/core/PhysicalGeometry.hpp>
#include <corsika/framework/geometry/Point.hpp>
#include <corsika/framework/geometry/Vector.hpp>

#include <cmath>

namespace corsika {

  class Helix {

  public:
    Helix(Point const& pR0, FrequencyType pOmegaC, VelocityVector const& pvPar,
          VelocityVector const& pvPerp)
        : r0_(pR0)
        , omegaC_(pOmegaC)
        , vPar_(pvPar)
        , vPerp_(pvPerp)
        , uPerp_(vPerp_.cross(vPar_.normalized()))
        , radius_(pvPar.getNorm() / abs(pOmegaC)) {}

    LengthType getRadius() const;

    Point getPosition(TimeType const t) const;

    VelocityVector getVelocity(TimeType const t) const;

    Point getPositionFromArclength(LengthType const l) const;

    LengthType getArcLength(TimeType const t1, TimeType const t2) const;

    TimeType getTimeFromArclength(LengthType const l) const;

  private:
    Point r0_;             
    FrequencyType omegaC_; 
    VelocityVector vPar_;  
    VelocityVector vPerp_, uPerp_;
    LengthType radius_;
  };

} // namespace corsika

#include <corsika/detail/framework/geometry/Helix.inl>