diff --git a/tofu/geom/_core_optics.py b/tofu/geom/_core_optics.py
index <HASH>..<HASH> 100644
--- a/tofu/geom/_core_optics.py
+++ b/tofu/geom/_core_optics.py
@@ -2801,6 +2801,14 @@ class CrystalBragg(utils.ToFuObject):
             )
 
         # ---------------
+        # refactor pts and lambok
+
+        indok = np.any(lambok, axis=0)
+        pts = pts[:, indok]
+        ptsXYZ = ptsXYZ[:, indok]
+        lambok = lambok[:, indok]
+
+        # ---------------
         # check strict
         if strict is True:
[#<I>] Cryst.get_plasmadomain_at_lamb() now faster and lighter thanks to refactoring, TODO: parallelizing?