|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
#include <iostream>
#include <sstream>
#include <fstream>
#include "CondFormats/EcalObjects/interface/EcalTPGLinearizationConst.h"
#include "CondTools/Ecal/interface/EcalTPGLinearizationConstXMLTranslator.h"
#include "CondTools/Ecal/interface/DOMHelperFunctions.h"
#include "CondTools/Ecal/interface/XMLTags.h"
#include "Utilities/Xerces/interface/Xerces.h"
#include "Utilities/Xerces/interface/XercesStrUtils.h"
using namespace XERCES_CPP_NAMESPACE;
using namespace xuti;
using namespace std;
int EcalTPGLinearizationConstXMLTranslator::writeXML(const std::string& filename,
const EcalCondHeader& header,
const EcalTPGLinearizationConst& record) {
cms::concurrency::xercesInitialize();
std::fstream fs(filename.c_str(), ios::out);
fs << dumpXML(header, record);
cms::concurrency::xercesTerminate();
return 0;
}
std::string EcalTPGLinearizationConstXMLTranslator::dumpXML(const EcalCondHeader& header,
const EcalTPGLinearizationConst& record) {
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(Linearization_tag.c_str()).ptr(), doctype);
DOMElement* root = doc->getDocumentElement();
xuti::writeHeader(root, header);
// open also a flat text file
std::ofstream fout;
fout.open("Linearization.txt");
for (int cellid = EBDetId::MIN_HASH; cellid < EBDetId::kSizeForDenseIndexing; ++cellid) {
uint32_t rawid = EBDetId::unhashIndex(cellid);
DOMElement* cellnode = writeCell(root, rawid);
float m12 = (record)[rawid].mult_x12;
float s12 = (record)[rawid].shift_x12;
float m6 = (record)[rawid].mult_x6;
float s6 = (record)[rawid].shift_x6;
float m1 = (record)[rawid].mult_x1;
float s1 = (record)[rawid].shift_x1;
WriteNodeWithValue(cellnode, Linearization_m12_tag, m12);
WriteNodeWithValue(cellnode, Linearization_m6_tag, m6);
WriteNodeWithValue(cellnode, Linearization_m1_tag, m1);
WriteNodeWithValue(cellnode, Linearization_s12_tag, s12);
WriteNodeWithValue(cellnode, Linearization_s6_tag, s6);
WriteNodeWithValue(cellnode, Linearization_s1_tag, s1);
fout << rawid << " " << m12 << " " << m6 << " " << m1 << " " << s12 << " " << s6 << " " << s1 << "\n";
}
for (int cellid = 0; cellid < EEDetId::kSizeForDenseIndexing; ++cellid) {
if (!EEDetId::validHashIndex(cellid))
continue;
uint32_t rawid = EEDetId::unhashIndex(cellid);
DOMElement* cellnode = writeCell(root, rawid);
float m12 = (record)[rawid].mult_x12;
float s12 = (record)[rawid].shift_x12;
float m6 = (record)[rawid].mult_x6;
float s6 = (record)[rawid].shift_x6;
float m1 = (record)[rawid].mult_x1;
float s1 = (record)[rawid].shift_x1;
WriteNodeWithValue(cellnode, Linearization_m12_tag, m12);
WriteNodeWithValue(cellnode, Linearization_m6_tag, m6);
WriteNodeWithValue(cellnode, Linearization_m1_tag, m1);
WriteNodeWithValue(cellnode, Linearization_s12_tag, s12);
WriteNodeWithValue(cellnode, Linearization_s6_tag, s6);
WriteNodeWithValue(cellnode, Linearization_s1_tag, s1);
fout << rawid << " " << m12 << " " << m6 << " " << m1 << " " << s12 << " " << s6 << " " << s1 << "\n";
}
fout.close();
std::string dump = cms::xerces::toString(writer->writeToString(root));
doc->release();
doctype->release();
writer->release();
return dump;
}
|