CMS 3D CMS Logo
EcalPulseShapesXMLTranslator.cc
Go to the documentation of this file.
1 #include <iostream>
2 #include <sstream>
3 #include <fstream>
4 #include <xercesc/dom/DOMNode.hpp>
5 #include <xercesc/dom/DOM.hpp>
6 #include <xercesc/parsers/XercesDOMParser.hpp>
7 #include "Utilities/Xerces/interface/Xerces.h"
8 #include "Utilities/Xerces/interface/XercesStrUtils.h"
9 #include <xercesc/util/XMLString.hpp>
10 #include <xercesc/sax/SAXException.hpp>
11 #include <xercesc/framework/LocalFileFormatTarget.hpp>
12 
13 #include "CondFormats/EcalObjects/interface/EcalPulseShapes.h"
14 #include "CondTools/Ecal/interface/EcalPulseShapesXMLTranslator.h"
15 #include "CondTools/Ecal/interface/DOMHelperFunctions.h"
16 #include "CondTools/Ecal/interface/XMLTags.h"
17 
18 using namespace XERCES_CPP_NAMESPACE;
19 using namespace xuti;
20 using namespace std;
21 
22 int EcalPulseShapesXMLTranslator::readXML(const std::string& filename,
23  EcalCondHeader& header,
24  EcalPulseShapes& record) {
25  cms::concurrency::xercesInitialize();
26 
27  XercesDOMParser* parser = new XercesDOMParser;
28  parser->setValidationScheme(XercesDOMParser::Val_Never);
29  parser->setDoNamespaces(false);
30  parser->setDoSchema(false);
31 
32  parser->parse(filename.c_str());
33 
34  DOMDocument* xmlDoc = parser->getDocument();
35  if (!xmlDoc) {
36  std::cout << "EcalPulseShapesXMLTranslator::Error parsing document" << std::endl;
37  return -1;
38  }
39 
40  DOMElement* elementRoot = xmlDoc->getDocumentElement();
41 
42  xuti::readHeader(elementRoot, header);
43 
44  DOMNode* cellnode = getChildNode(elementRoot, Cell_tag);
45 
46  int chan = 0;
47  while (cellnode) {
48  // std::cout << " Channel " << chan << std::endl;
49  float samples[EcalPulseShape::TEMPLATESAMPLES];
50 
51  DetId detid = readCellId(dynamic_cast<DOMElement*>(cellnode));
52  // std::cout << " readCell Id Channel " << chan << " tag " << mean12_tag << std::endl;
53 
54  DOMNode* sample0_node = getChildNode(cellnode, sample0_tag);
55  GetNodeData(sample0_node, samples[0]);
56  // std::cout << " tag " << sample0_tag << " sample0 " << sample0 << std::endl;
57 
58  DOMNode* sample1_node = getChildNode(cellnode, sample1_tag);
59  GetNodeData(sample1_node, samples[1]);
60  // std::cout << " tag " << sample1_tag << " sample1 " << sample1 << std::endl;
61 
62  DOMNode* sample2_node = getChildNode(cellnode, sample2_tag);
63  GetNodeData(sample2_node, samples[2]);
64  // std::cout << " tag " << sample2_tag << " sample2 " << sample2 << std::endl;
65 
66  DOMNode* sample3_node = getChildNode(cellnode, sample3_tag);
67  GetNodeData(sample3_node, samples[3]);
68  // std::cout << " tag " << sample3_tag << " sample3 " << sample3 << std::endl;
69 
70  DOMNode* sample4_node = getChildNode(cellnode, sample4_tag);
71  GetNodeData(sample4_node, samples[4]);
72  // std::cout << " tag " << sample4_tag << " sample4 " << sample4 << std::endl;
73 
74  DOMNode* sample5_node = getChildNode(cellnode, sample5_tag);
75  GetNodeData(sample5_node, samples[5]);
76  // std::cout << " tag " << sample5_tag << " sample5 " << sample5 << std::endl;
77 
78  DOMNode* sample6_node = getChildNode(cellnode, sample6_tag);
79  GetNodeData(sample6_node, samples[6]);
80  // std::cout << " tag " << sample6_tag << " sample6 " << sample6 << std::endl;
81 
82  DOMNode* sample7_node = getChildNode(cellnode, sample7_tag);
83  GetNodeData(sample7_node, samples[7]);
84  // std::cout << " tag " << sample7_tag << " sample7 " << sample7 << std::endl;
85 
86  DOMNode* sample8_node = getChildNode(cellnode, sample8_tag);
87  GetNodeData(sample8_node, samples[8]);
88  // std::cout << " tag " << sample8_tag << " sample8 " << sample8 << std::endl;
89 
90  DOMNode* sample9_node = getChildNode(cellnode, sample9_tag);
91  GetNodeData(sample9_node, samples[9]);
92  // std::cout << " tag " << sample9_tag << " sample9 " << sample9 << std::endl;
93 
94  DOMNode* sample10_node = getChildNode(cellnode, sample10_tag);
95  GetNodeData(sample10_node, samples[10]);
96  // std::cout << " tag " << sample10_tag << " sample10 " << sample10 << std::endl;
97 
98  DOMNode* sample11_node = getChildNode(cellnode, sample11_tag);
99  GetNodeData(sample11_node, samples[11]);
100  // std::cout << " tag " << sample11_tag << " sample11 " << sample11 << std::endl;
101 
102  for (int s = 0; s < EcalPulseShape::TEMPLATESAMPLES; ++s)
103  record[detid].pdfval[s] = samples[s];
104 
105  cellnode = cellnode->getNextSibling();
106 
107  while (cellnode && cellnode->getNodeType() != DOMNode::ELEMENT_NODE)
108  cellnode = cellnode->getNextSibling();
109  chan++;
110  }
111 
112  delete parser;
113  cms::concurrency::xercesTerminate();
114  std::cout << " nb of channels found in xml file " << chan << std::endl;
115  return 0;
116 }
117 
118 int EcalPulseShapesXMLTranslator::writeXML(const std::string& filename,
119  const EcalCondHeader& header,
120  const EcalPulseShapes& record) {
121  cms::concurrency::xercesInitialize();
122 
123  std::fstream fs(filename.c_str(), ios::out);
124  fs << dumpXML(header, record);
125 
126  cms::concurrency::xercesTerminate();
127 
128  return 0;
129 }
130 
131 std::string EcalPulseShapesXMLTranslator::dumpXML(const EcalCondHeader& header, const EcalPulseShapes& record) {
132  unique_ptr<DOMImplementation> impl(DOMImplementationRegistry::getDOMImplementation(cms::xerces::uStr("LS").ptr()));
133 
134  DOMLSSerializer* writer = impl->createLSSerializer();
135  if (writer->getDomConfig()->canSetParameter(XMLUni::fgDOMWRTFormatPrettyPrint, true))
136  writer->getDomConfig()->setParameter(XMLUni::fgDOMWRTFormatPrettyPrint, true);
137 
138  DOMDocumentType* doctype = impl->createDocumentType(cms::xerces::uStr("XML").ptr(), nullptr, nullptr);
139  DOMDocument* doc = impl->createDocument(nullptr, cms::xerces::uStr(PulseShapes_tag.c_str()).ptr(), doctype);
140 
141  DOMElement* root = doc->getDocumentElement();
142 
143  xuti::writeHeader(root, header);
144  if (record.barrelItems().empty())
145  return std::string();
146  for (int cellid = EBDetId::MIN_HASH; cellid < EBDetId::kSizeForDenseIndexing; ++cellid) {
147  uint32_t rawid = EBDetId::unhashIndex(cellid);
148 
149  if (record.find(rawid) == record.end())
150  continue;
151  if (!record[rawid].pdfval[5])
152  continue;
153 
154  DOMElement* cellnode = writeCell(root, rawid);
155 
156  WriteNodeWithValue(cellnode, sample0_tag, record[rawid].pdfval[0]);
157  WriteNodeWithValue(cellnode, sample1_tag, record[rawid].pdfval[1]);
158  WriteNodeWithValue(cellnode, sample2_tag, record[rawid].pdfval[2]);
159  WriteNodeWithValue(cellnode, sample3_tag, record[rawid].pdfval[3]);
160  WriteNodeWithValue(cellnode, sample4_tag, record[rawid].pdfval[4]);
161  WriteNodeWithValue(cellnode, sample5_tag, record[rawid].pdfval[5]);
162  WriteNodeWithValue(cellnode, sample6_tag, record[rawid].pdfval[6]);
163  WriteNodeWithValue(cellnode, sample7_tag, record[rawid].pdfval[7]);
164  WriteNodeWithValue(cellnode, sample8_tag, record[rawid].pdfval[8]);
165  WriteNodeWithValue(cellnode, sample9_tag, record[rawid].pdfval[9]);
166  WriteNodeWithValue(cellnode, sample10_tag, record[rawid].pdfval[10]);
167  WriteNodeWithValue(cellnode, sample11_tag, record[rawid].pdfval[11]);
168  }
169 
170  if (record.endcapItems().empty())
171  return std::string();
172  for (int cellid = 0; cellid < EEDetId::kSizeForDenseIndexing; ++cellid) {
173  if (!EEDetId::validHashIndex(cellid))
174  continue;
175 
176  uint32_t rawid = EEDetId::unhashIndex(cellid);
177 
178  if (record.find(rawid) == record.end())
179  continue;
180  if (!record[rawid].pdfval[5])
181  continue;
182 
183  DOMElement* cellnode = writeCell(root, rawid);
184 
185  WriteNodeWithValue(cellnode, sample0_tag, record[rawid].pdfval[0]);
186  WriteNodeWithValue(cellnode, sample1_tag, record[rawid].pdfval[1]);
187  WriteNodeWithValue(cellnode, sample2_tag, record[rawid].pdfval[2]);
188  WriteNodeWithValue(cellnode, sample3_tag, record[rawid].pdfval[3]);
189  WriteNodeWithValue(cellnode, sample4_tag, record[rawid].pdfval[4]);
190  WriteNodeWithValue(cellnode, sample5_tag, record[rawid].pdfval[5]);
191  WriteNodeWithValue(cellnode, sample6_tag, record[rawid].pdfval[6]);
192  WriteNodeWithValue(cellnode, sample7_tag, record[rawid].pdfval[7]);
193  WriteNodeWithValue(cellnode, sample8_tag, record[rawid].pdfval[8]);
194  WriteNodeWithValue(cellnode, sample9_tag, record[rawid].pdfval[9]);
195  WriteNodeWithValue(cellnode, sample10_tag, record[rawid].pdfval[10]);
196  WriteNodeWithValue(cellnode, sample11_tag, record[rawid].pdfval[11]);
197  }
198 
199  std::string dump = cms::xerces::toString(writer->writeToString(root));
200  doc->release();
201  doctype->release();
202  writer->release();
203 
204  return dump;
205 }
xuti::getChildNode
xercesc::DOMNode * getChildNode(xercesc::DOMNode *node, const std::string &nodename)
get the child of
XercesStrUtils.h
common_cff.doc
doc
Definition: common_cff.py:54
EcalCondTools.dumpXML
def dumpXML(db, tag, since, filename='dump.xml')
Definition: EcalCondTools.py:44
xuti::writeHeader
void writeHeader(xercesc::DOMNode *parentNode, const EcalCondHeader &header)
write
Definition: DOMHelperFunctions.cc:151
xuti::writeCell
xercesc::DOMElement * writeCell(xercesc::DOMNode *node, const DetId &detid)
Append a Cell node with attributes to.
Definition: DOMHelperFunctions.cc:71
xuti::sample4_tag
const std::string sample4_tag("sample_4")
xuti::sample1_tag
const std::string sample1_tag("sample_1")
EEDetId::unhashIndex
static EEDetId unhashIndex(int hi)
Definition: EEDetId.cc:65
gather_cfg.cout
cout
Definition: gather_cfg.py:144
EBDetId::unhashIndex
static EBDetId unhashIndex(int hi)
get a DetId from a compact index for arrays
Definition: EBDetId.h:110
EcalPulseShapesXMLTranslator::readXML
static int readXML(const std::string &filename, EcalCondHeader &header, EcalPulseShapes &record)
Definition: EcalPulseShapesXMLTranslator.cc:22
GlobalPosition_Frontier_DevDB_cff.record
record
Definition: GlobalPosition_Frontier_DevDB_cff.py:10
writedatasetfile.parser
parser
Definition: writedatasetfile.py:7
cscNeutronWriter_cfi.writer
writer
Definition: cscNeutronWriter_cfi.py:6
EcalCondObjectContainer
Definition: EcalCondObjectContainer.h:13
EgammaValidation_cff.samples
samples
Definition: EgammaValidation_cff.py:19
DetId
Definition: DetId.h:17
alignCSCRings.s
s
Definition: alignCSCRings.py:92
cms::xerces::toString
std::string toString(XMLCh const *toTranscode)
Definition: XercesStrUtils.h:34
xuti::sample6_tag
const std::string sample6_tag("sample_6")
EcalPulseShapesXMLTranslator::writeXML
static int writeXML(const std::string &filename, const EcalCondHeader &header, const EcalPulseShapes &record)
Definition: EcalPulseShapesXMLTranslator.cc:118
EEDetId::kSizeForDenseIndexing
Definition: EEDetId.h:329
xuti::WriteNodeWithValue
void WriteNodeWithValue(xercesc::DOMNode *parentNode, const std::string &tag, const T &value)
write a node with
Definition: DOMHelperFunctions.h:46
XERCES_CPP_NAMESPACE
Definition: XmlConfigParser.h:12
corrVsCorr.filename
filename
Definition: corrVsCorr.py:123
EcalPulseShapesXMLTranslator::dumpXML
static std::string dumpXML(const EcalCondHeader &header, const EcalPulseShapes &record)
Definition: EcalPulseShapesXMLTranslator.cc:131
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
xuti::sample7_tag
const std::string sample7_tag("sample_7")
EcalPulseShapes.h
xuti::sample11_tag
const std::string sample11_tag("sample_11")
DOMHelperFunctions.h
EcalCondHeader
Definition: EcalCondHeader.h:16
xuti::sample0_tag
const std::string sample0_tag("sample_0")
FrontierConditions_GlobalTag_cff.dump
dump
Definition: FrontierConditions_GlobalTag_cff.py:12
Xerces.h
cms::concurrency::xercesTerminate
void xercesTerminate()
Definition: Xerces.cc:23
xuti
Definition: DOMHelperFunctions.h:21
root
Definition: RooFitFunction.h:10
xuti::PulseShapes_tag
const std::string PulseShapes_tag("EcalPulseShapes")
EEDetId::validHashIndex
static bool validHashIndex(int i)
Definition: EEDetId.h:239
xuti::sample3_tag
const std::string sample3_tag("sample_3")
cms::xerces::uStr
ZStr< XMLCh > uStr(char const *str)
Definition: XercesStrUtils.h:30
xuti::sample9_tag
const std::string sample9_tag("sample_9")
xuti::sample5_tag
const std::string sample5_tag("sample_5")
EcalPulseShapesXMLTranslator.h
impl
Definition: trackAlgoPriorityOrder.h:18
std
Definition: JetResolutionObject.h:76
xuti::sample2_tag
const std::string sample2_tag("sample_2")
officialStyle.chan
chan
lumi = TPaveText(lowX+0.38, lowY+0.061, lowX+0.45, lowY+0.161, "NDC") lumi.SetBorderSize( 0 ) lumi....
Definition: officialStyle.py:106
EcalPulseShape::TEMPLATESAMPLES
static const int TEMPLATESAMPLES
Definition: EcalPulseShapes.h:10
EBDetId::kSizeForDenseIndexing
Definition: EBDetId.h:155
RecoTauValidation_cfi.header
header
Definition: RecoTauValidation_cfi.py:292
MillePedeFileConverter_cfg.out
out
Definition: MillePedeFileConverter_cfg.py:31
xuti::readCellId
const DetId readCellId(xercesc::DOMElement *node)
Assuming.
Definition: DOMHelperFunctions.cc:24
EBDetId::MIN_HASH
static const int MIN_HASH
Definition: EBDetId.h:149
XMLTags.h
xuti::sample8_tag
const std::string sample8_tag("sample_8")
cms::concurrency::xercesInitialize
void xercesInitialize()
Definition: Xerces.cc:18
xuti::readHeader
void readHeader(xercesc::DOMNode *parentNode, EcalCondHeader &header)
read header from
Definition: DOMHelperFunctions.cc:164
xuti::Cell_tag
const std::string Cell_tag("cell")
xuti::sample10_tag
const std::string sample10_tag("sample_10")
xuti::GetNodeData
void GetNodeData(xercesc::DOMNode *node, T &value)
get the node data
Definition: DOMHelperFunctions.h:38
Generated for CMSSW Reference Manual by doxygen 1.8.16