testSuperMieModel.cc

Go to the documentation of this file.

#include <cppunit/extensions/HelperMacros.h>

#include <det/Detector.h>
#include <fdet/FDetector.h>
#include <fdet/Eye.h>

#include <fwk/CoordinateSystemRegistry.h>

#include <atm/AerosolDB.h>
#include <atm/AttenuationResult.h>
#include <atm/ScatteringResult.h>

#include <fwk/CentralConfig.h>

#include <utl/TimeStamp.h>
#include <utl/UTCDateTime.h>
#include <utl/ErrorLogger.h>

#include <utl/Point.h>
#include <utl/TabulatedFunctionErrors.h>
#include <utl/ReferenceEllipsoid.h>

#include <tst/Verify.h>

#include <iostream>

using namespace std;
using namespace det;
using namespace fdet;
using namespace atm;
using namespace fwk;
using namespace utl;
using namespace tst;


class testSuperMieModel : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(testSuperMieModel);
  CPPUNIT_TEST(testGetAttFromDB);
  CPPUNIT_TEST(testGetAttFromXML);
  CPPUNIT_TEST(testGetScatFromDB);
  CPPUNIT_TEST(testGetScatFromXML);

  CPPUNIT_TEST(testUncertaintyDB);
  CPPUNIT_TEST(testUncertaintyXML);

  CPPUNIT_TEST_SUITE_END();

