VRayleighModel.cc

Go to the documentation of this file.

#include <utl/Point.h>
#include <utl/VRandomSampler.h>
#include <utl/RandomSamplerFromPDF.h>
#include <utl/PhysicalConstants.h>
#include <utl/TabulatedFunction.h>

#include <atm/VRayleighModel.h>

using namespace utl;


namespace atm {

  VRayleighModel::~VRayleighModel()
  {
    delete fAngularDistribution;
  }


  const utl::VRandomSampler&
  VRayleighModel::EvaluateScatteringAngleDistribution(const utl::Point& p, const double wLength)
    const 
  {
    const double angleMin = 0*deg;
    const double angleMax = 360*deg;
    const int n = 180; // 2deg steps
    const double delta = (angleMax - angleMin) / (n-1);  
    
    TabulatedFunction pdf;
    for (int i = 0; i < n; ++i) {
      const double angle = angleMin + delta * i;
      const double prob = EvaluateScatteringAngle(p, angle, wLength);
      pdf.PushBack(angle, prob);
    }
    delete fAngularDistribution;
    fAngularDistribution = new RandomSamplerFromPDF(pdf, RandomSamplerFromPDF::eLinear);
    return *fAngularDistribution;
  }

}