have fitcentroid correctly set position to be invalid when charge = 0. (

#215)
rat-pac · Jan 7, 2025 · 5d7039e · 5d7039e
1 parent df92ceb
commit 5d7039e
Showing 1 changed file with 7 additions and 2 deletions.
diff --git a/src/fit/src/FitCentroidProc.cc b/src/fit/src/FitCentroidProc.cc
@@ -50,9 +50,14 @@ Processor::Result FitCentroidProc::Event(DS::Root *ds, DS::EV *ev) {
     centroid += Qpow * pmtpos;
   }
 
-  centroid *= 1.0 / totalQ;
   DS::FitResult *fit = new DS::FitResult("FitCentroid");
-  fit->SetPosition(centroid);
+  fit->SetEnablePosition(true);
+  if (totalQ) {
+    centroid *= 1.0 / totalQ;
+    fit->SetPosition(centroid);
+  } else {
+    fit->SetValidPosition(false);
+  }
   ev->AddFitResult(fit);
 
   return Processor::OK;