|
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
|
#include "CondTools/Ecal/interface/EcalTPGFineGrainTowerfromFile.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSetfwd.h"
#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"
#include "Geometry/EcalMapping/interface/EcalMappingRcd.h"
#include <iostream>
#include <fstream>
#include <sstream>
popcon::EcalTPGFineGrainTowerfromFile::EcalTPGFineGrainTowerfromFile(const edm::ParameterSet& ps)
: m_name(ps.getUntrackedParameter<std::string>("name", "EcalTPGFineGrainTowerfromFile")) {
fname = ps.getParameter<std::string>("FileName");
std::cout << "EcalTPGFineGrainTowerfromFile constructor\n" << std::endl;
}
popcon::EcalTPGFineGrainTowerfromFile::~EcalTPGFineGrainTowerfromFile() {}
void popcon::EcalTPGFineGrainTowerfromFile::getNewObjects() {
std::cout << "------- Ecal -> getNewObjects\n";
edm::LogInfo("EcalTPGFineGrainTowerfromFile") << "Started GetNewObjects!!!";
int fileIOV;
std::cout << "LinPed which input IOV do you want " << std::endl;
std::cin >> fileIOV;
std::ifstream fLin;
std::ostringstream oss;
oss << fileIOV;
// std::string fname = "/afs/cern.ch/cms/ECAL/triggerTransp/TPG_beamv6_trans_" + oss.str() + "_spikekill.txt";
fLin.open(fname.c_str());
if (!fLin.is_open()) {
std::cout << "ERROR : can't open file '" << fname << std::endl;
return;
}
std::cout << " file " << fname << " opened" << std::endl;
/* structure of the file:
TOWER_EB 365224/375015 : 9792 lines 2448 towers 3 lines : 0, 0, 96 (LUTGroupId, FgGroupId, spike_killing_threshold)
empty line 375016
TOWER_EE 375017/379768 : 4752 lines 1584 towers 2 lines : 0, 0x0 (LUTGroupId, tower_lut_fg)
*/
std::string line;
for (int i = 0; i < 375016; i++)
getline(fLin, line);
char tower[8];
unsigned int towerId, LUTFg;
EcalTPGFineGrainTowerEE* fgrMap = new EcalTPGFineGrainTowerEE;
for (int itower = 0; itower < 1584; itower++) {
getline(fLin, line);
sscanf(line.c_str(), "%s %u", tower, &towerId);
if (itower < 10 || (itower > 1574 && itower < 1584))
std::cout << " string " << tower << " Id " << towerId;
getline(fLin, line); // LUTGroupId
getline(fLin, line); // tower_lut_fg
if (itower < 10 || (itower > 1574 && itower < 1584))
std::cout << " line " << line << std::endl;
sscanf(line.c_str(), "%x", &LUTFg);
// EcalTPGFineGrainTowerEE::Item item;
// item.lut = LUTFg;
// fgrMap->setValue(towerId, item);
fgrMap->setValue(towerId, LUTFg);
} // end loop over EE towers
fLin.close();
m_to_transfer.push_back(std::make_pair(fgrMap, fileIOV));
std::cout << "Ecal -> end of getNewObjects -----------\n";
}
|