Add Reset() to Container and EventHeader and add UndefValueDouble

core/Constants.hpp

@@ -30,6 +30,7 @@ typedef UInt_t UInteger_t;
 typedef Short_t ShortInt_t;
 typedef Long64_t PdgCode_t;
 
+constexpr double UndefValueDouble = -999.;
 constexpr Floating_t UndefValueFloat = -999.f;
 constexpr ShortInt_t UndefValueShort = -999;
 constexpr Integer_t UndefValueInt = -999;

core/Container.cpp

@@ -5,6 +5,8 @@
 
 #include "BranchConfig.hpp"
 
+#include <algorithm>
+
 namespace AnalysisTree {
 
 template<>
@@ -32,6 +34,13 @@ void Container::Init(const AnalysisTree::BranchConfig& branch) {
   bools_.resize(branch.GetSize<bool>());
 }
 
+void Container::Reset() {
+  std::fill(doubles_.begin(), doubles_.end(), UndefValueDouble);
+  std::fill(floats_.begin(), floats_.end(), UndefValueFloat);
+  std::fill(ints_.begin(), ints_.end(), UndefValueInt);
+  std::fill(bools_.begin(), bools_.end(), false);
+}
+
 void Container::Print() const noexcept {
   if (!ints_.empty()) {
     std::cout << "Integer fields: ";

core/Container.hpp

@@ -55,6 +55,12 @@ class Container : public IndexedObject {
   }
 
   void Init(const BranchConfig& branch);
+
+  /**
+   * @brief Sets all user-defined fields to undefined values (UndefValueDouble, UndefValueFloat, UndefValueInt, false)
+   */
+  virtual void Reset();
+
   virtual void Print() const noexcept;
 
  protected:

core/EventHeader.cpp

@@ -17,6 +17,11 @@ void EventHeader::SetVertexPosition3(const TVector3& pos) {
   vtx_pos_[Exyz::kZ] = pos.Pz();
 }
 
+void EventHeader::Reset() {
+  Container::Reset();
+  vtx_pos_.fill(UndefValueFloat);
+}
+
 const EventHeader& EventHeader::GetChannel(size_t i) const {
   if (i == 0) {
     return *this;

core/EventHeader.hpp

@@ -34,6 +34,11 @@ class EventHeader : public Container {
 
   void SetVertexPosition3(const TVector3& pos);
 
+  /**
+   * @brief Sets vertex position and all user-defined fields to undefined values
+   */
+  void Reset() override;
+
   template<typename T>
   T GetField(Integer_t iField) const {
     if (iField >= 0)