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

rat-pac · Dec 10, 2024 · 288fedd · 288fedd
1 parent 6ec3fd4
commit 288fedd
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;