VMieModel.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/VMieModel.h>

using namespace atm;
using namespace utl;


VMieModel::VMieModel() :
  fAngularDistribution(0)
{
}


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


const utl::VRandomSampler&
VMieModel::EvaluateScatteringAngleDistribution(const 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;
}


// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// compile-command: "make -C .. -k"
// End: