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
|
#include <iostream>
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "CondTools/Ecal/interface/EcalIntercalibConstantsXMLTranslator.h"
#include "CondTools/Ecal/interface/EcalCondHeader.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include <string>
#include <sstream>
#include <fstream>
using namespace std;
void usage() {
cout << endl;
cout << "array2xml [arrayfile] [xmlfile]" << endl;
cout << "Read coefficients from straight array [denseindex]" << endl;
cout << "and write in xml format" << endl;
}
int main(int argc, char* argv[]) {
if (argc != 3) {
usage();
exit(0);
}
string arrayfilename(argv[1]);
string xmlfilename(argv[2]);
fstream arrayfile(arrayfilename.c_str(), ios::in);
EcalIntercalibConstants rcd;
float c = 0;
int idx = 0;
while (arrayfile >> c) {
uint32_t id = EBDetId::unhashIndex(idx);
rcd[id] = c;
++idx;
}
cout << idx << endl;
for (int cellid = 0; cellid < EEDetId::kSizeForDenseIndexing; ++cellid) { // loop on EB cells
if (EEDetId::validHashIndex(cellid)) {
uint32_t rawid = EEDetId::unhashIndex(cellid);
rcd[rawid] = 0.0;
} // if
}
// write new format
EcalCondHeader h;
EcalIntercalibConstantsXMLTranslator::writeXML(xmlfilename, h, rcd);
}
|