sEPD_TreeGen: Cleanup and Refactor

- Remove saving to regular TTrees, instead only save to DST
- No need to save sepd channel phi as it can be reconstructed from the EpdGeom node that's saved to the DST by default
- Remove usage of std::format
- Remove usage of unique_ptr
- Add Print Method for verbosity and debug
- Ensure EventPlaneData has proper Reset Method to refresh buffer between events

Author: Steepspace

calibrations/sepd/sepd_eventplanecalib/EventPlaneData.cc

@@ -3,7 +3,40 @@
 EventPlaneData::EventPlaneData()
 {
   sepd_charge.fill(0);
-  sepd_phi.fill(0);
 }
 
-  
+void EventPlaneData::Reset()
+{
+  event_id = 0;
+  event_zvertex = std::numeric_limits<double>::quiet_NaN();
+  event_centrality = std::numeric_limits<double>::quiet_NaN();
+  sepd_totalcharge = std::numeric_limits<double>::quiet_NaN();
+  sepd_charge.fill(0);
+}
+
+void EventPlaneData::identify(std::ostream& os) const
+{
+  os << "--- EventPlaneData Identify ---" << std::endl;
+  os << "Event ID: " << event_id << std::endl;
+  os << "Z-Vertex: " << event_zvertex << std::endl;
+  os << "Centrality: " << event_centrality << std::endl;
+  os << "sEPD Total Charge: " << sepd_totalcharge << std::endl;
+  os << "-------------------------------" << std::endl;
+}
+
+int EventPlaneData::isValid() const
+{
+  // An object is considered invalid if the Z-vertex is still NaN
+  // or if the event ID hasn't been set (remains 0).
+  if (std::isnan(event_zvertex))
+  {
+    return 0;
+  }
+
+  if (event_id == 0)
+  {
+    return 0;
+  }
+
+  return 1;
+}

calibrations/sepd/sepd_eventplanecalib/EventPlaneData.h

@@ -1,9 +1,12 @@
 #ifndef SEPD_EVENTPLANECALIB_EVENTPLANEDATA_H
 #define SEPD_EVENTPLANECALIB_EVENTPLANEDATA_H
 
+#include "QVecDefs.h"
+
 #include <phool/PHObject.h>
 
 #include <array>
+#include <cmath>
 #include <limits>
 
 class EventPlaneData : public PHObject
@@ -12,9 +15,12 @@ class EventPlaneData : public PHObject
   EventPlaneData();
   ~EventPlaneData() override = default;
 
-  void Reset() override {*this = EventPlaneData();} // check if this works
-  // this should be in an sepd define (e.g. ../../../offline/packages/epd/EPDDefs.h)
-  static constexpr int SEPD_CHANNELS = 744;
+  EventPlaneData(const EventPlaneData&) = default;
+  EventPlaneData& operator=(const EventPlaneData&) = default;
+  EventPlaneData(EventPlaneData&&) = default;
+  EventPlaneData& operator=(EventPlaneData&&) = default;
+
+  void Reset() override;
   void set_event_id(int id) {event_id = id;}
   int get_event_id() const {return event_id;}
 
@@ -27,17 +33,19 @@ class EventPlaneData : public PHObject
   void set_sepd_charge(int channel, double chg) {sepd_charge[channel] = chg;}
   double get_sepd_charge(int channel) const {return sepd_charge[channel];}
 
-  void set_sepd_phi(int channel, double phi) {sepd_phi[channel] = phi;}
-  double get_sepd_phi(int channel) const {return sepd_phi[channel];}
-  
+  void set_event_centrality(double cent) { event_centrality = cent; }
+  double get_event_centrality() const { return event_centrality; }
+
+  void identify(std::ostream& os = std::cout) const override;
+  int isValid() const override;
+
  private:
   int event_id {0};
   double event_zvertex {std::numeric_limits<double>::quiet_NaN()};
   double event_centrality{std::numeric_limits<double>::quiet_NaN()};
   double sepd_totalcharge{std::numeric_limits<double>::quiet_NaN()};
 
-  std::array<double, SEPD_CHANNELS> sepd_charge {};
-  std::array<double, SEPD_CHANNELS> sepd_phi {};
+  std::array<double, QVecShared::SEPD_CHANNELS> sepd_charge {};
   ClassDefOverride(EventPlaneData, 1);
 };
 

calibrations/sepd/sepd_eventplanecalib/QVecDefs.h

@@ -10,7 +10,7 @@ namespace QVecShared
 {
   static constexpr size_t CENT_BINS = 8;
   static constexpr std::array<int, 3> HARMONICS = {2, 3, 4};
-  static constexpr int sepd_channels = 744;
+  static constexpr int SEPD_CHANNELS = 744;
 
   enum class ChannelStatus : int
   {