diff --git a/Device/Beam/Beam.cpp b/Device/Beam/Beam.cpp
index 2618de0965ff275a2e3cd8e4908229f005bdd5d4..a7ed2a4b0665b5a629c389fc6f929f7734a3202a 100644
--- a/Device/Beam/Beam.cpp
+++ b/Device/Beam/Beam.cpp
@@ -42,23 +42,6 @@ Beam::Beam(const Beam& other)
         m_shape_factor.reset(other.m_shape_factor->clone());
 }
 
-Beam& Beam::operator=(const Beam& other)
-{
-    if (&other == this)
-        return *this;
-    m_intensity = other.m_intensity;
-    m_wavelength = other.m_wavelength;
-    // m_direction = other.m_direction;
-    m_alpha = other.m_alpha;
-    m_phi = other.m_phi;
-    m_beamPolarization = other.m_beamPolarization;
-    if (other.m_shape_factor)
-        m_shape_factor.reset(other.m_shape_factor->clone());
-    else
-        m_shape_factor.release();
-    return *this;
-}
-
 Beam::~Beam() = default;
 
 Beam* Beam::clone() const
diff --git a/Device/Beam/Beam.h b/Device/Beam/Beam.h
index 1a71fbccc7cee127d1f5278f394fc6b7ba09b8a8..3e07ad28ad584374fb729010d47cfb55e55c0c01 100644
--- a/Device/Beam/Beam.h
+++ b/Device/Beam/Beam.h
@@ -28,7 +28,6 @@ class Beam : public INode {
 public:
     Beam(double intensity, double wavelength, const Direction& direction);
     Beam(const Beam& other);
-    Beam& operator=(const Beam& other);
     ~Beam() override;
 
     Beam* clone() const;
diff --git a/Device/Histo/SimulationResult.cpp b/Device/Histo/SimulationResult.cpp
index 244442f0d84252c1e445af3e89697179bb1616ab..116a75cc0cb23b7900f997961a211adae1abdb65 100644
--- a/Device/Histo/SimulationResult.cpp
+++ b/Device/Histo/SimulationResult.cpp
@@ -50,15 +50,6 @@ SimulationResult::SimulationResult(SimulationResult&& other)
 {
 }
 
-SimulationResult& SimulationResult::operator=(const SimulationResult& other)
-{
-    if (&other == this)
-        return *this;
-    m_data.reset(other.m_data->clone());
-    m_coordsys.reset(other.m_coordsys->clone());
-    return *this;
-}
-
 SimulationResult& SimulationResult::operator=(SimulationResult&& other)
 {
     m_data = std::move(other.m_data);
diff --git a/Device/Histo/SimulationResult.h b/Device/Histo/SimulationResult.h
index 6d6e434c4e3ec2572fbe13915e118f64d5cbc818..bc12d3f79d73d8b77ba862a1cf69dcc23b07a843 100644
--- a/Device/Histo/SimulationResult.h
+++ b/Device/Histo/SimulationResult.h
@@ -38,7 +38,6 @@ public:
     SimulationResult(const SimulationResult& other);
     SimulationResult(SimulationResult&& other);
 
-    SimulationResult& operator=(const SimulationResult& other);
     SimulationResult& operator=(SimulationResult&& other);
 
     std::vector<double> flatVector(Coords units = Coords::UNDEFINED) const;
diff --git a/Device/Mask/DetectorMask.cpp b/Device/Mask/DetectorMask.cpp
index 8dcb6debe69cbc4c5bfbd392f86ab2bcbc367cfd..2c6cc2e82f6c105d482e979b496695ad7e7866f9 100644
--- a/Device/Mask/DetectorMask.cpp
+++ b/Device/Mask/DetectorMask.cpp
@@ -50,16 +50,6 @@ DetectorMask::DetectorMask(const DetectorMask& other)
 {
 }
 
-DetectorMask& DetectorMask::operator=(const DetectorMask& other)
-{
-    if (this != &other) {
-        m_stack = other.m_stack;
-        m_masked = other.m_masked;
-        m_number_of_masked_channels = other.m_number_of_masked_channels;
-    }
-    return *this;
-}
-
 void DetectorMask::addMask(const IShape2D& shape, bool mask_value)
 {
     m_stack.emplace_back(new MaskPattern(shape.clone(), mask_value));
diff --git a/Device/Mask/DetectorMask.h b/Device/Mask/DetectorMask.h
index 4ee51de2e30aa293087c83f47b212089446f3967..51f7013c7ee6c930a8b6aa92aaa12159dc19a832 100644
--- a/Device/Mask/DetectorMask.h
+++ b/Device/Mask/DetectorMask.h
@@ -44,7 +44,6 @@ public:
     DetectorMask(const IAxis& xAxis, const IAxis& yAxis);
     ~DetectorMask();
     DetectorMask(const DetectorMask& other);
-    DetectorMask& operator=(const DetectorMask& other);
 
     //! Add mask to the stack of detector masks.
     //! The value "true" means that the area will be excluded from the analysis.
diff --git a/Param/Node/INode.h b/Param/Node/INode.h
index fef0e3244fb25b469e54c07e00a0ae3c3a12ed0a..5cd824a12960eea2a06294f53603a5d033b03278 100644
--- a/Param/Node/INode.h
+++ b/Param/Node/INode.h
@@ -49,7 +49,7 @@ public:
     //! Returns all children.
     virtual std::vector<const INode*> nodeChildren() const;
 
-    //! Returns all descendants.
+    //! Returns all descendants. // TODO sep22: rm if it remains UNUSED
     std::vector<const INode*> nodeOffspring() const;
 
     static void requestGt0(std::vector<std::string>& errs, const double& val,
diff --git a/Resample/Slice/ZLimits.cpp b/Resample/Slice/ZLimits.cpp
index eeece937dba194f673869b1e57b2756312668cd9..b99e18eb720054e9effdbe3b1a4ae07bd2a9698c 100644
--- a/Resample/Slice/ZLimits.cpp
+++ b/Resample/Slice/ZLimits.cpp
@@ -34,6 +34,8 @@ bool ZLimits::isFinite() const
     return !std::isinf(m_zmin) && !std::isinf(m_zmax);
 }
 
+// TODO sep22: rm all the following if they remain UNUSED
+
 ZLimits ZLimits::enclosingInterval(const ZLimits& left, const ZLimits& right)
 {
     return {std::min(left.zBottom(), right.zBottom()), std::max(left.zTop(), right.zTop())};
diff --git a/auto/Wrap/doxygenParam.i b/auto/Wrap/doxygenParam.i
index f64da0e9c20c85ab08e56958ab6171b9c5585da8..8c41b0c548016513678a43709a5027a739bfc276 100644
--- a/auto/Wrap/doxygenParam.i
+++ b/auto/Wrap/doxygenParam.i
@@ -505,7 +505,7 @@ Returns all children.
 
 %feature("docstring")  INode::nodeOffspring "std::vector< const INode * > INode::nodeOffspring() const
 INode::nodeOffspring
-Returns all descendants. 
+Returns all descendants. // TODO sep22: rm if it remains UNUSED. 
 ";
 
 %feature("docstring")  INode::jointError "std::string INode::jointError(const std::vector< std::string > errs) const
@@ -658,6 +658,16 @@ C++ includes: ParameterSample.h
 ParameterSample::ParameterSample";
 
 
+// File: classPDouble.xml
+%feature("docstring") PDouble "";
+
+%feature("docstring")  PDouble::PDouble "PDouble::PDouble(double val)
+PDouble::PDouble";
+
+%feature("docstring")  PDouble::PDouble "PDouble::PDouble(double *ptr)
+PDouble::PDouble";
+
+
 // File: classRangedDistributionCosine.xml
 %feature("docstring") RangedDistributionCosine "
 
@@ -866,6 +876,9 @@ NodeUtils::AllDescendantsOfType";
 // File: NodeUtils_8h.xml
 
 
+// File: PDouble_8h.xml
+
+
 // File: dir_dabc42ad3745509abd3a496944bb880e.xml
 
 
diff --git a/auto/Wrap/libBornAgainParam.py b/auto/Wrap/libBornAgainParam.py
index 6b9d0f6710dbdb5e142c7e81137512729f45d994..4f83de4d3bb24edd6c77ad768dab43de5f401cf8 100644
--- a/auto/Wrap/libBornAgainParam.py
+++ b/auto/Wrap/libBornAgainParam.py
@@ -2108,7 +2108,7 @@ class INode(object):
         nodeOffspring(INode self) -> swig_dummy_type_const_inode_vector
         std::vector< const INode * > INode::nodeOffspring() const
         INode::nodeOffspring
-        Returns all descendants. 
+        Returns all descendants. // TODO sep22: rm if it remains UNUSED. 
 
         """
         return _libBornAgainParam.INode_nodeOffspring(self)
diff --git a/auto/Wrap/libBornAgainParam_wrap.cpp b/auto/Wrap/libBornAgainParam_wrap.cpp
index 78d8cbc24852206ddf39064e2f3cee38894625cd..af1b47658813f9af125fe67e9bb65426e712f104 100644
--- a/auto/Wrap/libBornAgainParam_wrap.cpp
+++ b/auto/Wrap/libBornAgainParam_wrap.cpp
@@ -34992,7 +34992,7 @@ static PyMethodDef SwigMethods[] = {
 		"INode_nodeOffspring(INode self) -> swig_dummy_type_const_inode_vector\n"
 		"std::vector< const INode * > INode::nodeOffspring() const\n"
 		"INode::nodeOffspring\n"
-		"Returns all descendants. \n"
+		"Returns all descendants. // TODO sep22: rm if it remains UNUSED. \n"
 		"\n"
 		""},
 	 { "INode_requestGt0", _wrap_INode_requestGt0, METH_VARARGS, "INode_requestGt0(vector_string_t errs, double const & val, std::string const & name)"},