Merge pull request #589 from SBNSoftware/feature/rtriozzi_HWTriggerInfoToCAFs_v10_06_00_01

Add ICARUS HW trigger information to CAFs

Author: kjplows

sbncode/CAFMaker/CAFMaker_module.cc

@@ -1703,6 +1703,8 @@ void CAFMaker::produce(art::Event& evt) noexcept {
   caf::SRTrigger srtrigger;
   if (isValidTrigger) {
       FillTrigger(*extratrig_handle, trig_handle->at(0), srtrigger, triggerShift);
+      if (fDet == kICARUS) 
+        FillTriggerICARUS(*extratrig_handle, srtrigger);
   }
   // Fill trigger emulation information
   if (isValidEmulationTrigger) { 

sbncode/CAFMaker/FillTrigger.cxx

@@ -4,20 +4,16 @@
 namespace caf
 {
   void FillTrigger(const sbn::ExtraTriggerInfo& addltrig_info,
-		   const raw::Trigger& trig,
-		   caf::SRTrigger& triggerInfo,
-		   const double time_offset = 0.0)
+                   const raw::Trigger& trig,
+                   caf::SRTrigger& triggerInfo,
+                   const double time_offset = 0.0)
   {
     triggerInfo.global_trigger_time = addltrig_info.triggerTimestamp;
     triggerInfo.beam_gate_time_abs = addltrig_info.beamGateTimestamp;
     triggerInfo.beam_gate_det_time = trig.BeamGateTime() + time_offset;
     triggerInfo.global_trigger_det_time = trig.TriggerTime() + time_offset;
-	std::cout << "Gen triggerInfo.global_trigger_det_time: " << triggerInfo.global_trigger_det_time << std::endl;
-	std::cout << "Gen triggerInfo.beam_gate_det_time: " << triggerInfo.beam_gate_det_time << std::endl;
     double diff_ts = triggerInfo.global_trigger_det_time - triggerInfo.beam_gate_det_time;
-    std::cout << "diff_ts: " << diff_ts << std::endl;
     triggerInfo.trigger_within_gate = diff_ts;
-
     triggerInfo.prev_global_trigger_time = addltrig_info.previousTriggerTimestamp;
     triggerInfo.source_type = sbn::bits::value(addltrig_info.sourceType);
     triggerInfo.trigger_type = sbn::bits::value(addltrig_info.triggerType);
@@ -29,23 +25,34 @@ namespace caf
 
   }
 
-  void FillTriggerMC(double absolute_time, caf::SRTrigger& triggerInfo) {
+  void FillTriggerMC(double absolute_time, caf::SRTrigger& triggerInfo) 
+  {
     triggerInfo.global_trigger_time = absolute_time;
     triggerInfo.beam_gate_time_abs = absolute_time;
+    triggerInfo.trigger_within_gate = 0.; // Set this to 0 since the "MC" trigger is (for now) always at the spill time
+  }
 
-    // Set this to 0 since the "MC" trigger is (for now) always at the spill time
-    triggerInfo.trigger_within_gate = 0.;
+  void FillTriggerICARUS(const sbn::ExtraTriggerInfo& addltrig_info,
+                         caf::SRTrigger& triggerInfo) 
+  {
+    // Choose the cryostat that triggered first (if both are available)
+    int const cryo = addltrig_info.cryostats[sbn::ExtraTriggerInfo::EastCryostat].beamToTrigger < addltrig_info.cryostats[sbn::ExtraTriggerInfo::WestCryostat].beamToTrigger
+      ? sbn::ExtraTriggerInfo::EastCryostat 
+      : sbn::ExtraTriggerInfo::WestCryostat;
 
+    sbn::ExtraTriggerInfo::CryostatInfo const& cryoInfo = addltrig_info.cryostats[cryo];
+    if (cryoInfo.hasTrigger()) {
+      triggerInfo.trigger_cryo_source  = cryo;
+      triggerInfo.trigger_logic_bits   = cryoInfo.triggerLogicBits;
+      triggerInfo.gate_to_trigger_time = static_cast<float>(cryoInfo.beamToTrigger) / 1000.0f; // [us]
+    }  
   }
 
   void FillTriggerSBND(caf::SRSBNDTimingInfo& timingInfo, caf::SRTrigger& triggerInfo){
 
     triggerInfo.global_trigger_time = timingInfo.hltEtrig;
     triggerInfo.beam_gate_time_abs = timingInfo.hltBeamGate;
-    std::cout << "triggerInfo.global_trigger_det_time: " << triggerInfo.global_trigger_det_time << std::endl;
-	std::cout << "triggerInfo.beam_gate_det_time: " << triggerInfo.beam_gate_det_time << std::endl;
     double diff_ts = triggerInfo.global_trigger_det_time - triggerInfo.beam_gate_det_time;
-	std::cout << "diff_ts: " << diff_ts << std::endl;
     triggerInfo.trigger_within_gate = diff_ts;
   }
 

sbncode/CAFMaker/FillTrigger.h

@@ -8,18 +8,18 @@
 #include "lardataobj/RawData/TriggerData.h"
 #include "art/Framework/Principal/Handle.h"
 
-#include <vector>
-
 namespace caf
 {
-
   void FillTrigger(const sbn::ExtraTriggerInfo& addltrig_info,
                    const raw::Trigger& trig_info,
                    caf::SRTrigger& triggerInfo,
                    const double time_offset);
 
   void FillTriggerMC(double absolute_time, caf::SRTrigger& triggerInfo);
 
+  void FillTriggerICARUS(const sbn::ExtraTriggerInfo& addltrig_info,
+                         caf::SRTrigger& triggerInfo);            
+                                        
   void FillTriggerSBND(caf::SRSBNDTimingInfo& timingInfo, caf::SRTrigger& triggerInfo);
 
   void FillTriggerEmulation(art::Handle<std::vector<int>> const& monpulsesFlat,