MieModelOff.cc

Go to the documentation of this file.

#include <iostream>
#include <string>
#include <sstream>
#include <vector>

#include <det/Detector.h>
#include <atm/MieModelOff.h>
#include <atm/ScatteringResult.h>
#include <atm/AttenuationResult.h>
#include <utl/TabulatedFunctionErrors.h>

#include <fwk/CentralConfig.h>
#include <fwk/CoordinateSystemRegistry.h>

#include <utl/Point.h>
#include <utl/AugerUnits.h>
#include <utl/Reader.h>
#include <utl/ErrorLogger.h>
#include <utl/Vector.h>

#include <utl/ReferenceEllipsoid.h>

using namespace utl;
using namespace atm;
using namespace fwk;
using namespace std;


ScatteringResult 
MieModelOff::EvaluateMieScattering(const utl::Point& /*xA*/, const utl::Point& /*xB*/,
                                   const double angle, const double /*distance*/,
                                   const std::vector<double>& wLength)
  const
{
  utl::TabulatedFunctionErrors frac;
  for (std::vector<double>::const_iterator it = wLength.begin();
       it != wLength.end(); ++it)
    frac.PushBack(*it, 0,0,0);

  return ScatteringResult(frac, angle);
}


ScatteringResult 
MieModelOff::EvaluateMieScattering(const utl::Point& xA, const utl::Point& xB,
                                   const double angle, const double distance,
                                   const AttenuationResult& mieAttenuation)
  const
{
  // This method just delegates because the scattering disregards the
  // attenuation for "Off".
  const TabulatedFunctionErrors& attenuation = mieAttenuation.GetTransmissionFactor();
  vector<double> wLength(attenuation.GetNPoints());
  std::copy(attenuation.XBegin(), attenuation.XEnd(), wLength.begin());
  return EvaluateMieScattering(xA, xB, angle, distance, wLength);
}


AttenuationResult
MieModelOff::EvaluateMieAttenuation(const utl::Point& /*xInit*/, const utl::Point& /*xFinal*/,
                                    const std::vector<double>&  wLength)
  const
{
  utl::TabulatedFunctionErrors attenuation;
  for (std::vector<double>::const_iterator it = wLength.begin();
       it != wLength.end(); ++it)
    attenuation.PushBack(*it, 0,1,0);

  return AttenuationResult(attenuation);
}


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