FDsimG4XMLManager.cc

Go to the documentation of this file.

#include <boost/lambda/lambda.hpp>

#include "FDsimG4XMLManager.hh"
#include <fwk/CentralConfig.h>
#include <utl/ErrorLogger.h>
#include <utl/Reader.h>
#include <utl/AugerException.h>

using namespace TelescopeSimulatorLX;
using namespace std;
using namespace utl;
using namespace fwk;
using namespace boost::lambda;


void
FDsimG4XMLManager::Init(const string& configLink)
{
  FindConfig(configLink);
  //  FillMaps();
  fIsInitialized = true;
}


void
FDsimG4XMLManager::FindConfig(const string& configLink)
{
  if (fIsInitialized) {
    ostringstream warn;
    warn << "Attempt to initialize an already initialized manager: " << configLink;
    WARNING(warn);
    return;
  }

  fBranch = CentralConfig::GetInstance()->GetTopBranch(configLink);

  if (!fBranch) {
    ostringstream err;
    err << '\'' << configLink << "' manager requested in the manager configuration but "
           "the bootstrap file does not specify a configuration file for it.";
    ERROR(err);
    throw utl::XMLParseException(err.str());
  }
}


//Generic data getters
//--------------------
FDsimG4XMLManager::Status 
FDsimG4XMLManager::GetData(TabulatedFunction& tabulatedReturnData,
                           const string& componentProperty,
                           const string& componentName) 
  const
{
  // Find where requested data resides in DOM
  Branch dataBranch =
    FindBranch(componentProperty, componentName);

  if (!dataBranch)
    return eNotFound;

  // Since this is tabulated data, retrieve the ordinate and abscissa
  // NOTE THAT IT IS ASSUMED THAT THE ABSCISSA IS LABELED "x" AND 
  // THE ORDINATE IS LABELED "y"
  Branch abscissaBranch = dataBranch.GetChild("x");
  if (!abscissaBranch)
    return eNotFound;
  vector<double> abscissaVector;
  abscissaBranch.GetData(abscissaVector);

  Branch ordinateBranch = dataBranch.GetChild("y");
  if (!ordinateBranch)
    return eNotFound;
  vector<double> ordinateVector;
  ordinateBranch.GetData(ordinateVector);

  // As a convenience, scale all ordinate values if a scale factor is given
  double scaleY;
  Branch scaleYBranch = dataBranch.GetChild("scaleY");
  if (scaleYBranch) {
    scaleYBranch.GetData(scaleY);
    for_each(ordinateVector.begin(), ordinateVector.end(), _1 *= scaleY);
  }

  for (unsigned int i = 0; i < abscissaVector.size(); ++i)
    tabulatedReturnData.PushBack(abscissaVector[i], ordinateVector[i]);

  return eFound;
}


Branch
FDsimG4XMLManager::FindBranch(const string& componentProperty,
                              const string& componentName)
  const
{
  if (!fBranch)
    return fBranch;

  Branch telescopeB = fBranch.GetChild("telescopeConfig");
  if (!telescopeB)
    return telescopeB;

  Branch NameB = telescopeB.GetChild(componentName);
  if (!NameB)
    return NameB;

  Branch dataB = NameB.GetChild(componentProperty);

  return dataB;
}