Light Calorimetry: add SRLightCalo and SRTrueDeposit classes

sbnanaobj/StandardRecord/SRLightCalo.h

@@ -0,0 +1,27 @@
+/**
+ * @file   sbnanaobj/StandardRecord/SRLightCalo.h
+ * @brief  Defines data structures to store light calorimetry info (header only).
+ * @author Lynn Tung
+ * @date   November 19, 2025
+ */
+
+#ifndef SBNANAOBJ_STANDARDRECORD_SRLIGHTCALO_H
+#define SBNANAOBJ_STANDARDRECORD_SRLIGHTCALO_H
+
+#include "sbnanaobj/StandardRecord/SRConstants.h"
+
+namespace caf {
+
+  /**
+   * @brief A struct to store the reconstructed visible energy, light, and charge in a recob::Slice.
+   * 
+   */
+  struct SRLightCalo {
+    float charge  = caf::kSignalingNaN;      // the total charge in the slice (on bestplane) [# electrons]
+    float light   = caf::kSignalingNaN;      // the reconstructed total number of photons  [# photons]
+    float energy  = caf::kSignalingNaN;      // (charge+light)*scaling, [GeV]
+    int bestplane = caf::kUninitializedInt;  // plane with most number of hits 
+};
+}
+
+#endif

sbnanaobj/StandardRecord/SRSlice.h

@@ -19,6 +19,7 @@
 #include "sbnanaobj/StandardRecord/SRNuID.h"
 #include "sbnanaobj/StandardRecord/SRConstants.h"
 #include "sbnanaobj/StandardRecord/SRCVNScore.h"
+#include "sbnanaobj/StandardRecord/SRLightCalo.h"
 
 #include <climits>
 
@@ -54,6 +55,8 @@ namespace caf
 
       SRCorrectedOpFlash correctedOpFlash; //!< OpFlash corrected with using tpc information matched to this slice
 
+      SRLightCalo  lightcalo; // !< Reconstructed visible energy, light, and charge
+
       SRFakeReco fake_reco;
 
       bool is_clear_cosmic { false };         //!< Whether pandora marks the slice as a "clear" cosmic

sbnanaobj/StandardRecord/SRTrueDeposit.h

@@ -0,0 +1,26 @@
+/**
+ * @file   sbnanaobj/StandardRecord/SrTrueDeposit.h
+ * @brief  Defines data structures to store summed SimEnergyDeposits (header only).
+ * @author Lynn Tung
+ * @date   November 18, 2025
+ */
+
+#ifndef SBNANAOBJ_STANDARDRECORD_SRTRUEDEPOSIT_H
+#define SBNANAOBJ_STANDARDRECORD_SRTRUEDEPOSIT_H
+
+#include "sbnanaobj/StandardRecord/SRConstants.h"
+
+namespace caf {
+
+  /**
+   * @brief A struct to store true deposited energy, charge, and photons for each MCTruth interaction.
+   * 
+   */
+  struct SRTrueDeposit {
+    float electrons = caf::kSignalingNaN; ///< # of electrons at the site of energy deposition (after recombination, before drift)
+    float photons   = caf::kSignalingNaN; ///< # of photons at the site of energy deposition   (after recombination, before propagation) 
+    float energy    = caf::kSignalingNaN; ///< total **deposited energy** [GeV]
+  };
+}
+
+#endif

sbnanaobj/StandardRecord/SRTruthBranch.h

@@ -5,6 +5,7 @@
 #define SRTRUTHBRANCH_H
 
 #include "sbnanaobj/StandardRecord/SRTrueInteraction.h"
+#include "sbnanaobj/StandardRecord/SRTrueDeposit.h"
 #include "sbnanaobj/StandardRecord/SRMeVPrtl.h"
 
 #include <vector>
@@ -19,6 +20,7 @@ namespace caf
     ~SRTruthBranch();
 
     std::vector<SRTrueInteraction> nu;   ///< Vector of true nu or cosmic
+    std::vector<SRTrueDeposit>    dep;   ///< Vector of true energy deposits (summed), same size as nu
     size_t                        nnu;   ///< Number of true nu or cosmic
 
     std::vector<SRMeVPrtl>       prtl;   ///< If present -- information on decay of MeV "Portal" particle

sbnanaobj/StandardRecord/classes_def.xml

@@ -95,6 +95,7 @@
   </class>
 
   <class name="caf::SROpT0Finder" />
+  <class name="caf::SRLightCalo" />
 
   <class name="caf::SRSliceRecoBranch" ClassVersion="13">
    <version ClassVersion="13" checksum="4211357366"/>
@@ -174,6 +175,8 @@
    <version ClassVersion="10" checksum="4056415501"/>
   </class>
 
+  <class name="caf::SRTrueDeposit"/>
+
   <class name="caf::SRTrueParticlePlaneInfo" ClassVersion="11">
    <version ClassVersion="11" checksum="317847253"/>
    <version ClassVersion="10" checksum="1689803547"/>
@@ -215,7 +218,8 @@
    <version ClassVersion="10" checksum="1609299342"/>
   </class>
 
-  <class name="caf::SRTruthBranch" ClassVersion="11">
+  <class name="caf::SRTruthBranch" ClassVersion="12">
+   <version ClassVersion="12" checksum="3925767433"/>
    <version ClassVersion="11" checksum="1553456992"/>
    <version ClassVersion="10" checksum="2608785459"/>
   </class>
@@ -336,6 +340,7 @@
   <class name="std::vector<caf::SRTrueInteraction>" />
   <class name="std::vector<caf::SRMultiverse>" />
   <class name="std::vector<caf::SRTrueParticle>" />
+  <class name="std::vector<caf::SRTrueDeposit>" />
   <class name="std::vector<caf::SRParticleMatch>" />
   <class name="std::vector<caf::SRTrueCaloPoint>" />
   <class name="std::vector<caf::SRCaloPoint>" />