StraightTrajectory.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/geometry/PhysicalGeometry.hpp>
#include <corsika/framework/geometry/BaseTrajectory.hpp>

namespace corsika {

  class StraightTrajectory : public BaseTrajectory {

  public:
    StraightTrajectory() = delete;
    StraightTrajectory(StraightTrajectory const&) = default;
    StraightTrajectory(StraightTrajectory&&) = default;
    StraightTrajectory& operator=(StraightTrajectory const&) = delete;

    StraightTrajectory(Line const& theLine, TimeType timeLength)
        : line_(theLine)
        , timeLength_(timeLength)
        , timeStep_(timeLength)
        , initialVelocity_(theLine.getVelocity(TimeType::zero()))
        , finalVelocity_(theLine.getVelocity(timeLength)) {}

    StraightTrajectory(Line const& theLine,
                       TimeType const timeLength, // length of theLine (straight)
                       TimeType const timeStep,   // length of bend step (curved)
                       VelocityVector const& initialV, VelocityVector const& finalV)
        : line_(theLine)
        , timeLength_(timeLength)
        , timeStep_(timeStep)
        , initialVelocity_(initialV)
        , finalVelocity_(finalV) {}

    Line const& getLine() const { return line_; }

    Point getPosition(double const u) const { return line_.getPosition(timeLength_ * u); }

    VelocityVector getVelocity(double const u) const;

    DirectionVector getDirection(double const u) const {
      return getVelocity(u).normalized();
    }

    TimeType getDuration(double const u = 1) const;

    template <typename Particle>
    TimeType getTime(Particle const& particle, double const u) const;

    LengthType getLength(double const u = 1) const;

    void setLength(LengthType const limit);

    //   Scale other properties by "limit/timeLength_"
    void setDuration(TimeType const limit);

  protected:
    LengthType getDistance(double const u) const;

    void setFinalVelocity(VelocityVector const& v) { finalVelocity_ = v; }

  private:
    Line line_;
    TimeType timeLength_; 
    TimeType timeStep_;   
    VelocityVector initialVelocity_;
    VelocityVector finalVelocity_;
  };

} // namespace corsika

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