|
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
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
|
#include <memory>
// user include files
#include "EventFilter/EcalDigiToRaw/interface/BlockFormatter.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "DataFormats/EcalDetId/interface/EcalElectronicsId.h"
using namespace std;
BlockFormatter::BlockFormatter(Config const& iC, Params const& iP)
: plistDCCId_{iC.plistDCCId_},
orbit_number_{iP.orbit_number_},
bx_{iP.bx_},
lv1_{iP.lv1_},
runnumber_{iP.runnumber_},
debug_{iC.debug_},
doBarrel_{iC.doBarrel_},
doEndCap_{iC.doEndCap_},
doTCC_{iC.doTCC_},
doSR_{iC.doSR_},
doTower_{iC.doTower_} {}
void BlockFormatter::DigiToRaw(FEDRawDataCollection* productRawData) {
auto const run_number = runnumber_;
auto const bx = bx_;
auto const lv1 = lv1_;
if (debug_)
cout << "in BlockFormatter::DigiToRaw run_number orbit_number bx lv1 " << dec << run_number << " " << orbit_number_
<< " " << bx << " " << lv1 << endl;
for (int idcc = 1; idcc <= 54; idcc++) {
if ((!doBarrel_) && (idcc >= EcalElectronicsId::MIN_DCCID_EBM && idcc <= EcalElectronicsId::MAX_DCCID_EBP))
continue;
if ((!doEndCap_) && (idcc <= EcalElectronicsId::MAX_DCCID_EEM || idcc >= EcalElectronicsId::MIN_DCCID_EEP))
continue;
int FEDid = FEDNumbering::MINECALFEDID + idcc;
FEDRawData& rawdata = productRawData->FEDData(FEDid);
unsigned char* pData;
short int DCC_ERRORS = 0;
if (rawdata.size() == 0) {
rawdata.resize(8);
pData = rawdata.data();
Word64 word = 0x18 + ((FEDid & 0xFFF) << 8) + ((Word64)((Word64)bx & 0xFFF) << 20) +
((Word64)((Word64)lv1 & 0xFFFFFF) << 32) + (Word64)((Word64)0x51 << 56);
Word64* pw = reinterpret_cast<Word64*>(const_cast<unsigned char*>(pData));
*pw = word; // DAQ header
rawdata.resize(rawdata.size() + 8 * 8); // DCC header
pData = rawdata.data();
pData[11] = DCC_ERRORS & 0xFF;
pData[12] = run_number & 0xFF;
pData[13] = (run_number >> 8) & 0xFF;
pData[14] = (run_number >> 16) & 0xFF;
pData[15] = 0x01;
for (int i = 16; i <= 22; i++) {
pData[i] = 0; // to be filled for local data taking or calibration
}
pData[23] = 0x02;
pData[24] = orbit_number_ & 0xFF;
pData[25] = (orbit_number_ >> 8) & 0xFF;
pData[26] = (orbit_number_ >> 16) & 0xFF;
pData[27] = (orbit_number_ >> 24) & 0xFF;
int SRenable_ = 1;
int SR = SRenable_;
int ZS = 0;
int TZS = 0;
// int SR_CHSTATUS = 0;
pData[28] = (SR & 0x1) + ((ZS & 0x1) << 1) + ((TZS & 0x1) << 2);
pData[31] = 0x03;
for (int i = 0; i <= 4; i++) {
for (int j = 0; j < 7; j++) {
pData[32 + 8 * i + j] = 0;
}
pData[32 + 8 * i + 7] = 0x04;
}
} // endif rawdatasize == 0
} // loop on id
}
void BlockFormatter::print(FEDRawData& rawdata) {
int size = rawdata.size();
cout << "Print RawData size " << dec << size << endl;
unsigned char* pData = rawdata.data();
int n = size / 8;
for (int i = 0; i < n; i++) {
for (int j = 7; j >= 0; j--) {
if (8 * i + j <= size)
cout << hex << (int)pData[8 * i + j] << " ";
}
cout << endl;
}
}
void BlockFormatter::CleanUp(FEDRawDataCollection* productRawData, map<int, map<int, int> >* FEDorder) {
for (int id = 0; id < 36 + 18; id++) {
if ((!doBarrel_) && (id >= 9 && id <= 44))
continue;
if ((!doEndCap_) && (id <= 8 || id >= 45))
continue;
int FEDid = FEDNumbering::MINECALFEDID + id + 1;
FEDRawData& rawdata = productRawData->FEDData(FEDid);
// ---- if raw need not be made for a given fed, set its size to empty and return
if (find((*plistDCCId_).begin(), (*plistDCCId_).end(), (id + 1)) == (*plistDCCId_).end()) {
rawdata.resize(0);
continue;
}
// ---- Add the trailer word
int lastline = rawdata.size();
rawdata.resize(lastline + 8);
unsigned char* pData = rawdata.data();
int event_length = (lastline + 8) / 8; // in 64 bits words
pData[lastline + 7] = 0xa0;
// pData[lastline+4] = event_length & 0xFFFFFF;
pData[lastline + 4] = event_length & 0xFF;
pData[lastline + 5] = (event_length >> 8) & 0xFF;
pData[lastline + 6] = (event_length >> 16) & 0xFF;
int event_status = 0;
pData[lastline + 1] = event_status & 0x0F;
int tts = 0 << 4;
pData[lastline] = tts & 0xF0;
// ---- Write the event length in the DCC header
// pData[8] = event_length & 0xFFFFFF;
pData[8] = event_length & 0xFF;
pData[9] = (event_length >> 8) & 0xFF;
pData[10] = (event_length >> 16) & 0xFF;
// cout << " in BlockFormatter::CleanUp. FEDid = " << FEDid << " event_length*8 " << dec << event_length*8 << endl;
map<int, map<int, int> >::iterator fen = FEDorder->find(FEDid);
bool FED_has_data = true;
if (fen == FEDorder->end())
FED_has_data = false;
if (debug_ && (!FED_has_data))
cout << " FEDid is not in FEDorder ! " << endl;
if (!FED_has_data) {
int ch_status = 7;
for (int iFE = 1; iFE <= 68; iFE++) {
int irow = (iFE - 1) / 14;
int kval = ((iFE - 1) % 14) / 2;
if (iFE % 2 == 1)
pData[32 + 8 * irow + kval] |= ch_status & 0xFF;
else
pData[32 + 8 * irow + kval] |= ((ch_status << 4) & 0xFF);
}
}
if (FED_has_data) {
map<int, int>& FEorder = (*fen).second;
for (int iFE = 1; iFE <= 68; iFE++) {
map<int, int>::iterator fe = FEorder.find(iFE);
int ch_status = 0;
if (fe == FEorder.end()) // FE not present due to SRP, update CH_status
ch_status = 7; // CH_SUPPRESS
int irow = (iFE - 1) / 14;
int kval = ((iFE - 1) % 14) / 2;
if (iFE % 2 == 1)
pData[32 + 8 * irow + kval] |= ch_status & 0xFF;
else
pData[32 + 8 * irow + kval] |= ((ch_status << 4) & 0xFF);
}
}
}
}
void BlockFormatter::PrintSizes(FEDRawDataCollection* productRawData) {
for (int id = 0; id < 36 + 18; id++) {
// if ( (! doBarrel_) && (id >= 9 && id <= 44)) continue;
// if ( (! doEndCap_) && (id <= 8 || id >= 45)) continue;
int FEDid = FEDNumbering::MINECALFEDID + id;
FEDRawData& rawdata = productRawData->FEDData(FEDid);
if (rawdata.size() > 0)
cout << "Size of FED id " << dec << FEDid << " is : " << dec << rawdata.size() << endl;
}
}
|