Fix mask fraction corrected FSC

teamtomo · Sep 7, 2024 · 28e6385 · 28e6385
1 parent 0ba5c0d
commit 28e6385
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/src/ttfsc/_cli.py b/src/ttfsc/_cli.py
@@ -144,8 +144,9 @@ def ttfsc_cli(
                     rlnCorrectedFourierShellCorrelationPhaseRandomizedMaskedMaps=fsc_values_randomized[i]
                     if correct_for_masking
                     else nan,
-                    rlnFourierShellCorrelationParticleMaskFraction=1
-                    - ((1 - f) * data_general.rlnParticleBoxFractionSolventMask)
+                    rlnFourierShellCorrelationParticleMaskFraction=f
+                    / data_general.rlnParticleBoxFractionSolventMask
+                    / (1.0 + (1.0 / data_general.rlnParticleBoxFractionSolventMask - 1.0) * f.abs())
                     if mask != Masking.none
                     else nan,
                 )