X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fgui%2FStructureChooser.java;h=d63c94a81f70c1b918fe28f2b0ca57cb0b7c69e1;hb=744e4ade5413aa7e8006a8d4d29fdb718489bf08;hp=97ff4046126df89a79dddfc75187c6c74d729dbb;hpb=7534e6f358b2bba0db179cb2f496bd2e5a4da43d;p=jalview.git diff --git a/src/jalview/gui/StructureChooser.java b/src/jalview/gui/StructureChooser.java index 97ff404..d63c94a 100644 --- a/src/jalview/gui/StructureChooser.java +++ b/src/jalview/gui/StructureChooser.java @@ -21,44 +21,72 @@ package jalview.gui; +import java.awt.event.ActionEvent; +import java.awt.event.ActionListener; +import java.awt.event.ItemEvent; +import java.io.File; +import java.util.ArrayList; +import java.util.Collection; +import java.util.HashSet; +import java.util.LinkedHashMap; +import java.util.LinkedHashSet; +import java.util.List; +import java.util.Locale; +import java.util.Map; +import java.util.concurrent.Executors; + +import javax.swing.JCheckBox; +import javax.swing.JComboBox; +import javax.swing.JLabel; +import javax.swing.JMenuItem; +import javax.swing.JPopupMenu; +import javax.swing.JTable; +import javax.swing.SwingUtilities; +import javax.swing.table.AbstractTableModel; + +import com.stevesoft.pat.Regex; + +import jalview.analysis.AlignmentUtils; +import jalview.api.AlignmentViewPanel; import jalview.api.structures.JalviewStructureDisplayI; import jalview.bin.Cache; +import jalview.bin.Console; import jalview.bin.Jalview; -import jalview.datamodel.DBRefEntry; -import jalview.datamodel.DBRefSource; +import jalview.datamodel.AlignmentAnnotation; +import jalview.datamodel.AlignmentI; import jalview.datamodel.PDBEntry; +import jalview.datamodel.SequenceGroup; import jalview.datamodel.SequenceI; +import jalview.ext.jmol.JmolParser; import jalview.fts.api.FTSData; import jalview.fts.api.FTSDataColumnI; import jalview.fts.api.FTSRestClientI; +import jalview.fts.core.FTSDataColumnPreferences; import jalview.fts.core.FTSRestRequest; import jalview.fts.core.FTSRestResponse; import jalview.fts.service.pdb.PDBFTSRestClient; +import jalview.fts.service.threedbeacons.TDB_FTSData; +import jalview.gui.StructureViewer.ViewerType; +import jalview.gui.structurechooser.PDBStructureChooserQuerySource; +import jalview.gui.structurechooser.StructureChooserQuerySource; +import jalview.gui.structurechooser.ThreeDBStructureChooserQuerySource; import jalview.io.DataSourceType; +import jalview.io.JalviewFileChooser; +import jalview.io.JalviewFileView; +import jalview.jbgui.FilterOption; import jalview.jbgui.GStructureChooser; +import jalview.structure.StructureImportSettings.TFType; import jalview.structure.StructureMapping; import jalview.structure.StructureSelectionManager; import jalview.util.MessageManager; +import jalview.util.Platform; +import jalview.util.StringUtils; import jalview.ws.DBRefFetcher; +import jalview.ws.DBRefFetcher.FetchFinishedListenerI; +import jalview.ws.datamodel.alphafold.PAEContactMatrix; +import jalview.ws.seqfetcher.DbSourceProxy; import jalview.ws.sifts.SiftsSettings; -import java.awt.event.ItemEvent; -import java.util.ArrayList; -import java.util.Collection; -import java.util.HashSet; -import java.util.LinkedHashSet; -import java.util.List; -import java.util.Objects; -import java.util.Set; -import java.util.Vector; - -import javax.swing.JCheckBox; -import javax.swing.JComboBox; -import javax.swing.JLabel; -import javax.swing.JTable; -import javax.swing.SwingUtilities; -import javax.swing.table.AbstractTableModel; - /** * Provides the behaviors for the Structure chooser Panel * @@ -71,7 +99,10 @@ public class StructureChooser extends GStructureChooser { private static final String AUTOSUPERIMPOSE = "AUTOSUPERIMPOSE"; - private static int MAX_QLENGTH = 7820; + /** + * warn user if need to fetch more than this many uniprot records at once + */ + private static final int THRESHOLD_WARN_UNIPROT_FETCH_NEEDED = 20; private SequenceI selectedSequence; @@ -81,29 +112,105 @@ public class StructureChooser extends GStructureChooser private Collection discoveredStructuresSet; - private FTSRestRequest lastPdbRequest; + private StructureChooserQuerySource data; - private FTSRestClientI pdbRestCleint; + @Override + protected FTSDataColumnPreferences getFTSDocFieldPrefs() + { + return data.getDocFieldPrefs(); + } private String selectedPdbFileName; + private TFType localPdbTempfacType; + + private String localPdbPaeMatrixFileName; + private boolean isValidPBDEntry; private boolean cachedPDBExists; + private Collection lastDiscoveredStructuresSet; + + private boolean canQueryTDB = false; + + private boolean notQueriedTDBYet = true; + + List seqsWithoutSourceDBRef = null; + + private boolean showChooserGUI = true; + private static StructureViewer lastTargetedView = null; public StructureChooser(SequenceI[] selectedSeqs, SequenceI selectedSeq, AlignmentPanel ap) { + this(selectedSeqs, selectedSeq, ap, true); + } + + public StructureChooser(SequenceI[] selectedSeqs, SequenceI selectedSeq, + AlignmentPanel ap, boolean showGUI) + { + // which FTS engine to use + data = StructureChooserQuerySource.getQuerySourceFor(selectedSeqs); + initDialog(); + this.ap = ap; this.selectedSequence = selectedSeq; this.selectedSequences = selectedSeqs; this.progressIndicator = (ap == null) ? null : ap.alignFrame; + this.showChooserGUI = showGUI; init(); + } /** + * sets canQueryTDB if protein sequences without a canonical uniprot ref or at + * least one structure are discovered. + */ + private void populateSeqsWithoutSourceDBRef() + { + seqsWithoutSourceDBRef = new ArrayList(); + boolean needCanonical = false; + for (SequenceI seq : selectedSequences) + { + if (seq.isProtein()) + { + int dbRef = ThreeDBStructureChooserQuerySource + .checkUniprotRefs(seq.getDBRefs()); + if (dbRef < 0) + { + if (dbRef == -1) + { + // need to retrieve canonicals + needCanonical = true; + seqsWithoutSourceDBRef.add(seq); + } + else + { + // could be a sequence with pdb ref + if (seq.getAllPDBEntries() == null + || seq.getAllPDBEntries().size() == 0) + { + seqsWithoutSourceDBRef.add(seq); + } + } + } + } + } + // retrieve database refs for protein sequences + if (!seqsWithoutSourceDBRef.isEmpty()) + { + canQueryTDB = true; + if (needCanonical) + { + // triggers display of the 'Query TDB' button + notQueriedTDBYet = true; + } + } + }; + + /** * Initializes parameters used by the Structure Chooser Panel */ protected void init() @@ -114,31 +221,188 @@ public class StructureChooser extends GStructureChooser } chk_superpose.setSelected(Cache.getDefault(AUTOSUPERIMPOSE, true)); + btn_queryTDB.addActionListener(new ActionListener() + { + + @Override + public void actionPerformed(ActionEvent e) + { + promptForTDBFetch(false); + } + }); + + Executors.defaultThreadFactory().newThread(new Runnable() + { + @Override + public void run() + { + populateSeqsWithoutSourceDBRef(); + initialStructureDiscovery(); + } + + }).start(); + + } + + // called by init + private void initialStructureDiscovery() + { + // check which FTS engine to use + data = StructureChooserQuerySource.getQuerySourceFor(selectedSequences); // ensure a filter option is in force for search populateFilterComboBox(true, cachedPDBExists); - Thread discoverPDBStructuresThread = new Thread(new Runnable() + + // looks for any existing structures already loaded + // for the sequences (the cached ones) + // then queries the StructureChooserQuerySource to + // discover more structures. + // + // Possible optimisation is to only begin querying + // the structure chooser if there are no cached structures. + + long startTime = System.currentTimeMillis(); + updateProgressIndicator( + MessageManager.getString("status.loading_cached_pdb_entries"), + startTime); + loadLocalCachedPDBEntries(); + updateProgressIndicator(null, startTime); + updateProgressIndicator( + MessageManager.getString("status.searching_for_pdb_structures"), + startTime); + fetchStructuresMetaData(); + // revise filter options if no results were found + populateFilterComboBox(isStructuresDiscovered(), cachedPDBExists); + discoverStructureViews(); + updateProgressIndicator(null, startTime); + mainFrame.setVisible(showChooserGUI); + updateCurrentView(); + } + + /** + * raises dialog for Uniprot fetch followed by 3D beacons search + * + * @param ignoreGui + * - when true, don't ask, just fetch + */ + public void promptForTDBFetch(boolean ignoreGui) + { + final long progressId = System.currentTimeMillis(); + + // final action after prompting and discovering db refs + final Runnable strucDiscovery = new Runnable() { @Override public void run() { - long startTime = System.currentTimeMillis(); - updateProgressIndicator(MessageManager - .getString("status.loading_cached_pdb_entries"), startTime); - loadLocalCachedPDBEntries(); - updateProgressIndicator(null, startTime); - updateProgressIndicator(MessageManager.getString( - "status.searching_for_pdb_structures"), startTime); - fetchStructuresMetaData(); - // revise filter options if no results were found - populateFilterComboBox(isStructuresDiscovered(), cachedPDBExists); - discoverStructureViews(); - updateProgressIndicator(null, startTime); - mainFrame.setVisible(true); - updateCurrentView(); + mainFrame.setEnabled(false); + cmb_filterOption.setEnabled(false); + progressBar.setProgressBar( + MessageManager.getString("status.searching_3d_beacons"), + progressId); + btn_queryTDB.setEnabled(false); + // TODO: warn if no accessions discovered + populateSeqsWithoutSourceDBRef(); + // redo initial discovery - this time with 3d beacons + // Executors. + previousWantedFields = null; + lastSelected = (FilterOption) cmb_filterOption.getSelectedItem(); + cmb_filterOption.setSelectedItem(null); + cachedPDBExists = false; // reset to initial + initialStructureDiscovery(); + if (!isStructuresDiscovered()) + { + progressBar.setProgressBar(MessageManager.getString( + "status.no_structures_discovered_from_3d_beacons"), + progressId); + btn_queryTDB.setToolTipText(MessageManager.getString( + "status.no_structures_discovered_from_3d_beacons")); + btn_queryTDB.setEnabled(false); + pnl_queryTDB.setVisible(false); + } + else + { + cmb_filterOption.setSelectedIndex(0); // select 'best' + btn_queryTDB.setVisible(false); + pnl_queryTDB.setVisible(false); + progressBar.setProgressBar(null, progressId); + } + mainFrame.setEnabled(true); + cmb_filterOption.setEnabled(true); } - }); - discoverPDBStructuresThread.start(); + }; + + final FetchFinishedListenerI afterDbRefFetch = new FetchFinishedListenerI() + { + + @Override + public void finished() + { + // filter has been selected, so we set flag to remove ourselves + notQueriedTDBYet = false; + // new thread to discover structures - via 3d beacons + Executors.defaultThreadFactory().newThread(strucDiscovery).start(); + + } + }; + + // fetch db refs if OK pressed + final Runnable discoverCanonicalDBrefs = () -> { + btn_queryTDB.setEnabled(false); + populateSeqsWithoutSourceDBRef(); + + final int y = seqsWithoutSourceDBRef.size(); + if (y > 0) + { + final SequenceI[] seqWithoutSrcDBRef = seqsWithoutSourceDBRef + .toArray(new SequenceI[y]); + DBRefFetcher dbRefFetcher = new DBRefFetcher(seqWithoutSrcDBRef, + progressBar, new DbSourceProxy[] + { new jalview.ws.dbsources.Uniprot() }, null, false); + dbRefFetcher.addListener(afterDbRefFetch); + // ideally this would also gracefully run with callbacks + + dbRefFetcher.fetchDBRefs(true); + } + else + { + // call finished action directly + afterDbRefFetch.finished(); + } + }; + final Runnable revertview = () -> { + if (lastSelected != null) + { + cmb_filterOption.setSelectedItem(lastSelected); + } + }; + int threshold = Cache.getDefault("UNIPROT_AUTOFETCH_THRESHOLD", + THRESHOLD_WARN_UNIPROT_FETCH_NEEDED); + Console.debug("Using Uniprot fetch threshold of " + threshold); + if (ignoreGui || seqsWithoutSourceDBRef.size() < threshold) + { + Executors.newSingleThreadExecutor().submit(discoverCanonicalDBrefs); + return; + } + // need cancel and no to result in the discoverPDB action - mocked is + // 'cancel' TODO: mock should be OK + + StructureChooser thisSC = this; + JvOptionPane.newOptionDialog(thisSC.getFrame()) + .setResponseHandler(JvOptionPane.OK_OPTION, + discoverCanonicalDBrefs) + .setResponseHandler(JvOptionPane.CANCEL_OPTION, revertview) + .setResponseHandler(JvOptionPane.NO_OPTION, revertview) + .showDialog( + MessageManager.formatMessage( + "label.fetch_references_for_3dbeacons", + seqsWithoutSourceDBRef.size()), + MessageManager.getString("label.3dbeacons"), + JvOptionPane.YES_NO_OPTION, JvOptionPane.PLAIN_MESSAGE, + null, new Object[] + { MessageManager.getString("action.ok"), + MessageManager.getString("action.cancel") }, + MessageManager.getString("action.ok"), false); } /** @@ -165,8 +429,7 @@ public class StructureChooser extends GStructureChooser if (view.isLinkedWith(ap)) { - targetView.insertItemAt(viewHandler, - linkedViewsAt++); + targetView.insertItemAt(viewHandler, linkedViewsAt++); } else { @@ -217,36 +480,36 @@ public class StructureChooser extends GStructureChooser void fetchStructuresMetaData() { long startTime = System.currentTimeMillis(); - pdbRestCleint = PDBFTSRestClient.getInstance(); - Collection wantedFields = pdbDocFieldPrefs + Collection wantedFields = data.getDocFieldPrefs() .getStructureSummaryFields(); discoveredStructuresSet = new LinkedHashSet<>(); HashSet errors = new HashSet<>(); + + FilterOption selectedFilterOpt = ((FilterOption) cmb_filterOption + .getSelectedItem()); + for (SequenceI seq : selectedSequences) { - FTSRestRequest pdbRequest = new FTSRestRequest(); - pdbRequest.setAllowEmptySeq(false); - pdbRequest.setResponseSize(500); - pdbRequest.setFieldToSearchBy("("); - FilterOption selectedFilterOpt = ((FilterOption) cmb_filterOption - .getSelectedItem()); - pdbRequest.setFieldToSortBy(selectedFilterOpt.getValue(), - !chk_invertFilter.isSelected()); - pdbRequest.setWantedFields(wantedFields); - pdbRequest.setSearchTerm(buildQuery(seq) + ")"); - pdbRequest.setAssociatedSequence(seq); + FTSRestResponse resultList; try { - resultList = pdbRestCleint.executeRequest(pdbRequest); + resultList = data.fetchStructuresMetaData(seq, wantedFields, + selectedFilterOpt, !chk_invertFilter.isSelected()); + // null response means the FTSengine didn't yield a query for this + // consider designing a special exception if we really wanted to be + // OOCrazy + if (resultList == null) + { + continue; + } } catch (Exception e) { e.printStackTrace(); errors.add(e.getMessage()); continue; } - lastPdbRequest = pdbRequest; if (resultList.getSearchSummary() != null && !resultList.getSearchSummary().isEmpty()) { @@ -260,9 +523,11 @@ public class StructureChooser extends GStructureChooser if (discoveredStructuresSet != null && !discoveredStructuresSet.isEmpty()) { - getResultTable().setModel(FTSRestResponse - .getTableModel(lastPdbRequest, discoveredStructuresSet)); + getResultTable() + .setModel(data.getTableModel(discoveredStructuresSet)); + noOfStructuresFound = discoveredStructuresSet.size(); + lastDiscoveredStructuresSet = discoveredStructuresSet; mainFrame.setTitle(MessageManager.formatMessage( "label.structure_chooser_no_of_structures", noOfStructuresFound, totalTime)); @@ -309,145 +574,6 @@ public class StructureChooser extends GStructureChooser } /** - * Builds a query string for a given sequences using its DBRef entries - * - * @param seq - * the sequences to build a query for - * @return the built query string - */ - - static String buildQuery(SequenceI seq) - { - boolean isPDBRefsFound = false; - boolean isUniProtRefsFound = false; - StringBuilder queryBuilder = new StringBuilder(); - Set seqRefs = new LinkedHashSet<>(); - - if (seq.getAllPDBEntries() != null - && queryBuilder.length() < MAX_QLENGTH) - { - for (PDBEntry entry : seq.getAllPDBEntries()) - { - if (isValidSeqName(entry.getId())) - { - queryBuilder.append("pdb_id:").append(entry.getId().toLowerCase()) - .append(" OR "); - isPDBRefsFound = true; - } - } - } - - if (seq.getDBRefs() != null && seq.getDBRefs().length != 0) - { - for (DBRefEntry dbRef : seq.getDBRefs()) - { - if (isValidSeqName(getDBRefId(dbRef)) - && queryBuilder.length() < MAX_QLENGTH) - { - if (dbRef.getSource().equalsIgnoreCase(DBRefSource.UNIPROT)) - { - queryBuilder.append("uniprot_accession:") - .append(getDBRefId(dbRef)).append(" OR "); - queryBuilder.append("uniprot_id:").append(getDBRefId(dbRef)) - .append(" OR "); - isUniProtRefsFound = true; - } - else if (dbRef.getSource().equalsIgnoreCase(DBRefSource.PDB)) - { - - queryBuilder.append("pdb_id:") - .append(getDBRefId(dbRef).toLowerCase()).append(" OR "); - isPDBRefsFound = true; - } - else - { - seqRefs.add(getDBRefId(dbRef)); - } - } - } - } - - if (!isPDBRefsFound && !isUniProtRefsFound) - { - String seqName = seq.getName(); - seqName = sanitizeSeqName(seqName); - String[] names = seqName.toLowerCase().split("\\|"); - for (String name : names) - { - // System.out.println("Found name : " + name); - name.trim(); - if (isValidSeqName(name)) - { - seqRefs.add(name); - } - } - - for (String seqRef : seqRefs) - { - queryBuilder.append("text:").append(seqRef).append(" OR "); - } - } - - int endIndex = queryBuilder.lastIndexOf(" OR "); - if (queryBuilder.toString().length() < 6) - { - return null; - } - String query = queryBuilder.toString().substring(0, endIndex); - return query; - } - - /** - * Remove the following special characters from input string +, -, &, !, (, ), - * {, }, [, ], ^, ", ~, *, ?, :, \ - * - * @param seqName - * @return - */ - static String sanitizeSeqName(String seqName) - { - Objects.requireNonNull(seqName); - return seqName.replaceAll("\\[\\d*\\]", "") - .replaceAll("[^\\dA-Za-z|_]", "").replaceAll("\\s+", "+"); - } - - /** - * Ensures sequence ref names are not less than 3 characters and does not - * contain a database name - * - * @param seqName - * @return - */ - static boolean isValidSeqName(String seqName) - { - // System.out.println("seqName : " + seqName); - String ignoreList = "pdb,uniprot,swiss-prot"; - if (seqName.length() < 3) - { - return false; - } - if (seqName.contains(":")) - { - return false; - } - seqName = seqName.toLowerCase(); - for (String ignoredEntry : ignoreList.split(",")) - { - if (seqName.contains(ignoredEntry)) - { - return false; - } - } - return true; - } - - static String getDBRefId(DBRefEntry dbRef) - { - String ref = dbRef.getAccessionId().replaceAll("GO:", ""); - return ref; - } - - /** * Filters a given list of discovered structures based on supplied argument * * @param fieldToFilterBy @@ -457,54 +583,33 @@ public class StructureChooser extends GStructureChooser { Thread filterThread = new Thread(new Runnable() { + @Override public void run() { long startTime = System.currentTimeMillis(); - pdbRestCleint = PDBFTSRestClient.getInstance(); lbl_loading.setVisible(true); - Collection wantedFields = pdbDocFieldPrefs + Collection wantedFields = data.getDocFieldPrefs() .getStructureSummaryFields(); Collection filteredResponse = new HashSet<>(); HashSet errors = new HashSet<>(); for (SequenceI seq : selectedSequences) { - FTSRestRequest pdbRequest = new FTSRestRequest(); - if (fieldToFilterBy.equalsIgnoreCase("uniprot_coverage")) - { - pdbRequest.setAllowEmptySeq(false); - pdbRequest.setResponseSize(1); - pdbRequest.setFieldToSearchBy("("); - pdbRequest.setSearchTerm(buildQuery(seq) + ")"); - pdbRequest.setWantedFields(wantedFields); - pdbRequest.setAssociatedSequence(seq); - pdbRequest.setFacet(true); - pdbRequest.setFacetPivot(fieldToFilterBy + ",entry_entity"); - pdbRequest.setFacetPivotMinCount(1); - } - else - { - pdbRequest.setAllowEmptySeq(false); - pdbRequest.setResponseSize(1); - pdbRequest.setFieldToSearchBy("("); - pdbRequest.setFieldToSortBy(fieldToFilterBy, - !chk_invertFilter.isSelected()); - pdbRequest.setSearchTerm(buildQuery(seq) + ")"); - pdbRequest.setWantedFields(wantedFields); - pdbRequest.setAssociatedSequence(seq); - } + FTSRestResponse resultList; try { - resultList = pdbRestCleint.executeRequest(pdbRequest); + resultList = data.selectFirstRankedQuery(seq, + discoveredStructuresSet, wantedFields, fieldToFilterBy, + !chk_invertFilter.isSelected()); + } catch (Exception e) { e.printStackTrace(); errors.add(e.getMessage()); continue; } - lastPdbRequest = pdbRequest; if (resultList.getSearchSummary() != null && !resultList.getSearchSummary().isEmpty()) { @@ -520,8 +625,8 @@ public class StructureChooser extends GStructureChooser Collection reorderedStructuresSet = new LinkedHashSet<>(); reorderedStructuresSet.addAll(filteredResponse); reorderedStructuresSet.addAll(discoveredStructuresSet); - getResultTable().setModel(FTSRestResponse - .getTableModel(lastPdbRequest, reorderedStructuresSet)); + getResultTable() + .setModel(data.getTableModel(reorderedStructuresSet)); FTSRestResponse.configureTableColumn(getResultTable(), wantedFields, tempUserPrefs); @@ -565,11 +670,12 @@ public class StructureChooser extends GStructureChooser @Override protected void pdbFromFile_actionPerformed() { - // TODO: JAL-3048 not needed for Jalview-JS until JSmol dep and StructureChooser + // TODO: JAL-3048 not needed for Jalview-JS until JSmol dep and + // StructureChooser // works - jalview.io.JalviewFileChooser chooser = new jalview.io.JalviewFileChooser( - jalview.bin.Cache.getProperty("LAST_DIRECTORY")); - chooser.setFileView(new jalview.io.JalviewFileView()); + JalviewFileChooser chooser = new JalviewFileChooser( + Cache.getProperty("LAST_DIRECTORY")); + chooser.setFileView(new JalviewFileView()); chooser.setDialogTitle( MessageManager.formatMessage("label.select_pdb_file_for", selectedSequence.getDisplayId(false))); @@ -578,61 +684,180 @@ public class StructureChooser extends GStructureChooser selectedSequence.getDisplayId(false))); int value = chooser.showOpenDialog(null); - if (value == jalview.io.JalviewFileChooser.APPROVE_OPTION) + if (value == JalviewFileChooser.APPROVE_OPTION) { selectedPdbFileName = chooser.getSelectedFile().getPath(); - jalview.bin.Cache.setProperty("LAST_DIRECTORY", selectedPdbFileName); + Cache.setProperty("LAST_DIRECTORY", selectedPdbFileName); + boolean guessTFType = localPdbPaeMatrixFileName == null; + localPdbPaeMatrixFileName = guessPAEFilename(); + guessTFType |= localPdbPaeMatrixFileName != null; + Regex alphaFold = JmolParser.getNewAlphafoldValidator(); + if (guessTFType + && alphaFold.search(new File(selectedPdbFileName).getName()) + && !tempFacAsChanged) + { + // localPdbPaeMatrixFileName was null and now isn't and filename could + // well be AlphaFold and user hasn't adjusted the tempFacType + combo_tempFacAs.setSelectedItem(TFType.PLDDT); + } validateSelections(); } } /** + * Handles action event for btn_paeMatrixFile + */ + @Override + protected void paeMatrixFile_actionPerformed() + { + File pdbFile = new File(selectedPdbFileName); + String setFile = Cache.getProperty("LAST_DIRECTORY"); + if (localPdbPaeMatrixFileName != null) + { + File paeFile = new File(localPdbPaeMatrixFileName); + if (paeFile.exists()) + setFile = paeFile.getAbsolutePath(); + else if (paeFile.getParentFile().exists()) + setFile = paeFile.getParentFile().getAbsolutePath(); + } + else + { + String guess = guessPAEFilename(); + if (guess != null) + setFile = guess; + } + JalviewFileChooser chooser = new JalviewFileChooser(setFile); + chooser.setFileView(new JalviewFileView()); + chooser.setDialogTitle(MessageManager.formatMessage( + "label.select_pae_matrix_file_for", pdbFile.getName())); + chooser.setToolTipText(MessageManager.formatMessage( + "label.load_pae_matrix_file_associate_with_structure", + pdbFile.getName())); + + // TODO convert to Callable/Promise + int value = chooser.showOpenDialog(null); + if (value == JalviewFileChooser.APPROVE_OPTION) + { + String fileName = chooser.getSelectedFile().getPath(); + try + { + PAEContactMatrix.validateContactMatrixFile(fileName); + } catch (Exception thr) + { + JvOptionPane.showInternalMessageDialog(this, MessageManager + .formatMessage("label.couldnt_load_file", new Object[] + { fileName }) + "
" + thr.getLocalizedMessage(), + MessageManager.getString("label.error_loading_file"), + JvOptionPane.WARNING_MESSAGE); + Console.error("Couldn't import " + fileName + " as a PAE matrix", + thr); + return; + } + localPdbPaeMatrixFileName = fileName; + Cache.setProperty("LAST_DIRECTORY", localPdbPaeMatrixFileName); + } + validateAssociationFromFile(); + } + + private String guessPAEFilename() + { + if (selectedPdbFileName.toLowerCase(Locale.ROOT).endsWith(".pdb") + || selectedPdbFileName.toLowerCase(Locale.ROOT) + .endsWith(".cif")) + { + String jsonExt = selectedPdbFileName.substring(0, + selectedPdbFileName.length() - 4) + ".json"; + // AlphaFold naming scheme + String guessFile1 = StringUtils.replaceLast(jsonExt, "model", + "predicted_aligned_error"); + // nf-core mode naming scheme + String guessFile2 = StringUtils.replaceLast(jsonExt, ".json", + "_scores.json"); + if (new File(guessFile1).exists()) + { + return guessFile1; + } + else if (new File(jsonExt).exists()) + { + return jsonExt; + } + else if (new File(guessFile2).exists()) + { + return guessFile2; + } + } + return null; + } + + /** * Populates the filter combo-box options dynamically depending on discovered * structures */ protected void populateFilterComboBox(boolean haveData, boolean cachedPDBExist) { + populateFilterComboBox(haveData, cachedPDBExist, null); + } + + /** + * Populates the filter combo-box options dynamically depending on discovered + * structures + */ + protected void populateFilterComboBox(boolean haveData, + boolean cachedPDBExist, FilterOption lastSel) + { + /* * temporarily suspend the change listener behaviour */ cmb_filterOption.removeItemListener(this); - + int selSet = -1; cmb_filterOption.removeAllItems(); if (haveData) { - cmb_filterOption.addItem(new FilterOption( - MessageManager.getString("label.best_quality"), - "overall_quality", VIEWS_FILTER, false)); - cmb_filterOption.addItem(new FilterOption( - MessageManager.getString("label.best_resolution"), - "resolution", VIEWS_FILTER, false)); - cmb_filterOption.addItem(new FilterOption( - MessageManager.getString("label.most_protein_chain"), - "number_of_protein_chains", VIEWS_FILTER, false)); - cmb_filterOption.addItem(new FilterOption( - MessageManager.getString("label.most_bound_molecules"), - "number_of_bound_molecules", VIEWS_FILTER, false)); - cmb_filterOption.addItem(new FilterOption( - MessageManager.getString("label.most_polymer_residues"), - "number_of_polymer_residues", VIEWS_FILTER, true)); + List filters = data + .getAvailableFilterOptions(VIEWS_FILTER); + data.updateAvailableFilterOptions(VIEWS_FILTER, filters, + lastDiscoveredStructuresSet); + int p = 0; + for (FilterOption filter : filters) + { + if (lastSel != null && filter.equals(lastSel)) + { + selSet = p; + } + p++; + cmb_filterOption.addItem(filter); + } } + cmb_filterOption.addItem( new FilterOption(MessageManager.getString("label.enter_pdb_id"), - "-", VIEWS_ENTER_ID, false)); + "-", VIEWS_ENTER_ID, false, null)); cmb_filterOption.addItem( new FilterOption(MessageManager.getString("label.from_file"), - "-", VIEWS_FROM_FILE, false)); + "-", VIEWS_FROM_FILE, false, null)); + if (canQueryTDB && notQueriedTDBYet) + { + btn_queryTDB.setVisible(true); + pnl_queryTDB.setVisible(true); + } if (cachedPDBExist) { FilterOption cachedOption = new FilterOption( - MessageManager.getString("label.cached_structures"), - "-", VIEWS_LOCAL_PDB, false); + MessageManager.getString("label.cached_structures"), "-", + VIEWS_LOCAL_PDB, false, null); cmb_filterOption.addItem(cachedOption); - cmb_filterOption.setSelectedItem(cachedOption); + if (selSet == -1) + { + cmb_filterOption.setSelectedItem(cachedOption); + } + } + if (selSet > -1) + { + cmb_filterOption.setSelectedIndex(selSet); } - cmb_filterOption.addItemListener(this); } @@ -643,16 +868,41 @@ public class StructureChooser extends GStructureChooser { FilterOption selectedFilterOpt = ((FilterOption) cmb_filterOption .getSelectedItem()); + + if (lastSelected == selectedFilterOpt) + { + // don't need to do anything, probably + return; + } + // otherwise, record selection + // and update the layout and dialog accordingly + lastSelected = selectedFilterOpt; + layout_switchableViews.show(pnl_switchableViews, selectedFilterOpt.getView()); String filterTitle = mainFrame.getTitle(); mainFrame.setTitle(frameTitle); chk_invertFilter.setVisible(false); + if (selectedFilterOpt.getView() == VIEWS_FILTER) { mainFrame.setTitle(filterTitle); - chk_invertFilter.setVisible(true); - filterResultSet(selectedFilterOpt.getValue()); + // TDB Query has no invert as yet + chk_invertFilter.setVisible(selectedFilterOpt + .getQuerySource() instanceof PDBStructureChooserQuerySource); + + if (data != selectedFilterOpt.getQuerySource() + || data.needsRefetch(selectedFilterOpt)) + { + data = selectedFilterOpt.getQuerySource(); + // rebuild the views completely, since prefs will also change + tabRefresh(); + return; + } + else + { + filterResultSet(selectedFilterOpt.getValue()); + } } else if (selectedFilterOpt.getView() == VIEWS_ENTER_ID || selectedFilterOpt.getView() == VIEWS_FROM_FILE) @@ -718,6 +968,46 @@ public class StructureChooser extends GStructureChooser .setEnabled(selectedCount > 1 || targetView.getItemCount() > 0); } + @Override + protected boolean showPopupFor(int selectedRow, int x, int y) + { + FilterOption selectedFilterOpt = ((FilterOption) cmb_filterOption + .getSelectedItem()); + String currentView = selectedFilterOpt.getView(); + + if (currentView == VIEWS_FILTER + && data instanceof ThreeDBStructureChooserQuerySource) + { + + TDB_FTSData row = ((ThreeDBStructureChooserQuerySource) data) + .getFTSDataFor(getResultTable(), selectedRow, + discoveredStructuresSet); + String pageUrl = row.getModelViewUrl(); + JPopupMenu popup = new JPopupMenu("3D Beacons"); + JMenuItem viewUrl = new JMenuItem("View model web page"); + viewUrl.addActionListener(new ActionListener() + { + @Override + public void actionPerformed(ActionEvent e) + { + Desktop.showUrl(pageUrl); + } + }); + popup.add(viewUrl); + SwingUtilities.invokeLater(new Runnable() + { + @Override + public void run() + { + popup.show(getResultTable(), x, y); + } + }); + return true; + } + // event not handled by us + return false; + } + /** * Validates inputs from the Manual PDB entry panel */ @@ -766,7 +1056,9 @@ public class StructureChooser extends GStructureChooser { AssociateSeqOptions assSeqOpt = (AssociateSeqOptions) fileChooserAssSeqPanel .getCmb_assSeq().getSelectedItem(); - lbl_fromFileStatus.setIcon(errorImage); + // lbl_fromFileStatus.setIcon(errorImage); + String pdbFileString = ""; + String pdbFileTooltip = ""; if (selectedSequences.length == 1 || (assSeqOpt != null && !assSeqOpt .getName().equalsIgnoreCase("-Select Associated Seq-"))) { @@ -774,14 +1066,46 @@ public class StructureChooser extends GStructureChooser if (selectedPdbFileName != null && selectedPdbFileName.length() > 0) { btn_add.setEnabled(true); - lbl_fromFileStatus.setIcon(goodImage); + // lbl_fromFileStatus.setIcon(goodImage); + pdbFileString = new File(selectedPdbFileName).getName(); + pdbFileTooltip = new File(selectedPdbFileName).getAbsolutePath(); + setPdbOptionsEnabled(true); + } + else + { + pdbFileString = MessageManager.getString("label.none"); + pdbFileTooltip = MessageManager.getString("label.nothing_selected"); + setPdbOptionsEnabled(false); } } else { btn_pdbFromFile.setEnabled(false); - lbl_fromFileStatus.setIcon(errorImage); + setPdbOptionsEnabled(false); + // lbl_fromFileStatus.setIcon(errorImage); + pdbFileString = MessageManager.getString("label.none"); + pdbFileTooltip = MessageManager.getString("label.nothing_selected"); } + lbl_pdbFile.setText(pdbFileString); + lbl_pdbFile.setToolTipText(pdbFileTooltip); + + // PAE file choice + String paeFileString = ""; + String paeFileTooltip = ""; + if (localPdbPaeMatrixFileName != null + && localPdbPaeMatrixFileName.length() > 0) + { + paeFileString = new File(localPdbPaeMatrixFileName).getName(); + paeFileTooltip = new File(localPdbPaeMatrixFileName) + .getAbsolutePath(); + } + else + { + paeFileString = MessageManager.getString("label.none"); + paeFileTooltip = MessageManager.getString("label.nothing_selected"); + } + lbl_paeFile.setText(paeFileString); + lbl_paeFile.setToolTipText(paeFileTooltip); } @Override @@ -790,6 +1114,8 @@ public class StructureChooser extends GStructureChooser validateSelections(); } + private FilterOption lastSelected = null; + /** * Handles the state change event for the 'filter' combo-box and 'invert' * check-box @@ -849,7 +1175,7 @@ public class StructureChooser extends GStructureChooser } return found; } - + /** * Handles the 'New View' action */ @@ -880,7 +1206,14 @@ public class StructureChooser extends GStructureChooser final StructureSelectionManager ssm = ap.getStructureSelectionManager(); final int preferredHeight = pnl_filter.getHeight(); + btn_add.setEnabled(false); + btn_newView.setEnabled(false); + btn_cancel.setEnabled(false); + actionsPanel.setEnabled(false); + final String progress = MessageManager + .getString("label.working_ellipsis"); + setProgressBar(progress, progress.hashCode()); Runnable viewStruc = new Runnable() { @Override @@ -894,37 +1227,12 @@ public class StructureChooser extends GStructureChooser if (currentView == VIEWS_FILTER) { - int pdbIdColIndex = restable.getColumn("PDB Id") - .getModelIndex(); - int refSeqColIndex = restable.getColumn("Ref Sequence") - .getModelIndex(); int[] selectedRows = restable.getSelectedRows(); PDBEntry[] pdbEntriesToView = new PDBEntry[selectedRows.length]; - int count = 0; List selectedSeqsToView = new ArrayList<>(); - for (int row : selectedRows) - { - String pdbIdStr = restable - .getValueAt(row, pdbIdColIndex).toString(); - SequenceI selectedSeq = (SequenceI) restable - .getValueAt(row, refSeqColIndex); - selectedSeqsToView.add(selectedSeq); - PDBEntry pdbEntry = selectedSeq.getPDBEntry(pdbIdStr); - if (pdbEntry == null) - { - pdbEntry = getFindEntry(pdbIdStr, - selectedSeq.getAllPDBEntries()); - } + pdbEntriesToView = data.collectSelectedRows(restable, + selectedRows, selectedSeqsToView); - if (pdbEntry == null) - { - pdbEntry = new PDBEntry(); - pdbEntry.setId(pdbIdStr); - pdbEntry.setType(PDBEntry.Type.PDB); - selectedSeq.getDatasetSequence().addPDBId(pdbEntry); - } - pdbEntriesToView[count++] = pdbEntry; - } SequenceI[] selectedSeqs = selectedSeqsToView .toArray(new SequenceI[selectedSeqsToView.size()]); sViewer = launchStructureViewer(ssm, pdbEntriesToView, ap, @@ -942,8 +1250,9 @@ public class StructureChooser extends GStructureChooser List selectedSeqsToView = new ArrayList<>(); for (int row : selectedRows) { - PDBEntry pdbEntry = (PDBEntry) tbl_local_pdb.getValueAt(row, - pdbIdColIndex); + PDBEntry pdbEntry = ((PDBEntryTableModel) tbl_local_pdb + .getModel()).getPDBEntryAt(row).getPdbEntry(); + pdbEntriesToView[count++] = pdbEntry; SequenceI selectedSeq = (SequenceI) tbl_local_pdb .getValueAt(row, refSeqColIndex); @@ -970,7 +1279,8 @@ public class StructureChooser extends GStructureChooser if (pdbIdStr.split(":").length > 1) { pdbEntry.setId(pdbIdStr.split(":")[0]); - pdbEntry.setChainCode(pdbIdStr.split(":")[1].toUpperCase()); + pdbEntry.setChainCode( + pdbIdStr.split(":")[1].toUpperCase(Locale.ROOT)); } else { @@ -987,28 +1297,28 @@ public class StructureChooser extends GStructureChooser } else if (currentView == VIEWS_FROM_FILE) { - SequenceI userSelectedSeq = ((AssociateSeqOptions) fileChooserAssSeqPanel - .getCmb_assSeq().getSelectedItem()).getSequence(); + StructureChooser sc = StructureChooser.this; + TFType tft = (TFType) sc.combo_tempFacAs.getSelectedItem(); + String paeFilename = sc.localPdbPaeMatrixFileName; + AssociateSeqOptions assSeqOpt = (AssociateSeqOptions) fileChooserAssSeqPanel + .getCmb_assSeq().getSelectedItem(); + SequenceI userSelectedSeq = assSeqOpt.getSequence(); if (userSelectedSeq != null) { selectedSequence = userSelectedSeq; } - PDBEntry fileEntry = new AssociatePdbFileWithSeq() - .associatePdbWithSeq(selectedPdbFileName, - DataSourceType.FILE, selectedSequence, true, - Desktop.instance); - - sViewer = launchStructureViewer( - ssm, new PDBEntry[] - { fileEntry }, ap, - new SequenceI[] - { selectedSequence }); + String pdbFilename = selectedPdbFileName; + + StructureChooser.openStructureFileForSequence(ssm, sc, ap, + selectedSequence, true, pdbFilename, tft, paeFilename, + true); } SwingUtilities.invokeLater(new Runnable() { @Override public void run() { + setProgressBar("Complete.", progress.hashCode()); closeAction(preferredHeight); mainFrame.dispose(); } @@ -1034,21 +1344,6 @@ public class StructureChooser extends GStructureChooser } } - private PDBEntry getFindEntry(String id, Vector pdbEntries) - { - Objects.requireNonNull(id); - Objects.requireNonNull(pdbEntries); - PDBEntry foundEntry = null; - for (PDBEntry entry : pdbEntries) - { - if (entry.getId().equalsIgnoreCase(id)) - { - return entry; - } - } - return foundEntry; - } - /** * Answers a structure viewer (new or existing) configured to superimpose * added structures or not according to the user's choice @@ -1056,8 +1351,7 @@ public class StructureChooser extends GStructureChooser * @param ssm * @return */ - StructureViewer getTargetedStructureViewer( - StructureSelectionManager ssm) + StructureViewer getTargetedStructureViewer(StructureSelectionManager ssm) { Object sv = targetView.getSelectedItem(); @@ -1074,10 +1368,18 @@ public class StructureChooser extends GStructureChooser * @return */ private StructureViewer launchStructureViewer( - StructureSelectionManager ssm, - final PDBEntry[] pdbEntriesToView, + StructureSelectionManager ssm, final PDBEntry[] pdbEntriesToView, final AlignmentPanel alignPanel, SequenceI[] sequences) { + return launchStructureViewer(ssm, pdbEntriesToView, alignPanel, + sequences, null); + } + + private StructureViewer launchStructureViewer( + StructureSelectionManager ssm, final PDBEntry[] pdbEntriesToView, + final AlignmentPanel alignPanel, SequenceI[] sequences, + ViewerType viewerType) + { long progressId = sequences.hashCode(); setProgressBar(MessageManager .getString("status.launching_3d_structure_viewer"), progressId); @@ -1138,17 +1440,21 @@ public class StructureChooser extends GStructureChooser } if (pdbEntriesToView.length > 1) { - setProgressBar(MessageManager.getString( - "status.fetching_3d_structures_for_selected_entries"), + setProgressBar( + MessageManager.getString( + "status.fetching_3d_structures_for_selected_entries"), progressId); - theViewer.viewStructures(pdbEntriesToView, sequences, alignPanel); + theViewer.viewStructures(pdbEntriesToView, sequences, alignPanel, + viewerType); } else { setProgressBar(MessageManager.formatMessage( "status.fetching_3d_structures_for", - pdbEntriesToView[0].getId()),progressId); - theViewer.viewStructures(pdbEntriesToView[0], sequences, alignPanel); + pdbEntriesToView[0].getId()), progressId); + // Can we pass a pre-computeMappinged pdbFile? + theViewer.viewStructures(pdbEntriesToView[0], sequences, alignPanel, + viewerType); } setProgressBar(null, progressId); // remember the last viewer we used... @@ -1192,52 +1498,63 @@ public class StructureChooser extends GStructureChooser && !discoveredStructuresSet.isEmpty(); } + protected int PDB_ID_MIN = 3;// or: (Jalview.isJS() ? 3 : 1); // Bob proposes + // this. + // Doing a search for "1" or "1c" is valuable? + // Those work but are enormously slow. + @Override protected void txt_search_ActionPerformed() { - new Thread() - { - @Override - public void run() + String text = txt_search.getText().trim(); + if (text.length() >= PDB_ID_MIN) + new Thread() { - errorWarning.setLength(0); - isValidPBDEntry = false; - if (txt_search.getText().length() > 0) + + @Override + public void run() { - String searchTerm = txt_search.getText().toLowerCase(); - searchTerm = searchTerm.split(":")[0]; - // System.out.println(">>>>> search term : " + searchTerm); - List wantedFields = new ArrayList<>(); - FTSRestRequest pdbRequest = new FTSRestRequest(); - pdbRequest.setAllowEmptySeq(false); - pdbRequest.setResponseSize(1); - pdbRequest.setFieldToSearchBy("(pdb_id:"); - pdbRequest.setWantedFields(wantedFields); - pdbRequest.setSearchTerm(searchTerm + ")"); - pdbRequest.setAssociatedSequence(selectedSequence); - pdbRestCleint = PDBFTSRestClient.getInstance(); - wantedFields.add(pdbRestCleint.getPrimaryKeyColumn()); - FTSRestResponse resultList; - try + errorWarning.setLength(0); + isValidPBDEntry = false; + if (text.length() > 0) { - resultList = pdbRestCleint.executeRequest(pdbRequest); - } catch (Exception e) - { - errorWarning.append(e.getMessage()); - return; - } finally - { - validateSelections(); - } - if (resultList.getSearchSummary() != null - && resultList.getSearchSummary().size() > 0) - { - isValidPBDEntry = true; + // TODO move this pdb id search into the PDB specific + // FTSSearchEngine + // for moment, it will work fine as is because it is self-contained + String searchTerm = text.toLowerCase(Locale.ROOT); + searchTerm = searchTerm.split(":")[0]; + // System.out.println(">>>>> search term : " + searchTerm); + List wantedFields = new ArrayList<>(); + FTSRestRequest pdbRequest = new FTSRestRequest(); + pdbRequest.setAllowEmptySeq(false); + pdbRequest.setResponseSize(1); + pdbRequest.setFieldToSearchBy("(pdb_id:"); + pdbRequest.setWantedFields(wantedFields); + pdbRequest.setSearchTerm(searchTerm + ")"); + pdbRequest.setAssociatedSequence(selectedSequence); + FTSRestClientI pdbRestClient = PDBFTSRestClient.getInstance(); + wantedFields.add(pdbRestClient.getPrimaryKeyColumn()); + FTSRestResponse resultList; + try + { + resultList = pdbRestClient.executeRequest(pdbRequest); + } catch (Exception e) + { + errorWarning.append(e.getMessage()); + return; + } finally + { + validateSelections(); + } + if (resultList.getSearchSummary() != null + && resultList.getSearchSummary().size() > 0) + { + isValidPBDEntry = true; + } } + validateSelections(); } - validateSelections(); - } - }.start(); + }.start(); } @Override @@ -1245,15 +1562,19 @@ public class StructureChooser extends GStructureChooser { if (selectedSequences != null) { + lbl_loading.setVisible(true); Thread refreshThread = new Thread(new Runnable() { @Override public void run() { fetchStructuresMetaData(); + // populateFilterComboBox(true, cachedPDBExists); + filterResultSet( ((FilterOption) cmb_filterOption.getSelectedItem()) .getValue()); + lbl_loading.setVisible(false); } }); refreshThread.start(); @@ -1307,7 +1628,7 @@ public class StructureChooser extends GStructureChooser value = entry.getSequence(); break; case 1: - value = entry.getPdbEntry(); + value = entry.getQualifiedId(); break; case 2: value = entry.getPdbEntry().getChainCode() == null ? "_" @@ -1348,6 +1669,15 @@ public class StructureChooser extends GStructureChooser this.pdbEntry = pdbEntry; } + public String getQualifiedId() + { + if (pdbEntry.hasProvider()) + { + return pdbEntry.getProvider() + ":" + pdbEntry.getId(); + } + return pdbEntry.toString(); + } + public SequenceI getSequence() { return sequence; @@ -1365,23 +1695,136 @@ public class StructureChooser extends GStructureChooser @Override public void setProgressBar(String message, long id) { - progressBar.setProgressBar(message, id); + if (!Platform.isHeadless() && progressBar != null) + progressBar.setProgressBar(message, id); } @Override public void registerHandler(long id, IProgressIndicatorHandler handler) { - progressBar.registerHandler(id, handler); + if (progressBar != null) + progressBar.registerHandler(id, handler); } @Override public boolean operationInProgress() { - return progressBar.operationInProgress(); + return progressBar == null ? false : progressBar.operationInProgress(); } public JalviewStructureDisplayI getOpenedStructureViewer() { return sViewer == null ? null : sViewer.sview; } + + @Override + protected void setFTSDocFieldPrefs(FTSDataColumnPreferences newPrefs) + { + data.setDocFieldPrefs(newPrefs); + + } + + /** + * + * @return true when all initialisation threads have finished and dialog is + * visible + */ + public boolean isDialogVisible() + { + return mainFrame != null && data != null && cmb_filterOption != null + && mainFrame.isVisible() + && cmb_filterOption.getSelectedItem() != null; + } + + /** + * + * @return true if the 3D-Beacons query button will/has been displayed + */ + public boolean isCanQueryTDB() + { + return canQueryTDB; + } + + public boolean isNotQueriedTDBYet() + { + return notQueriedTDBYet; + } + + /** + * Open a single structure file for a given sequence + */ + public static void openStructureFileForSequence( + StructureSelectionManager ssm, StructureChooser sc, + AlignmentPanel ap, SequenceI seq, boolean prompt, + String sFilename, TFType tft, String paeFilename, + boolean doXferSettings) + { + openStructureFileForSequence(ssm, sc, ap, seq, prompt, sFilename, tft, + paeFilename, false, true, doXferSettings, null); + } + + public static StructureViewer openStructureFileForSequence( + StructureSelectionManager ssm, StructureChooser sc, + AlignmentPanel ap, SequenceI seq, boolean prompt, + String sFilename, TFType tft, String paeFilename, + boolean forceHeadless, boolean showRefAnnotations, + boolean doXferSettings, ViewerType viewerType) + { + StructureViewer sv = null; + boolean headless = forceHeadless; + if (sc == null) + { + // headless = true; + prompt = false; + sc = new StructureChooser(new SequenceI[] { seq }, seq, ap, false); + } + if (ssm == null) + { + ssm = ap.getStructureSelectionManager(); + } + + PDBEntry fileEntry = new AssociatePdbFileWithSeq().associatePdbWithSeq( + sFilename, DataSourceType.FILE, seq, prompt, Desktop.instance, + tft, paeFilename, doXferSettings); + + // if headless, "false" in the sc constructor above will avoid GUI behaviour + // in sc.launchStructureViewer() + if (!headless && !(viewerType == null)) + { + sv = sc.launchStructureViewer(ssm, new PDBEntry[] { fileEntry }, ap, + new SequenceI[] + { seq }, viewerType); + } + + sc.mainFrame.dispose(); + + if (showRefAnnotations) + { + showReferenceAnnotationsForSequence(ap.alignFrame, seq); + } + + return sv; + } + + public static void showReferenceAnnotationsForSequence(AlignFrame af, + SequenceI sequence) + { + AlignViewport av = af.getCurrentView(); + AlignmentI al = av.getAlignment(); + + List forSequences = new ArrayList<>(); + forSequences.add(sequence); + final Map> candidates = new LinkedHashMap<>(); + AlignmentUtils.findAddableReferenceAnnotations(forSequences, null, + candidates, al); + final SequenceGroup selectionGroup = av.getSelectionGroup(); + AlignmentUtils.addReferenceAnnotations(candidates, al, selectionGroup); + for (AlignmentViewPanel ap : af.getAlignPanels()) + { + // required to readjust the height and position of the PAE + // annotation + ap.adjustAnnotationHeight(); + } + + } }