Skip to content

Commit a7d4f4c

Browse files
bear-is-asleepbear-is-asleep
authored
Merge pull request #490 from SBNSoftware/feature/trj_ntb_decoder
NTB Decoder and data product
2 parents fd10e7e + f193775 commit a7d4f4c

7 files changed

Lines changed: 242 additions & 0 deletions

File tree

‎sbndcode/Decoders/TPC/CMakeLists.txt‎

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -19,6 +19,19 @@ cet_build_plugin( SBNDTPCDecoder art::module
1919
sbndcode_ChannelMaps_TPC_TPCChannelMapService_service
2020
)
2121

22+
cet_build_plugin( NTBDecoder art::module
23+
SOURCE NTBDecoder_module.cc
24+
LIBRARIES
25+
sbndaq_artdaq_core::sbndaq-artdaq-core_Overlays_SBND
26+
artdaq_core::artdaq-core_Utilities
27+
art::Utilities
28+
fhiclcpp::fhiclcpp
29+
messagefacility::MF_MessageLogger
30+
art::Framework_Core
31+
ROOT::Core
32+
ROOT::Tree
33+
)
34+
2235
install_headers()
2336
install_fhicl()
2437
install_source()

‎sbndcode/Decoders/TPC/NTBDecoderDefaults.fcl‎

Lines changed: 12 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,12 @@
1+
BEGIN_PROLOG
2+
3+
NTBDecoderDefaults: {
4+
module_type: NTBDecoder
5+
ntb_data_label: "daq:NEVISTB"
6+
ntb_container_data_label: "daq:ContainerNEVISTB"
7+
OutputInstance: ""
8+
DebugLevel: 0
9+
}
10+
11+
END_PROLOG
12+

‎sbndcode/Decoders/TPC/NTBDecoder_module.cc‎

Lines changed: 142 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,142 @@
1+
////////////////////////////////////////////////////////////////////////
2+
// Class: NTBDecoder
3+
// Plugin Type: producer (Unknown Unknown)
4+
// File: NTBDecoder_module.cc
5+
//
6+
// Generated at Thu Aug 1 11:38:21 2024 by Thomas Junk using cetskelgen
7+
// from cetlib version 3.18.02.
8+
////////////////////////////////////////////////////////////////////////
9+
10+
#include "art/Framework/Core/EDProducer.h"
11+
#include "art/Framework/Core/ModuleMacros.h"
12+
#include "art/Framework/Principal/Event.h"
13+
#include "art/Framework/Principal/Handle.h"
14+
#include "art/Framework/Principal/Run.h"
15+
#include "art/Framework/Principal/SubRun.h"
16+
#include "canvas/Utilities/InputTag.h"
17+
#include "fhiclcpp/ParameterSet.h"
18+
#include "messagefacility/MessageLogger/MessageLogger.h"
19+
20+
#include <memory>
21+
22+
#include "artdaq-core/Data/ContainerFragment.hh"
23+
#include "sbndntb.h"
24+
#include "sbndaq-artdaq-core/Overlays/SBND/NevisTPCFragment.hh"
25+
#include "sbndaq-artdaq-core/Overlays/SBND/NevisTBFragment.hh"
26+
#include "sbndaq-artdaq-core/Overlays/SBND/NevisTB_dataFormat.hh"
27+
28+
class NTBDecoder;
29+
30+
31+
class NTBDecoder : public art::EDProducer {
32+
public:
33+
explicit NTBDecoder(fhicl::ParameterSet const& p);
34+
// The compiler-generated destructor is fine for non-base
35+
// classes without bare pointers or other resource use.
36+
37+
// Plugins should not be copied or assigned.
38+
NTBDecoder(NTBDecoder const&) = delete;
39+
NTBDecoder(NTBDecoder&&) = delete;
40+
NTBDecoder& operator=(NTBDecoder const&) = delete;
41+
NTBDecoder& operator=(NTBDecoder&&) = delete;
42+
43+
// Required functions.
44+
void produce(art::Event& e) override;
45+
46+
private:
47+
48+
std::string fTag;
49+
std::string fContainerTag;
50+
std::string fOutputInstance;
51+
int fDebugLevel;
52+
53+
void process_NTB_AUX(const artdaq::Fragment& frag, std::vector<raw::ntb::sbndntb> &sbndntbs);
54+
};
55+
56+
57+
NTBDecoder::NTBDecoder(fhicl::ParameterSet const& p)
58+
: EDProducer{p}
59+
{
60+
fTag = p.get<std::string>("ntb_data_label","daq:NEVISTB");
61+
fContainerTag = p.get<std::string>("ntb_container_data_label","daq:ContainerNEVISTB");
62+
fOutputInstance = p.get<std::string>("OutputInstance","");
63+
fDebugLevel = p.get<int>("DebugLevel",0);
64+
65+
consumes<artdaq::Fragments>(fTag);
66+
consumes<artdaq::Fragments>(fContainerTag);
67+
produces<std::vector<raw::ntb::sbndntb>>(fOutputInstance);
68+
}
69+
70+
void NTBDecoder::produce(art::Event& e)
71+
{
72+
// Implementation of required member function here.
73+
74+
// look first for container fragments and then non-container fragments
75+
76+
std::vector<raw::ntb::sbndntb> sbndntbs;
77+
78+
art::InputTag itag1(fContainerTag);
79+
auto cont_frags = e.getHandle<artdaq::Fragments>(itag1);
80+
if (cont_frags)
81+
{
82+
for (auto const& cont : *cont_frags)
83+
{
84+
artdaq::ContainerFragment cont_frag(cont);
85+
for (size_t ii = 0; ii < cont_frag.block_count(); ++ii)
86+
{
87+
process_NTB_AUX(*cont_frag[ii], sbndntbs);
88+
}
89+
}
90+
}
91+
92+
art::InputTag itag2(fTag);
93+
auto frags = e.getHandle<artdaq::Fragments>(itag2);
94+
if (frags)
95+
{
96+
for(auto const& frag: *frags)
97+
{
98+
process_NTB_AUX(frag, sbndntbs);
99+
}
100+
}
101+
102+
e.put(std::make_unique<std::vector<raw::ntb::sbndntb>>(std::move(sbndntbs)),fOutputInstance);
103+
}
104+
105+
void NTBDecoder::process_NTB_AUX(const artdaq::Fragment& frag, std::vector<raw::ntb::sbndntb> &sbndntbs)
106+
{
107+
if (fDebugLevel > 0)
108+
{
109+
std::cout << "NTBDecoder_module: got into aux" << std::endl;
110+
}
111+
sbndaq::NevisTBFragment tbfrag(frag);
112+
const auto header = *tbfrag.header();
113+
const auto fdata = *tbfrag.data();
114+
const auto trailer = *tbfrag.trailer();
115+
116+
raw::ntb::sbndntb ostruct;
117+
ostruct.boardreader_timestamp = frag.timestamp();
118+
ostruct.event_number = tbfrag.metadata()->EventNumber();
119+
ostruct.frame_number = tbfrag.metadata()->FrameNumber();
120+
ostruct.sample_number = tbfrag.metadata()->SampleNumber();
121+
ostruct.busy = header.busy;
122+
ostruct.sixteen_mhz_remainder = header.remainder;
123+
ostruct.two_mhz_sample = header.sample;
124+
ostruct.frame = header.getFrame();
125+
ostruct.ntrig = header.getTriggerNumber();
126+
ostruct.pmt_trig_data = fdata.pmt_trig_data;
127+
ostruct.pc = fdata.pc;
128+
ostruct.external = fdata.external;
129+
ostruct.active = fdata.active;
130+
ostruct.gate2 = fdata.gate2;
131+
ostruct.gate1 = fdata.gate1;
132+
ostruct.veto = fdata.veto;
133+
ostruct.calib = fdata.calib;
134+
ostruct.phase = fdata.getPhase();
135+
ostruct.gatefake = fdata.gatefake;
136+
ostruct.beamfake = fdata.beamfake;
137+
ostruct.spare1 = fdata.spare1;
138+
ostruct.trailer = trailer.getTrailerWord();
139+
sbndntbs.push_back(ostruct);
140+
}
141+
142+
DEFINE_ART_MODULE(NTBDecoder)

