ReadCloudLL.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/TimeStamp.h>
using namespace utl;

#include "ReadCloudLL.h"
using namespace ReadCloudLLNS;
using namespace fwk;


ReadCloudLL::ReadCloudLL()
{
}


ReadCloudLL::~ReadCloudLL()
{
}


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


VModule::ResultFlag
ReadCloudLL::Run(evt::Event& /*event*/)
{
  INFO("ReadCloudLL::Run()");

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

  try {
    // Get an fdet::Pixel and output its coverage (same syntax for fevt::Pixel).
    const Eye& eye = det.GetFDetector().GetEye(1);
    const Pixel& pixel = eye.GetTelescope(1).GetPixel(2);

    if ( ! pixel.HasCloudFraction() ) {
      ostringstream err;
      err << " no cloud fraction for pixel "
          << pixel.GetId()
          << " eye: " << pixel.GetEyeId()
          << " tel:" << pixel.GetTelescopeId();
      ERROR(err.str());
      return eFailure;
    }

    float coverage = pixel.GetCloudFraction();
    cout << "\nCLOUD FRACTION = "  << coverage << "%"
         << " for Pixel " << pixel.GetId()
         << " in Eye "     << pixel.GetEyeId()
         << ", Telescope " << pixel.GetTelescopeId()
         << endl;

    // Output the coverage data for a full site.
    int eyeId = eye.GetId();
    cout << "\nEYE: " << eyeId << endl;

    for (unsigned int telId = eye.GetFirstTelescopeId(); telId <=
         eye.GetLastTelescopeId(); ++telId) {

      const Telescope& telescope = eye.GetTelescope(telId);
      cout << "\nTELESCOPE: " << telId << endl;

      for (unsigned int pixId = telescope.GetFirstPixelId();
           pixId <= telescope.GetLastPixelId(); ++pixId) {

        const Pixel& detPixel = telescope.GetPixel(pixId);
        if ( detPixel.HasCloudFraction() ) {

          coverage = detPixel.GetCloudFraction();

          cout << setw(3) << coverage;
          if (pixId % 20 == 0)
            cout << endl;

        }
      }
    }
  }
  catch (const AugerException& e) {
    ERROR(e.GetMessage());
  }

  return eSuccess;
}


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