|
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
100
101
102
103
104
105
106
107
108
109
110
111
|
#include "CondTools/Ecal/interface/EcalTPGLinPed.h"
#include "CondTools/Ecal/interface/EcalTPGPedestalsHandler.h"
#include "OnlineDB/EcalCondDB/interface/EcalLogicID.h"
#include "OnlineDB/EcalCondDB/interface/LMFSextuple.h"
#include "FWCore/ParameterSet/interface/ParameterSetfwd.h"
#include "DataFormats/Provenance/interface/Timestamp.h"
#include "OnlineDB/EcalCondDB/interface/Tm.h"
#include <iostream>
#include <iomanip>
#include <sstream>
popcon::EcalTPGLinPed::EcalTPGLinPed(const edm::ParameterSet& ps)
: m_name(ps.getUntrackedParameter<std::string>("name", "EcalTPGLinPedHandler")) {
std::cout << "EcalTPGLinPed constructor\n" << std::endl;
}
popcon::EcalTPGLinPed::~EcalTPGLinPed() {
// do nothing
}
void popcon::EcalTPGLinPed::getNewObjects() {
std::cout << "------- Ecal -> getNewObjects\n";
unsigned long long max_since = 1;
Ref payload = lastPayload();
// here popcon tells us which is the last since of the last object in the offline DB
max_since = tagInfo().lastInterval.since;
Tm max_since_tm(max_since);
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;
std::string line;
for (int i = 0; i < 85; i++)
getline(fLin, line);
char cryst[10];
uint32_t ped[kGains], mult[kGains], shift[kGains];
uint32_t id;
EcalTPGLinearizationConst* linC = new EcalTPGLinearizationConst;
EcalTPGPedestals* peds = new EcalTPGPedestals;
for (int iChannel = 0; iChannel < kEBChannels; iChannel++) {
getline(fLin, line);
sscanf(line.c_str(), "%s %u", cryst, &id);
for (int gain = 0; gain < kGains; gain++) {
getline(fLin, line);
sscanf(line.c_str(), "%X %X %X", &ped[gain], &mult[gain], &shift[gain]);
}
EcalTPGLinearizationConst::Item item;
item.mult_x1 = mult[2];
item.mult_x6 = mult[1];
item.mult_x12 = mult[0];
item.shift_x1 = shift[2];
item.shift_x6 = shift[1];
item.shift_x12 = shift[0];
EcalTPGPedestals::Item itemPed;
itemPed.mean_x1 = ped[2];
itemPed.mean_x6 = ped[1];
itemPed.mean_x12 = ped[0];
linC->insert(std::make_pair(id, item));
peds->insert(std::make_pair(id, itemPed));
} // end loop over EB channels
getline(fLin, line); // comment before EE crystals
std::cout << " comment line " << line << std::endl;
for (int iChannel = 0; iChannel < kEEChannels; iChannel++) {
getline(fLin, line);
// std::cout << " line " << line << std::endl;
sscanf(line.c_str(), "%s %u", cryst, &id);
// std::cout << cryst << " id " << id << std::endl;
for (int gain = 0; gain < kGains; gain++) {
getline(fLin, line);
// std::cout << " line g " << line << std::endl;
sscanf(line.c_str(), "%X %X %X", &ped[gain], &mult[gain], &shift[gain]);
// std::cout << " gain " << gain << " ped " << ped[gain] << " mult " << mult[gain] << " shift " << shift[gain]<< std::endl;
}
EcalTPGLinearizationConst::Item item;
item.mult_x1 = mult[2];
item.mult_x6 = mult[1];
item.mult_x12 = mult[0];
item.shift_x1 = shift[2];
item.shift_x6 = shift[1];
item.shift_x12 = shift[0];
EcalTPGPedestals::Item itemPed;
itemPed.mean_x1 = ped[2];
itemPed.mean_x6 = ped[1];
itemPed.mean_x12 = ped[0];
linC->insert(std::make_pair(id, item));
peds->insert(std::make_pair(id, itemPed));
} // end loop over EE channels
fLin.close();
// for the time beeing just transfer pedestal
// m_to_transfer.push_back(std::make_pair(linC, fileIOV));
m_to_transfer.push_back(std::make_pair(peds, fileIOV));
std::cout << "Ecal -> end of getNewObjects -----------\n";
}
|