ReadFDCalib.cc

Go to the documentation of this file.

#include <sstream>

#include "ReadFDCalib.h"
#include <utl/ErrorLogger.h>

#include <evt/Event.h>
#include <det/Detector.h>

#include <utl/TimeStamp.h>
#include <utl/UTCDateTime.h>
#include <utl/TabulatedFunction.h>
#include <utl/AugerUnits.h>
#include <utl/AugerException.h>

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

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

using namespace ReadFDCalibNS;

ReadFDCalib::ReadFDCalib(){}

ReadFDCalib::~ReadFDCalib(){}

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

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

  // Set detector time
  Detector::GetInstance().Update(UTCDateTime(2012,1,1,0,0).GetTimeStamp());

  const FDetector& fDet = Detector::GetInstance().GetFDetector();

  // iterate over all eyes
   for (FDetector::EyeIterator eyeIt = fDet.EyesBegin() ;
        eyeIt != fDet.EyesEnd() ; ++eyeIt) {


    cout << "eye : " << eyeIt->GetName() << endl;
    cout << "---------------------------" << endl;

    // iterate over all telescopes in the eye
    for (Eye::TelescopeIterator telIt = eyeIt->TelescopesBegin();
         telIt != eyeIt->TelescopesEnd(); ++telIt) {

      cout << "  telescope : " << telIt->GetId() << endl;

      // loop on all pixels

      for (unsigned int iPixId = telIt->GetFirstPixelId() ;
           iPixId <= telIt->GetLastPixelId() ; ++iPixId) {

        cout << iPixId << ": ";

        // Get calibration constants as function of wavelength
        TabulatedFunction calVsWavelength;
        try {
          calVsWavelength = telIt->GetPixel(iPixId).GetEndToEndCalibrationConstant();
        } catch (DataNotFoundInDBException& ed) {
          cout << "could not retrieve calib consts for eye : "
               << eyeIt->GetName() << ", telescope " << telIt->GetId()
               << " " << endl;
          break;
        }

        // Iterate through calibration constants for each wavelength
        for (TabulatedFunction::ConstIterator calIt = calVsWavelength.Begin() ;
             calIt != calVsWavelength.End() ; ++calIt) {
          cout << "(" << calIt->X()/nanometer << " nm, " << calIt->Y() << ") ";
        }
        cout << ", ";
      }
      cout << endl;
    }
  }
  return eSuccess;
}

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