SStationListXMLManager.cc

Go to the documentation of this file.

#include <iostream>
#include <algorithm>

#include <sdet/SStationListXMLManager.h>
#include <utl/ErrorLogger.h>
#include <utl/Reader.h>
#include <sdet/SManagerRegister.h>

using namespace std;
using namespace sdet;
using namespace utl;
using namespace det;


REGISTER_S_MANAGER("SStationListXMLManager", SStationListXMLManager);


void
SStationListXMLManager::Init(const std::string& configLink)
{
  // first do the normal Init()
  VManager::Init(configLink);

  for (auto cB = fBranch.GetFirstChild(); cB; cB = cB.GetNextSibling()) {

    if (cB.GetName() != "station")
      continue;

    SStationListManager::StationData station;

    station.fId = FindComponent<int>("id", cB.GetAttributes());

    cB.GetChild("northing").GetData(station.fNorthing);
    cB.GetChild("easting").GetData(station.fEasting);
    cB.GetChild("altitude").GetData(station.fAltitude);
    station.fName = cB.GetChild("name").GetDataString();
    station.fCommissionTime = cB.GetChild("commission").GetDataString();
    station.fDecommissionTime = cB.GetChild("decommission").GetDataString();

    /* If the isUUB property does not exist (as is the case for real
       measurements and may be the case for pre-upgrade simulation
       station list XML), set isUUB to 0. Same holds for small PMT and
       scintillator. */

    auto isUUBB = cB.GetChild("isUUB");
    if (isUUBB)
      isUUBB.GetData(station.fIsUUB);
    else
      station.fIsUUB = false;

    auto hasSmallPMTB = cB.GetChild("hasSmallPMT");
    if (hasSmallPMTB)
      hasSmallPMTB.GetData(station.fHasSmallPMT);
    else
      station.fHasSmallPMT = false;

    auto hasScintillatorB = cB.GetChild("hasScintillator");
    if (hasScintillatorB)
      hasScintillatorB.GetData(station.fHasScintillator);
    else
      station.fHasScintillator = false;

    /* If the commissionUUB property does not exist (as is the case
       for simulations and uncommissioned stations), return the
       decommission time in all cases except when isUUB property
       exists (which is the case for simulations). In this case,
       return the station commission time. */

    auto commissionUUBB = cB.GetChild("commissionUUB");
    if (commissionUUBB)
      commissionUUBB.GetData(station.fUUBCommissionTime);
    else
      station.fUUBCommissionTime =
        station.fIsUUB ? station.fCommissionTime : station.fDecommissionTime;

    // interpret inGrid as indices where vector<bool> is true
    const auto inGrid = cB.GetChild("inGrid").Get<vector<int>>();
    if (inGrid.empty())
      station.fInGrid.clear();
    else {
      const int maxGrid = *max_element(inGrid.begin(), inGrid.end());
      station.fInGrid = vector<bool>(maxGrid);
      for (const auto ig : inGrid) {
        // skip 0, map 1 -> 1000, 2 -> 0100 etc
        if (ig <= 0)
          continue;
        station.fInGrid[ig - 1] = true;
      }
    }

    cB.GetChild("groupId").GetData(station.fGroupId);

    const auto axes = cB.GetChild("axes").Get<vector<double>>();
    station.fAxis1 = axes[0];
    station.fAxis2 = axes[1];

    station.fEllipsoid = cB.GetChild("ellipsoid").GetDataString();
    cB.GetChild("zone").GetData(station.fZone);
    station.fBand = cB.GetChild("band").GetDataString();

    const bool ok = fStationManager.ConditionalAddStationData(station);

    if (!ok) {
      ostringstream err;
      err << "multiple instances of station id=" << station.fId
          << " found in XML configuration";
      ERROR(err);
      throw utl::XMLParseException(err.str());
    }

  }
}