MolecularDB.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 <utl/AugerException.h>
#include <atm/Atmosphere.h>
#include <atm/MolecularZone.h>
#include <atm/MolecularDB.h>

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


const string MolecularDB::fgProfileNameLookup[MolecularIds::eUSStdA+1] = {
  "weather station",
  "radiosonde data",
  "Malargue January model",
  "Malargue February model",
  "Malargue March model",
  "Malargue April model",
  "Malargue May model",
  "Malargue June model",
  "Malargue July model",
  "Malargue August model",
  "Malargue September model",
  "Malargue October model",
  "Malargue November model",
  "Malargue December model",
  "Malargue GDAS data",
  "US-StdA"
};


MolecularDB::MolecularDB(const MolecularIds::ProfileId id) :
  fZoneMap(),
  fProfileId(id)
{
  //Detector::GetInstance().GetTime();
  const VManager& manager =
    Detector::GetInstance().GetAManagerRegister().GetManager("AMolecularSQLManager");

  // find the profile_type_id
  VManager::IndexMap profileIndexMap;
  profileIndexMap["profile_type_name"] = string("'") + fgProfileNameLookup[id] + string("'");
  string profileTypeIdString;
  const VManager::Status stat1 =
    manager.GetData(profileTypeIdString, "profile_type", "profile_type_id", profileIndexMap);
  if (stat1 == VManager::eNotFound) {
    ostringstream msg;
    msg << "Request for table 'profile_type' with profile_type_name '"
        << profileIndexMap["profile_type_name"] << "' not found.";
    WARNING(msg);
    return; // return MolecularDB with 0 zones
  }

  VManager::IndexMap molecularIndexMap;
  molecularIndexMap["profile_type_id"] = profileTypeIdString;
  //  molecularIndexMap["profile_type_name"] = fgProfileNameLookup[id];


  const VManager::Status stat2 =
    manager.GetData(fMolecularIdStrings, "molecular", "molecular_id", molecularIndexMap);
  if (stat2 == VManager::eNotFound) {
    ostringstream msg;
    msg << "Request for table 'molecular' with profile_type_id '"
        << molecularIndexMap["profile_type_id"] << "' turned up with no data "
           "for time: " << Detector::GetInstance().GetTime() << ". "
           "Returning 0 zones.";
    WARNING(msg);
    return; // return MolecularDB with 0 zones
  }

  // for (vector<string>::iterator sIt = fMolecularIdStrings.begin() ;
  //      sIt != fMolecularIdStrings.end(); ++sIt) {
  //   cout << "    " << *sIt << endl;
  // }

  CacheZones();
}


void
MolecularDB::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("AMolecularSQLManager");

  // Make an MolecularZone object for each zone valid at the current
  // detector time (and for specified software version, etc.)
  //

  for (unsigned int i = 0; i < fMolecularIdStrings.size(); ++i) {
    string zoneIdString;
    VManager::IndexMap indexMap;
    indexMap["molecular_id"] = fMolecularIdStrings[i];

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

    VManager::IndexMap zoneNameIndexMap;
    zoneNameIndexMap["zone_id"] = zoneIdString;

    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 MolecularZone* const aZone =
      new MolecularZone(zoneIdString, fMolecularIdStrings[i],
                        zoneName, zoneNorthing, zoneEasting);

    fZoneMap[zoneName] = aZone;
  }
}


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