Added changes to read 192b words

Author: maxdubnowski

sbndcode/Decoders/CMakeLists.txt

@@ -3,5 +3,6 @@ add_subdirectory(PTB)
 add_subdirectory(SPECTDC)
 add_subdirectory(PMT)
 add_subdirectory(CRT)
+add_subdirectory(PTBAna)
 
 install_fhicl()

sbndcode/Decoders/PTB/SBNDPTBDecoder_module.cc

@@ -76,9 +76,13 @@ SBNDPTBDecoder::SBNDPTBDecoder(fhicl::ParameterSet const & p)
   produces<std::vector<raw::ptb::sbndptb> >(fOutputInstance);
 }
 
+int _run; int _subrun; int _event;
 void SBNDPTBDecoder::produce(art::Event & evt)
 {
-
+  _run = evt.id().run();
+  _subrun = evt.id().subRun();
+  _event = evt.id().event();
+   std::cout << "Run: " << _run  << "  SubRun: " << _subrun  << "  Event: " << _event  << std::endl; 
 
   // look first for container fragments and then non-container fragments
 
@@ -117,6 +121,61 @@ void SBNDPTBDecoder::produce(art::Event & evt)
   evt.put(std::make_unique<std::vector<raw::ptb::sbndptb>>(std::move(sbndptbs)),fOutputInstance);
 }
 
+
+void print_fragment_words(const artdaq::Fragment& frag, size_t wordNum, size_t bitsPerWord ) {
+  
+  const __uint8_t* data_ptr = reinterpret_cast<const __uint8_t*>(frag.dataBegin());
+  size_t wordCount = frag.dataSizeBytes() / (bitsPerWord / 8);
+  for (size_t w = 0; w < wordCount; w++) {
+    // Check if the current word index matches the specified word index
+    if (w == wordNum) {
+        // Print the bits for the specified n-bit word
+      for (int i = (bitsPerWord-1); i >= 0; --i) {
+            // Calculate the byte index and bit index
+	size_t byteIndex = w * (bitsPerWord / 8) + (i / 8);
+	int bitIndex = (i % 8); // Get the correct bit position in the byte
+	// Print the bit
+	std::cout << static_cast<int>((data_ptr[byteIndex] >> bitIndex) & 1);
+	if (i % 8 == 0) {
+	  std::cout << " "; 
+	}
+      }
+      std::cout << std::endl; // New line after printing the word
+    }
+  }      
+}
+
+
+//Copied from Lines 1382-1406 from     sbndaq-artdaq/sbndaq-artdaq/ArtModules/SBND/EventAna_module.cc
+
+ /*******************************************
+  * Note: Below for the Timestamp conversion,
+  *       The PTB TS is in UTC seconds since the Unix epoch in 50MHz clock ticks.
+  *       This means to recover seconds since the Unix epoch use: sec_since_epoch = TS / 50e6
+  *       and of course nanosec_since_epoch = (TS / 50e6) * 1e9 = TS * 20
+  *
+  * Note: The `CTBFragment` constructor grabs the chunk of memory the artDAQ fragment occupies.
+  *       Since we know a priori the number of PTB words (from the TCP header) and the size (128b for all words)
+  *       we can loop over the memory addresses, casting each 128b into the correct word type given by each word's `type`.
+  *
+  * There are 5 word types with the following bit arrangement:
+  *       Feeback        = | 3b Word Type | 61b Payload | 64b Timestamp | (here Payload is split into code, source, payload1, payload2)
+  *       LLT, HLT       = | 3b Word Type | 61b Payload | 64b Timestamp |
+  *       Channel Status = | 3b Word Type | 64b Payload | 61b Timestamp | (Larger Payload to fit all input channels)
+  *
+  * Word Type:
+  *       - 0x0 = Feedback/Error Word -> Errors from the firmware, should abort the run
+  *       - 0x1 = Low Level Trigger (LLT) -> Holds a record of any asserted LLTs
+  *       - 0x2 = High Level Trigger (HLT) -> Holds a record of any asserted HLTs
+  *       - 0x3 = Channel status -> Holds a record of any asserted inputs
+  *       - 0x7 = Timestamp Word -> No payload just a timestamp, these are periodic
+  *
+  * Note: Payload AND'd with mask the size of the expected number of bits just
+  *       to make sure we don't get any unexpected garbage bits since we aren't using
+  *       the full bits of the variable type e.g. uint64_t, unint16_t..
+  **************************************/
+
+
 void SBNDPTBDecoder::_process_PTB_AUX(const artdaq::Fragment& frag, ptbsv_t &sout)
 {
   sbndaq::CTBFragment ctbfrag(frag);   // somehow the name CTBFragment stuck
@@ -145,8 +204,24 @@ void SBNDPTBDecoder::_process_PTB_AUX(const artdaq::Fragment& frag, ptbsv_t &sou
 	  tstruct.timestamp = ctbfrag.Trigger(iword)->timestamp;
 	  if (ctbfrag.Trigger(iword)->IsHLT()) 
 	    {
+	      
+	      if (*(frag.metadata<int>()) == 2){ //If data is taken with 192b PTB words
+		tstruct.prev_timestamp = ctbfrag.PTBWord(iword)->prevTS;
+		tstruct.gate_counter = ctbfrag.Trigger(iword)->gate_counter;
+              }
+              else if (*(frag.metadata<int>()) != 2){ //If data is taken with 128b PTB words
+		tstruct.prev_timestamp = 0;
+		tstruct.gate_counter = 0;
+	      }
+
 	      ix = sout.HLTrigs.size();
 	      sout.HLTrigs.push_back(tstruct);
+
+	      //Make sure HLT words are printing correctly
+	      print_fragment_words(frag, iword, 192 ); 
+	      std::cout <<"Gate Counter: " << ctbfrag.Trigger(iword)->gate_counter  << "  |  Trigger Payload: " << tstruct.trigger_word << "  | Previous Timestamp: " << ctbfrag.PTBWord(iword)->prevTS*20 << std::endl;
+
+
 	      if (fDebugLevel > 0)
 		{
 		  std::cout << "SBNDPTBDecoder_module: found HLT: " << wt << " " << ix << std::endl;

sbndcode/Decoders/PTB/sbndptb.h

@@ -18,12 +18,15 @@ namespace raw {
   namespace ptb {
 
     struct Trigger {
+      ULong64_t prev_timestamp;
       ULong64_t timestamp;
       ULong64_t trigger_word;
+      uint32_t gate_counter;
       uint32_t word_type;
     };
 
     struct ChStatus {
+      ULong64_t prev_timestamp;
       ULong64_t timestamp;
       uint32_t beam;
       uint32_t crt;
@@ -35,6 +38,7 @@ namespace raw {
     };
 
     struct Feedback {
+      ULong64_t prev_timestamp;
       ULong64_t timestamp;
       uint32_t code;
       uint32_t source;
@@ -43,6 +47,7 @@ namespace raw {
     };
 
     struct Misc {
+      ULong64_t prev_timestamp;
       ULong64_t timestamp;
       ULong64_t payload;
       uint32_t word_type;

sbndcode/Decoders/PTBAna/CMakeLists.txt

@@ -0,0 +1,10 @@
+art_make(
+         MODULE_LIBRARIES
+	 ROOT::Tree
+         larsim::Utils
+	 sbndcode_GeoWrappers
+)
+
+
+
+install_fhicl()