AerosolDB.cc

Go to the documentation of this file.

#include <sstream>
#include <iostream>

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

#include <utl/AugerException.h>

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


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

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

  const VManager::Status stat =
    manager.GetData(fAerosolIdString, "aerosol", "aerosol_id", emptyIndexMap);

  if (stat == VManager::eNotFound) {
    /*ostringstream msg;
    msg << "Request for AerosolDB, but no aerosol data available for current time: ";
        << Detector::GetInstance().GetTime() << ". Returning 0 zones.";
    WARNING(msg.str());*/
    return; // return AerosolDB with 0 zones
  }

  CacheZones();
}


void
AerosolDB::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");

  //fZoneMap = new InternalZoneMap;

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

  // Make an AerosolZone object for each zone valid at the current
  // detector time (and for specified software version, etc.)
  //
  vector<string> zoneIdStrings;
  det::VManager::IndexMap indexMap;
  indexMap["aerosol_id"] = fAerosolIdString;
  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) {

    det::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 AerosolZone* const aZone =
       new AerosolZone(*it, fAerosolIdString, zoneName, zoneNorthing, zoneEasting);

    fZoneMap[zoneName] = aZone;
  }
}

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