LidarDB.cc

Go to the documentation of this file.

#include <sstream>
#include <iostream>
#include <boost/lexical_cast.hpp>

#include <det/Detector.h>
#include <utl/ErrorLogger.h>
#include <atm/Atmosphere.h>
#include <atm/LidarDB.h>

#include <utl/AugerException.h>

using namespace det;
using namespace std;
using namespace atm;
using namespace utl;


LidarDB::LidarDB()
{
  // Find the id for this AerosolDB
  // (this is just the primary key in the aerosol table for the given
  //  time, database version, etc.
  VManager::IndexMap emptyIndexMap;

  const VManager& manager =
    Detector::GetInstance().GetAManagerRegister().GetManager("ALidarSQLManager");

  const VManager::Status stat =
    manager.GetData(fLidarIdString, "lidar", "lidar_id", emptyIndexMap);

  if (stat == VManager::eNotFound)
    return; // return LidarDB with 0 zones

  CacheZones();
}


void
LidarDB::CacheZones()
{
  // CacheZones should only be invoked when fZoneMap is 0
  if (!fZoneMap.empty())
    ERROR("Call to CacheZones with non-empty fZoneMap. This is a bug .. please report it");

  const VManager& manager =
    Detector::GetInstance().GetAManagerRegister().GetManager("ALidarSQLManager");

  // Make an LidarZone object for each zone valid at the current
  // detector time (and for specified software version, etc.)
  vector<string> zoneIdStrings;
  VManager::IndexMap indexMap;
  indexMap["lidar_id"] = fLidarIdString;
  const VManager::Status stat =
    manager.GetData(zoneIdStrings, "zone", "zone_id", indexMap);

  if (stat == VManager::eNotFound)
    throw DataNotFoundInDBException("zone", "zone_id", indexMap,
                                    boost::lexical_cast<string>(Detector::GetInstance().GetTime()));

  for (vector<string>::iterator it = zoneIdStrings.begin();
       it != zoneIdStrings.end(); ++it) {

    VManager::IndexMap zoneNameIndexMap;
    zoneNameIndexMap["zone_id"] = *it;

    string zoneName;
    manager.GetData(zoneName, "zone", "name", zoneNameIndexMap);

    double zoneNorthing;
    manager.GetData(zoneNorthing, "zone", "northing", zoneNameIndexMap);

    double zoneEasting;
    manager.GetData(zoneEasting, "zone", "easting", zoneNameIndexMap);

    const LidarZone* const aZone =
       new LidarZone(*it, fLidarIdString, zoneName, zoneNorthing, zoneEasting);

    fZoneMap[zoneName] = aZone;
  }

}

LidarDB::~LidarDB()
{
  for (InternalZoneIterator it = fZoneMap.begin();
       it != fZoneMap.end(); ++it)
    delete it->second;
}