Merge pull request #42 from DUNE/feature/flash

Feature/flash

Author: sindhu-ku

duneanaobj/StandardRecord/SREnums.h

@@ -9,6 +9,7 @@
 #define DUNEANAOBJ_SRENUMS_H
 
 #include <cstddef>
+#include <limits>
 
 namespace caf
 {
@@ -100,6 +101,20 @@ namespace caf
       PartType type = kUnknown; ///< Which of the particle collections this particle lives in
       int      part = -1;       ///< Index of SRParticle in the SRInteraction
   };
+  
+  class FlashMatch
+  {
+    private:
+
+      static constexpr float NaN = std::numeric_limits<float>::signaling_NaN();
+
+    public:
+    
+      int   id            = -1;  ///< id of the matched flash in SROpticalFlash
+      float time          = NaN; ///< time of the matched flash
+      float total_pe      = NaN; ///< total pe of the matched flash
+      float hypothesis_pe = NaN; ///< hypothesis pe from reconstruction for this interaction
+  };
 
   /// Which reconstruction toolkit was used to reconstruct this FD event?
   enum FD_RECO_STACK

duneanaobj/StandardRecord/SRMeta.h

@@ -18,10 +18,10 @@ namespace caf
     public:
       bool enabled = false;     ///< Does this detector have data present in this event?
 
-      unsigned int run    = 0;
-      unsigned int subrun = 0;
-      unsigned int event  = 0;
-      unsigned int subevt = 0;
+      int run    = -1;
+      int subrun = -1;
+      int event  = -1;
+      int subevt = -1;
 
       /// detector-dependent trigger type for the relevant readout window
       int triggertype = -1;

duneanaobj/StandardRecord/SRNDLAr.h

@@ -9,6 +9,7 @@
 
 #include "duneanaobj/StandardRecord/SRTrack.h"
 #include "duneanaobj/StandardRecord/SRShower.h"
+#include "duneanaobj/StandardRecord/SROpticalFlash.h"
 
 namespace caf
 {
@@ -21,6 +22,9 @@ namespace caf
 
       std::vector<SRShower> showers;
       std::size_t           nshowers = 0;
+      
+      std::vector<FlashMatch> flash;        ///< Attributes of matched flash
+      
   };
 
   /// The information needed to uniquely identify an ND-LAr reco object
@@ -40,7 +44,8 @@ namespace caf
       std::size_t ndlp = 0;
       std::vector<SRNDLArInt> pandora;   ///< Reconstructed interactions from Pandora
       std::size_t npandora = 0;
-
+      std::vector<SROpticalFlash> flashes;      ///< Collection of flashes
+      std::size_t                 nflashes = 0;
       /// Convenience function for use mainly with SRNDTrackAssn.
       /// Given a specific reco pathway (specified with a SRNDLAr::RECO_STACK value),
       /// an interaction index, and a track index, return the associated reco object

duneanaobj/StandardRecord/SROpticalFlash.h

@@ -0,0 +1,31 @@
+////////////////////////////////////////////////////////////////////////
+/// \file    SROpticalFlash.h
+/// \brief   A general optical flash container
+/// \author  S. Kumaran <s.kumaran@uci.edu>
+
+
+#ifndef DUNEANAOBJ_SROPTICALFLASH_H
+#define DUNEANAOBJ_SROPTICALFLASH_H
+
+
+namespace caf
+{
+  /// \brief optical flash candidate
+  class SROpticalFlash
+  {
+    private:
+      // make the uses of it below more readable
+      static constexpr float NaN = std::numeric_limits<float>::signaling_NaN();
+    
+    public:
+      int    id            = -1;     ///< id of the flash>
+      int    tpc_id        = -1;     ///< id of the TPC>
+      double time          = NaN;    ///< time of the flash w.r.t. trigger time in s
+      double time_width    = NaN;    ///< width of the flash in s
+      float  total_pe      = NaN;    ///< total photoelectrons recorded by all the channels
+      
+  };
+
+} // caf
+
+#endif //DUNEANAOBJ_SROPTICALFLASH_H

duneanaobj/StandardRecord/classes_def.xml

@@ -42,7 +42,8 @@
    <version ClassVersion="10" checksum="175116589"/>
   </class>
 
-  <class name="caf::SRDetectorMeta" ClassVersion="11">
+  <class name="caf::SRDetectorMeta" ClassVersion="12">
+   <version ClassVersion="12" checksum="989983909"/>
    <version ClassVersion="11" checksum="810717581"/>
    <version ClassVersion="10" checksum="2557124789"/>
   </class>
@@ -196,12 +197,16 @@
   <class name="std::vector<caf::SRNDTrackAssn>">
   </class>
 
-  <class name="caf::SRNDLAr" ClassVersion="11">
+  <class name="caf::SRNDLAr" ClassVersion="12">
+   <version ClassVersion="12" checksum="620685221"/>
    <version ClassVersion="11" checksum="1183768120"/>
    <version ClassVersion="10" checksum="3845592869"/>
   </class>
 
-  <class name="caf::SRNDLArInt" ClassVersion="10">
+  <class name="caf::SRNDLArInt" ClassVersion="13">
+   <version ClassVersion="13" checksum="1379574528"/>
+   <version ClassVersion="12" checksum="947356810"/>
+   <version ClassVersion="11" checksum="1873045675"/>
    <version ClassVersion="10" checksum="947356810"/>
   </class>