RotationIndexer: Refine detector angles

image_analysis/IndexAndRefine.cpp

@@ -20,6 +20,8 @@ IndexAndRefine::IndexAndRefine(const DiffractionExperiment &x, IndexerThreadPool
         rotation_indexer = std::make_unique<RotationIndexer>(x, *indexer);
 }
 
+#include <iostream>
+
 IndexAndRefine::IndexingOutcome IndexAndRefine::DetermineLatticeAndSymmetry(DataMessage &msg) {
     IndexingOutcome outcome(experiment);
 
@@ -30,8 +32,9 @@ IndexAndRefine::IndexingOutcome IndexAndRefine::DetermineLatticeAndSymmetry(Data
 
             outcome.experiment.BeamX_pxl(result->geom.GetBeamX_pxl())
                                   .BeamY_pxl(result->geom.GetBeamY_pxl())
-                                  .DetectorDistance_mm(result->geom.GetDetectorDistance_mm());
-
+                                  .DetectorDistance_mm(result->geom.GetDetectorDistance_mm())
+            .PoniRot1_rad(result->geom.GetPoniRot1_rad())
+            .PoniRot2_rad(result->geom.GetPoniRot2_rad());
             outcome.symmetry.centering = result->search_result.centering;
             outcome.symmetry.niggli_class = result->search_result.niggli_class;
             outcome.symmetry.crystal_system = result->search_result.system;

image_analysis/RotationIndexer.cpp

@@ -85,6 +85,7 @@ void RotationIndexer::TryIndex() {
             .min_spots = experiment.GetIndexingSettings().GetViableCellMinSpots(),
             .refine_beam_center = true,
             .refine_distance_mm =  false,
+            .refine_detector_angles = true,
             .index_ice_rings = experiment.GetIndexingSettings().GetIndexIceRings(),
             .axis = axis_
         };