|
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
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
|
#include "CondTools/Ecal/interface/EcalLaser_weekly_Linearization.h"
#include "OnlineDB/EcalCondDB/interface/EcalLogicID.h"
#include "OnlineDB/EcalCondDB/interface/LMFSextuple.h"
#include "FWCore/ParameterSet/interface/ParameterSetfwd.h"
#include "CondCore/CondDB/interface/Time.h"
#include "DataFormats/Provenance/interface/Timestamp.h"
#include "OnlineDB/EcalCondDB/interface/Tm.h"
#include <iostream>
#include <iomanip>
#include <sstream>
popcon::EcalLaser_weekly_Linearization::EcalLaser_weekly_Linearization(const edm::ParameterSet &ps)
: m_name(ps.getUntrackedParameter<std::string>("name", "EcalLaser_weekly_Handler")) {
std::cout << "EcalLaser Source handler constructor\n" << std::endl;
}
popcon::EcalLaser_weekly_Linearization::~EcalLaser_weekly_Linearization() {
// do nothing
}
void popcon::EcalLaser_weekly_Linearization::getNewObjects() {
// int file[1] = {190708};
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);
// const EcalLaserAPDPNRatios::EcalLaserAPDPNRatiosMap& laserRatiosMap =
// payload->getLaserMap();
// std::cout << "payload->getLaserMap(): OK " << std::endl;
// std::cout << "Its size is " << laserRatiosMap.size() << std::endl;
// const EcalLaserAPDPNRatios::EcalLaserTimeStampMap& laserTimeMap =
// payload->getTimeMap();
// std::cout << "payload->getTimeMap(): OK " << std::endl;
// std::cout << "Last Object in Offline DB has SINCE = " << max_since
// << " -> " << max_since_tm.cmsNanoSeconds()
// << " (" << max_since_tm << ")"
// << " and SIZE = " << tagInfo().size
// << std::endl;
for (int week = 0; week < 1; week++) {
int fileIOV;
std::cout << " which input IOV do you want " << std::endl;
std::cin >> fileIOV;
std::ifstream fWeek;
std::ostringstream oss;
oss << fileIOV;
std::string fname = "/afs/cern.ch/cms/ECAL/triggerTransp/TPG_beamv6_trans_" + oss.str() + "_spikekill.txt";
/*
oss << file[week];
// std::string fname = "/afs/cern.ch/cms/ECAL/triggerTransp/TPG_beamv5_trans_" + oss.str() + "_";
std::string fname = "/afs/cern.ch/cms/ECAL/triggerTransp/TPG_beamv6_trans_" + oss.str() + "_";
oss.str("");
// if(week == 32) oss << 200000;
// else oss << file[week + 1] - 1;
oss << 200000;
fname += oss.str() + ".txt";
*/
fWeek.open(fname.c_str());
if (!fWeek.is_open()) {
std::cout << "ERROR : can't open file '" << fname << std::endl;
break;
}
std::cout << " file " << fname << " opened" << std::endl;
// int rawId;
// float corrp;
std::string line;
for (int i = 0; i < 85; i++)
getline(fWeek, line);
char cryst[10];
uint32_t ped[kGains], mult[kGains], shift[kGains];
uint32_t id;
EcalTPGLinearizationConst *linC = new EcalTPGLinearizationConst;
for (int iChannel = 0; iChannel < kEBChannels; iChannel++) {
getline(fWeek, line);
// std::cout << " line " << line << std::endl;
// fWeek >> cryst >> id;
sscanf(line.c_str(), "%s %u", cryst, &id);
// std::cout << cryst << " id " << id << std::endl;
// EBDetId ebId = DetId(id);
for (int gain = 0; gain < kGains; gain++) {
// fWeek >> std::hex >> ped[gain] >> mult[gain] >> shift[gain];
getline(fWeek, 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];
linC->insert(std::make_pair(id, item));
// corr.p1=corr.p2=corr.p3 = corrp;
// corrSet->setValue((int)ebId, corr );
// cryst ="";
} // end loop over EB channels
getline(fWeek, line); // cmment before EE crystals
std::cout << " comment line " << line << std::endl;
for (int iChannel = 0; iChannel < kEEChannels; iChannel++) {
getline(fWeek, line);
// std::cout << " line " << line << std::endl;
sscanf(line.c_str(), "%s %u", cryst, &id);
// std::cout << cryst << " id " << id << std::endl;
// EEDetId eeId = DetId(id);
for (int gain = 0; gain < kGains; gain++) {
getline(fWeek, 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;
} // corr.p1=corr.p2=corr.p3 = corrp;
// corrSet->setValue((int)eeId, corr );
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];
linC->insert(std::make_pair(id, item));
} // end loop over EE channels
fWeek.close();
// special tag for Stephanie
// m_to_transfer.push_back(std::make_pair((EcalTPGLinearizationConst*)linC, file[week]));
m_to_transfer.push_back(std::make_pair(linC, fileIOV));
// end special
} // end loop over week
std::cout << "Ecal -> end of getNewObjects -----------\n";
}
|