79{
81 LArDigit * scDigit=0 ;
82
83
84 uint32_t nData = robFrag->rod_ndata();
86 if (nData < 3) {
88 << "Empty ROD block (less than 3 words) for source ID " << robFrag->rod_source_id() );
89 return;
90 }
91 else
92
93 robFrag->rod_source_id();
94
95
96 unsigned int n = robFrag->payload_size_word();
97
98
99 std::vector<std::string> string;
102
103 for(
unsigned int s=1;
s<
n-13;
s+=6) {
104
106 for (
unsigned int k=0;
k<6;++
k){
108
109 std::ostringstream oss;
110 oss << std::hex <<
x;
111 if (oss.str().size() == 7){
112
113
115
116 }
118
119 }
120
121 string.push_back(std::move(
ss));
122 }
123
124
125
126
127constexpr int nFibers = 20;
128constexpr int nchannels = 8;
129constexpr int Nabbachannels = nFibers * nchannels;
130
131constexpr int fiberSeq[16] = {0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1};
132constexpr int adcSeq[16] = {1,0,1,0,3,2,3,2,5,4,5,4,7,6,7,6};
133
134std::vector< std::vector< std::vector<int> > > abbasamples;
135abbasamples.resize(nFibers);
136
137std::vector<std::vector<short>> abba_channel_samples;
138abba_channel_samples.resize(Nabbachannels);
139
140for (int fib=0;fib < nFibers; ++fib){
141 abbasamples[fib].resize(nchannels);
142 for (
int ch = 0;
ch < nchannels; ++
ch){
143 abbasamples[fib][
ch].resize(nsamples);
144 int Nchannel = fib*8 +
ch;
145 abba_channel_samples[Nchannel].resize(nsamples);
146 }
147 }
148
149
150
151 for (int fib=0;fib < nFibers; fib+=2){
152 for (
int sam=0; sam <
nsamples; sam++){
153 int zx = sam + (fib/2)*nsamples;
154 const std::string &wx = string[zx];
156 for (
int m = 0;
m < 16; ++
m){
157 int fiber_number = fiberSeq[
m];
158 int adc_number = adcSeq[
m];
160 std::stringstream
ss;
161 char mysample[3];
162
163 for (
int j = idx; j <
idx+3; ++j){
165
167 }
169
170 short mysample_int;
171 ss >> std::hex >> mysample_int;
172 int fb = fiber_number+fib;
174 int abbachannel = fb*8+
ch;
176
177
178
179
180
181 abbasamples[fb][
ch][sam] = mysample_int;
182 abba_channel_samples[abbachannel][sam] = mysample_int;
183
184
185 }
186 }
187 }
188
189
190std::map<int, HWIdentifier> abba_mapping;
191
192
193int iphi = 0;
194switch (robFrag->rod_source_id())
195{
196case 0x410511:
197 iphi=19;
198 break;
200
201case 0x410512:
202 iphi=18;
203 break;
204case 0x410521:
205 iphi=20;
206 break;
208case 0x410522:
209 iphi=21;
210 break;
211
212default:
213msg(MSG::ERROR) <<
"ROD_source_id not found" <<
endmsg;
214}
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
281
282
283for (std::size_t index_abba_channel_samples = 0; index_abba_channel_samples < abba_channel_samples.size(); ++index_abba_channel_samples) {
284 auto abba_mapping_search = abba_mapping.find(index_abba_channel_samples);
285 if (abba_mapping_search == abba_mapping.end()) {
286 msg(MSG::WARNING) <<
"No mapping for index_abba_channel_samples: "<<std::dec<<index_abba_channel_samples <<
endmsg;
287 continue;
288 }
289
290 HWIdentifier channelID = (*abba_mapping_search).second;
291
292 std::vector<short> adcValues = abba_channel_samples[index_abba_channel_samples];
293
294 scDigit=new LArDigit(channelID,dummyGain,std::move(adcValues));
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
335
336 }
337
338abbasamples.clear();
339
340 return;
341
342
343}
static void fill(std::map< int, HWIdentifier > *toFill, int iphi)
value_type push_back(value_type pElem)
Add an element to the end of the collection.