MultipleScatterer.h

Go to the documentation of this file.

#ifndef _MultipleScatterer_h_
#define _MultipleScatterer_h_


#include <utl/Photon.h>
#include <utl/Point.h>
#include <utl/Vector.h>
#include <utl/RandomEngine.h>

#include <vector>

using namespace std;



class MultipleScatterer {


  public:

    MultipleScatterer();

enum EAerosolModel {
      eLongtin,
      eHG01,
      eHG05,
      eHG09
    };

    void AddPhotons(const utl::Point& PointOfOrigin,
                        const utl::Photon& PhotonDirect,
                        const utl::Point& PointOfDetection,
                        std::vector<utl::Photon>& photonVector);

    void AddPhotons(const utl::Point& PointOfOrigin,
                        const utl::Photon& PhotonDirect,
                        const double Transmission,
                        std::vector<utl::Photon>& photonVector);


  private:

    double MaxZeta;
    double MaxTime;
    //double zeta; // unused. LN.
    EAerosolModel fAerosolModel;
    utl::RandomEngine* fRandomEngine;
    bool GetRandomZeta(double& zeta, double& weightfactor, double Time, double OptDistance, double height, double MaxZeta);
    double IntegralA(double A, double B, double C, double zeta);
    double IntegralB(double D, double E, double F, double G, double H, double zeta);
};


#endif