1 #ifndef mu2e_artdaq_core_Overlays_ArtFragmentReader_hh
2 #define mu2e_artdaq_core_Overlays_ArtFragmentReader_hh
12 #include "artdaq-core/Data/Fragment.hh"
13 #include "mu2e-artdaq-core/Overlays/ArtFragment.hh"
19 class ArtFragmentReader;
21 std::ostream &operator<<(std::ostream &, ArtFragmentReader
const &);
31 adc_t DBH_ByteCount(adc_t
const *pos);
32 bool DBH_Valid(adc_t
const *pos);
33 adc_t DBH_ROCID(adc_t
const *pos);
34 adc_t DBH_PacketType(adc_t
const *pos);
35 adc_t DBH_PacketCount(adc_t
const *pos);
36 uint64_t DBH_Timestamp(adc_t
const *pos);
37 adc_t DBH_TimestampLow(adc_t
const *pos);
38 adc_t DBH_TimestampMedium(adc_t
const *pos);
39 adc_t DBH_TimestampHigh(adc_t
const *pos);
40 adc_t DBH_Status(adc_t
const *pos);
41 adc_t DBH_FormatVersion(adc_t
const *pos);
42 adc_t DBH_EVBMode(adc_t
const *pos);
43 adc_t DBH_SubsystemID(adc_t
const *pos);
44 adc_t DBH_DTCID(adc_t
const *pos);
47 adc_t DBT_StrawIndex(adc_t
const *pos);
48 adc_t DBT_TDC0(adc_t
const *pos);
49 adc_t DBT_TDC1(adc_t
const *pos);
50 adc_t DBT_TOT0(adc_t
const *pos);
51 adc_t DBT_TOT1(adc_t
const *pos);
52 std::array<adc_t, 15> DBT_Waveform(adc_t
const *pos);
53 adc_t DBT_Flags(adc_t
const *pos);
56 adc_t DBCR_NumHits(adc_t
const *pos);
58 adc_t DBC_ChannelID(adc_t
const *pos,
size_t hitIdx);
59 adc_t DBC_BoardID(adc_t
const *pos,
size_t hitIdx);
60 adc_t DBC_Time(adc_t
const *pos,
size_t hitIdx);
61 adc_t DBC_NumSamples(adc_t
const *pos,
size_t hitIdx);
62 adc_t DBC_PeakSampleIdx(adc_t
const *pos,
size_t hitIdx);
63 std::vector<int> DBC_Waveform(adc_t
const *pos,
size_t hitIdx);
65 adc_t DBC_DIRACOutputA(adc_t
const *pos,
size_t hitIdx);
66 adc_t DBC_DIRACOutputB(adc_t
const *pos,
size_t hitIdx);
67 adc_t DBC_ErrorFlags(adc_t
const *pos,
size_t hitIdx);
70 adc_t DBVR_ControllerID(adc_t
const *pos);
71 adc_t DBVR_PacketType(adc_t
const *pos);
72 adc_t DBVR_EventWordCount(adc_t
const *pos);
73 adc_t DBVR_ActiveFEBFlags0(adc_t
const *pos);
74 adc_t DBVR_ActiveFEBFlags1(adc_t
const *pos);
75 adc_t DBVR_ActiveFEBFlags2(adc_t
const *pos);
76 adc_t DBVR_TriggerCount(adc_t
const *pos);
77 adc_t DBVR_EventType(adc_t
const *pos);
78 adc_t DBVR_ErrorFlags(adc_t
const *pos);
80 adc_t DBVR_NumHits(adc_t
const *pos);
83 adc_t DBV_sipmID(adc_t
const *pos,
size_t hitIdx);
85 std::array<unsigned int, 8> DBV_ADCs(adc_t
const *pos,
size_t hitIdx);
86 adc_t DBV_startTDC(adc_t
const *pos,
size_t hitIdx);
93 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_ByteCount(adc_t
const *pos)
98 bool mu2e::ArtFragmentReader::DBH_Valid(adc_t
const *pos)
100 return (*(pos + 1) >> 15) & 0x0001;
103 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_ROCID(adc_t
const *pos)
105 return *(pos + 1) & 0x000F;
108 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_PacketType(adc_t
const *pos)
110 return (*(pos + 1) >> 4) & 0x000F;
113 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_PacketCount(adc_t
const *pos)
115 return *(pos + 2) & 0x07FF;
118 uint64_t mu2e::ArtFragmentReader::DBH_Timestamp(adc_t
const *pos)
120 return uint64_t(*(pos + 3)) + (uint64_t(*(pos + 4)) << 16) + (uint64_t(*(pos + 5)) << 32);
123 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_TimestampLow(adc_t
const *pos)
128 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_TimestampMedium(adc_t
const *pos)
133 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_TimestampHigh(adc_t
const *pos)
138 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_Status(adc_t
const *pos)
140 return *(pos + 6) & 0x00FF;
143 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_FormatVersion(adc_t
const *pos)
145 return *(pos + 6) >> 8;
148 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_EVBMode(adc_t
const *pos)
150 return *(pos + 7) >> 8;
153 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_SubsystemID(adc_t
const *pos)
155 return (*(pos + 7) >> 6) & 0x0003;
158 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBH_DTCID(adc_t
const *pos)
160 return *(pos + 7) & 0x003F;
166 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBT_StrawIndex(adc_t
const *pos)
168 return *(pos + 8 + 0);
171 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBT_TDC0(adc_t
const *pos)
173 return (uint32_t(*(pos + 8 + 1)) & 0xFFFF);
176 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBT_TDC1(adc_t
const *pos)
178 return (uint32_t(*(pos + 8 + 2)) & 0xFFFF);
181 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBT_TOT0(adc_t
const *pos)
183 return (uint32_t(*(pos + 8 + 3)) & 0x00FF);
186 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBT_TOT1(adc_t
const *pos)
188 return ((uint32_t(*(pos + 8 + 3)) >> 8) & 0x00FF);
191 std::array<mu2e::ArtFragmentReader::adc_t, 15> mu2e::ArtFragmentReader::DBT_Waveform(adc_t
const *pos)
193 std::array<adc_t, 15> waveform;
198 for (
size_t i = 0; i < 4; i += 1) {
199 waveform[0 + i * 4] = *(pos + 8 + 4 + i * 3) & 0x0FFF;
200 waveform[1 + i * 4] = ((*(pos + 8 + 4 + i * 3 + 1) & 0x00FF) << 4) | (*(pos + 8 + 4 + i * 3) >> 12);
201 waveform[2 + i * 4] = ((*(pos + 8 + 4 + i * 3 + 2) & 0x000F) << 8) | (*(pos + 8 + 4 + i * 3 + 1) >> 8);
203 waveform[3 + i * 4] = (*(pos + 8 + 4 + i * 3 + 2) >> 4);
210 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBT_Flags(adc_t
const *pos)
212 return (*(pos + 8 + 15) >> 8);
219 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBCR_NumHits(adc_t
const *pos)
221 return *(pos + 8 + 0);
224 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBC_ChannelID(adc_t
const *pos,
size_t hitIdx)
226 return *(pos + 8 + *(pos + 8 + 1 + hitIdx) + 0) >> 10;
229 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBC_BoardID(adc_t
const *pos,
size_t hitIdx)
231 return *(pos + 8 + *(pos + 8 + 1 + hitIdx) + 0) & 0x03FF;
235 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBC_Time(adc_t
const *pos,
size_t hitIdx)
237 return *(pos + 8 + *(pos + 8 + 1 + hitIdx) + 1);
240 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBC_NumSamples(adc_t
const *pos,
size_t hitIdx)
242 return *(pos + 8 + *(pos + 8 + 1 + hitIdx) + 2) & 0x00FF;
245 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBC_PeakSampleIdx(adc_t
const *pos,
size_t hitIdx)
247 return *(pos + 8 + *(pos + 8 + 1 + hitIdx) + 2) >> 8;
250 std::vector<int> mu2e::ArtFragmentReader::DBC_Waveform(adc_t
const *pos,
size_t hitIdx)
252 std::vector<int> waveform(DBC_NumSamples(pos,hitIdx));
253 for (
size_t i = 0; i < waveform.size(); i++) {
254 waveform[i] = *(pos + 8 + *(pos + 8 + 1 + hitIdx) + 3 + i);
259 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBC_DIRACOutputA(adc_t
const *pos,
size_t hitIdx)
261 return *(pos + 8 + *(pos + 8 + 1 + hitIdx) + 3 + DBC_NumSamples(pos,hitIdx) + 0);
264 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBC_DIRACOutputB(adc_t
const *pos,
size_t hitIdx)
266 return *(pos + 8 + *(pos + 8 + 1 + hitIdx) + 3 + DBC_NumSamples(pos,hitIdx) + 1);
269 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBC_ErrorFlags(adc_t
const *pos,
size_t hitIdx)
271 return *(pos + 8 + *(pos + 8 + 1 + hitIdx) + 3 + DBC_NumSamples(pos,hitIdx) + 2);
278 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_ControllerID(adc_t
const *pos)
280 return *(pos + 8 + 0) >> 8;
283 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_PacketType(adc_t
const *pos)
285 return *(pos + 8 + 0) & 0x00FF;
288 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_EventWordCount(adc_t
const *pos)
290 return *(pos + 8 + 1);
293 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_ActiveFEBFlags0(adc_t
const *pos)
295 return *(pos + 8 + 3) & 0x00FF;
298 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_ActiveFEBFlags1(adc_t
const *pos)
300 return *(pos + 8 + 3) >> 8;
303 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_ActiveFEBFlags2(adc_t
const *pos)
305 return *(pos + 8 + 2) & 0x00FF;
308 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_TriggerCount(adc_t
const *pos)
310 return *(pos + 8 + 5);
313 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_EventType(adc_t
const *pos)
315 return *(pos + 8 + 7) >> 8;
318 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_ErrorFlags(adc_t
const *pos)
320 return *(pos + 8 + 7) & 0x00FF;
323 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBVR_NumHits(adc_t
const *pos)
326 adc_t numHits = (*(pos + 8 + 1) - 16) / 12;
330 if(*(pos + 8 + 8 + 6*(numHits-1) + 1) == 0) {
331 numHits = numHits - 1;
345 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBV_sipmID(adc_t
const *pos,
size_t hitIdx)
347 return *(pos + 16 + hitIdx * 6 + 0);
355 std::array<unsigned int, 8> mu2e::ArtFragmentReader::DBV_ADCs(adc_t
const *pos,
size_t hitIdx)
357 std::array<unsigned int, 8> ADCs;
358 for (
size_t i = 0; i < 8; i += 2) {
359 ADCs[i] = *(pos + 16 + hitIdx * 6 + 2 + i / 2) & 0x00FF;
360 ADCs[i + 1] = *(pos + 16 + hitIdx * 6 + 2 + i / 2) >> 8;
365 mu2e::ArtFragmentReader::adc_t mu2e::ArtFragmentReader::DBV_startTDC(adc_t
const *pos,
size_t hitIdx)
367 return *(pos + 16 + hitIdx * 6 + 1) & 0x00FF;