MolecularZone.cc

Go to the documentation of this file.

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

#include <atm/MolecularZone.h>
#include <atm/MolecularLayer.h>
#include <atm/Atmosphere.h>
#include <utl/ErrorLogger.h>

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


MolecularZone::MolecularZone(const string& zoneId, const string& molecularId,
                             const string& zoneName,
                             const double northing, const double easting) :
  VZone(zoneId, zoneName, northing, easting),
  fLayerVector(0),
  fMolecularIdString(molecularId)
{
  // Find the id for this MolecularZone object.
  // This is just the primary key in the molecular_zone table for the
  // given molecular_id and zone_id
  //
  VManager::IndexMap indexMap;
  indexMap["molecular_id"] = fMolecularIdString;
  indexMap["zone_id"] = fZoneIdString;

  const VManager& manager = Detector::GetInstance().GetAManagerRegister();
  const VManager::Status stat =
    manager.GetData(fMolecularZoneIdString, "molecular_zone", "molecular_zone_id", indexMap);
  if (stat  == VManager::eNotFound)
    throw DataNotFoundInDBException("molecular_zone", "molecular_zone_id", indexMap,
                                    boost::lexical_cast<string>(Detector::GetInstance().GetTime()));
}


MolecularZone::LayerIterator
MolecularZone::LayersBegin()
  const
{
  if (!fLayerVector)
    CacheLayers(fLayerVector, fMolecularZoneIdString);
  return LayerIterator(fLayerVector->begin());
}


MolecularZone::LayerIterator
MolecularZone::LayersEnd()
  const
{
  if (!fLayerVector)
    CacheLayers(fLayerVector, fMolecularZoneIdString);
  return LayerIterator(fLayerVector->end());
}


MolecularZone::~MolecularZone()
{
  if (fLayerVector) {
    for (InternalLayerIterator sIt = fLayerVector->begin();
         sIt != fLayerVector->end(); ++sIt)
      delete *sIt;
    delete fLayerVector;
  }
}


void
MolecularZone::CacheLayers(vector<const MolecularLayer*>*& layerVector,
                           const string& zoneIdString)
  const
{
  if (layerVector)
    return;

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

  layerVector = new vector<const MolecularLayer*>;

  VManager::IndexMap indexMap;
  indexMap["molecular_zone_id"] = zoneIdString;

  vector<string> layerIdStrings;
  const VManager::Status stat =
    manager.GetData(layerIdStrings, "molecular_layer", "molecular_layer_id", indexMap);

  if (stat == VManager::eNotFound) {
    ostringstream msg;
    msg << "Did not find requested data, header table: molecular_layer, "
           " zoneId: " << zoneIdString;
    ERROR(msg);
    throw DataNotFoundInDBException("molecularLayer", "zone_id", indexMap,
                                    boost::lexical_cast<string>(Detector::GetInstance().GetTime()));
  }

  for (vector<string>::iterator it = layerIdStrings.begin();
       it != layerIdStrings.end(); ++it) {
    MolecularLayer* layer = new MolecularLayer(*it);
    layerVector->push_back(layer);
  }
}


// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k"
// End: