fix skills rules

Geant4's G4EXPATShim creates EXPAT::EXPAT (uppercase) with
IMPORTED_LOCATION set to ${EXPAT_LIBRARY}, which is empty when EXPAT
is found via conda's config-mode package (expat::expat, lowercase).

After find_package(Geant4), patch EXPAT::EXPAT with the real library
path taken from expat::expat IMPORTED_LOCATION_NOCONFIG, falling back
to find_library if needed.

Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>

.agents/rules/context_exclude_vtk.md

@@ -0,0 +1,9 @@
+---
+trigger: manual
+---
+# Context Restriction: No VTK
+When this rule is active, restrict the operational context to libraries excluding the VTK layer.
+- **Exclude Path**: `src/Vtk`
+- **Include Paths**: `src/Core`, `src/Math`, `src/HEP`, `src/Root`, `src/Python`, `src/utils`
+- **Focus**: Tomographic reconstruction algorithms, data structures in `Core`, and physical modeling in `HEP`. 
+- **Constraint**: Avoid referencing `Puppet`, `Viewport`, or any VTK-specific headers unless the user overrides this restriction.

.agents/rules/context_gcompose.md

@@ -0,0 +1,9 @@
+---
+trigger: manual
+---
+# Context Focus: gcompose Application
+When this rule is active, prioritize the `gcompose` GUI application.
+- **Primary Path**: `app/gcompose`
+- **Focus**: `MainPanel`, `ContextPanel`, `PropertiesPanel`, and `ViewportPane`.
+- **Integration**: Wiring of Qt signals/slots between the `uLib` core model and the GUI widgets.
+- **Dependency**: Reference `src/Vtk` and `src/Core` as the underlying framework for the application.

.agents/rules/context_vtk.md

@@ -0,0 +1,9 @@
+---
+trigger: always_on
+---
+# Context Inclusion: VTK
+When this rule is active, include the VTK visualization layer in the operational context.
+- **Priority Path**: `src/Vtk`
+- **Focus**: `Puppet` hierarchy, `Viewport` management, and the synchronization between domain objects and VTK props.
+- **Key Classes**: `vtkViewport`, `vtkQViewport`, `vtkObjectsContext`, and all classes in `src/Vtk/HEP/Geant`.
+- **Logic**: Ensure transformations (TRS) applied to domain objects are correctly mirrored in the visualization layer and vice versa.

.agents/rules/micromamba.md

@@ -1,7 +0,0 @@
----
-trigger: always_on
----
-
-build in build directory using always micromamba "mutom" env.
-build with make flag -j$(nproc).
-

.agents/skills/micromamba_build.md

@@ -0,0 +1,39 @@
+# Skill: Build uLib with Micromamba
+
+This skill provides instructions for building the uLib project using the micromamba environment.
+
+## Context
+- **Environment**: micromamba `uLib`
+- **Output Directory**: `build`
+- **CPU Usage**: All available cores
+
+## Instructions
+
+1. **Environment Setup**:
+   Ensure micromamba is properly initialized and the `uLib` environment is active.
+   ```bash
+   export MAMBA_EXE="/home/share/micromamba/bin/micromamba"
+   export MAMBA_ROOT_PREFIX="/home/share/micromamba"
+   eval "$(/home/share/micromamba/bin/micromamba shell hook --shell bash)"
+   micromamba activate uLib
+   ```
+
+2. **Full Rebuild (if needed)**:
+   If the `build` directory does not exist or a full reconfiguration is required:
+   ```bash
+   conan profile detect --force
+   conan install . --output-folder=build --build=missing
+   cmake --preset conan-release
+   ```
+
+3. **Incremental Build**:
+   Run the build command from the root directory, pointing to the `build` folder and using all cores.
+   ```bash
+   cmake --build build -j$(nproc)
+   ```
+
+4. **Specific Target Build**:
+   To build a specific target (e.g., gcompose):
+   ```bash
+   cmake --build build --target gcompose -j$(nproc)
+   ```

.gitignore

@@ -18,3 +18,5 @@ test_props.xml
 test_props2.xml
 test_boost.cpp
 .claude/settings.json
+build_output.log
+configure_output.log

CLAUDE.md

@@ -10,7 +10,7 @@ This file provides guidance to Claude Code (claude.ai/code) when working with co
 export MAMBA_EXE="/home/share/micromamba/bin/micromamba"
 export MAMBA_ROOT_PREFIX="/home/share/micromamba"
 eval "$(/home/share/micromamba/bin/micromamba shell hook --shell bash)"
-micromamba activate mutom
+micromamba activate uLib
 
 # Configure (from repo root, using Conan preset — uses Ninja + ccache)
 cmake --preset conan-release

CMakeLists.txt

@@ -15,11 +15,20 @@ if(POLICY CMP0167)
     cmake_policy(SET CMP0167 NEW)
 endif()
 
-
 ## -------------------------------------------------------------------------- ##
 
 project(uLib)
 
+# Applica la flag SOLO se il compilatore è GCC
+if(CMAKE_CXX_COMPILER_ID STREQUAL "GNU")
+    add_compile_options(-fno-merge-constants)
+endif()
+
+# Disabilita il warning se il compilatore è Clang (o AppleClang)
+if(CMAKE_CXX_COMPILER_ID MATCHES "Clang")
+    add_compile_options(-Wno-ignored-optimization-argument)
+endif()
+
 # CUDA Toolkit seems to be missing locally. Toggle ON if nvcc is made available.
 option(USE_CUDA "Enable CUDA support" OFF)
 if(USE_CUDA)
@@ -124,6 +133,8 @@ find_package(Eigen3 CONFIG REQUIRED)
 get_target_property(EIGEN3_INCLUDE_DIRS Eigen3::Eigen INTERFACE_INCLUDE_DIRECTORIES)
 include_directories(${EIGEN3_INCLUDE_DIRS})
 
+find_package(OpenMP)
+
 find_package(ROOT CONFIG REQUIRED)
 include(${ROOT_USE_FILE})
 
@@ -145,6 +156,8 @@ else()
                                 IOXML
                                 IOXMLParser
                                 ImagingCore
+                                ImagingHybrid
+                                ImagingSources
                                 InteractionStyle
                                 InteractionWidgets
                                 RenderingAnnotation

README.md

@@ -30,13 +30,13 @@ You can create and activate the environment using either `micromamba` or `conda`
 **Using Micromamba:**
 ```bash
 micromamba env create -f condaenv.yml
-micromamba activate mutom
+micromamba activate uLib
 ```
 
 **Using Conda:**
 ```bash
 conda env create -f condaenv.yml
-conda activate mutom
+conda activate uLib
 ```
 
 ### Configure and Build
@@ -65,12 +65,24 @@ cmake --preset conan-release
 cmake --build build -j$(nproc)
 ```
 
+5. **Clean build (wipe and rebuild everything):**
+```bash
+cmake --build build --clean-first -j$(nproc)
+```
+
+6. **Run tests:**
+```bash
+cmake --build build --target test -j$(nproc)
+# or equivalently:
+ctest --test-dir build --output-on-failure -j$(nproc)
+```
+
 #### LLVM/Clang build (clang + lld + ccache — fastest)
 
 A `fast` conan profile is provided that uses **clang**, **lld** (LLVM linker), and **ccache**. Install them into your environment first:
 
 ```bash
-micromamba install -n mutom -y clang clangxx lld -c conda-forge
+micromamba install -n uLib -y clang clangxx lld -c conda-forge
 ```
 
 Then build using the `fast` profile:
@@ -91,6 +103,6 @@ The `fast` profile is defined at `~/.conan2/profiles/fast` and sets:
 ### Make python package
 
 ```bash
-micromamba run -n mutom env USE_CUDA=ON poetry install
+micromamba run -n uLib env USE_CUDA=ON poetry install
 ```
 

conanfile.txt

@@ -1,6 +1,6 @@
 [requires]
 eigen/3.4.0
-boost/1.83.0
+boost/1.86.0
 # pybind11/3.0.2
 hdf5/1.14.3
 

condaenv.yml

@@ -1,4 +1,4 @@
-name: mutom
+name: uLib
 channels:
   - conda-forge
 dependencies:
@@ -7,4 +7,13 @@ dependencies:
   - cmake
   - conan
   - root
-  - vtk
+  - vtk=9.4 # VTK 9.4
+  - pybind11
+#  - boost=1.86.0 # requested by VTK 9.4
+  - ninja
+  - clang
+  - clangxx
+  - lld
+  - ccache
+  - OpenMP
+  - Geant4

src/Core/Archives.h

@@ -80,6 +80,9 @@ template <class ArchiveImplementation> class polymorphic_iarchive_route;
 namespace boost {
 namespace serialization {
 template <typename T> struct hrp;
+template <typename T> struct hrp_val;
+template <typename T> struct hrp_enum;
+template <typename T> struct hrp_enum_val;
 }
 } // namespace boost
 
@@ -177,6 +180,24 @@ public:
     return *this->This();
   }
 
+  template <class T>
+  Archive &operator>>(const boost::serialization::hrp_val<T> &t) {
+    this->This()->load_override(const_cast<boost::serialization::hrp_val<T> &>(t));
+    return *this->This();
+  }
+
+  template <class T>
+  Archive &operator>>(const boost::serialization::hrp_enum<T> &t) {
+    this->This()->load_override(const_cast<boost::serialization::hrp_enum<T> &>(t));
+    return *this->This();
+  }
+
+  template <class T>
+  Archive &operator>>(const boost::serialization::hrp_enum_val<T> &t) {
+    this->This()->load_override(const_cast<boost::serialization::hrp_enum_val<T> &>(t));
+    return *this->This();
+  }
+
   // the & operator
   template <class T> Archive &operator&(T &t) { return *(this->This()) >> t; }
 
@@ -190,6 +211,21 @@ public:
     return *(this->This()) >> t;
   }
 
+  template <class T>
+  Archive &operator&(const boost::serialization::hrp_val<T> &t) {
+    return *(this->This()) >> t;
+  }
+
+  template <class T>
+  Archive &operator&(const boost::serialization::hrp_enum<T> &t) {
+    return *(this->This()) >> t;
+  }
+
+  template <class T>
+  Archive &operator&(const boost::serialization::hrp_enum_val<T> &t) {
+    return *(this->This()) >> t;
+  }
+
   // the == operator
   template <class T> Archive &operator==(T &t) { return this->operator&(t); }
 
@@ -230,11 +266,61 @@ public:
     return *this->This();
   }
 
+  template <class T> Archive &operator<<(const boost::serialization::hrp<T> &t) {
+    this->This()->save_override(t);
+    return *this->This();
+  }
+
+  template <class T> Archive &operator<<(const boost::serialization::hrp_val<T> &t) {
+    this->This()->save_override(t);
+    return *this->This();
+  }
+
+  template <class T> Archive &operator<<(const boost::serialization::hrp_enum<T> &t) {
+    this->This()->save_override(t);
+    return *this->This();
+  }
+
+  template <class T> Archive &operator<<(const boost::serialization::hrp_enum_val<T> &t) {
+    this->This()->save_override(t);
+    return *this->This();
+  }
+
+  template <class T> Archive &operator<<(const boost::serialization::nvp<T> &t) {
+    this->This()->save_override(t);
+    return *this->This();
+  }
+ 
   // the & operator
   template <class T> Archive &operator&(const T &t) {
     return *this->This() << t;
   }
 
+  template <class T>
+  Archive &operator&(const boost::serialization::hrp<T> &t) {
+    return *this->This() << t;
+  }
+
+  template <class T>
+  Archive &operator&(const boost::serialization::hrp_val<T> &t) {
+    return *this->This() << t;
+  }
+
+  template <class T>
+  Archive &operator&(const boost::serialization::hrp_enum<T> &t) {
+    return *this->This() << t;
+  }
+
+  template <class T>
+  Archive &operator&(const boost::serialization::hrp_enum_val<T> &t) {
+    return *this->This() << t;
+  }
+
+  template <class T>
+  Archive &operator&(const boost::serialization::nvp<T> &t) {
+    return *this->This() << t;
+  }
+
   // the == operator
   template <class T> Archive &operator==(T &t) { return this->operator&(t); }
 

src/Core/CMakeLists.txt

@@ -21,6 +21,7 @@ set(HEADERS
     StringReader.h 
     Threads.h
     Monitor.h
+    Property.h
     Types.h 
     Uuid.h 
     Vector.h

src/Core/Object.cpp

@@ -65,10 +65,19 @@ public:
   std::vector<Slot> slov;
   std::vector<PropertyBase*> m_Properties;
   std::vector<PropertyBase*> m_DynamicProperties;
+  std::vector<PropertyBase*> m_DisplayProperties;
   bool m_SignalsBlocked;
 };
 
 // Implementations of Property methods
+void Object::RegisterDisplayProperty(PropertyBase* prop) {
+    if (prop) d->m_DisplayProperties.push_back(prop);
+}
+
+const std::vector<PropertyBase*>& Object::GetDisplayProperties() const {
+    return d->m_DisplayProperties;
+}
+
 void Object::RegisterProperty(PropertyBase* prop) {
     if (prop) {
         d->m_Properties.push_back(prop);
@@ -104,24 +113,9 @@ void Object::NotifyPropertiesUpdated() {
     for (auto* p : d->m_DynamicProperties) p->Updated();
 }
 
-// In Object.h, the template serialize needs to be updated to call property serialization.
-// However, since Object::serialize is a template in the header, we might need a helper here.
-
-template <class ArchiveT>
-void Object::serialize(ArchiveT &ar, const unsigned int version) {
-    ar & boost::serialization::make_nvp("InstanceName", d->m_InstanceName);
-    for (auto* prop : d->m_Properties) {
-        prop->serialize(ar, version);
-    }
-}
-
 void Object::Updated() { ULIB_SIGNAL_EMIT(Object::Updated); }
 void Object::PropertyUpdated() { ULIB_SIGNAL_EMIT(Object::PropertyUpdated); }
 
-template <class ArchiveT>
-void Object::save_override(ArchiveT &ar, const unsigned int version) {}
-
-
 ////////////////////////////////////////////////////////////////////////////////
 ////////////////////////////////////////////////////////////////////////////////
 ////////////////////////////////////////////////////////////////////////////////

src/Core/Object.h

@@ -92,9 +92,11 @@ public:
 
   ////////////////////////////////////////////////////////////////////////////
   // PROPERTIES //
-  void RegisterProperty(PropertyBase* prop);
+  virtual void RegisterProperty(PropertyBase* property);
-  void RegisterDynamicProperty(PropertyBase* prop);
+  virtual void RegisterDynamicProperty(PropertyBase* property);
-  const std::vector<PropertyBase*>& GetProperties() const;
+  virtual void RegisterDisplayProperty(PropertyBase* property);
+  virtual const std::vector<PropertyBase*>& GetProperties() const;
+  virtual const std::vector<PropertyBase*>& GetDisplayProperties() const;
   PropertyBase* GetProperty(const std::string& name) const;
 
   /** @brief Sends an Updated signal for all properties of this object. useful for real-time UI refresh. */
@@ -124,7 +126,7 @@ public:
   virtual void serialize(Archive::log_archive & ar, const unsigned int version) {}
 
   template <class ArchiveT>
-  void save_override(ArchiveT &ar, const unsigned int version);
+  void save_override(ArchiveT &ar, const unsigned int version) {}
 
   void SaveConfig(std::ostream &os, int version = 0);
   void LoadConfig(std::istream &is, int version = 0);

src/Core/Property.h

@@ -5,19 +5,23 @@
 #include <vector>
 #include <sstream>
 #include <typeinfo>
-#include <typeindex> // Added
+#include <typeindex>
 #include <boost/serialization/nvp.hpp>
 #include <boost/lexical_cast.hpp>
-#include <vector>
 #include <boost/type_traits/is_class.hpp>
 #include <boost/mpl/bool.hpp>
-#include <boost/serialization/serialization.hpp>
+#include <set>
+#include <boost/type_traits/is_base_of.hpp>
 #include "Core/Archives.h"
 #include "Core/Signal.h"
 #include "Core/Object.h"
 
 namespace uLib {
 
+namespace Archive {
+class property_register_archive;
+}
+
 /**
  * @brief Base class for properties to allow runtime listing and identification.
  */
@@ -27,7 +31,7 @@ public:
   virtual const std::string& GetName() const = 0;
   virtual const char* GetTypeName() const = 0;
   virtual std::string GetValueAsString() const = 0;
-  virtual std::type_index GetTypeIndex() const = 0; // Added
+  virtual std::type_index GetTypeIndex() const = 0;
   virtual const std::string& GetUnits() const = 0;
   virtual void SetUnits(const std::string& units) = 0;
   virtual const std::vector<std::string>& GetEnumLabels() const {
@@ -61,42 +65,29 @@ public:
   virtual void serialize(Archive::hrt_oarchive & ar, const unsigned int version) override = 0;
   virtual void serialize(Archive::hrt_iarchive & ar, const unsigned int version) override = 0;
   virtual void serialize(Archive::log_archive & ar, const unsigned int version) override = 0;
+  
+  virtual void serialize(Archive::property_register_archive & ar, const unsigned int v) = 0;
 };
 
-
-
-
-
-
-
-
 /**
  * @brief Template class for typed properties.
  */
 template <typename T>
 class Property : public PropertyBase {
 public:
-  // PROXY: Use an existing variable as back-end storage
   Property(Object* owner, const std::string& name, T* valuePtr, const std::string& units = "", const std::string& group = "")
     : m_owner(owner), m_name(name), m_units(units), m_group(group), m_value(valuePtr), m_own(false),
       m_HasRange(false), m_HasDefault(false), m_ReadOnly(false) {
-      if (m_owner) {
+      if (m_owner) m_owner->RegisterProperty(this);
-          m_owner->RegisterProperty(this);
-      }
   }
 
-  // MANAGED: Create and own internal storage
   Property(Object* owner, const std::string& name, const T& defaultValue = T(), const std::string& units = "", const std::string& group = "")
     : m_owner(owner), m_name(name), m_units(units), m_group(group), m_value(new T(defaultValue)), m_own(true),
       m_HasRange(false), m_HasDefault(true), m_Default(defaultValue), m_ReadOnly(false) {
-      if (m_owner) {
+      if (m_owner) m_owner->RegisterProperty(this);
-          m_owner->RegisterProperty(this);
-      }
   }
 
-  virtual ~Property() {
+  virtual ~Property() { if (m_own) delete m_value; }
-      if (m_own) delete m_value;
-  }
 
   // Identification
   virtual const std::string& GetName() const override { return m_name; }
@@ -107,36 +98,18 @@ public:
   virtual const std::string& GetGroup() const override { return m_group; }
   virtual void SetGroup(const std::string& group) override { m_group = group; }
 
-  
   std::string GetValueAsString() const override {
-      try {
+      try { return boost::lexical_cast<std::string>(*m_value); }
-          return boost::lexical_cast<std::string>(*m_value);
+      catch (...) { std::stringstream ss; ss << *m_value; return ss.str(); }
-      } catch (const boost::bad_lexical_cast&) {
-          std::stringstream ss;
-          ss << *m_value;
-          return ss.str();
-      }
   }
 
   // Accessors
   const T& Get() const { return *m_value; }
-  template<typename U = T>
-  typename std::enable_if<std::is_arithmetic<U>::value, void>::type
-  ValidateT(T& val) {
-    if (m_HasRange) {
-        if (val < m_Min) val = m_Min;
-        if (val > m_Max) val = m_Max;
-    }
-  }
-
-  template<typename U = T>
-  typename std::enable_if<!std::is_arithmetic<U>::value, void>::type
-  ValidateT(T& val) {
-  }
-
   void Set(const T& value) {
     T val = value;
-    ValidateT<T>(val);
+    if constexpr (std::is_arithmetic<T>::value) {
+      if (m_HasRange) { if (val < m_Min) val = m_Min; if (val > m_Max) val = m_Max; }
+    }
     if (*m_value != val) {
       *m_value = val;
       ULIB_SIGNAL_EMIT(Property<T>::PropertyChanged);
@@ -150,65 +123,55 @@ public:
   void SetReadOnly(bool ro) { m_ReadOnly = ro; }
 
   virtual bool IsReadOnly() const override { return m_ReadOnly; }
-
   virtual bool HasRange() const override { return m_HasRange; }
+  virtual double GetMin() const override { return m_HasRange ? convert_to_double(m_Min) : 0.0; }
+  virtual double GetMax() const override { return m_HasRange ? convert_to_double(m_Max) : 0.0; }
   
-  template<typename U = T>
+  const T& GetMinTyped() const { return m_Min; }
-  typename std::enable_if<std::is_arithmetic<U>::value, double>::type
+  const T& GetMaxTyped() const { return m_Max; }
-  GetMinT() const { return (double)m_Min; }
-
-  template<typename U = T>
-  typename std::enable_if<!std::is_arithmetic<U>::value, double>::type
-  GetMinT() const { return 0.0; }
-
-  template<typename U = T>
-  typename std::enable_if<std::is_arithmetic<U>::value, double>::type
-  GetMaxT() const { return (double)m_Max; }
-
-  template<typename U = T>
-  typename std::enable_if<!std::is_arithmetic<U>::value, double>::type
-  GetMaxT() const { return 0.0; }
-
-  virtual double GetMin() const override { return GetMinT<T>(); }
-  virtual double GetMax() const override { return GetMaxT<T>(); }
 
   virtual bool HasDefault() const override { return m_HasDefault; }
   virtual std::string GetDefaultValueAsString() const override {
-      try { return boost::lexical_cast<std::string>(m_Default); } 
+      try { return boost::lexical_cast<std::string>(m_Default); } catch (...) { return ""; }
-      catch (...) { return ""; }
   }
 
-  // Operators for seamless usage
+  // Operators
   operator const T&() const { return *m_value; }
-  Property& operator=(const T& value) {
+  Property& operator=(const T& value) { Set(value); return *this; }
-    Set(value);
-    return *this;
-  }
 
-  // Signals
   signals:
   virtual void PropertyChanged() { ULIB_SIGNAL_EMIT(Property<T>::PropertyChanged); }
 
+private:
+  template <typename U>
+  static double convert_to_double(const U& val) {
+    return convert_to_double_impl(val, typename std::is_arithmetic<U>::type());
+  }
+  template <typename U>
+  static double convert_to_double_impl(const U& val, std::true_type) { return (double)val; }
+  template <typename U>
+  static double convert_to_double_impl(const U& val, std::false_type) { return 0.0; }
+
+public:
   // Serialization
   template <class ArchiveT>
-  void serialize_impl(ArchiveT & ar, const unsigned int version) {
+  void serialize_helper(ArchiveT & ar, const unsigned int version) {
-    ar & boost::serialization::make_nvp(m_name.c_str(), *m_value);
+    ar & boost::serialization::make_hrp(m_name.c_str(), *m_value, m_units.c_str());
   }
 
-  void serialize(Archive::xml_oarchive & ar, const unsigned int v) override { serialize_impl(ar, v); }
+  virtual void serialize(Archive::xml_oarchive & ar, const unsigned int v) override { serialize_helper(ar, v); }
-  void serialize(Archive::xml_iarchive & ar, const unsigned int v) override { serialize_impl(ar, v); }
+  virtual void serialize(Archive::xml_iarchive & ar, const unsigned int v) override { serialize_helper(ar, v); }
-  void serialize(Archive::text_oarchive & ar, const unsigned int v) override { serialize_impl(ar, v); }
+  virtual void serialize(Archive::text_oarchive & ar, const unsigned int v) override { serialize_helper(ar, v); }
-  void serialize(Archive::text_iarchive & ar, const unsigned int v) override { serialize_impl(ar, v); }
+  virtual void serialize(Archive::text_iarchive & ar, const unsigned int v) override { serialize_helper(ar, v); }
-  void serialize(Archive::hrt_oarchive & ar, const unsigned int v) override { serialize_impl(ar, v); }
+  virtual void serialize(Archive::hrt_oarchive & ar, const unsigned int v) override { serialize_helper(ar, v); }
-  void serialize(Archive::hrt_iarchive & ar, const unsigned int v) override { serialize_impl(ar, v); }
+  virtual void serialize(Archive::hrt_iarchive & ar, const unsigned int v) override { serialize_helper(ar, v); }
-  void serialize(Archive::log_archive & ar, const unsigned int v) override { serialize_impl(ar, v); }
+  virtual void serialize(Archive::log_archive & ar, const unsigned int v) override { serialize_helper(ar, v); }
   
-  virtual void Updated() override {
+  virtual void serialize(Archive::property_register_archive & ar, const unsigned int v) override;
-      PropertyBase::Updated();
-      this->PropertyChanged();
-  }
 
-private:
+  virtual void Updated() override { PropertyBase::Updated(); this->PropertyChanged(); }
+
+protected:
   std::string m_name;
   std::string m_units;
   std::string m_group;
@@ -224,20 +187,7 @@ private:
 };
 
 /**
- * @brief Conveninent typedefs for common property types.
+ * @brief Property specialized for enumerations.
- */
-typedef Property<std::string> StringProperty;
-typedef Property<int>           IntProperty;
-typedef Property<unsigned int>  UIntProperty;
-typedef Property<long>          LongProperty;
-typedef Property<unsigned long> ULongProperty;
-typedef Property<float>         FloatProperty;
-typedef Property<double>        DoubleProperty;
-typedef Property<Bool_t>        BoolProperty;
-
-
-/**
- * @brief Property specialized for enumerations, providing labels for GUI representations.
  */
 class EnumProperty : public Property<int> {
 public:
@@ -248,68 +198,49 @@ public:
   const char* GetTypeName() const override { return "Enum"; }
   virtual std::type_index GetTypeIndex() const override { return std::type_index(typeid(EnumProperty)); }
 
+  template <class ArchiveT>
+  void serialize_enum_helper(ArchiveT & ar, const unsigned int version) {
+    ar & boost::serialization::make_hrp_enum(m_name.c_str(), *m_value, m_Labels, m_units.c_str());
+  }
+
+  virtual void serialize(Archive::xml_oarchive & ar, const unsigned int v) override { serialize_enum_helper(ar, v); }
+  virtual void serialize(Archive::xml_iarchive & ar, const unsigned int v) override { serialize_enum_helper(ar, v); }
+  virtual void serialize(Archive::text_oarchive & ar, const unsigned int v) override { serialize_enum_helper(ar, v); }
+  virtual void serialize(Archive::text_iarchive & ar, const unsigned int v) override { serialize_enum_helper(ar, v); }
+  virtual void serialize(Archive::hrt_oarchive & ar, const unsigned int v) override { serialize_enum_helper(ar, v); }
+  virtual void serialize(Archive::hrt_iarchive & ar, const unsigned int v) override { serialize_enum_helper(ar, v); }
+  virtual void serialize(Archive::log_archive & ar, const unsigned int v) override { serialize_enum_helper(ar, v); }
+  
+  virtual void serialize(Archive::property_register_archive & ar, const unsigned int v) override;
+
 private:
   std::vector<std::string> m_Labels;
 };
 
-
-
-
-
-
-
-
-
-
-
-/**
- * @brief Macro to simplify property declaration within a class.
- * Usage: ULIB_PROPERTY(float, Width, 1.0f)
- * It creates a raw member m_Width and a Property proxy Width.
- */
-#define ULIB_PROPERTY(type, name, defaultValue) \
-    type m_##name = defaultValue; \
-    Property<type> name = Property<type>(this, #name, &m_##name);
-
-
 } // namespace uLib
 
 namespace uLib {
 namespace Archive {
-class property_register_archive;
-} // namespace Archive
-} // namespace uLib
 
-namespace boost {
+class property_register_archive
-namespace archive {
+    : public boost::archive::detail::common_oarchive<property_register_archive> {
-namespace detail {
+protected:
-template <>
-class interface_oarchive<uLib::Archive::property_register_archive>
-    : public uLib_interface_oarchive<uLib::Archive::property_register_archive> {};
-} // namespace detail
-} // namespace archive
-} // namespace boost
-
-namespace uLib {
-namespace Archive {
-
-/**
- * @brief A special archive that creates and registers Property proxies
- * for any member it encounters wrapped in HRP().
- */
-class property_register_archive : 
-    public boost::archive::detail::common_oarchive<property_register_archive>
-{
     Object* m_Object;
+    bool m_DisplayOnly;
 public:
     friend class boost::archive::detail::interface_oarchive<property_register_archive>;
     friend class boost::archive::save_access;
 
-    typedef boost::archive::detail::common_oarchive<property_register_archive> detail_common_oarchive;
+    using boost::archive::detail::common_oarchive<property_register_archive>::save_override;
 
-    property_register_archive(Object* obj) : 
+    property_register_archive(Object* obj, bool displayOnly = false) : 
         boost::archive::detail::common_oarchive<property_register_archive>(boost::archive::no_header),
-        m_Object(obj) {}
+        m_Object(obj), m_DisplayOnly(displayOnly) {
+        if (obj) m_Visited.insert(dynamic_cast<const void*>(obj));
+    }
+
+    template<class T> property_register_archive &operator&(const T &t) { this->save_override(t); return *this; }
+    template<class T> property_register_archive &operator<<(const T &t) { this->save_override(t); return *this; }
     
     std::string GetCurrentGroup() const {
         std::string group;
@@ -320,77 +251,91 @@ public:
         return group;
     }
 
-    // Core logic: encounter HRP -> Create Dynamic Property
+    template<class T> void register_property(Property<T>& p) {
-    template<class T>
+        save_property_impl(p.GetName().c_str(), const_cast<T&>(p.Get()), p.GetUnits().c_str(), 
-    void save_override(const boost::serialization::hrp<T> &t) {
+                          p.HasRange(), p.GetMinTyped(), p.GetMaxTyped(), p.IsReadOnly());
-        if (m_Object) {
+    }
-            Property<T>* p = new Property<T>(m_Object, t.name(), &const_cast<boost::serialization::hrp<T>&>(t).value(), t.units() ? t.units() : "", GetCurrentGroup());
+
-            if (t.has_range()) p->SetRange(t.min_val(), t.max_val());
+    void register_enum_property(EnumProperty& p) {
-            if (t.has_default()) p->SetDefault(t.default_val());
+        if (!m_Object) return;
-            p->SetReadOnly(t.is_read_only());
+        EnumProperty* newP = new EnumProperty(m_Object, p.GetName(), const_cast<int*>(&p.Get()), p.GetEnumLabels(), p.GetUnits(), GetCurrentGroup());
-            m_Object->RegisterDynamicProperty(p);
+        newP->SetReadOnly(p.IsReadOnly());
+        if (m_DisplayOnly) {
+            m_Object->RegisterDisplayProperty(newP);
+            Object* obj = m_Object;
+            Object::connect(newP, &PropertyBase::Updated, [obj]() { obj->Updated(); });
+        } else {
+            m_Object->RegisterDynamicProperty(newP);
         }
     }
 
-    template<class T>
+    template<class T> void save_property_impl(const char* name, T& val, const char* units, bool hasRange, const T& minVal, const T& maxVal, bool isReadOnly) {
-    void save_override(const boost::serialization::hrp_val<T> &t) {
         if (m_Object) {
-            // Note: hrp_val stores by value. Property usually points to existing data.
+            Property<T>* p = new Property<T>(m_Object, name, &val, units ? units : "", GetCurrentGroup());
-            // But here we are registering properties from HRP wrappers.
+            set_range_helper(p, hasRange, minVal, maxVal, typename std::is_arithmetic<T>::type());
-            // If it's hrp_val, it means it's an rvalue from a getter. 
+            p->SetReadOnly(isReadOnly);
-            // The hrp_val wrapper itself owns the value.
+            if (m_DisplayOnly) {
-            // However, the property_register_archive is temporary.
+                m_Object->RegisterDisplayProperty(p);
-            // This is a bit tricky. Usually HRP(rvalue) is meant for read-only display.
+                Object* obj = m_Object;
-            // Let's use the address of the value in the wrapper, but mark it read-only.
+                Object::connect(p, &PropertyBase::Updated, [obj]() { obj->Updated(); });
-            Property<T>* p = new Property<T>(m_Object, t.name(), &const_cast<boost::serialization::hrp_val<T>&>(t).value(), t.units() ? t.units() : "", GetCurrentGroup());
+            } else {
-            if (t.has_range()) p->SetRange(t.min_val(), t.max_val());
+                m_Object->RegisterDynamicProperty(p);
-            if (t.has_default()) p->SetDefault(t.default_val());
+            }
-            p->SetReadOnly(t.is_read_only());
-            m_Object->RegisterDynamicProperty(p);
         }
     }
 
-    template<class T>
+    template<class U> static void set_range_helper(Property<U>* p, bool hasRange, const U& minVal, const U& maxVal, std::true_type) { if (hasRange) p->SetRange(minVal, maxVal); }
-    void save_override(const boost::serialization::hrp_enum<T> &t) {
+    template<class U> static void set_range_helper(Property<U>* p, bool hasRange, const U& minVal, const U& maxVal, std::false_type) {}
+
+    template<class T> void save_override(const boost::serialization::hrp<T> &t) {
+        // To handle T correctly without deduction issues, we assume T can be passed to save_property_impl
+        T dummy = T(); // Ensure we can construct T
+        save_property_impl(t.name(), const_cast<boost::serialization::hrp<T>&>(t).value(), t.units(), t.has_range(), t.has_range() ? t.min_val() : dummy, t.has_range() ? t.max_val() : dummy, t.is_read_only());
+    }
+    template<class T> void save_override(const boost::serialization::hrp_val<T> &t) {
+        T dummy = T();
+        save_property_impl(t.name(), const_cast<boost::serialization::hrp_val<T>&>(t).value(), t.units(), t.has_range(), t.has_range() ? t.min_val() : dummy, t.has_range() ? t.max_val() : dummy, t.is_read_only());
+    }
+    template<class T> void save_override(const boost::serialization::hrp_enum<T> &t) {
         if (m_Object) {
             EnumProperty* p = new EnumProperty(m_Object, t.name(), (int*)&const_cast<boost::serialization::hrp_enum<T>&>(t).value(), t.labels(), t.units() ? t.units() : "", GetCurrentGroup());
             p->SetReadOnly(t.is_read_only());
-            m_Object->RegisterDynamicProperty(p);
+            if (m_DisplayOnly) { m_Object->RegisterDisplayProperty(p); Object* obj = m_Object; Object::connect(p, &PropertyBase::Updated, [obj]() { obj->Updated(); }); }
+            else { m_Object->RegisterDynamicProperty(p); }
         }
     }
-
+    template<class T> void save_override(const boost::serialization::hrp_enum_val<T> &t) {
-    template<class T>
-    void save_override(const boost::serialization::hrp_enum_val<T> &t) {
         if (m_Object) {
             EnumProperty* p = new EnumProperty(m_Object, t.name(), (int*)&const_cast<boost::serialization::hrp_enum_val<T>&>(t).value(), t.labels(), t.units() ? t.units() : "", GetCurrentGroup());
             p->SetReadOnly(t.is_read_only());
-            m_Object->RegisterDynamicProperty(p);
+            if (m_DisplayOnly) { m_Object->RegisterDisplayProperty(p); Object* obj = m_Object; Object::connect(p, &PropertyBase::Updated, [obj]() { obj->Updated(); }); }
+            else { m_Object->RegisterDynamicProperty(p); }
         }
     }
 
-    // Handle standard NVPs by recursing (important for base classes)
+    template<class T> void save_override(const boost::serialization::nvp<T> &t) {
-    template<class T>
-    void save_override(const boost::serialization::nvp<T> &t) {
         if (t.name()) m_GroupStack.push_back(t.name());
         this->save_helper(t.const_value(), typename boost::is_class<T>::type());
         if (t.name()) m_GroupStack.pop_back();
     }
 
-    // Recursion for nested classes, ignore primitives
+    void save_override(const std::string &t) {}
-    template<class T>
+    template<class T> void save_override(T * const & t) {
-    void save_override(const T &t) {
+        if (!t) return;
-        this->save_helper(t, typename boost::is_class<T>::type());
+        this->save_pointer_helper(t, typename boost::is_base_of<Object, T>::type());
     }
-
+    template<class T> void save_pointer_helper(T* t, boost::mpl::true_) {
-    template<class T>
+        const void* ptr = dynamic_cast<const void*>(t);
-    void save_helper(const T &t, boost::mpl::true_) {
+        if (m_Visited.find(ptr) != m_Visited.end()) return;
-        boost::serialization::serialize_adl(*this, const_cast<T&>(t), 0);
+        m_Visited.insert(ptr);
+        this->save_override(*t);
     }
+    template<class T> void save_pointer_helper(T* t, boost::mpl::false_) {}
+    template<class T> void save_override(const T &t) { this->save_helper(t, typename boost::is_class<T>::type()); }
+    template<class T> void save_helper(const T &t, boost::mpl::true_) { boost::serialization::serialize_adl(*this, const_cast<T&>(t), 0); }
+    void save_helper(const std::string &t, boost::mpl::true_) {}
+    template<class T> void save_helper(const T &t, boost::mpl::false_) {}
 
-    template<class T>
-    void save_helper(const T &t, boost::mpl::false_) {}
-
-    // Required attribute overrides for common_oarchive
     void save_override(const boost::archive::object_id_type & t) {}
     void save_override(const boost::archive::object_reference_type & t) {}
     void save_override(const boost::archive::version_type & t) {}
@@ -400,30 +345,30 @@ public:
     void save_override(const boost::archive::class_name_type & t) {}
     void save_override(const boost::archive::tracking_type & t) {}
 
-private:
+protected:
     std::vector<std::string> m_GroupStack;
+    std::set<const void*> m_Visited;
 };
 
+} // namespace Archive
+} // namespace uLib
 
+namespace uLib {
+
+template <typename T>
+inline void Property<T>::serialize(Archive::property_register_archive & ar, const unsigned int v) {
+  ar.register_property(*this);
+}
+
+inline void EnumProperty::serialize(Archive::property_register_archive & ar, const unsigned int v) {
+  ar.register_enum_property(*this);
+}
+
+namespace Archive {
 
-/**
- * @brief Convenience macro to automatically activate and register all HRP members
- * as uLib properties. Usage: ULIB_ACTIVATE_PROPERTIES(obj)
- */
 #define ULIB_ACTIVATE_PROPERTIES(obj) \
     { uLib::Archive::property_register_archive _ar_tmp(&(obj)); _ar_tmp & (obj); }
 
-/**
- * @brief Declares a private member that automatically calls ULIB_ACTIVATE_PROPERTIES
- * in every constructor of the class. Place this macro as the last declaration
- * inside the class body (before the closing brace).
- *
- * Usage: ULIB_DECLARE_PROPERTIES(ClassName)
- *
- * This replaces per-constructor ULIB_ACTIVATE_PROPERTIES(*this) calls.
- * RegisterDynamicProperty deduplicates by qualified name, so re-registration
- * from inherited activators in a hierarchy is safe.
- */
 #define ULIB_DECLARE_PROPERTIES(SelfType) \
 private: \
     struct _PropActivator { \
@@ -434,6 +379,31 @@ private: \
     } _prop_activator{this};
 
 } // namespace Archive
+
+// Convenience macro: declares a named Property<T> member with a default value.
+// Usage inside a class body (requires 'this' to be available, so use in-class initializer):
+//   ULIB_PROPERTY(int, MyProp, 42)
+#define ULIB_PROPERTY(type, name, defaultVal) \
+    ::uLib::Property<type> name{this, #name, (type)(defaultVal)};
+
+// Common property type aliases
+typedef Property<bool>        BoolProperty;
+typedef Property<int>         IntProperty;
+typedef Property<float>       FloatProperty;
+typedef Property<double>      DoubleProperty;
+typedef Property<std::string> StringProperty;
+
+template <class ArchiveT>
+void serialize_properties_helper(ArchiveT &ar, const std::vector<PropertyBase*> &props, unsigned int version) {
+    for (auto* prop : props) prop->serialize(ar, version);
+}
+
+template <class ArchiveT>
+void Object::serialize(ArchiveT &ar, const unsigned int version) {
+    ar & boost::serialization::make_nvp("InstanceName", this->GetInstanceName());
+    serialize_properties_helper(ar, this->GetProperties(), version);
+}
+
 } // namespace uLib
 
 #endif // U_CORE_PROPERTY_H

src/Core/Types.h

@@ -182,7 +182,7 @@ typedef bool Bool_t; // Boolean (0=false, 1=true) (bool)
                                                                                \
 public:                                                                        \
   typedef type_info::BaseClass BaseClass;                                      \
-  virtual const char *type_name() const override { return type_info::name; }            \
+  virtual const char *type_name() const override { return type_info::name; }    \
   /**/
 
 /**

src/Core/testing/VectorTest.cpp

@@ -38,7 +38,7 @@ inline const unsigned long VectorSplice(const _Tp &_it, const _Tp &_end,
 
   _Tp it = _it;
   _Tp end = _end - 1;
-  for (it; it != end;) {
+  for (; it != end;) {
     if (_comp(*it, value))
       it++;
     else if (_comp(*end, value)) {

src/HEP/Geant/EmitterPrimary.hh

@@ -32,7 +32,7 @@ class EmitterPrimary : public G4VUserPrimaryGeneratorAction, public AffineTransf
     virtual ~EmitterPrimary();
 
     // Metodo principale chiamato all'inizio di ogni evento
-    virtual void GeneratePrimaries(G4Event*);         
+    virtual void GeneratePrimaries(G4Event*) override;
 
     virtual void Updated() override { ULIB_SIGNAL_EMIT(EmitterPrimary::Updated); }
     
@@ -51,7 +51,7 @@ class SkyPlaneEmitterPrimary : public EmitterPrimary
     SkyPlaneEmitterPrimary();
     virtual ~SkyPlaneEmitterPrimary();
 
-    virtual void GeneratePrimaries(G4Event*);
+    virtual void GeneratePrimaries(G4Event*) override;
 
     void SetPlane(const uLib::Vector3f& p0, const uLib::Vector3f& normal);
     void SetSkySize(const uLib::Vector2f& size);
@@ -72,7 +72,7 @@ class CylinderEmitterPrimary : public EmitterPrimary
     CylinderEmitterPrimary();
     virtual ~CylinderEmitterPrimary();
 
-    virtual void GeneratePrimaries(G4Event*);
+    virtual void GeneratePrimaries(G4Event*) override;
 
     void SetRadius(float r);
     float GetRadius() const { return m_Radius; }
@@ -101,7 +101,7 @@ class QuadMeshEmitterPrimary : public EmitterPrimary
     virtual ~QuadMeshEmitterPrimary();
 
     // Metodo principale chiamato all'inizio di ogni evento
-    virtual void GeneratePrimaries(G4Event*);
+    virtual void GeneratePrimaries(G4Event*) override;
 
     void SetMesh(uLib::QuadMesh* mesh);
     

src/HEP/Geant/Solid.cpp

@@ -105,6 +105,8 @@ void Solid::SetTransform(Matrix4f transform) {
         m_Physical->SetTranslation(*m_Position);
         m_Physical->SetRotation(m_Rotation);
     }
+
+    std::cout << "Solid " << GetName() << " position: " << pos << " rotation: " << m << std::endl;
 }
 
 void Solid::SetParent(Solid *parent) {

src/Math/Units.h

@@ -13,41 +13,41 @@ namespace uLib {
   using namespace CLHEP;
 
   inline namespace literals {
-    constexpr double operator"" _m(long double v) { return static_cast<double>(v) * CLHEP::meter; }
+    constexpr double operator""_m(long double v) { return static_cast<double>(v) * CLHEP::meter; }
-    constexpr double operator"" _cm(long double v) { return static_cast<double>(v) * CLHEP::centimeter; }
+    constexpr double operator""_cm(long double v) { return static_cast<double>(v) * CLHEP::centimeter; }
-    constexpr double operator"" _mm(long double v) { return static_cast<double>(v) * CLHEP::millimeter; }
+    constexpr double operator""_mm(long double v) { return static_cast<double>(v) * CLHEP::millimeter; }
-    constexpr double operator"" _um(long double v) { return static_cast<double>(v) * CLHEP::micrometer; }
+    constexpr double operator""_um(long double v) { return static_cast<double>(v) * CLHEP::micrometer; }
-    constexpr double operator"" _nm(long double v) { return static_cast<double>(v) * CLHEP::nanometer; }
+    constexpr double operator""_nm(long double v) { return static_cast<double>(v) * CLHEP::nanometer; }
-    constexpr double operator"" _km(long double v) { return static_cast<double>(v) * CLHEP::kilometer; }
+    constexpr double operator""_km(long double v) { return static_cast<double>(v) * CLHEP::kilometer; }
 
-    constexpr double operator"" _m(unsigned long long v) { return static_cast<double>(v) * CLHEP::meter; }
+    constexpr double operator""_m(unsigned long long v) { return static_cast<double>(v) * CLHEP::meter; }
-    constexpr double operator"" _cm(unsigned long long v) { return static_cast<double>(v) * CLHEP::centimeter; }
+    constexpr double operator""_cm(unsigned long long v) { return static_cast<double>(v) * CLHEP::centimeter; }
-    constexpr double operator"" _mm(unsigned long long v) { return static_cast<double>(v) * CLHEP::millimeter; }
+    constexpr double operator""_mm(unsigned long long v) { return static_cast<double>(v) * CLHEP::millimeter; }
-    constexpr double operator"" _um(unsigned long long v) { return static_cast<double>(v) * CLHEP::micrometer; }
+    constexpr double operator""_um(unsigned long long v) { return static_cast<double>(v) * CLHEP::micrometer; }
-    constexpr double operator"" _nm(unsigned long long v) { return static_cast<double>(v) * CLHEP::nanometer; }
+    constexpr double operator""_nm(unsigned long long v) { return static_cast<double>(v) * CLHEP::nanometer; }
-    constexpr double operator"" _km(unsigned long long v) { return static_cast<double>(v) * CLHEP::kilometer; }
+    constexpr double operator""_km(unsigned long long v) { return static_cast<double>(v) * CLHEP::kilometer; }
 
-    constexpr double operator"" _deg(long double v) { return static_cast<double>(v) * CLHEP::degree; }
+    constexpr double operator""_deg(long double v) { return static_cast<double>(v) * CLHEP::degree; }
-    constexpr double operator"" _rad(long double v) { return static_cast<double>(v) * CLHEP::radian; }
+    constexpr double operator""_rad(long double v) { return static_cast<double>(v) * CLHEP::radian; }
-    constexpr double operator"" _deg(unsigned long long v) { return static_cast<double>(v) * CLHEP::degree; }
+    constexpr double operator""_deg(unsigned long long v) { return static_cast<double>(v) * CLHEP::degree; }
-    constexpr double operator"" _rad(unsigned long long v) { return static_cast<double>(v) * CLHEP::radian; }
+    constexpr double operator""_rad(unsigned long long v) { return static_cast<double>(v) * CLHEP::radian; }
 
-    constexpr double operator"" _ns(long double v) { return static_cast<double>(v) * CLHEP::nanosecond; }
+    constexpr double operator""_ns(long double v) { return static_cast<double>(v) * CLHEP::nanosecond; }
-    constexpr double operator"" _s(long double v) { return static_cast<double>(v) * CLHEP::second; }
+    constexpr double operator""_s(long double v) { return static_cast<double>(v) * CLHEP::second; }
-    constexpr double operator"" _ms(long double v) { return static_cast<double>(v) * CLHEP::millisecond; }
+    constexpr double operator""_ms(long double v) { return static_cast<double>(v) * CLHEP::millisecond; }
-    constexpr double operator"" _ns(unsigned long long v) { return static_cast<double>(v) * CLHEP::nanosecond; }
+    constexpr double operator""_ns(unsigned long long v) { return static_cast<double>(v) * CLHEP::nanosecond; }
-    constexpr double operator"" _s(unsigned long long v) { return static_cast<double>(v) * CLHEP::second; }
+    constexpr double operator""_s(unsigned long long v) { return static_cast<double>(v) * CLHEP::second; }
-    constexpr double operator"" _ms(unsigned long long v) { return static_cast<double>(v) * CLHEP::millisecond; }
+    constexpr double operator""_ms(unsigned long long v) { return static_cast<double>(v) * CLHEP::millisecond; }
 
-    constexpr double operator"" _MeV(long double v) { return static_cast<double>(v) * CLHEP::megaelectronvolt; }
+    constexpr double operator""_MeV(long double v) { return static_cast<double>(v) * CLHEP::megaelectronvolt; }
-    constexpr double operator"" _eV(long double v) { return static_cast<double>(v) * CLHEP::electronvolt; }
+    constexpr double operator""_eV(long double v) { return static_cast<double>(v) * CLHEP::electronvolt; }
-    constexpr double operator"" _keV(long double v) { return static_cast<double>(v) * CLHEP::kiloelectronvolt; }
+    constexpr double operator""_keV(long double v) { return static_cast<double>(v) * CLHEP::kiloelectronvolt; }
-    constexpr double operator"" _GeV(long double v) { return static_cast<double>(v) * CLHEP::gigaelectronvolt; }
+    constexpr double operator""_GeV(long double v) { return static_cast<double>(v) * CLHEP::gigaelectronvolt; }
-    constexpr double operator"" _TeV(long double v) { return static_cast<double>(v) * CLHEP::teraelectronvolt; }
+    constexpr double operator""_TeV(long double v) { return static_cast<double>(v) * CLHEP::teraelectronvolt; }
-    constexpr double operator"" _MeV(unsigned long long v) { return static_cast<double>(v) * CLHEP::megaelectronvolt; }
+    constexpr double operator""_MeV(unsigned long long v) { return static_cast<double>(v) * CLHEP::megaelectronvolt; }
-    constexpr double operator"" _eV(unsigned long long v) { return static_cast<double>(v) * CLHEP::electronvolt; }
+    constexpr double operator""_eV(unsigned long long v) { return static_cast<double>(v) * CLHEP::electronvolt; }
-    constexpr double operator"" _keV(unsigned long long v) { return static_cast<double>(v) * CLHEP::kiloelectronvolt; }
+    constexpr double operator""_keV(unsigned long long v) { return static_cast<double>(v) * CLHEP::kiloelectronvolt; }
-    constexpr double operator"" _GeV(unsigned long long v) { return static_cast<double>(v) * CLHEP::gigaelectronvolt; }
+    constexpr double operator""_GeV(unsigned long long v) { return static_cast<double>(v) * CLHEP::gigaelectronvolt; }
   }
 }
 

src/Vtk/CMakeLists.txt

@@ -37,6 +37,8 @@ list(APPEND HEADERS ${HEP_GEANT_HEADERS})
 set(LIBRARIES Eigen3::Eigen
               ${ROOT_LIBRARIES}
               ${VTK_LIBRARIES}
+              VTK::ImagingHybrid
+              VTK::ImagingSources
               ${PACKAGE_LIBPREFIX}Math
               ${PACKAGE_LIBPREFIX}Detectors
               ${PACKAGE_LIBPREFIX}Geant)
@@ -56,7 +58,7 @@ set_target_properties(${libname} PROPERTIES
                        AUTOMOC ON
                        AUTOUIC ON
                        AUTORCC ON)
-target_link_libraries(${libname} ${LIBRARIES} Qt6::Widgets)
+target_link_libraries(${libname} PUBLIC ${LIBRARIES} Qt6::Widgets)
 
 install(TARGETS ${libname}
         EXPORT "uLibTargets"

src/Vtk/HEP/Geant/testing/vtkSolidsTest.cpp

@@ -61,8 +61,8 @@ int main(int argc, char** argv) {
     vtkTess.AddToViewer(viewer);
 
     // Color them differently
-    vtkActor::SafeDownCast(vtkBox.GetProp())->GetProperty()->SetColor(0.8, 0.2, 0.2); // Redish box
+    vtkBox.SetColor(0.8, 0.2, 0.2);  // Redish box
-    vtkActor::SafeDownCast(vtkTess.GetProp())->GetProperty()->SetColor(0.2, 0.8, 0.2); // Greenish tess
+    vtkTess.SetColor(0.2, 0.8, 0.2); // Greenish tess
     
     // Position tessellated solid away from box
     Matrix4f trans = Matrix4f::Identity();

src/Vtk/HEP/Geant/vtkBoxSolid.cpp

@@ -10,47 +10,74 @@
 //////////////////////////////////////////////////////////////////////////////*/
 
 #include "vtkBoxSolid.h"
+#include "Core/Monitor.h"
 #include <vtkCubeSource.h>
 #include <vtkPolyDataMapper.h>
 #include <vtkActor.h>
+#include <vtkProperty.h>
+#include <vtkAssembly.h>
+#include <vtkTransform.h>
+#include <vtkMatrix4x4.h>
+#include <Geant4/G4VPhysicalVolume.hh>
+#include "Vtk/Math/vtkDense.h"
 
 namespace uLib {
 namespace Vtk {
 
 vtkBoxSolid::vtkBoxSolid(Geant::BoxSolid *content)
     : vtkGeantSolid(content), m_BoxContent(content) {
-  // Re-run Update for box-specific pipe
+  this->InstallPipe();
+ 
+  // Connect the model's Updated event to updateTransform to ensure VTK sync
+  m_UpdateConnection = Object::connect(m_BoxContent, &uLib::Object::Updated, this, &vtkBoxSolid::Update);
+
+  // Initial sync
   this->Update();
 }
 
-vtkBoxSolid::~vtkBoxSolid() {}
+vtkBoxSolid::~vtkBoxSolid() {
+}
 
 void vtkBoxSolid::Update() {
+  ConnectionBlock blocker(m_UpdateConnection);
   this->UpdateGeometry();
-  this->UpdateTransform();
+  // Ensure base Puppet properties (color, opacity, etc) and transform are applied
+  this->Puppet::Update();
+}
+
+void vtkBoxSolid::SyncFromVtk() {
+  this->Puppet::SyncFromVtk();
+  if (auto* proxy = vtkProp3D::SafeDownCast(this->GetProxyProp())) {
+      if (vtkMatrix4x4* mat = proxy->GetUserMatrix()) {
+          m_BoxContent->SetTransform(VtkToMatrix4f(mat));
+      }
+  }
 }
 
 void vtkBoxSolid::UpdateGeometry() {
-  if (!m_BoxContent || !m_BoxContent->GetObject()) {
+  // Sync geometry from G4VSolid provided by vtkGeantSolid (tessellation)
-    // Fallback to base tessellation if no model object
+  vtkGeantSolid::UpdateGeometry();
-    vtkGeantSolid::UpdateGeometry();
+}
-    return;
+
+void vtkBoxSolid::UpdateTransform() {
+  // Take transform from Puppet base (which uses GetContent() -> ContainerBox TRS)
+  this->Puppet::Update();
+}
+
+void vtkBoxSolid::serialize_display(uLib::Archive::display_properties_archive &ar,
+                                    const unsigned int version) {
+  // Expose Geant solid properties and underlying Box/TRS properties
+  this->Puppet::serialize_display(ar, version);
+  if (m_BoxContent) {
+    ar & NVP("Box", *m_BoxContent);
+    if (m_BoxContent->GetObject()) {
+      ar & NVP("Container", *m_BoxContent->GetObject());
+    }
   }
+}
 
-  // Use the underlying ContainerBox for precise geometry
+void vtkBoxSolid::InstallPipe() {
-  Vector3f size = m_BoxContent->GetObject()->GetSize();
+  vtkGeantSolid::InstallPipe();
-  
-  vtkNew<vtkCubeSource> cube;
-  cube->SetXLength(size(0));
-  cube->SetYLength(size(1));
-  cube->SetZLength(size(2));
-  cube->Update();
-
-  vtkPolyData *poly = GetPolyData();
-  if (poly) {
-      poly->ShallowCopy(cube->GetOutput());
-      poly->Modified();
-  }
 }
 
 } // namespace Vtk

src/Vtk/HEP/Geant/vtkBoxSolid.h

@@ -26,8 +26,15 @@
 #ifndef U_VTKBOXSOLID_H
 #define U_VTKBOXSOLID_H
 
+#include "Core/Types.h"
+#include "Core/Property.h"
+#include "Core/Serializable.h"
+
 #include "vtkGeantSolid.h"
 
+class vtkCubeSource;
+class vtkActor;
+
 namespace uLib {
 namespace Vtk {
 
@@ -35,15 +42,36 @@ namespace Vtk {
  * @brief VTK Puppet for visualizing a Geant::BoxSolid.
  */
 class vtkBoxSolid : public vtkGeantSolid {
+  uLibTypeMacro(vtkBoxSolid, uLib::Vtk::vtkGeantSolid)
+
 public:
   vtkBoxSolid(Geant::BoxSolid *content);
   virtual ~vtkBoxSolid();
 
   virtual void Update() override;
   virtual void UpdateGeometry() override;
+  virtual void UpdateTransform() override;
+  virtual void SyncFromVtk() override;
+
+  virtual uLib::Object *GetContent() const override { 
+    return m_BoxContent ? (::uLib::Object*)m_BoxContent->GetObject() : nullptr; 
+  }
+
+  virtual void serialize_display(uLib::Archive::display_properties_archive &ar,
+                                 const unsigned int version = 0) override;
+
+  template <typename Ar>
+  void serialize(Ar &ar, const unsigned int version) {
+    ar & NVP("BoxSolid", *m_BoxContent);
+  }
 
 protected:
+  virtual void InstallPipe() override;
+
   Geant::BoxSolid *m_BoxContent;
+  uLib::Connection m_UpdateConnection;
+
+  ULIB_DECLARE_PROPERTIES(vtkBoxSolid)
 };
 
 } // namespace Vtk

src/Vtk/HEP/Geant/vtkGeantEvent.h

@@ -43,7 +43,7 @@ public:
 
   virtual class vtkPolyData *GetPolyData() const override;
 
-  virtual void Update();
+  virtual void Update() override;
 
 protected:
   virtual void InstallPipe();

src/Vtk/HEP/Geant/vtkGeantSolid.cpp

@@ -82,6 +82,7 @@ void vtkGeantSolid::UpdateGeometry() {
   int nVertices = polyhedron->GetNoVertices();
   for (int i = 1; i <= nVertices; ++i) {
     G4Point3D vtx = polyhedron->GetVertex(i);
+    
     points->InsertNextPoint(vtx.x(), vtx.y(), vtx.z());
   }
 

src/Vtk/Math/vtkContainerBox.cpp

@@ -54,6 +54,7 @@ struct ContainerBoxData {
   vtkSmartPointer<vtkMatrix4x4> m_Affine;
   uLib::Connection              m_UpdateSignal;
 
+
   ContainerBoxData() : m_Cube(vtkSmartPointer<vtkActor>::New()), 
                        m_Axes(vtkSmartPointer<vtkActor>::New()), 
                        m_VtkAsm(vtkSmartPointer<vtkAssembly>::New()),

src/Vtk/Math/vtkContainerBox.h

@@ -38,13 +38,15 @@ namespace Vtk {
 struct ContainerBoxData;
 
 class vtkContainerBox : public Puppet, public Polydata {
+  
+  uLibTypeMacro(vtkContainerBox, Puppet, Polydata)
   typedef ContainerBox Content;
   
 public:
   vtkContainerBox(Content *content);
   ~vtkContainerBox();
 
-  virtual class vtkPolyData *GetPolyData() const;
+  virtual class vtkPolyData *GetPolyData() const override;
 
   /**
    * @brief Updates the VTK representation from the internal state.
@@ -62,9 +64,9 @@ protected:
   virtual void InstallPipe();
 
   struct ContainerBoxData *d;
-  Content                 *m_Content;
+  ContainerBox            *m_Content;
-  bool                     m_BlockUpdate = false;
   
+  ULIB_DECLARE_PROPERTIES(vtkContainerBox)
 };
 
 } // namespace Vtk

src/Vtk/Math/vtkCylinder.cpp

@@ -72,8 +72,8 @@ void vtkCylinder::Update() {
 
           // Apply Axis alignment
           int axis = m_Content->GetAxis();
-          if (axis == 0) alignment->RotateZ(-90); // Y -> X
+          if (axis == 0) alignment->RotateZ(-90);      // Y -> X
-          else if (axis == 1) ;                   // Y -> Y (identity)
+          else if (axis == 1) ;                        // Y -> Y (identity)
           else if (axis == 2) alignment->RotateX(90);  // Y -> Z
       }
 
@@ -93,18 +93,10 @@ void vtkCylinder::SyncFromVtk() {
 
   // VTK -> Model: Extract new world TRS from proxy
   vtkMatrix4x4* rootMat = root->GetUserMatrix();  
-  if (rootMat) {
-      std::cout << "[vtkCylinder::SyncFromVtk] Read Proxy UserMatrix:" << std::endl;
-      rootMat->Print(std::cout);
-  }
-  
   Matrix4f vtkWorld = VtkToMatrix4f(rootMat);
   
   // Directly sync model from the world matrix
   m_Content->FromMatrix(vtkWorld);
-
-  std::cout << "[vtkCylinder::SyncFromVtk] New Model WorldMatrix:" << std::endl << m_Content->GetWorldMatrix() << std::endl;
-
   m_Content->Updated(); 
 }
 

src/Vtk/uLibVtkInterface.cxx

@@ -183,7 +183,7 @@ public:
             vtkPolyData* polydata = nullptr;
             if (vtkActor *actor = vtkActor::SafeDownCast(m_Prop)) {
                 if (actor->GetMapper()) {
-                    polydata = vtkPolyData::SafeDownCast(actor->GetMapper()->GetDataSetInput());
+                    polydata = vtkPolyData::SafeDownCast(actor->GetMapper()->GetInput());
                 }
             } else if (vtkAssembly *asm_p = vtkAssembly::SafeDownCast(m_Prop)) {
                 vtkPropCollection *parts = asm_p->GetParts();
@@ -192,7 +192,7 @@ public:
                     for (int i = 0; i < parts->GetNumberOfItems(); ++i) {
                         vtkActor *a = vtkActor::SafeDownCast(parts->GetNextProp());
                         if (a && a->GetMapper()) {
-                            polydata = vtkPolyData::SafeDownCast(a->GetMapper()->GetDataSetInput());
+                            polydata = vtkPolyData::SafeDownCast(a->GetMapper()->GetInput());
                             if (polydata) break;
                         }
                     }
@@ -560,28 +560,55 @@ bool Puppet::IsSelected() const
     return pd->m_Selected;
 }
 
+void Puppet::ApplyPuppetTransform(vtkProp3D* prop)
+{
+    if (!prop) return;
+    if (auto* content = this->GetContent()) {
+        if (auto* tr = dynamic_cast<uLib::TRS*>(content)) {
+            vtkNew<vtkMatrix4x4> m;
+            Matrix4fToVtk(tr->GetMatrix(), m);
+            prop->SetUserMatrix(m);
+            prop->Modified();
+        }
+    }
+}
+
+void Puppet::SyncFromVtk()
+{
+    if (auto* content = this->GetContent()) {
+        if (auto* tr = dynamic_cast<uLib::TRS*>(content)) {
+            if (auto* proxy = this->GetProxyProp()) {
+                if (vtkMatrix4x4* mat = proxy->GetUserMatrix()) {
+                    tr->FromMatrix(VtkToMatrix4f(mat));
+                    content->Updated();
+                }
+            }
+        }
+    }
+}
+
 void Puppet::Update()
 {
-    // Derived classes should have updated the transform if they override Update()
+    // Apply content transform via virtual GetProp() / ApplyPuppetTransform(),
-    // or we can apply base transform if it's default:
+    // so all derived classes benefit without duplicating the matrix code.
-    // pd->ApplyTransform(pd->m_Prop);
+    this->ApplyPuppetTransform(vtkProp3D::SafeDownCast(this->GetProp()));
 
-    pd->ApplyAppearance(pd->m_Prop);
+    // Use virtual GetProp() for appearance so overriders (e.g. vtkVoxImage)
+    // that never call SetProp() are handled correctly.
+    pd->ApplyAppearance(this->GetProp());
 
     if (pd->m_Selected) {
         pd->UpdateHighlight();
     }
 
-    if (pd->m_Prop) {
+    if (auto* prop = this->GetProp()) {
-        if (pd->m_ShowBoundingBox) {
+        if (pd->m_ShowBoundingBox && pd->m_OutlineSource) {
-            double* bounds = pd->m_Prop->GetBounds();
+            double* bounds = prop->GetBounds();
             pd->m_OutlineSource->SetBounds(bounds);
             pd->m_OutlineSource->Update();
         }
-        
+        if (pd->m_ShowScaleMeasures && pd->m_CubeAxesActor) {
-        if (pd->m_ShowScaleMeasures) {
+            pd->m_CubeAxesActor->SetBounds(prop->GetBounds());
-            double* bounds = pd->m_Prop->GetBounds();
-            pd->m_CubeAxesActor->SetBounds(bounds);
         }
     }
 
@@ -598,6 +625,7 @@ void Puppet::Update()
 }
 
 
+
 void Puppet::ConnectInteractor(vtkRenderWindowInteractor *interactor)
 {
 }

src/Vtk/uLibVtkInterface.h

@@ -35,6 +35,8 @@
 #include <iomanip>
 #include <ostream>
 #include <vector>
+#include <set>
+#include <boost/type_traits/is_base_of.hpp>
 
 // vtk classes forward declaration //
 class vtkProp;
@@ -106,7 +108,7 @@ uLibTypeMacro(Puppet, uLib::Object)
    * This method should be called when the VTK representation has been modified
    * (e.g., via a gizmo) and the changes need to be pushed back to the model.
    */
-  virtual void SyncFromVtk() {}
+  virtual void SyncFromVtk();
 
   enum Representation {
     Points = 0,
@@ -127,10 +129,10 @@ uLibTypeMacro(Puppet, uLib::Object)
 
   vtkRendererCollection *GetRenderers() const;
 
-  const std::vector<uLib::PropertyBase *> &GetDisplayProperties() const {
+  const std::vector<uLib::PropertyBase *> &GetDisplayProperties() const override {
     return m_DisplayProperties;
   }
-  void RegisterDisplayProperty(uLib::PropertyBase *prop) {
+  void RegisterDisplayProperty(uLib::PropertyBase *prop) override {
     m_DisplayProperties.push_back(prop);
   }
 
@@ -149,6 +151,7 @@ protected:
 
   void ApplyAppearance(vtkProp *prop);
   void ApplyTransform(vtkProp3D *p3d);
+  void ApplyPuppetTransform(vtkProp3D *p3d);
 
   std::vector<uLib::PropertyBase *> m_DisplayProperties;
   mutable uLib::RecursiveMutex m_UpdateMutex;
@@ -179,10 +182,17 @@ class display_properties_archive
     : public boost::archive::detail::common_oarchive<
           display_properties_archive> {
 public:
-  display_properties_archive(Vtk::Puppet *puppet)
+  friend class boost::archive::detail::interface_oarchive<display_properties_archive>;
+  friend class boost::archive::save_access;
+
+  using boost::archive::detail::common_oarchive<display_properties_archive>::save_override;
+  display_properties_archive(Vtk::Puppet *p)
       : boost::archive::detail::common_oarchive<display_properties_archive>(
             boost::archive::no_header),
-        m_Puppet(puppet) {}
+        m_Puppet(p) {
+    if (p)
+      m_Visited.insert(dynamic_cast<const void *>(p));
+  }
 
   std::string GetCurrentGroup() const {
     std::string group;
@@ -234,6 +244,24 @@ public:
       m_GroupStack.pop_back();
   }
 
+  // Follow pointers to discover properties in child objects
+  template<class T>
+  void save_override(T * const & t) {
+      if (!t) return;
+      this->save_pointer_helper(t, typename boost::is_base_of<uLib::Object, T>::type());
+  }
+
+  template<class T>
+  void save_pointer_helper(T* t, boost::mpl::true_) {
+      const void* ptr = dynamic_cast<const void*>(t);
+      if (m_Visited.find(ptr) != m_Visited.end()) return;
+      m_Visited.insert(ptr);
+      this->save_override(*t);
+  }
+
+  template<class T>
+  void save_pointer_helper(T* t, boost::mpl::false_) {}
+
   // Recursion for nested classes, ignore primitives
   template <class T> void save_override(const T &t) {
     this->save_helper(t, typename boost::is_class<T>::type());
@@ -243,6 +271,8 @@ public:
     boost::serialization::serialize_adl(*this, const_cast<T &>(t), 0);
   }
 
+  void save_helper(const std::string &t, boost::mpl::true_) {}
+
   template <class T> void save_helper(const T &t, boost::mpl::false_) {}
 
   void save_override(const boost::archive::object_id_type &t) {}
@@ -254,9 +284,25 @@ public:
   void save_override(const boost::archive::class_name_type &t) {}
   void save_override(const boost::archive::tracking_type &t) {}
 
+  // Called by Property<T>::serialize() and EnumProperty::serialize() to
+  // directly register an existing property object as a display property.
+  void register_property(uLib::PropertyBase &p) {
+    if (m_Puppet) {
+      m_Puppet->RegisterDisplayProperty(&p);
+      Vtk::Puppet *puppet = m_Puppet;
+      uLib::Object::connect(&p, &uLib::PropertyBase::Updated,
+                            [puppet]() { puppet->Update(); });
+    }
+  }
+
+  void register_enum_property(uLib::EnumProperty &p) {
+    register_property(p);
+  }
+
 private:
   Vtk::Puppet *m_Puppet;
   std::vector<std::string> m_GroupStack;
+  std::set<const void *> m_Visited;
 };
 
 } // namespace Archive