TabulatedFlatAtmospherePropagator.hpp

/*
 * (c) Copyright 2023 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/PhysicalConstants.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp>
#include <corsika/modules/radio/propagators/RadioPropagator.hpp>

namespace corsika {

  template <typename TEnvironment>
  class TabulatedFlatAtmospherePropagator final
      : public RadioPropagator<TabulatedFlatAtmospherePropagator<TEnvironment>,
                               TEnvironment> {

    using Base =
        RadioPropagator<TabulatedFlatAtmospherePropagator<TEnvironment>, TEnvironment>;
    using SignalPathCollection = typename Base::SignalPathCollection;

  public:
    TabulatedFlatAtmospherePropagator(TEnvironment const& env, Point const& upperLimit,
                                      Point const& lowerLimit, LengthType const step);

    template <typename Particle>
    SignalPathCollection propagate(Particle const& particle, Point const& source,
                                   Point const& destination);

  private:
    Point const upperLimit_; 
    Point const lowerLimit_; 
    LengthType const step_;  
    InverseLengthType const
        inverseStep_;            
    LengthType const maxHeight_; 
    LengthType const minHeight_; 
    std::vector<double> refractivityTable_; 
    std::vector<LengthType>
        heightTable_; 
    std::vector<double>
        integratedRefractivityTable_; 
    double slopeRefrLower_; 
    double slopeIntRefrLower_; 
    double slopeRefrUpper_; 
    double slopeIntRefrUpper_; 
    double lastElement_;       
    std::deque<Point> points;  
    std::vector<double>
        rindex; 

  }; // End: FlatEarthPropagator

  template <typename TEnvironment>
  TabulatedFlatAtmospherePropagator<TEnvironment>
  make_tabulated_flat_atmosphere_radio_propagator(TEnvironment const& env,
                                                  Point const& upperLimit,
                                                  Point const& lowerLimit,
                                                  LengthType const step) {
    return TabulatedFlatAtmospherePropagator<TEnvironment>(env, upperLimit, lowerLimit,
                                                           step);
  }

} // namespace corsika

#include <corsika/detail/modules/radio/propagators/TabulatedFlatAtmospherePropagator.inl>