Beam timing integration and optical calibration nutples update

fcl/reco/Definitions/stage0_icarus_defs.fcl

@@ -9,6 +9,7 @@
 #include "recowire_icarus.fcl"
 #include "hitfindermodules_icarus.fcl"
 #include "timing_icarus.fcl"
+#include "timing_beam.fcl"
 #include "icarus_ophitfinder.fcl"
 #include "icarus_flashfinder.fcl"
 #include "trigger_emulation_icarus.fcl"
@@ -115,6 +116,9 @@ icarus_stage0_producers:
 
   daqPMTonbeam:                   @local::copyPMTonBeam
 
+  ### Beam timing
+  beamTiming:                     @local::icarus_beam_signal_extractor
+
   ### Purity monitoring
   purityana0:                     { module_type: "ICARUSPurityDQM" }
   purityana1:                     { module_type: "ICARUSPurityDQM" }
@@ -243,6 +247,7 @@ icarus_stage0_PMT:                  [ triggerconfig,
                                       daqTrigger,
                                       pmtconfig,
                                       daqPMT,
+                                      beamTiming,
                                       pmtconfigbaselines,
                                       pmtthr,
                                       pmtbaselines,

fcl/reco/Stage1/Run2/stage1_run2_icarus_MC.fcl

@@ -50,7 +50,7 @@ physics.producers: {
 physics.reco: [ 
                 @sequence::icarus_reco_Gauss_CryoE ,
                 @sequence::icarus_reco_Gauss_CryoW ,
-		            @sequence::icarus_reco_fm,
+                @sequence::icarus_reco_fm,
                 @sequence::icarus_tpcpmtbarycentermatch,
                 caloskimCalorimetryCryoE, caloskimCalorimetryCryoW,
                 mcassociationsGausCryoE,  mcassociationsGausCryoW
@@ -82,6 +82,11 @@ physics.analyzers.caloskimW.RawDigitproducers: ["MCDecodeTPCROI:PHYSCRATEDATATPC
 physics.producers.mcassociationsGausCryoE.HitParticleAssociations.HitModuleLabelVec: ["cluster3DCryoE"]
 physics.producers.mcassociationsGausCryoW.HitParticleAssociations.HitModuleLabelVec: ["cluster3DCryoW"]
 
+# Remove missing products in MC
+physics.analyzers.simpleLightAna.TriggerLabel: ""
+physics.analyzers.simpleLightAna.RWMLabel: ""
+physics.analyzers.simpleLightAna.OpDetWaveformLabels: ["opdaq"]
+
 services.message.destinations :
 {
   STDCOUT:

icaruscode/IcarusObj/PMTBeamSignal.h

@@ -0,0 +1,69 @@
+/**
+ * @file   icaruscode/IcarusObj/PMTBeamSignal.h
+ * @brief  Holds the event-by-event RWM or EW times
+ * @author Matteo Vicenzi (mvicenzi@bnl.gov)
+ * @date   March 14 2024
+ */
+
+#ifndef ICARUSCODE_ICARUSOBJ_PMTBEAMSIGNAL_H
+#define ICARUSCODE_ICARUSOBJ_PMTBEAMSIGNAL_H
+
+// C/C++ standard libraries
+#include <limits>
+#include <string>
+#include <cstddef>
+
+namespace icarus::timing
+{
+
+  /// Special value to denote no special channel information.
+  static constexpr auto NoChannel = std::numeric_limits<unsigned int>::max();
+  /// Special value to denote no time channel information.
+  static constexpr double NoTime = std::numeric_limits<double>::max();
+  // Special value to denote no sample information.
+  static constexpr std::size_t NoSample = 0;
+
+  /**
+   * @brief Beam time as seen by a PMT readout board.
+   *
+   * This could either be an early warning (EW) or a resistive wall monitor (RWM) time.
+   * These signals are delivered via fibers and digitized in special PMT channels.
+   *
+   * Both the time in @ref DetectorClocksElectronicsTime "electronics time scale"
+   * and time time relative to the hardware trigger are included.
+   *
+   * The information in this object may be missing: its validity should
+   * always be checked in advance with `isValid()`.
+   */
+
+  struct PMTBeamSignal
+  {
+
+    /// Special channel this time was extracted from.
+    /// These are defined in `CAEN_V1730_setup_icarus.fcl`.
+    unsigned int specialChannel = NoChannel;
+
+    /// Board on which the special channel is on (e.g: WW-TOP-A).
+    /// Should match the same format as `icarusDB::PMTChannelInfo_t::digitizerLabel`.
+    std::string digitizerLabel = "";
+
+    /// Crate this time applies to (e.g.: WW-TOP).
+    /// Corresponds to the first part of `digitizerLabel`.
+    std::string crate = "";
+
+    /// Sample within the waveform where the reference signal is found.
+    std::size_t sample = NoSample;
+
+    /// Start time in electronics time [us].
+    double startTimeAbs = NoTime;
+
+    /// Start time relative to trigger time [us].
+    double startTime = NoTime;
+
+    /// Returns whether the time is valid.
+    bool isValid() const { return (sample != NoSample); }
+  };
+
+} // namespace icarus::timing
+
+#endif // ICARUSCODE_ICARUSOBJ_PMTBEAMSIGNAL_H

icaruscode/PMT/OpReco/CMakeLists.txt

@@ -16,9 +16,10 @@ cet_build_plugin(ICARUSOpHitFinder art::module
 
 
 set(	MODULE_LIBRARIES
-                        icarusalg::Utilities
+            icarusalg::Utilities
 			sbnobj::Common_Trigger
-			larcorealg::Geometry
+    	    icaruscode::Decode_DataProducts
+            larcorealg::Geometry
 			larcore::Geometry_Geometry_service
 			lardataobj::RecoBase
 			lardataobj::Simulation
@@ -48,7 +49,7 @@ set(	MODULE_LIBRARIES
 			ROOT::Gdml
 			ROOT::FFTW
 			ROOT::Core
-                        ROOT::Tree
+            ROOT::Tree
 )
 cet_build_plugin(FakeFlash art::module LIBRARIES ${MODULE_LIBRARIES})
 cet_build_plugin(FakePhotoS art::module LIBRARIES ${MODULE_LIBRARIES})
@@ -61,6 +62,7 @@ cet_build_plugin(ICARUSOpHitAna art::module LIBRARIES ${MODULE_LIBRARIES})
 cet_build_plugin(ICARUSOpHitTuple art::module LIBRARIES ${MODULE_LIBRARIES})
 cet_build_plugin(ICARUSParticleAna art::module LIBRARIES ${MODULE_LIBRARIES})
 cet_build_plugin(TPCPMTBarycenterMatchProducer art::module LIBRARIES ${MODULE_LIBRARIES})
+cet_build_plugin(ICARUSBeamStructureAna art::module LIBRARIES ${MODULE_LIBRARIES})
 
 
 install_headers()