public:
  void
  setUp()
  {
    CentralConfig::GetInstance(BOOTSTRAPFILE);
    ErrorLogger::GetInstance().SetVerbosity(Verbosity::eVerbose);
  }

  void tearDown() { }

  void testGetAttFromDB()
  {
    // set a time which should be valid for the test database.
    //
    Detector::GetInstance().Update(UTCDateTime(2005,1,1,0,0,0).GetTimeStamp());
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

    vector<double> wlength;
    wlength.push_back(300.*nanometer);
    wlength.push_back(355.*nanometer);
    wlength.push_back(400.*nanometer);

    // Los Leones
    //
    CoordinateSystemPtr
      CSll(Detector::GetInstance().GetFDetector().GetEye("Los Leones").GetLocalCoordinateSystem() );

    Point All(0.,      0.,   0.,        CSll);
    Point Bll(0.,      0.,   9000.*m,   CSll);

    AttenuationResult AttABll = theAtm.EvaluateMieAttenuation(All, Bll, wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(AttABll.GetTransmissionFactor().Y(355*nanometer),0.979985, 1.e-5));
  }

  void
  testGetAttFromXML()
  {
    // Set time for which there is no data in the test database.
    // The smart model should use the parametric model in this case
    //
    Detector::GetInstance().Update(TimeStamp(111));
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

    vector<double> wlength;
    wlength.push_back(300.*nanometer);
    wlength.push_back(355.*nanometer);
    wlength.push_back(400.*nanometer);

    // Los Leones
    //
    CoordinateSystemPtr
      CSll(Detector::GetInstance().GetFDetector().GetEye("Los Leones").GetLocalCoordinateSystem() );

    Point All(0.,      0.,   0.,        CSll);
    Point Bll(0.,      0.,   2600.*m,   CSll);

    AttenuationResult AttABll = theAtm.EvaluateMieAttenuation(All, Bll, wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(AttABll.GetTransmissionFactor().Y(300.*nanometer), 0.915745));
  }

  void
  testGetScatFromDB()
  {
    // set a time which should be valid for the test database.
    //
    Detector::GetInstance().Update(UTCDateTime(2005,1,1,0,0,0).GetTimeStamp());
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();
    CoordinateSystemPtr
      CSll(Detector::GetInstance().GetFDetector().GetEye("Los Leones").GetLocalCoordinateSystem() );

    vector<double> wlength;
    wlength.push_back(300*nanometer);
    wlength.push_back(355*nanometer);
    wlength.push_back(400*nanometer);

    Point pointA(0., 0.    ,  0.,      CSll);
    Point pointB(0., 0.    ,  1400.*m, CSll);

    ScatteringResult scatAB1 =
      theAtm.EvaluateMieScattering(pointA, pointB,
                                   20.*degree, 1.*km, wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(1.e8*scatAB1.GetScatteringFactor().Y(400.*nanometer), 0.856146, 1.e-5));
  }

  void
  testGetScatFromXML()
  {
    // set to time for which no data exist in the test database
    //
    Detector::GetInstance().Update(TimeStamp(111));
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

    vector<double> wlength;
    wlength.push_back(300*nanometer);
    wlength.push_back(350*nanometer);
    wlength.push_back(400*nanometer);

    CoordinateSystemPtr
      CSll(Detector::GetInstance().GetFDetector().GetEye("Los Leones").GetLocalCoordinateSystem() );
    Point pointA(0., 0.    ,  0.,      CSll);
    Point pointB(0., 0.    ,  1400.*m, CSll);

    ScatteringResult scatAB1 =
      theAtm.EvaluateMieScattering(pointA, pointB,
                                   20.*degree, 1.*km, wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(1.e8*scatAB1.GetScatteringFactor().Y(400.*nanometer), 2.70631, 1.e-5));
  }

  void testUncertaintyDB() {
    Detector::GetInstance().Update(UTCDateTime(2005,1,1,0,0,0).GetTimeStamp());

    CoordinateSystemPtr
      CSll(Detector::GetInstance().GetFDetector().GetEye("Los Leones").GetLocalCoordinateSystem() );

    ReferenceEllipsoid wgs84(ReferenceEllipsoid::Get(ReferenceEllipsoid::eWGS84));

    Point All(0.,      0.,   0.,        CSll);
    Point Bll(0.,      0.,   2600.*m,   CSll);

    vector<double> wlength;
    wlength.push_back(300.*nanometer);
    wlength.push_back(355.*nanometer);
    wlength.push_back(400.*nanometer);

    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();
    AttenuationResult attABll = theAtm.EvaluateMieAttenuation(All,
                                                              Bll,
                                                              wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(attABll.GetTransmissionFactor().Y(355*nanometer), 0.9814559317));
    theAtm.SetUncertaintyBound(Atmosphere::eMie, 1.);
    attABll = theAtm.EvaluateMieAttenuation(All,
                                            Bll,
                                            wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(attABll.GetTransmissionFactor().Y(355*nanometer), 0.9796204908));
    theAtm.SetUncertaintyBound(Atmosphere::eMie, -1.);
    attABll = theAtm.EvaluateMieAttenuation(All,
                                            Bll,
                                            wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(attABll.GetTransmissionFactor().Y(355*nanometer), 0.9832947194));
  }

  void testUncertaintyXML() {
    // Set time for which there is no data in the test database.
    // The smart model should use the parametric model in this case
    //
    Detector::GetInstance().Update(TimeStamp(111));
    const Atmosphere& theAtm = Detector::GetInstance().GetAtmosphere();

    vector<double> wlength;
    wlength.push_back(300.*nanometer);
    wlength.push_back(355.*nanometer);
    wlength.push_back(400.*nanometer);

    // Los Leones
    //
    CoordinateSystemPtr
      CSll(Detector::GetInstance().GetFDetector().GetEye("Los Leones").GetLocalCoordinateSystem() );

    Point All(0.,      0.,   0.,        CSll);
    Point Bll(0.,      0.,   2600.*m,   CSll);

    AttenuationResult attABll = theAtm.EvaluateMieAttenuation(All,
                                                              Bll,
                                                              wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(attABll.GetTransmissionFactor().Y(355*nanometer), 0.928318));

    theAtm.SetUncertaintyBound(Atmosphere::eMie, 1.);
    attABll = theAtm.EvaluateMieAttenuation(All,
                                            Bll,
                                            wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(attABll.GetTransmissionFactor().Y(355*nanometer), 0.928318));
    theAtm.SetUncertaintyBound(Atmosphere::eMie, -1.);
    attABll = theAtm.EvaluateMieAttenuation(All,
                                            Bll,
                                            wlength);
    CPPUNIT_ASSERT(Verify<CloseTo>(attABll.GetTransmissionFactor().Y(355*nanometer), 0.928318));

  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(testSuperMieModel);


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