Calculate in_fiducial flag correctly in LIController

Harvard-Neutrino · Mar 6, 2024 · ebc150f · ebc150f
1 parent 384002b
commit ebc150f
Showing 1 changed file with 1 addition and 3 deletions.
diff --git a/python/LIController.py b/python/LIController.py
@@ -417,11 +417,9 @@ def SaveEvents(self, filename, fill_tables_at_exit=True, hdf5=True, parquet=True
 
                 if self.fid_vol is not None:
                     pos = _math.Vector3D(datasets["vertex"][-1][-1])
-                    geo_pos = self.detector_model.DetPositionToGeoPosition(pos)
                     dir = _math.Vector3D(datasets["primary_momentum"][-1][-1][1:])
                     dir.normalize()
-                    geo_dir = self.detector_model.DetDirectionToGeoDirection(dir)
-                    datasets["in_fiducial"][-1].append(self.fid_vol.IsInside(geo_pos.get(),geo_dir.get()))
+                    datasets["in_fiducial"][-1].append(self.fid_vol.IsInside(pos,dir))
                 else:
                     datasets["in_fiducial"][-1].append(False)