Penalize far apart pairs not more than as if they were two single ends

ksahlin · Aug 26, 2023 · e5689c2 · e5689c2
1 parent 4893484
commit e5689c2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/aln.cpp b/src/aln.cpp
@@ -910,7 +910,7 @@ inline void align_PE(
         if (r1_r2 || r2_r1) {
             // Treat a1/a2 as a pair
             float x = std::abs(a1.ref_start - a2.ref_start);
-            combined_score = (double)a1.score + (double)a2.score + log(normal_pdf(x, mu, sigma));
+            combined_score = (double)a1.score + (double)a2.score + std::max(-20.0f, log(normal_pdf(x, mu, sigma)));
             //* (1 - s2 / s1) * min_matches * log(s1);
         } else {
             // Treat a1/a2 as two single-end reads