EcalPulseSymmCovariancesXMLTranslator Class Reference

#include <EcalPulseSymmCovariancesXMLTranslator.h>

Static Public Member Functions
static int	readXML (const std::string &filename, EcalCondHeader &header, EcalPulseSymmCovariances &record)
static int	writeXML (const std::string &filename, const EcalCondHeader &header, const EcalPulseSymmCovariances &record)

Static Private Member Functions
static std::string	dumpXML (const EcalCondHeader &header, const EcalPulseSymmCovariances &record)

Detailed Description

Translates a EcalGainRatio record to XML and vice versa

Author

Francesco RUBBO

Version

Id

EcalPulseSymmCovariancesXMLTranslator.h,v 0 2010/04/26 fay Exp

Date

26 Apr 2010

Definition at line 15 of file EcalPulseSymmCovariancesXMLTranslator.h.

Member Function Documentation

dumpXML()

std::string EcalPulseSymmCovariancesXMLTranslator::dumpXML ( const EcalCondHeader &  header, const EcalPulseSymmCovariances &  record  )

staticprivate

Definition at line 105 of file EcalPulseSymmCovariancesXMLTranslator.cc.

References hcalRecHitTable_cff::doc, GCP_Ntuples_cfg::dump, RecoTauValidation_cfi::header, mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, EBDetId::kSizeForDenseIndexing, EEDetId::kSizeForDenseIndexing, EBDetId::MIN_HASH, conifer::pow(), xuti::PulseSymmCovariances_tag(), AlCaHarvesting_cff::record, AlCaHLTBitMon_QueryRunRegistry::string, EcalPulseShape::TEMPLATESAMPLES, cms::xerces::toString(), EBDetId::unhashIndex(), EEDetId::unhashIndex(), cms::xerces::uStr(), EEDetId::validHashIndex(), xuti::writeCell(), xuti::writeHeader(), xuti::WriteNodeWithValue(), and convertToRaw::writer.

                                                                                                   {
  unique_ptr<DOMImplementation> impl(DOMImplementationRegistry::getDOMImplementation(cms::xerces::uStr("LS").ptr()));

  DOMLSSerializer* writer = impl->createLSSerializer();
  if (writer->getDomConfig()->canSetParameter(XMLUni::fgDOMWRTFormatPrettyPrint, true))
    writer->getDomConfig()->setParameter(XMLUni::fgDOMWRTFormatPrettyPrint, true);

  DOMDocumentType* doctype = impl->createDocumentType(cms::xerces::uStr("XML").ptr(), nullptr, nullptr);
  DOMDocument* doc = impl->createDocument(nullptr, cms::xerces::uStr(PulseSymmCovariances_tag.c_str()).ptr(), doctype);

  DOMElement* root = doc->getDocumentElement();

  std::vector<std::string> pulsecov_tag(static_cast<size_t>(std::pow(EcalPulseShape::TEMPLATESAMPLES, 2)));
  for (int k = 0; k < std::pow(EcalPulseShape::TEMPLATESAMPLES, 2); ++k) {
    int i = k / EcalPulseShape::TEMPLATESAMPLES;
    int j = k % EcalPulseShape::TEMPLATESAMPLES;
    int linK = -1;
    if (j >= i)
      linK = j + (EcalPulseShape::TEMPLATESAMPLES - 1) * i;
    else
      linK = i + (EcalPulseShape::TEMPLATESAMPLES - 1) * j;
    pulsecov_tag[linK] = Form("samplecov_%d_%d", i, j);
  }

  xuti::writeHeader(root, header);
  if (record.barrelItems().empty())
    return std::string();
  for (int cellid = EBDetId::MIN_HASH; cellid < EBDetId::kSizeForDenseIndexing; ++cellid) {
    uint32_t rawid = EBDetId::unhashIndex(cellid);

    if (record.find(rawid) == record.end())
      continue;

    DOMElement* cellnode = writeCell(root, rawid);

    for (int k = 0; k < std::pow(EcalPulseShape::TEMPLATESAMPLES, 2); ++k) {
      WriteNodeWithValue(cellnode, pulsecov_tag[k], record[rawid].covval[k]);
    }
  }
  if (record.endcapItems().empty())
    return std::string();
  for (int cellid = 0; cellid < EEDetId::kSizeForDenseIndexing; ++cellid) {
    if (!EEDetId::validHashIndex(cellid))
      continue;

    uint32_t rawid = EEDetId::unhashIndex(cellid);

    if (record.find(rawid) == record.end())
      continue;

    DOMElement* cellnode = writeCell(root, rawid);

    for (int k = 0; k < std::pow(EcalPulseShape::TEMPLATESAMPLES, 2); ++k) {
      WriteNodeWithValue(cellnode, pulsecov_tag[k], record[rawid].covval[k]);
    }
  }

  std::string dump = cms::xerces::toString(writer->writeToString(root));
  doc->release();
  doctype->release();
  writer->release();

  return dump;
}

