Record recombinations at run time in ExpManager_7

src/libaevol/7/ExpManager_7.cpp

@@ -37,6 +37,7 @@
 #include "Rna_7.h"
 #include "List_Metadata.h"
 #include "7/fuzzy/FuzzyFactory_7.h"
+#include "mutation/Recombination.h"
 
 #include <algorithm>
 #include <err.h>
@@ -420,6 +421,15 @@ void ExpManager_7::do_mutation(int indiv_id) {
                                                                   exp_m_->world()->grid(x, y)->mut_prng(),
                                                                   exp_m_->exp_s()->mut_params()->max_tries(),
                                                                   exp_m_->exp_s()->mut_params()->min_align_score());
+
+  if (exp_manager->record_tree()) {
+    Recombination* recomb = new Recombination(current_individuals[indiv_id]->recomb_posA_,
+                                            current_individuals[indiv_id]->recomb_posB_,
+                                            current_individuals[indiv_id]->recomb_score_);
+    exp_manager->tree()->report_by_index(AeTime::time(), indiv_id)->dna_replic_report().add_recomb(recomb);
+    delete recomb;
+  }
+
   dna_length = current_individuals[indiv_id]->annotated_chromosome_[A]->dna_->length();
   #else
   dna_length = previous_individuals[next_generation_reproducer_[indiv_id]]->annotated_chromosome_[A]->dna_->length();