JAL-2794 most method parameters/class variables made final
[jalview.git] / src / jalview / ext / archaeopteryx / ArchaeopteryxTreeConverter.java
index 53cc2c2..0bdd4e2 100644 (file)
@@ -2,28 +2,43 @@ package jalview.ext.archaeopteryx;
 
 import jalview.analysis.TreeBuilder;
 import jalview.datamodel.SequenceI;
+import jalview.ext.forester.ForesterConversions;
 import jalview.ext.forester.ForesterMatrix;
 import jalview.util.MessageManager;
 
+import java.util.HashMap;
+import java.util.Map;
+
 import org.forester.evoinference.matrix.distance.DistanceMatrix;
 import org.forester.phylogeny.Phylogeny;
 import org.forester.phylogeny.PhylogenyNode;
-
+import org.forester.phylogeny.data.NodeData;
+import org.forester.phylogeny.data.Sequence;
+
+/**
+ * Note that this currently demands a 1:1 relationship between nodes and
+ * sequences
+ * 
+ * @author kjvanderheide
+ *
+ */
 public class ArchaeopteryxTreeConverter
 {
   protected final SequenceI[] sequences;
 
-  private Phylogeny aptxTree;
-
-  private PhylogenyNode rootNode;
-
   protected final DistanceMatrix distances;
   
   protected final TreeBuilder jalviewTree;
   
   public String treeTitle;
 
+  private final Phylogeny aptxTree;
+
+  private PhylogenyNode rootNode;
+
+  private final Map<SequenceI, PhylogenyNode> alignmentBoundNodes;
 
+  private final Map<PhylogenyNode, SequenceI> nodesBoundAlignment;
 
   public ArchaeopteryxTreeConverter(final TreeBuilder calculatedTree)
   {
@@ -33,13 +48,19 @@ public class ArchaeopteryxTreeConverter
             jalviewTree.getDistances(), sequences);
     aptxTree = new Phylogeny();
     rootNode = new PhylogenyNode();
+    alignmentBoundNodes = new HashMap<>(distances.getSize());
+    nodesBoundAlignment = new HashMap<>(distances.getSize());
+
 
   }
 
   public Phylogeny buildAptxTree(final PhylogenyNode treeRoot)
   {
 
-    rootNode = treeRoot;
+    if (treeRoot != null)
+    {
+      rootNode = treeRoot;
+    }
 
     buildAptxTree();
 
@@ -47,13 +68,26 @@ public class ArchaeopteryxTreeConverter
 
   }
 
+
   public Phylogeny buildAptxTree()
   {
-    // NeighborJoiningF foresterClustering = NeighborJoiningF
-    // .createInstance();
-    // aptxTree = foresterClustering.execute(distances);
+
+    for (SequenceI sequence : sequences)
+    {
+      Sequence seq = ForesterConversions
+              .createForesterSequence(sequence, true);
+      PhylogenyNode sequenceNode = new PhylogenyNode(sequence.getName());
+      NodeData nodeData = sequenceNode.getNodeData();
+      nodeData.setSequence(seq);
+
+      nodesBoundAlignment.put(sequenceNode, sequence);
+      alignmentBoundNodes.put(sequence, sequenceNode);
+      rootNode.addAsChild(sequenceNode);
+    }
 
 
+    aptxTree.setRoot(rootNode);
+
     treeTitle = generateTreeName();
     aptxTree.setName(treeTitle);
 
@@ -61,13 +95,21 @@ public class ArchaeopteryxTreeConverter
 
   }
 
-  private Phylogeny clusterNodes()
+  public Map<SequenceI, PhylogenyNode> getAlignmentBoundNodes()
   {
-    return aptxTree;
+    return alignmentBoundNodes;
+  }
 
+  public Map<PhylogenyNode, SequenceI> getNodesBoundAlignment()
+  {
+    return nodesBoundAlignment;
   }
 
+  private Phylogeny clusterNodes()
+  {
+    return aptxTree;
 
+  }
   /**
    * Formats a localised title for the tree panel, like
    * <p>
@@ -77,7 +119,8 @@ public class ArchaeopteryxTreeConverter
    * 
    * @return
    */
-  public String generateTreeName()
+  public String generateTreeName() // Move this and add selection region to the
+                                   // title when applicable
   {
     if (treeTitle != null) // will currently never happen, loaded tree file will
                            // take a different path