InclinedAtmosphericProfile.h

Go to the documentation of this file.

#ifndef _atm_InclinedAtmosphericProfile_h_
#define _atm_InclinedAtmosphericProfile_h_

#include <vector>
#include <utl/AugerException.h>


namespace utl {
  class Point;
  class Vector;
}

namespace atm {

  class ProfileResult;

  class InclinedAtmosphericProfile {

  public:
    InclinedAtmosphericProfile(const utl::Point& core,
                               const utl::Vector& axis,
                               const double deltaDepth);

    ~InclinedAtmosphericProfile();

    const atm::ProfileResult& EvaluateSlantDepthVsDistance() const;

    const atm::ProfileResult& EvaluateDistanceVsSlantDepth() const;

    const atm::ProfileResult& EvaluateHeightVsSlantDepth() const;

    const atm::ProfileResult& EvaluateHeightVsDistance() const;

    class InclinedAtmosphereModelException;

    static double IntegratedGrammage(const utl::Point& pStart,
                                     const utl::Point& pStop,
                                     const double delta);

  private:
    atm::ProfileResult* fSlantDepthVsDistance = nullptr;
    atm::ProfileResult* fDistanceVsSlantDepth = nullptr;
    atm::ProfileResult* fHeightVsSlantDepth = nullptr;
    atm::ProfileResult* fHeightVsDistance = nullptr;

    class UTMZoneException;

  };


  class InclinedAtmosphericProfile::InclinedAtmosphereModelException : public utl::AugerException {
  public:
    InclinedAtmosphereModelException(const std::string& message = ""):
      AugerException(message) { }

    virtual std::string GetExceptionName() const override
    { return "Problem in inclined atmosphere model"; }
  };

}


#endif