ParametricXMLRayleighModel.h

Go to the documentation of this file.

#ifndef _atm_ParametricXMLRayleighModel_h_
#define _atm_ParametricXMLRayleighModel_h_


#include <atm/VRayleighModel.h>
#include <vector>

namespace atm {

  class ScatteringResult;
  class AttenuationResult;
  
  class ParametricXMLRayleighModel : public atm::VRayleighModel {

    public:

      ParametricXMLRayleighModel();

      virtual ~ParametricXMLRayleighModel() { }

      atm::ScatteringResult EvaluateRayleighScattering(
                                      const utl::Point& xA,
                                      const utl::Point&  xB,
                                      const double angle,
                                      const double distance,
                                      const std::vector<double>& wLength) const;
      
      atm::ScatteringResult EvaluateRayleighScattering(
                                      const utl::Point& xA,
                                      const utl::Point&  xB,
                                      const double angle,
                                      const double distance,
                                      const atm::AttenuationResult& raylAtt) const;
      
      atm::AttenuationResult EvaluateRayleighAttenuation(
                                      const utl::Point& xInit, 
                                      const utl::Point& xFinal,  
                                      const std::vector<double>& wLength) const;

      double EvaluateRayleighScattering(const utl::Point&  xA,
                                        const utl::Point&  xB,
                                        const double angle,
                                        const double distance,
                                        double wLength) const;
      
      double EvaluateRayleighScattering(const utl::Point&  xA,
                                        const utl::Point&  xB,
                                        const double angle,
                                        const double distance,
                                        double wLength,
                                        const double raylAtt) const;
      
      double EvaluateRayleighAttenuation(const utl::Point& xInit, 
                                         const utl::Point& xFinal,  
                                         double wLength) const;

      double GetAttenuationLength(const utl::Point& p, 
                                  const double wLength) const;

      double EvaluateScatteringAngle(const utl::Point& p,
                                     const double angle,
                                     const double wLength) const;
      
      bool HasData() const {return true;};

      void Init();

    private:

      bool fIntegratedGDistance;
      double fIntegrationStepWidth;

      // Depolarisation factor, Buchholtz Applied Optics 34, 2765 (1995)
      double fRhoN; 
      double GDistance(const utl::Point& xA, const utl::Point& xB) const;

      REGISTER_MODEL(atm::VRayleighModel, "ParametricXML", ParametricXMLRayleighModel);
  };

} // atm


#endif // _atm_ParametricXMLRayleighModel_h_