ShowerAxis.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/media/Environment.hpp>
#include <corsika/framework/geometry/Point.hpp>
#include <corsika/framework/geometry/Vector.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp>

#include <cassert>
#include <cstdlib>
#include <fstream>
#include <functional>
#include <iterator>
#include <memory>
#include <stdexcept>
#include <vector>

#include <iostream>

#include <boost/math/quadrature/gauss_kronrod.hpp>

namespace corsika {

  class ShowerAxis {
  public:
    template <typename TEnvModel>
    ShowerAxis(Point const& pStart, Vector<length_d> const& length,
               Environment<TEnvModel> const& env, bool const doThrow = false,
               int const steps = 10'000);

    LengthType getSteplength() const;

    GrammageType getMaximumX() const;

    GrammageType getMinimumX() const;

    GrammageType getProjectedX(Point const& p) const;

    GrammageType getX(LengthType) const;

    DirectionVector const& getDirection() const;

    Point const& getStart() const;

  private:
    Point const pointStart_;
    Vector<length_d> const length_;
    bool throw_ = false;
    LengthType const max_length_, steplength_;
    DirectionVector const axis_normalized_;
    std::vector<GrammageType> X_;

    // for storing the lengths corresponding to equidistant X values
    GrammageType const X_binning_ = 1_g / 1_cm / 1_cm;
    std::vector<LengthType> d_;
  };
} // namespace corsika

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