Use recalculated or best cell

src/dials/util/reindex.py

@@ -6,6 +6,7 @@
 from cctbx import sgtbx
 from rstbx.symmetry.constraints import parameter_reduction
 
+from dials.algorithms.scaling.scaling_library import determine_best_unit_cell
 from dials.algorithms.symmetry.reindex_to_reference import (
     determine_reindex_operator_against_reference,
 )
@@ -74,6 +75,9 @@ def change_of_basis_op_against_reference(
             ), "No 'intensity.sum.value' in reflections"
             table.set_flags(flex.bool(table.size(), True), table.flags.integrated_sum)
     # Make miller array of the datasets
+    cell = determine_best_unit_cell(experiments)
+    for expt in experiments:
+        expt.crystal.set_unit_cell(cell)
     filter_logger = logging.getLogger("dials.util.filter_reflections")
     filter_logger.disabled = True
     test_miller_sets = filtered_arrays_from_experiments_reflections(