readXML()

int EcalPulseSymmCovariancesXMLTranslator::readXML ( const std::string &  filename, EcalCondHeader &  header, EcalPulseSymmCovariances &  record  )

static

Definition at line 24 of file EcalPulseSymmCovariancesXMLTranslator.cc.

References xuti::Cell_tag(), officialStyle::chan, gather_cfg::cout, corrVsCorr::filename, xuti::getChildNode(), xuti::GetNodeData(), RecoTauValidation_cfi::header, mps_fire::i, createfilelist::int, dqmiolumiharvest::j, dqmdumpme::k, writedatasetfile::parser, xuti::readCellId(), xuti::readHeader(), AlCaHarvesting_cff::record, EcalPulseShape::TEMPLATESAMPLES, cms::concurrency::xercesInitialize(), and cms::concurrency::xercesTerminate().

                                                                                     {
  cms::concurrency::xercesInitialize();

  XercesDOMParser* parser = new XercesDOMParser;
  parser->setValidationScheme(XercesDOMParser::Val_Never);
  parser->setDoNamespaces(false);
  parser->setDoSchema(false);

  parser->parse(filename.c_str());

  DOMDocument* xmlDoc = parser->getDocument();
  if (!xmlDoc) {
    std::cout << "EcalPulseSymmCovariancesXMLTranslator::Error parsing document" << std::endl;
    return -1;
  }

  DOMElement* elementRoot = xmlDoc->getDocumentElement();

  xuti::readHeader(elementRoot, header);

  DOMNode* cellnode = getChildNode(elementRoot, Cell_tag);

  int chan = 0;
  while (cellnode) {
    //    std::cout << " Channel " << chan << std::endl;
    float covs[EcalPulseShape::TEMPLATESAMPLES * (EcalPulseShape::TEMPLATESAMPLES + 1) / 2];

    DetId detid = readCellId(dynamic_cast<DOMElement*>(cellnode));
    //    std::cout << " readCell Id Channel " << chan << " tag " << mean12_tag << std::endl;

    std::vector<std::string> pulsecov_tag(
        static_cast<size_t>(EcalPulseShape::TEMPLATESAMPLES * (EcalPulseShape::TEMPLATESAMPLES + 1) / 2));
    for (int i = 0; i < EcalPulseShape::TEMPLATESAMPLES; ++i) {
      for (int j = 0; j < EcalPulseShape::TEMPLATESAMPLES; ++j) {
        int k = -1;
        if (j >= i)
          k = j + (EcalPulseShape::TEMPLATESAMPLES - 1) * i;
        else
          k = i + (EcalPulseShape::TEMPLATESAMPLES - 1) * j;
        pulsecov_tag[k] = Form("samplecov_%d_%d", i, j);
      }
    }

    DOMNode** covs_node =
        new DOMNode*[int(EcalPulseShape::TEMPLATESAMPLES * (EcalPulseShape::TEMPLATESAMPLES + 1) / 2)];
    for (int k = 0; k < EcalPulseShape::TEMPLATESAMPLES * (EcalPulseShape::TEMPLATESAMPLES + 1) / 2; ++k) {
      covs_node[k] = getChildNode(cellnode, pulsecov_tag[k]);
      GetNodeData(covs_node[k], covs[k]);
    }

    for (int k = 0; k < EcalPulseShape::TEMPLATESAMPLES * (EcalPulseShape::TEMPLATESAMPLES + 1) / 2; ++k)
      record[detid].covval[k] = covs[k];

    cellnode = cellnode->getNextSibling();

    while (cellnode && cellnode->getNodeType() != DOMNode::ELEMENT_NODE)
      cellnode = cellnode->getNextSibling();
    chan++;
  }

  delete parser;
  cms::concurrency::xercesTerminate();
  std::cout << " nb of channels found in xml file " << chan << std::endl;
  return 0;
}

writeXML()

int EcalPulseSymmCovariancesXMLTranslator::writeXML ( const std::string &  filename, const EcalCondHeader &  header, const EcalPulseSymmCovariances &  record  )

static

Definition at line 92 of file EcalPulseSymmCovariancesXMLTranslator.cc.

References EcalCondTools::dumpXML(), corrVsCorr::filename, compareTotals::fs, RecoTauValidation_cfi::header, MillePedeFileConverter_cfg::out, AlCaHarvesting_cff::record, cms::concurrency::xercesInitialize(), and cms::concurrency::xercesTerminate().

                                                                                            {
  cms::concurrency::xercesInitialize();

  std::fstream fs(filename.c_str(), ios::out);
  fs << dumpXML(header, record);

  cms::concurrency::xercesTerminate();

  return 0;
}