MiscalibReaderXML.cc

Go to the documentation of this file.

#include "CalibCalorimetry/CaloMiscalibTools/interface/MiscalibReaderFromXML.h"
#include "CalibCalorimetry/CaloMiscalibTools/interface/MiscalibReaderFromXMLDomUtils.h"
#include "CalibCalorimetry/CaloMiscalibTools/interface/CaloMiscalibMap.h"
#include "FWCore/Concurrency/interface/Xerces.h"
#include <stdio.h>
#include <stdlib.h>
#include <string>
#include <vector>

using namespace xercesc;

int MiscalibReaderFromXML::s_numberOfInstances = 0; //to check that there is only 1 instance

inline std::string _toString(const XMLCh *toTranscode){
    std::string tmp(XMLString::transcode(toTranscode));
    return tmp;
}

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



MiscalibReaderFromXML::MiscalibReaderFromXML(CaloMiscalibMap & caloMap):caloMap_(caloMap){
      
         
    try { 
        //std::cout << "Xerces-c initialization Number "
        //<< s_numberOfInstances<<std::endl;
        if (s_numberOfInstances==0) 
        cms::concurrency::xercesInitialize();  
    }
    catch (const XMLException& e) {
        std::cout << "Xerces-c error in initialization \n"
        << "Exception message is:  \n"
        << _toString(e.getMessage()) <<std::endl;
        // throw an exception here
    }
 
    ++s_numberOfInstances;
    
    
    
    
}


int MiscalibReaderFromXML::getIntAttribute(DOMNamedNodeMap * attribute, const std::string & attribute_name)
{
bool well_formed_string;
int retval = MiscalibReaderFromXMLDomUtils::getIntAttribute(attribute,attribute_name,well_formed_string);
if(!well_formed_string) std::cout << "MiscalibReaderFromXML::getIntAttribute PROBLEMS ...!!!" << std::endl;

  return retval;
}


double MiscalibReaderFromXML::getScalingFactor(XERCES_CPP_NAMESPACE::DOMNamedNodeMap *attribute)
{
return MiscalibReaderFromXML::getFloatAttribute(attribute,"scale_factor");
}


double MiscalibReaderFromXML::getFloatAttribute(DOMNamedNodeMap * attribute, const std::string & attribute_name)
{
bool well_formed_string;
double retval = MiscalibReaderFromXMLDomUtils::getFloatAttribute(attribute,attribute_name,well_formed_string);
if(!well_formed_string) std::cout << "MiscalibReaderFromXML::getFloatAttribute PROBLEMS ...!!!" << std::endl;
  
  return retval;
}


bool MiscalibReaderFromXML::parseXMLMiscalibFile(std::string configFile){

    XercesDOMParser* parser = new XercesDOMParser;     
    parser->setValidationScheme(XercesDOMParser::Val_Auto);
    parser->setDoNamespaces(false);
    parser->parse(configFile.c_str()); 
    DOMDocument* doc = parser->getDocument();
    assert(doc);

        unsigned int linkTagsNum = doc->getElementsByTagName(_toDOMS("Cell"))->getLength();
        // The following should be on LogInfo
        //std::cout << "Read number of Cells = " << linkTagsNum << std::endl;

        if(linkTagsNum==0) std::cout <<"Number of Cells in file is 0 - probably bad file format"<<std::endl;

        int count=0;    
    for (unsigned int i=0; i<linkTagsNum; i++){
            
        
        DOMNode* linkNode = doc->getElementsByTagName(_toDOMS("Cell"))->item(i);
        if (! linkNode){
            std::cout<<"Node LINK does not exist, i="<<i<<std::endl;
            return true;
        }
        DOMElement* linkElement = static_cast<DOMElement *>(linkNode);          
        if (! linkElement){
            std::cout<<"Element LINK does not exist, i="<<i<<std::endl;
            return true;         
        }


                DOMNamedNodeMap *attributes = linkNode->getAttributes();
        double scalingfactor= getScalingFactor(attributes);
        
        
        
        DetId cell = parseCellEntry(attributes);
        
        if(cell!= DetId(0)) 
        { 
        count++;
        caloMap_.addCell(cell,scalingfactor);
        } else  
        {
          //        std::cout << "Null received" << std::endl;
        }
        
    }
 
        // The following should be on LogInfo
        // std::cout << "Number of good Cells = " << count << std::endl;
    return false;

}