ReadCloudHL.cc

Go to the documentation of this file.

#include <iomanip>
#include <sstream>
using namespace std;

#include <evt/Event.h>
using namespace evt;

#include <det/Detector.h>
using namespace det;

#include <fdet/FDetector.h>
#include <fdet/Eye.h>
#include <fdet/Telescope.h>
#include <fdet/Pixel.h>
using namespace fdet;

#include <utl/ErrorLogger.h>
#include <utl/Point.h>
#include <utl/TimeStamp.h>
#include <utl/Vector.h>
using namespace utl;

#include <atm/CloudResult.h>
using namespace atm;

#include "ReadCloudHL.h"
using namespace ReadCloudHLNS;
using namespace fwk;

#include <atm/LidarDB.h>
#include <atm/LidarZone.h>

ReadCloudHL::ReadCloudHL()
{
}

ReadCloudHL::~ReadCloudHL()
{
}

VModule::ResultFlag
ReadCloudHL::Init()
{
  INFO("ReadCloudHL::Init()");
  return eSuccess;
}

VModule::ResultFlag
ReadCloudHL::Run(evt::Event& event)
{
  INFO("ReadCloudHL::Run()");

  // Set detector time to something .
  Detector& detector = Detector::GetInstance();
  detector.Update(TimeStamp(832384850) );

  const Atmosphere& atmos = detector.GetAtmosphere();

  try {
    // Define points in the sky and determine whether a cloud obscures the
    // view of these points from a given FD pixel.
    //
    const fdet::Eye& detLL = detector.GetFDetector().GetEye("Los Leones");
    CoordinateSystemPtr csLL = detLL.GetEyeCoordinateSystem();

    // Get an fdet::Pixel and two points it looks at.
    const fdet::Pixel& pixel = detLL.GetTelescope(1).GetPixel(22);

    double scale;
    const Vector& pixelDir = pixel.GetDirection();
    const Point origin(0., 0., 0., csLL);

    scale = 1. * km;
    Point A = scale * pixelDir + origin;

    scale = 10. * km;
    Point B = scale * pixelDir + origin;

    CloudResult cResult;

    // Is the fdet::Pixel obscured at point A?  (Same syntax for fevt::Pixel.)
    cResult = atmos.EvaluateCloudCoverage(pixel, A);
    cout << "Cloud status in Pixel view to Point A:\n"
         << cResult.GetCoverageString() << endl;

    // Is the fdet::Pixel obscured at point B?  (Same syntax for fevt::Pixel.)
    cResult = atmos.EvaluateCloudCoverage(pixel, B);
    cout << "Cloud status in Pixel view to Point B:\n"
         << cResult.GetCoverageString() << endl;
  }
  catch (const AugerException& e) {
    ERROR(e.GetMessage());
  }
  return eSuccess;
}

VModule::ResultFlag
ReadCloudHL::Finish()
{
  INFO("ReadCloudHL::Finish()");
  return eSuccess;
}