2 * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
3 * Copyright (C) $$Year-Rel$$ The Jalview Authors
5 * This file is part of Jalview.
7 * Jalview is free software: you can redistribute it and/or
8 * modify it under the terms of the GNU General Public License
9 * as published by the Free Software Foundation, either version 3
10 * of the License, or (at your option) any later version.
12 * Jalview is distributed in the hope that it will be useful, but
13 * WITHOUT ANY WARRANTY; without even the implied warranty
14 * of MERCHANTABILITY or FITNESS FOR A PARTICULAR
15 * PURPOSE. See the GNU General Public License for more details.
17 * You should have received a copy of the GNU General Public License
18 * along with Jalview. If not, see <http://www.gnu.org/licenses/>.
19 * The Jalview Authors are detailed in the 'AUTHORS' file.
21 package jalview.workers;
23 import jalview.analysis.StructureFrequency;
24 import jalview.api.AlignViewportI;
25 import jalview.api.AlignmentViewPanel;
26 import jalview.datamodel.AlignmentAnnotation;
27 import jalview.datamodel.AlignmentI;
28 import jalview.datamodel.Annotation;
29 import jalview.datamodel.SequenceI;
31 import java.util.Hashtable;
33 public class StrucConsensusThread extends AlignCalcWorker
35 public StrucConsensusThread(AlignViewportI alignViewport,
36 AlignmentViewPanel alignPanel)
38 super(alignViewport, alignPanel);
41 AlignmentAnnotation strucConsensus;
43 Hashtable[] hStrucConsensus;
45 private long nseq = -1;
50 if (alignViewport.isClosed())
55 AlignmentI alignment = alignViewport.getAlignment();
59 if (alignment == null || (aWidth = alignment.getWidth()) < 0)
63 strucConsensus = alignViewport.getAlignmentStrucConsensusAnnotation();
64 hStrucConsensus = alignViewport.getRnaStructureConsensusHash();
65 strucConsensus.annotations = null;
66 strucConsensus.annotations = new Annotation[aWidth];
68 hStrucConsensus = new Hashtable[aWidth];
70 AlignmentAnnotation[] aa = alignViewport.getAlignment()
71 .getAlignmentAnnotation();
72 AlignmentAnnotation rnaStruc = null;
73 // select rna struct to use for calculation
76 for (int i = 0; i < aa.length; i++)
78 if (aa[i].visible && aa[i].isRNA() && aa[i].isValidStruc())
85 // check to see if its valid
87 if (rnaStruc == null || !rnaStruc.isValidStruc())
94 final SequenceI[] arr = alignment.getSequencesArray();
96 jalview.analysis.StructureFrequency.calculate(arr, 0,
97 alignment.getWidth(), hStrucConsensus, true, rnaStruc);
98 } catch (ArrayIndexOutOfBoundsException x)
102 alignViewport.setRnaStructureConsensusHash(hStrucConsensus);
103 // TODO AlignmentAnnotation rnaStruc!!!
104 updateResultAnnotation(true);
109 * update the consensus annotation from the sequence profile data using
110 * current visualization settings.
113 public void updateAnnotation()
115 updateResultAnnotation(false);
118 public void updateResultAnnotation(boolean immediate)
120 if (immediate || !calcMan.isWorking(this) && strucConsensus != null
121 && hStrucConsensus != null)
123 StructureFrequency.completeConsensus(strucConsensus, hStrucConsensus,
124 0, hStrucConsensus.length,
125 alignViewport.isIgnoreGapsConsensus(),
126 alignViewport.isShowSequenceLogo(), nseq);