XMLConfigReader.cc

Go to the documentation of this file.

#include <iostream>
#include <cmath>
#include <algorithm>
#include <utility>
#include <array>

#include "L1Trigger/L1TMuonOverlap/interface/XMLConfigReader.h"
#include "L1Trigger/L1TMuonOverlap/interface/GoldenPattern.h"
#include "L1Trigger/L1TMuonOverlap/interface/OMTFinput.h"

#include "CondFormats/L1TObjects/interface/L1TMuonOverlapParams.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "xercesc/framework/StdOutFormatTarget.hpp"
#include "xercesc/framework/LocalFileFormatTarget.hpp"
#include "xercesc/parsers/XercesDOMParser.hpp"
#include "xercesc/dom/DOM.hpp"
#include "xercesc/dom/DOMException.hpp"
#include "xercesc/dom/DOMImplementation.hpp"
#include "xercesc/sax/HandlerBase.hpp"
#include "xercesc/util/XMLString.hpp"
#include "xercesc/util/PlatformUtils.hpp"
#include "xercesc/util/XercesDefs.hpp"
XERCES_CPP_NAMESPACE_USE

#include "L1Trigger/RPCTrigger/interface/RPCConst.h"
// XMLConfigReader
inline std::string _toString(XMLCh const* toTranscode) {
std::string tmp(xercesc::XMLString::transcode(toTranscode));
return tmp;
}

inline XMLCh*  _toDOMS(std::string temp) {
  XMLCh* buff = XMLString::transcode(temp.c_str());
  return  buff;
}
XMLConfigReader::XMLConfigReader(){

  //XMLPlatformUtils::Initialize();
  
  //parser = new XercesDOMParser(); 
  //parser->setValidationScheme(XercesDOMParser::Val_Auto);
  //parser->setDoNamespaces(false);

  //doc = 0;  
}

