RHardwareProfileXMLManager.cc

Go to the documentation of this file.

#include <boost/lambda/lambda.hpp>

#include <rdet/RHardwareProfileXMLManager.h>
#include <utl/ErrorLogger.h>
#include <utl/TabulatedFunction.h>
#include <utl/TabulatedFunctionComplexLgAmpPhase.h>
#include <rdet/RDetector.h>
#include <rdet/RManagerRegister.h> 

using namespace rdet;
using namespace std;
using namespace utl;
using namespace det;
using namespace xercesc;
using namespace boost::lambda;

REGISTER_R_MANAGER("RHardwareProfileXMLManager", RHardwareProfileXMLManager)


void
RHardwareProfileXMLManager::Init(const string& configLink)
{
  VManager::Init(configLink);
}


// Meant for finding GetData("flatDemo", "Response")
Branch RHardwareProfileXMLManager::FindHardwareProfileProperty(const string& hardwareType,
                                                               const string& hardwareProp)
  const
{
  Branch hardwareB = FindHardwareProfile(hardwareType);
  
  // in case FindHardwareProfile could not find the desired branch
  if (!hardwareB)
    return hardwareB;
  
  Branch hardwarePropB = hardwareB.GetFirstChild();
  for ( ; hardwarePropB; hardwarePropB = hardwarePropB.GetNextSibling())
    if (hardwarePropB.GetName() == hardwareProp)
      return hardwarePropB;

  if (IsReportingErrors()) {
    ostringstream err;
    err << "Cannot find " + hardwareProp + " while parsing for hardwareType "+ hardwareType +" in XML." ;
    ERROR(err);
  }

  return hardwarePropB;  
}


// Find single hardware types ("flatDemo"):
Branch RHardwareProfileXMLManager::FindHardwareProfile(const string& hardwareType) 
  const
{
  Branch hardwareB = fBranch.GetFirstChild();
  for ( ; hardwareB; hardwareB = hardwareB.GetNextSibling())
    if (hardwareB.GetName() == "RResponse" &&
        hardwareB.GetAttributes().begin()->second == hardwareType)
      return hardwareB;

  if (IsReportingErrors()) {
    ostringstream err;
    err << "Cannot find specifications for '" + hardwareType + "' while XML parsing.";
    ERROR(err);
  }

  return hardwareB;
}


// propIndex not used, but to obey interface
VManager::Status
RHardwareProfileXMLManager::GetInternalData(TabulatedFunctionComplexLgAmpPhase& returnData,
                           const string& hardwareProp,
                           const string& hardwareType,
                           const IndexMap& /*propIndex*/)
  const
{

  // Find where requested data resides in DOM
  Branch dataBranch =
    FindHardwareProfileProperty(hardwareType, hardwareProp);  
  if (!dataBranch)
    return eNotFound;

  // Since this is tabulated data, retrieve the ordinate and abscissa
  // NOTE: it is assumed that the x, usually the frequency, is labelled "x" 
  // amplification can be given in one of three units: "LgAmp", "LgPow", "dB"
  // LgAmp means log10 of the amplification with regard to voltage,
  // LgPow means log10 of the amplification with regard to power,
  // dB means amplification in dB
  // the third column is the phase given in either degrees or radians
  Branch xB = dataBranch.GetChild("x");
  if (!xB)
    return eNotFound;

  vector<double> x;
  xB.GetData(x);

  // the branch AmpB has to be filled with vlues in log10(Amplitude); check for the three options and convert accordingly
  vector<double> lgAmp;
  Branch AmpB = dataBranch.GetChild("dB");
  if (AmpB) {
    AmpB.GetData(lgAmp);    
    for_each(lgAmp.begin(), lgAmp.end(), _1 *= 0.05); // convert dB to LgAmp
  } else {
    AmpB = dataBranch.GetChild("LgAmp");
    if (AmpB) {
      AmpB.GetData(lgAmp); // no conversion necessary
    } else {
      AmpB = dataBranch.GetChild("LgPow");
      if (AmpB) {
        AmpB.GetData(lgAmp);    
        for_each(lgAmp.begin(), lgAmp.end(), _1 *= 0.5); // convert LgPow to LgAmp
      } else      
        return eNotFound; // none of "dB", "LgAmp", "lgPow"
    }
  }

  Branch phaseB = dataBranch.GetChild("Phase");
  if (!phaseB)
    return eNotFound;

  vector<double> phase;
  phaseB.GetData(phase);

  const unsigned int n = x.size();
  returnData = TabulatedFunctionComplexLgAmpPhase();
  for (unsigned int i = 0; i < n; ++i)
    returnData.PushBack(x[i], ComplexLgAmpPhase(lgAmp[i],phase[i]));

  return eFound;
}