diff --git a/include/Ecal/Event/EcalCluster.h b/include/Ecal/Event/EcalCluster.h
index 19950c0..13eebaa 100644
--- a/include/Ecal/Event/EcalCluster.h
+++ b/include/Ecal/Event/EcalCluster.h
@@ -84,6 +84,19 @@ class EcalCluster {
     centroidZ_ = z;
   }
 
+  void setRMSXYZ(double x, double y, double z) {
+    rmsX_ = x;
+    rmsY_ = y;
+    rmsZ_ = z;
+  }
+  void setDXDZ(double x) { DXDZ_=x; }
+
+  void setDYDZ(double x) { DYDZ_=x; }
+
+  void setEDXDZ(double x) { errDXDZ_=x; }
+
+  void setEDYDZ(double x) { errDYDZ_=x; }
+
   double getEnergy() const { return energy_; }
 
   int getNHits() const { return nHits_; }
@@ -94,6 +107,20 @@ class EcalCluster {
 
   double getCentroidZ() const { return centroidZ_; }
 
+  double getRMSX() const { return rmsX_; }
+
+  double getRMSY() const { return rmsY_; }
+
+  double getRMSZ() const { return rmsZ_; }
+
+  double getDXDZ() const { return DXDZ_; }
+
+  double getDYDZ() const { return DYDZ_; }
+
+  double getEDXDZ() const { return errDXDZ_; }
+
+  double getEDYDZ() const { return errDYDZ_; }
+
   const std::vector<unsigned int>& getHitIDs() const { return hitIDs_; }
 
   bool operator<(const EcalCluster& rhs) const {
@@ -107,6 +134,13 @@ class EcalCluster {
   double centroidX_{0};
   double centroidY_{0};
   double centroidZ_{0};
+  double rmsX_{0};
+  double rmsY_{0};
+  double rmsZ_{0};
+  double DXDZ_{0};
+  double DYDZ_{0};
+  double errDXDZ_{0};
+  double errDYDZ_{0};
 
   ClassDef(EcalCluster, 1);
 };
diff --git a/src/Ecal/Event/EcalCluster.cxx b/src/Ecal/Event/EcalCluster.cxx
index 3c8ca92..683caf5 100644
--- a/src/Ecal/Event/EcalCluster.cxx
+++ b/src/Ecal/Event/EcalCluster.cxx
@@ -27,6 +27,13 @@ ClassImp(ldmx::EcalCluster)
     centroidX_ = 0;
     centroidY_ = 0;
     centroidZ_ = 0;
+    rmsX_ = 0;
+    rmsY_ = 0;
+    rmsZ_ = 0;
+    DXDZ_ = 0;
+    DYDZ_ = 0;
+    errDXDZ_ = 0;
+    errDYDZ_ = 0;
   }
 
   void EcalCluster::addHits(const std::vector<const EcalHit *> hitsVec) {