more comments

deconvolution/evaluate.py

@@ -405,8 +405,8 @@ def verify_graph_edges(d2_component):
 
         # Construct the molecule(s) that this udg really 'reflects'
         # i.e. the udg has a central node and two cliques
-        # that central nodes is the result of merging of several molecules
-        # ideally, only one of them is connected to the cliques
+        # that central node is the result of merging of several molecules
+        # ideally, only one of those molecules is connected to the molecules of the cliques
         # (there could be more than one though; in that case the udg is 'ambiguous')
         # udg_molecules aims to reflect the molecule(s) underlying this udg 
         udg_molecules = set()
@@ -431,7 +431,7 @@ def verify_graph_edges(d2_component):
 
             nexts.sort(key=lambda x: x)
             quality = sum([1.0/x if mol_idx+x in nexts else 0 for x in range(1,6)]) / sum([1.0/x for x in range(1,6)])
-            if quality > 0.6:
+            if quality > 0.6: # eyeballed but still arbitrary
                 udg_molecules.add(mol_idx)
             #print("mol",mol_idx,molecule_idxs,"quality",quality,"nexts",nexts)