XMLConfigReader::~XMLConfigReader()
{
  //  delete parser;
  //XMLPlatformUtils::Terminate();
}
void XMLConfigReader::readLUT(l1t::LUT *lut,const L1TMuonOverlapParams & aConfig, const std::string & type){

  std::stringstream strStream;
  int totalInWidth = 7;//Number of bits used to address LUT
  int outWidth = 6;//Number of bits used to store LUT value

  if(type=="iCharge") outWidth = 1;
  if(type=="iEta") outWidth = 2;
  if(type=="iPt") outWidth = 9;
  if(type=="meanDistPhi"){
    outWidth = 11;
    totalInWidth = 14;
  }
  if(type=="pdf"){
    outWidth = 6;
    totalInWidth = 21;
  }
  
  strStream <<"#<header> V1 "<<totalInWidth<<" "<<outWidth<<" </header> "<<std::endl;

  auto const & aGPs = readPatterns(aConfig);

  unsigned int in = 0;
  int out = 0;
  for(auto it: aGPs){
    if(type=="iCharge") out = it->key().theCharge==-1 ? 0:1;
    if(type=="iEta") out = it->key().theEtaCode;
    if(type=="iPt") out = it->key().thePtCode;
    if(type=="meanDistPhi"){
      for(unsigned int iLayer = 0;iLayer<(unsigned) aConfig.nLayers();++iLayer){
    for(unsigned int iRefLayer=0;iRefLayer<(unsigned) aConfig.nRefLayers();++iRefLayer){
      out = (1<<(outWidth-1)) + it->meanDistPhiValue(iLayer,iRefLayer);
      strStream<<in<<" "<<out<<std::endl;
      ++in;
    }
      }
    }
    if(type=="pdf"){
      for(unsigned int iLayer = 0;iLayer<(unsigned)aConfig.nLayers();++iLayer){
    for(unsigned int iRefLayer=0;iRefLayer<(unsigned)aConfig.nRefLayers();++iRefLayer){
      for(unsigned int iPdf=0;iPdf<exp2(aConfig.nPdfAddrBits());++iPdf){
        out = it->pdfValue(iLayer,iRefLayer,iPdf);
        strStream<<in<<" "<<out<<std::endl;
        ++in;
      }
    }
      }
    }
    if(type!="meanDistPhi" && type!="pdf"){
      strStream<<in<<" "<<out<<std::endl;
      ++in;
    }
  }

  lut->read(strStream);
}
unsigned int XMLConfigReader::getPatternsVersion() const{

  if(!patternsFile.size()) return 0;

  unsigned int version=0;
  XMLPlatformUtils::Initialize();
  {
    XercesDOMParser parser;
    parser.setValidationScheme(XercesDOMParser::Val_Auto);
    parser.setDoNamespaces(false);
    
    parser.parse(patternsFile.c_str()); 
    xercesc::DOMDocument* doc = parser.getDocument();
    assert(doc);

    XMLCh *xmlOmtf=_toDOMS("OMTF");
    XMLCh *xmlVersion= _toDOMS("version");
    DOMNode *aNode = doc->getElementsByTagName(xmlOmtf)->item(0);
    DOMElement* aOMTFElement = static_cast<DOMElement *>(aNode);
    
    version = std::stoul(_toString(aOMTFElement->getAttribute(xmlVersion)), nullptr, 16);
    XMLString::release(&xmlOmtf);
    XMLString::release(&xmlVersion);
    parser.resetDocumentPool();
  }
  XMLPlatformUtils::Terminate();
  
  return version;
}
std::vector<std::shared_ptr<GoldenPattern>> XMLConfigReader::readPatterns(const L1TMuonOverlapParams & aConfig){

  aGPs.clear();
  
  XMLPlatformUtils::Initialize();

  XMLCh *xmlGP= _toDOMS("GP");
  std::array<XMLCh *,4> xmliPt= {{_toDOMS("iPt1"),_toDOMS("iPt2"),_toDOMS("iPt3"),_toDOMS("iPt4") }};

  {
    XercesDOMParser parser;
    parser.setValidationScheme(XercesDOMParser::Val_Auto);
    parser.setDoNamespaces(false);
    
    
    parser.parse(patternsFile.c_str()); 
    xercesc::DOMDocument* doc = parser.getDocument();
    assert(doc);
    
    unsigned int nElem = doc->getElementsByTagName(xmlGP)->getLength();
    if(nElem<1){
      edm::LogError("critical")<<"Problem parsing XML file "<<patternsFile<<std::endl;
      edm::LogError("critical")<<"No GoldenPattern items: GP found"<<std::endl;
      return aGPs;
    }
    
    DOMNode *aNode = 0;
    DOMElement* aGPElement = 0;
    unsigned int iGPNumber=0;
    
    for(unsigned int iItem=0;iItem<nElem;++iItem){
      aNode = doc->getElementsByTagName(xmlGP)->item(iItem);
      aGPElement = static_cast<DOMElement *>(aNode);
      
      std::unique_ptr<GoldenPattern> aGP;
      for(unsigned int index = 1;index<5;++index){
    if(aGPElement->getAttributeNode(xmliPt[index-1])) {
      aGP = buildGP(aGPElement, aConfig, index, iGPNumber);
      if(aGP){    
        aGPs.emplace_back(std::move(aGP));
        iGPNumber++;
      }
    }
    else{
      aGP = buildGP(aGPElement, aConfig);
      if(aGP){
        aGPs.emplace_back(std::move(aGP));
        iGPNumber++;
      }
      break;
    }
      }
    }
    
    // Reset the documents vector pool and release all the associated memory back to the system.
  //parser->resetDocumentPool();
    parser.resetDocumentPool();
  }
  XMLString::release(&xmlGP);
  XMLString::release(&xmliPt[0]);
  XMLString::release(&xmliPt[1]);
  XMLString::release(&xmliPt[2]);
  XMLString::release(&xmliPt[3]);


  XMLPlatformUtils::Terminate();

  return aGPs;
}
std::unique_ptr<GoldenPattern> XMLConfigReader::buildGP(DOMElement* aGPElement,
                     const L1TMuonOverlapParams & aConfig,
                     unsigned int index,
                     unsigned int aGPNumber){


  XMLCh *xmliEta= _toDOMS("iEta");
  //index 0 means no number at the end
  std::ostringstream stringStr;
  if (index>0) stringStr<<"iPt"<<index;
  else stringStr.str("iPt");
  XMLCh *xmliPt=_toDOMS(stringStr.str().c_str());
  stringStr.str("");
  if (index>0) stringStr<<"value"<<index;
  else stringStr.str("value");
  XMLCh *xmlValue=_toDOMS(stringStr.str().c_str());
  
  XMLCh *xmliCharge= _toDOMS("iCharge");
  XMLCh *xmlLayer= _toDOMS("Layer");
  XMLCh *xmlRefLayer= _toDOMS("RefLayer");
  XMLCh *xmlmeanDistPhi= _toDOMS("meanDistPhi");
  XMLCh *xmlPDF= _toDOMS("PDF");
  
  unsigned int iPt = std::atoi(_toString(aGPElement->getAttribute(xmliPt)).c_str());  
  int iEta = std::atoi(_toString(aGPElement->getAttribute(xmliEta)).c_str());
  int iCharge = std::atoi(_toString(aGPElement->getAttribute(xmliCharge)).c_str());
  int val = 0;
  unsigned int nLayers = aGPElement->getElementsByTagName(xmlLayer)->getLength();
  assert(nLayers==(unsigned) aConfig.nLayers());

  DOMNode *aNode = 0;
  DOMElement* aLayerElement = 0;
  DOMElement* aItemElement = 0;
  GoldenPattern::vector2D meanDistPhi2D(nLayers);
  GoldenPattern::vector1D pdf1D(exp2(aConfig.nPdfAddrBits()));
  GoldenPattern::vector3D pdf3D(aConfig.nLayers());
  GoldenPattern::vector2D pdf2D(aConfig.nRefLayers());

  if(iPt==0){
    GoldenPattern::vector1D meanDistPhi1D(aConfig.nRefLayers());
    meanDistPhi2D.assign(aConfig.nLayers(),meanDistPhi1D);
    pdf1D.assign(exp2(aConfig.nPdfAddrBits()),0);
    pdf2D.assign(aConfig.nRefLayers(),pdf1D);
    pdf3D.assign(aConfig.nLayers(),pdf2D);

    Key aKey(iEta,iPt,iCharge, aGPNumber);
    auto aGP = std::make_unique<GoldenPattern>(aKey,static_cast<const OMTFConfiguration*>(nullptr));
    aGP->setMeanDistPhi(meanDistPhi2D);
    aGP->setPdf(pdf3D);
    return aGP;
  }
  
  for(unsigned int iLayer=0;iLayer<nLayers;++iLayer){
    aNode = aGPElement->getElementsByTagName(xmlLayer)->item(iLayer);
    aLayerElement = static_cast<DOMElement *>(aNode); 
    unsigned int nItems = aLayerElement->getElementsByTagName(xmlRefLayer)->getLength();
    assert(nItems==(unsigned) aConfig.nRefLayers());
    GoldenPattern::vector1D meanDistPhi1D(nItems);
    for(unsigned int iItem=0;iItem<nItems;++iItem){
      aNode = aLayerElement->getElementsByTagName(xmlRefLayer)->item(iItem);
      aItemElement = static_cast<DOMElement *>(aNode); 
      val = std::atoi(_toString(aItemElement->getAttribute(xmlmeanDistPhi)).c_str());
      meanDistPhi1D[iItem] = val;
    }
    meanDistPhi2D[iLayer] = meanDistPhi1D;

    nItems = aLayerElement->getElementsByTagName(xmlPDF)->getLength();
    assert(nItems==aConfig.nRefLayers()*exp2(aConfig.nPdfAddrBits()));
    for(unsigned int iRefLayer=0;iRefLayer<(unsigned) aConfig.nRefLayers();++iRefLayer){
      pdf1D.assign(exp2(aConfig.nPdfAddrBits()),0);
      for(unsigned int iPdf=0;iPdf<exp2(aConfig.nPdfAddrBits());++iPdf){
    aNode = aLayerElement->getElementsByTagName(xmlPDF)->item(iRefLayer*exp2(aConfig.nPdfAddrBits())+iPdf);
    aItemElement = static_cast<DOMElement *>(aNode);
    val = std::atoi(_toString(aItemElement->getAttribute(xmlValue)).c_str());
    pdf1D[iPdf] = val;
      }
      pdf2D[iRefLayer] = pdf1D;
    }
    pdf3D[iLayer] = pdf2D;
  }

  Key aKey(iEta,iPt,iCharge, aGPNumber);
  auto aGP = std::make_unique<GoldenPattern>(aKey,static_cast<const OMTFConfiguration*>(nullptr));
  aGP->setMeanDistPhi(meanDistPhi2D);
  aGP->setPdf(pdf3D);

  XMLString::release(&xmliEta);
  XMLString::release(&xmliPt);
  XMLString::release(&xmliCharge);
  XMLString::release(&xmlLayer);
  XMLString::release(&xmlRefLayer);
  XMLString::release(&xmlmeanDistPhi);
  XMLString::release(&xmlPDF);
  XMLString::release(&xmlValue);

  return aGP;
}
std::vector<std::vector<int> > XMLConfigReader::readEvent(unsigned int iEvent,
                              unsigned int iProcessor,
                              bool readEta){

  return OMTFinput::vector2D();
 
}
void XMLConfigReader::readConfig(L1TMuonOverlapParams *aConfig) const{

  XMLPlatformUtils::Initialize();
  {
    XercesDOMParser parser;
    parser.setValidationScheme(XercesDOMParser::Val_Auto);
    parser.setDoNamespaces(false);

    XMLCh *xmlOMTF= _toDOMS("OMTF");
    XMLCh *xmlversion= _toDOMS("version");
    XMLCh *xmlGlobalData= _toDOMS("GlobalData");
    XMLCh *xmlnPdfAddrBits= _toDOMS("nPdfAddrBits");
    XMLCh *xmlnPdfValBits= _toDOMS("nPdfValBits");
    XMLCh *xmlnPhiBits= _toDOMS("nPhiBits");
    XMLCh *xmlnPhiBins= _toDOMS("nPhiBins");
    XMLCh *xmlnProcessors = _toDOMS("nProcessors");
    XMLCh *xmlnLogicRegions = _toDOMS("nLogicRegions");
    XMLCh *xmlnInputs= _toDOMS("nInputs");
    XMLCh *xmlnLayers= _toDOMS("nLayers");
    XMLCh *xmlnRefLayers= _toDOMS("nRefLayers");
    XMLCh *xmliProcessor= _toDOMS("iProcessor");
    XMLCh *xmlbarrelMin= _toDOMS("barrelMin");
    XMLCh *xmlbarrelMax= _toDOMS("barrelMax");
    XMLCh *xmlendcap10DegMin= _toDOMS("endcap10DegMin");
    XMLCh *xmlendcap10DegMax= _toDOMS("endcap10DegMax");
    XMLCh *xmlendcap20DegMin= _toDOMS("endcap20DegMin");
    XMLCh *xmlendcap20DegMax= _toDOMS("endcap20DegMax");
    XMLCh *xmlLayerMap = _toDOMS("LayerMap");
    XMLCh *xmlhwNumber = _toDOMS("hwNumber");
    XMLCh *xmllogicNumber = _toDOMS("logicNumber");
    XMLCh *xmlbendingLayer = _toDOMS("bendingLayer");
    XMLCh *xmlconnectedToLayer = _toDOMS("connectedToLayer");
    XMLCh *xmlRefLayerMap = _toDOMS("RefLayerMap");
    XMLCh *xmlrefLayer = _toDOMS("refLayer");
    XMLCh *xmlProcessor = _toDOMS("Processor");
    XMLCh *xmlRefLayer = _toDOMS("RefLayer");
    XMLCh *xmliRefLayer = _toDOMS("iRefLayer");
    XMLCh *xmliGlobalPhiStart = _toDOMS("iGlobalPhiStart");
    XMLCh *xmlRefHit = _toDOMS("RefHit");
    XMLCh *xmliRefHit = _toDOMS("iRefHit");
    XMLCh *xmliPhiMin = _toDOMS("iPhiMin");
    XMLCh *xmliPhiMax = _toDOMS("iPhiMax");
    XMLCh *xmliInput = _toDOMS("iInput");
    XMLCh *xmliRegion = _toDOMS("iRegion");
    XMLCh *xmlLogicRegion = _toDOMS("LogicRegion");
    XMLCh *xmlLayer = _toDOMS("Layer");
    XMLCh *xmliLayer = _toDOMS("iLayer");
    XMLCh *xmliFirstInput = _toDOMS("iFirstInput");
    XMLCh *xmlnHitsPerLayer = _toDOMS("nHitsPerLayer");
    XMLCh *xmlnRefHits = _toDOMS("nRefHits");
    XMLCh *xmlnTestRefHits = _toDOMS("nTestRefHits");
    XMLCh *xmlnGoldenPatterns = _toDOMS("nGoldenPatterns");
    XMLCh *xmlConnectionMap = _toDOMS("ConnectionMap");
    parser.parse(configFile.c_str()); 
    xercesc::DOMDocument* doc = parser.getDocument();
    assert(doc);
    unsigned int nElem = doc->getElementsByTagName(xmlOMTF)->getLength();
    if(nElem!=1){
      edm::LogError("critical")<<"Problem parsing XML file "<<configFile<<std::endl;
      assert(nElem==1);
    }
    DOMNode *aNode = doc->getElementsByTagName(xmlOMTF)->item(0);
    DOMElement* aOMTFElement = static_cast<DOMElement *>(aNode);
    
    unsigned int version = std::stoul(_toString(aOMTFElement->getAttribute(xmlversion)), nullptr, 16);
    aConfig->setFwVersion(version);
    
    nElem = aOMTFElement->getElementsByTagName(xmlGlobalData)->getLength();
    assert(nElem==1);
    aNode = aOMTFElement->getElementsByTagName(xmlGlobalData)->item(0);
    DOMElement* aElement = static_cast<DOMElement *>(aNode); 
    
    unsigned int nPdfAddrBits = std::atoi(_toString(aElement->getAttribute(xmlnPdfAddrBits)).c_str()); 
    unsigned int nPdfValBits =  std::atoi(_toString(aElement->getAttribute(xmlnPdfValBits)).c_str()); 
    unsigned int nHitsPerLayer =  std::atoi(_toString(aElement->getAttribute(xmlnHitsPerLayer)).c_str()); 
    unsigned int nPhiBits =  std::atoi(_toString(aElement->getAttribute(xmlnPhiBits)).c_str()); 
    unsigned int nPhiBins =  std::atoi(_toString(aElement->getAttribute(xmlnPhiBins)).c_str()); 

    unsigned int nRefHits =  std::atoi(_toString(aElement->getAttribute(xmlnRefHits)).c_str()); 
    unsigned int nTestRefHits =  std::atoi(_toString(aElement->getAttribute(xmlnTestRefHits)).c_str());
    unsigned int nProcessors =  std::atoi(_toString(aElement->getAttribute(xmlnProcessors)).c_str());
    unsigned int nLogicRegions =  std::atoi(_toString(aElement->getAttribute(xmlnLogicRegions)).c_str());
    unsigned int nInputs =  std::atoi(_toString(aElement->getAttribute(xmlnInputs)).c_str());
    unsigned int nLayers =  std::atoi(_toString(aElement->getAttribute(xmlnLayers)).c_str());
    unsigned int nRefLayers =  std::atoi(_toString(aElement->getAttribute(xmlnRefLayers)).c_str());
    unsigned int nGoldenPatterns =  std::atoi(_toString(aElement->getAttribute(xmlnGoldenPatterns)).c_str());
    
    std::vector<int> paramsVec(L1TMuonOverlapParams::GENERAL_NCONFIG);
    paramsVec[L1TMuonOverlapParams::GENERAL_ADDRBITS] = nPdfAddrBits;
    paramsVec[L1TMuonOverlapParams::GENERAL_VALBITS] = nPdfValBits;
    paramsVec[L1TMuonOverlapParams::GENERAL_HITSPERLAYER] = nHitsPerLayer;
    paramsVec[L1TMuonOverlapParams::GENERAL_PHIBITS] = nPhiBits;
    paramsVec[L1TMuonOverlapParams::GENERAL_PHIBINS] = nPhiBins;
    paramsVec[L1TMuonOverlapParams::GENERAL_NREFHITS] = nRefHits;
    paramsVec[L1TMuonOverlapParams::GENERAL_NTESTREFHITS] = nTestRefHits;
    paramsVec[L1TMuonOverlapParams::GENERAL_NPROCESSORS] = nProcessors;
    paramsVec[L1TMuonOverlapParams::GENERAL_NLOGIC_REGIONS] = nLogicRegions;
    paramsVec[L1TMuonOverlapParams::GENERAL_NINPUTS] = nInputs;
    paramsVec[L1TMuonOverlapParams::GENERAL_NLAYERS] = nLayers;
    paramsVec[L1TMuonOverlapParams::GENERAL_NREFLAYERS] = nRefLayers;
    paramsVec[L1TMuonOverlapParams::GENERAL_NGOLDENPATTERNS] = nGoldenPatterns;
    aConfig->setGeneralParams(paramsVec);
    
    std::vector<int> sectorsStart(3*nProcessors), sectorsEnd(3*nProcessors);
    nElem = aOMTFElement->getElementsByTagName(xmlConnectionMap)->getLength();
    DOMElement* aConnectionElement = 0;
    for(unsigned int i=0;i<nElem;++i){
      aNode = aOMTFElement->getElementsByTagName(xmlConnectionMap)->item(i);
      aConnectionElement = static_cast<DOMElement *>(aNode);
      unsigned int iProcessor = std::atoi(_toString(aConnectionElement->getAttribute(xmliProcessor)).c_str());
      unsigned int barrelMin = std::atoi(_toString(aConnectionElement->getAttribute(xmlbarrelMin)).c_str());
      unsigned int barrelMax = std::atoi(_toString(aConnectionElement->getAttribute(xmlbarrelMax)).c_str());
      unsigned int endcap10DegMin = std::atoi(_toString(aConnectionElement->getAttribute(xmlendcap10DegMin)).c_str());
      unsigned int endcap10DegMax = std::atoi(_toString(aConnectionElement->getAttribute(xmlendcap10DegMax)).c_str());
      unsigned int endcap20DegMin = std::atoi(_toString(aConnectionElement->getAttribute(xmlendcap20DegMin)).c_str());
      unsigned int endcap20DegMax = std::atoi(_toString(aConnectionElement->getAttribute(xmlendcap20DegMax)).c_str());
      
      sectorsStart[iProcessor] = barrelMin;
      sectorsStart[iProcessor + nProcessors] = endcap10DegMin;
      sectorsStart[iProcessor  + 2*nProcessors] = endcap20DegMin;
      
      sectorsEnd[iProcessor] = barrelMax;
      sectorsEnd[iProcessor + nProcessors] = endcap10DegMax;
      sectorsEnd[iProcessor + 2*nProcessors] = endcap20DegMax;       
    }  
    aConfig->setConnectedSectorsStart(sectorsStart);
    aConfig->setConnectedSectorsEnd(sectorsEnd);
    
    
    std::vector<L1TMuonOverlapParams::LayerMapNode> aLayerMapVec;
    L1TMuonOverlapParams::LayerMapNode aLayerMapNode;
    
    nElem = aOMTFElement->getElementsByTagName(xmlLayerMap)->getLength();
    DOMElement* aLayerElement = 0;
    for(unsigned int i=0;i<nElem;++i){
      aNode = aOMTFElement->getElementsByTagName(xmlLayerMap)->item(i);
      aLayerElement = static_cast<DOMElement *>(aNode); 
      unsigned int hwNumber = std::atoi(_toString(aLayerElement->getAttribute(xmlhwNumber)).c_str());
      unsigned int logicNumber = std::atoi(_toString(aLayerElement->getAttribute(xmllogicNumber)).c_str());
      unsigned int isBendingLayer = std::atoi(_toString(aLayerElement->getAttribute(xmlbendingLayer)).c_str());
      unsigned int iConnectedLayer = std::atoi(_toString(aLayerElement->getAttribute(xmlconnectedToLayer)).c_str());
      aLayerMapNode.logicNumber = logicNumber;
      aLayerMapNode.hwNumber = hwNumber;
      aLayerMapNode.connectedToLayer = iConnectedLayer;
      aLayerMapNode.bendingLayer = isBendingLayer;
      aLayerMapVec.push_back(aLayerMapNode);
    }
    aConfig->setLayerMap(aLayerMapVec);
    
    std::vector<L1TMuonOverlapParams::RefLayerMapNode> aRefLayerMapVec;
    L1TMuonOverlapParams::RefLayerMapNode aRefLayerNode;
    
    nElem = aOMTFElement->getElementsByTagName(xmlRefLayerMap)->getLength();
    DOMElement* aRefLayerElement = 0;
    for(unsigned int i=0;i<nElem;++i){
      aNode = aOMTFElement->getElementsByTagName(xmlRefLayerMap)->item(i);
      aRefLayerElement = static_cast<DOMElement *>(aNode); 
      unsigned int refLayer = std::atoi(_toString(aRefLayerElement->getAttribute(xmlrefLayer)).c_str());
      unsigned int logicNumber = std::atoi(_toString(aRefLayerElement->getAttribute(xmllogicNumber)).c_str());
      aRefLayerNode.refLayer = refLayer;
      aRefLayerNode.logicNumber = logicNumber;
      aRefLayerMapVec.push_back(aRefLayerNode);
    }
    aConfig->setRefLayerMap(aRefLayerMapVec);
    
    std::vector<int> aGlobalPhiStartVec(nProcessors*nRefLayers);
    
    std::vector<L1TMuonOverlapParams::RefHitNode> aRefHitMapVec(nProcessors*nRefHits);
    L1TMuonOverlapParams::RefHitNode aRefHitNode;
    
    std::vector<L1TMuonOverlapParams::LayerInputNode> aLayerInputMapVec(nProcessors*nLogicRegions*nLayers);  
    L1TMuonOverlapParams::LayerInputNode aLayerInputNode;
    
    nElem = aOMTFElement->getElementsByTagName(xmlProcessor)->getLength();
    assert(nElem==nProcessors);
    DOMElement* aProcessorElement = 0;
    for(unsigned int i=0;i<nElem;++i){
      aNode = aOMTFElement->getElementsByTagName(xmlProcessor)->item(i);
      aProcessorElement = static_cast<DOMElement *>(aNode); 
      unsigned int iProcessor = std::atoi(_toString(aProcessorElement->getAttribute(xmliProcessor)).c_str());
      unsigned int nElem1 = aProcessorElement->getElementsByTagName(xmlRefLayer)->getLength();
      assert(nElem1==nRefLayers);
      DOMElement* aRefLayerElement = 0;
      for(unsigned int ii=0;ii<nElem1;++ii){
    aNode = aProcessorElement->getElementsByTagName(xmlRefLayer)->item(ii);
    aRefLayerElement = static_cast<DOMElement *>(aNode); 
    unsigned int iRefLayer = std::atoi(_toString(aRefLayerElement->getAttribute(xmliRefLayer)).c_str());
    int iPhi = std::atoi(_toString(aRefLayerElement->getAttribute(xmliGlobalPhiStart)).c_str());
    aGlobalPhiStartVec[iRefLayer + iProcessor*nRefLayers] = iPhi;
      }
      nElem1 = aProcessorElement->getElementsByTagName(xmlRefHit)->getLength();
      assert( (iProcessor==0 && nElem1==nRefHits) || (iProcessor!=0 && nElem1==0) );
      DOMElement* aRefHitElement = 0;
      for(unsigned int ii=0;ii<nElem1;++ii){
    aNode = aProcessorElement->getElementsByTagName(xmlRefHit)->item(ii);
    aRefHitElement = static_cast<DOMElement *>(aNode); 
    unsigned int iRefHit = std::atoi(_toString(aRefHitElement->getAttribute(xmliRefHit)).c_str());
    int iPhiMin = std::atoi(_toString(aRefHitElement->getAttribute(xmliPhiMin)).c_str());
    int iPhiMax = std::atoi(_toString(aRefHitElement->getAttribute(xmliPhiMax)).c_str());
    unsigned int iInput = std::atoi(_toString(aRefHitElement->getAttribute(xmliInput)).c_str());
    unsigned int iRegion = std::atoi(_toString(aRefHitElement->getAttribute(xmliRegion)).c_str());
    unsigned int iRefLayer = std::atoi(_toString(aRefHitElement->getAttribute(xmliRefLayer)).c_str());
    
    aRefHitNode.iRefHit = iRefHit;
    aRefHitNode.iPhiMin = iPhiMin;
    aRefHitNode.iPhiMax = iPhiMax;
    aRefHitNode.iInput = iInput;
    aRefHitNode.iRegion = iRegion;
    aRefHitNode.iRefLayer = iRefLayer;
    for (unsigned int iProcessor=0; iProcessor<nProcessors; iProcessor++) aRefHitMapVec[iRefHit + iProcessor*nRefHits] = aRefHitNode;
      }
      unsigned int nElem2 = aProcessorElement->getElementsByTagName(xmlLogicRegion)->getLength();
      assert( (iProcessor==0 && nElem2==nLogicRegions) || (iProcessor!=0 && nElem2==0) );
      DOMElement* aRegionElement = 0;
      for(unsigned int ii=0;ii<nElem2;++ii){
    aNode = aProcessorElement->getElementsByTagName(xmlLogicRegion)->item(ii);
    aRegionElement = static_cast<DOMElement *>(aNode); 
    unsigned int iRegion = std::atoi(_toString(aRegionElement->getAttribute(xmliRegion)).c_str());
    unsigned int nElem3 = aRegionElement->getElementsByTagName(xmlLayer)->getLength();
    assert(nElem3==nLayers);
    DOMElement* aLayerElement = 0;
    for(unsigned int iii=0;iii<nElem3;++iii){
      aNode = aRegionElement->getElementsByTagName(xmlLayer)->item(iii);
      aLayerElement = static_cast<DOMElement *>(aNode); 
      unsigned int iLayer = std::atoi(_toString(aLayerElement->getAttribute(xmliLayer)).c_str());
      unsigned int iFirstInput = std::atoi(_toString(aLayerElement->getAttribute(xmliFirstInput)).c_str());
      unsigned int nInputs = std::atoi(_toString(aLayerElement->getAttribute(xmlnInputs)).c_str());
      aLayerInputNode.iLayer = iLayer;
      aLayerInputNode.iFirstInput = iFirstInput;
      aLayerInputNode.nInputs = nInputs;
      for (unsigned int iProcessor=0; iProcessor<nProcessors; ++iProcessor) aLayerInputMapVec[iLayer + iRegion*nLayers + iProcessor*nLayers*nLogicRegions] = aLayerInputNode;
    }
      }   
    }
    
    aConfig->setGlobalPhiStartMap(aGlobalPhiStartVec);
    aConfig->setLayerInputMap(aLayerInputMapVec);
    aConfig->setRefHitMap(aRefHitMapVec);
    
    // Reset the documents vector pool and release all the associated memory back to the system.
    parser.resetDocumentPool();
  

  XMLString::release(&xmlOMTF);
  XMLString::release(&xmlversion);
  XMLString::release(&xmlGlobalData);
  XMLString::release(&xmlnPdfAddrBits);
  XMLString::release(&xmlnPdfValBits);
  XMLString::release(&xmlnPhiBits);
  XMLString::release(&xmlnPhiBins);
  XMLString::release(&xmlnProcessors);
  XMLString::release(&xmlnLogicRegions);
  XMLString::release(&xmlnInputs);
  XMLString::release(&xmlnLayers);
  XMLString::release(&xmlnRefLayers);
  XMLString::release(&xmliProcessor);
  XMLString::release(&xmlbarrelMin);
  XMLString::release(&xmlbarrelMax);
  XMLString::release(&xmlendcap10DegMin);
  XMLString::release(&xmlendcap10DegMax);
  XMLString::release(&xmlendcap20DegMin);
  XMLString::release(&xmlendcap20DegMax);
  XMLString::release(&xmlLayerMap);
  XMLString::release(&xmlhwNumber);
  XMLString::release(&xmllogicNumber);
  XMLString::release(&xmlbendingLayer);
  XMLString::release(&xmlconnectedToLayer);
  XMLString::release(&xmlRefLayerMap);
  XMLString::release(&xmlrefLayer);
  XMLString::release(&xmlProcessor);
  XMLString::release(&xmlRefLayer);
  XMLString::release(&xmliRefLayer);
  XMLString::release(&xmliGlobalPhiStart);
  XMLString::release(&xmlRefHit);
  XMLString::release(&xmliRefHit);
  XMLString::release(&xmliPhiMin);
  XMLString::release(&xmliPhiMax);
  XMLString::release(&xmliInput);
  XMLString::release(&xmliRegion);
  XMLString::release(&xmlLogicRegion);
  XMLString::release(&xmlLayer);
  XMLString::release(&xmliLayer);
  XMLString::release(&xmliFirstInput);
  XMLString::release(&xmlnHitsPerLayer);
  XMLString::release(&xmlnRefHits);
  XMLString::release(&xmlnTestRefHits);
  XMLString::release(&xmlnGoldenPatterns);
  XMLString::release(&xmlConnectionMap);
  }
  XMLPlatformUtils::Terminate();
}

  //  xercesc::XercesDOMParser *parser;
  //  xercesc::DOMDocument* doc;