BaseExponential.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/ParticleProperties.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp>
#include <corsika/framework/geometry/Line.hpp>
#include <corsika/framework/geometry/Point.hpp>
#include <corsika/framework/geometry/BaseTrajectory.hpp>
#include <limits>

namespace corsika {

  template <typename TDerived>
  class BaseExponential {

  public:
    BaseExponential(Point const& point, LengthType const referenceHeight,
                    MassDensityType rho0, LengthType lambda);

    Point const& getAnchorPoint() const { return point_; }

  protected:
    auto const& getImplementation() const;

    MassDensityType getMassDensity(LengthType const height) const;

    // clang-format off
    // clang-format on
    GrammageType getIntegratedGrammage(BaseTrajectory const& line,
                                       DirectionVector const& axis) const;

    // clang-format off
    // clang-format on
    LengthType getArclengthFromGrammage(BaseTrajectory const& line, GrammageType grammage,
                                        DirectionVector const& axis) const;

  private:
    MassDensityType const rho0_;
    LengthType const lambda_;
    InverseLengthType const invLambda_;
    Point const point_;
    LengthType const referenceHeight_;

  }; // class BaseExponential

} // namespace corsika

#include <corsika/detail/media/BaseExponential.inl>