JAL-885; Multi-helix bug is solved; One should think about a more

[jalview.git] / src / jalview / analysis / StructureFrequency.java

diff --git a/src/jalview/analysis/StructureFrequency.java b/src/jalview/analysis/StructureFrequency.java
index 29c0a0b..a5dd030 100644 (file)
--- a/src/jalview/analysis/StructureFrequency.java
+++ b/src/jalview/analysis/StructureFrequency.java
@@ -181,6 +181,15 @@ public class StructureFrequency
     }
   }
   
+  public static int findPair(SequenceFeature[] pairs,int indice){
+         for(int i=0; i<pairs.length; i++){
+                 if(pairs[i].getBegin()==indice){
+                         return pairs[i].getEnd();
+                 }
+         }
+         return -1;
+  }
+  
   /**
    * Method to calculate a 'base pair consensus row', very similar 
    * to nucleotide consensus but takes into account a given structure
@@ -193,7 +202,7 @@ public class StructureFrequency
    */
   public static final void calculate(SequenceI[] sequences, int start,
                  int end, Hashtable[] result, boolean profile, AlignmentAnnotation rnaStruc){
-
+         System.out.println("StructureFrequency.calculate");
          Hashtable residueHash;
          String maxResidue;
          char[] seq, struc=rnaStruc.getRNAStruc().toCharArray();
@@ -203,10 +212,15 @@ public class StructureFrequency
          int[] values = new int[255];
          float percentage;
          
+         /*for(int x=(rna.length-1); x>=0; x--){
+                 System.out.println("BP-"+((rna.length-1)-x)+" "+rna[x].getBegin()+","+rna[x].getEnd()+" "+rna[x].getFeatureGroup());
+         }*/
+         
+         
          for (i = start; i < end; i++) //foreach column
          {
                  residueHash = new Hashtable();
-                 maxResidue="";
+                 maxResidue="-";
                  values = new int[255];
                  bpEnd=-1;
 
@@ -216,7 +230,7 @@ public class StructureFrequency
                          s = '-';
                  }
 
-                 if(s == '-'){
+                 if(s != '('){
                          values['-']++;
                  }
                  else
@@ -245,35 +259,31 @@ public class StructureFrequency
                                                  values['-']++;
                                                  continue;
                                          }
-                                         if(s == '('){
-                                                 bpEnd=rna[(rna.length-1-nonGap)].getEnd();              
+                                         //if(s == '('){
+                                                 //bpEnd=rna[(rna.length-1-nonGap)].getEnd();
+                                                 bpEnd=findPair(rna,i);
                                                  cEnd=seq[bpEnd];
+                                                 //System.out.println(i+","+bpEnd);
                                                  if(checkBpType(c,cEnd)){
                                                          values['H']++; //H means it's a helix (structured)
                                                  }
-                                         }
+                                                 maxResidue="H";
+                                 //}
                                  }
                          }
                          nonGap++;
                  }
-                 /*UPDATE this for new values
-             if (profile)
-             {
-               residueHash.put(PROFILE, new int[][]
-               { values, new int[]
-               { jSize, nongap } });
-             }
-                  */
+                 //UPDATE this for new values
+             //if (profile)
+             //{
+            //   residueHash.put(PROFILE, new int[][]
+             //  { values, new int[]
+            //   { jSize, nongap } });
+            // }
+                  
 
                  count=values['H'];
-                 
-                 if (count == 0)
-                 {
-                         maxResidue = "-";
-                 }else{
-                         maxResidue="H";
-                 }
-                 
+                                 
                  residueHash.put(MAXCOUNT, new Integer(count));
              residueHash.put(MAXRESIDUE, maxResidue);
              
@@ -288,7 +298,6 @@ public class StructureFrequency
                  if(bpEnd>0){
                          result[bpEnd]=residueHash;
                  }
-
          }
   }