‎sbndcode/Decoders/TPC/classes.h‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,7 @@
11
#include "canvas/Persistency/Common/Wrapper.h"
22
#include <vector>
33
#include "sbndcode/Decoders/TPC/TPCDecodeAna.h"
4+
#include "sbndcode/Decoders/TPC/sbndntb.h"
45

56

67
namespace {

‎sbndcode/Decoders/TPC/classes_def.xml‎

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -7,4 +7,12 @@
77
<class name="std::vector<tpcAnalysis::TPCDecodeAna>"/>
88
<class name="art::Wrapper<tpcAnalysis::TPCDecodeAna>"/>
99
<class name="art::Wrapper<std::vector<tpcAnalysis::TPCDecodeAna>>"/>
10+
11+
<class name="raw::ntb::sbndntb" ClassVersion="10">
12+
<version ClassVersion="10" checksum="4110489055"/>
13+
</class>
14+
<class name="std::vector<raw::ntb::sbndntb>"/>
15+
<class name="art::Wrapper<raw::ntb::sbndntb>"/>
16+
<class name="art::Wrapper<std::vector<raw::ntb::sbndntb>>"/>
17+
1018
</lcgdict>

‎sbndcode/Decoders/TPC/runntbdecoder.fcl‎

Lines changed: 31 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,31 @@
1+
#include "NTBDecoderDefaults.fcl"
2+
3+
physics:
4+
{
5+
// now also have something produce the digits and headers
6+
producers:
7+
{
8+
daq: @local::NTBDecoderDefaults
9+
}
10+
11+
analyzers: {}
12+
13+
my_producer_modules: [daq]
14+
trigger_paths: [my_producer_modules]
15+
a: [commissioning]
16+
end_paths: [a]
17+
}
18+
19+
source: {}
20+
21+
outputs: {
22+
commissioning: {
23+
module_type: "RootOutput"
24+
fileName: "%ifb_ntbdecode.root"
25+
dataTier: "decoded-raw"
26+
}
27+
}
28+
29+
30+
process_name: NTBDECODER
31+

‎sbndcode/Decoders/TPC/sbndntb.h‎

Lines changed: 35 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,35 @@
1+
#ifndef RAW_NTB_SBNDNTB_HEADER_GUARD
2+
#define RAW_NTB_SBNDNTB_HEADER_GUARD
3+
4+
#include <stdlib.h>
5+
6+
namespace raw {
7+
namespace ntb {
8+
struct sbndntb {
9+
uint64_t boardreader_timestamp; // from artdaq::fragment
10+
uint32_t event_number; // from fragment metadata
11+
uint32_t frame_number; // from fragment metadata
12+
uint16_t sample_number; // from fragment metadata
13+
uint16_t busy; // from NTB header words
14+
uint16_t sixteen_mhz_remainder; // from NTB header words
15+
uint16_t two_mhz_sample; // from NTB header words
16+
uint32_t frame; // from NTB header words
17+
uint32_t ntrig; // from NTB header words
18+
uint16_t pmt_trig_data; // from NTB trigger data words
19+
uint16_t pc; // from NTB trigger data words
20+
uint16_t external; // from NTB trigger data words
21+
uint16_t active; // from NTB trigger data words
22+
uint16_t gate2; // from NTB trigger data words
23+
uint16_t gate1; // from NTB trigger data words
24+
uint16_t veto; // from NTB trigger data words
25+
uint16_t calib; // from NTB trigger data words
26+
uint16_t phase; // from NTB trigger data words
27+
uint16_t gatefake; // from NTB trigger data words
28+
uint16_t beamfake; // from NTB trigger data words
29+
uint16_t spare1; // from NTB trigger data words
30+
uint32_t trailer; // from NTB trailer words
31+
};
32+
}
33+
}
34+
35+
#endif

0 commit comments

Comments
 (0)