Warning
This document is for an old release of Galaxy. You can alternatively view this page in the latest release if it exists or view the top of the latest release's documentation.
Source code for galaxy.tools.data
"""
Manage tool data tables, which store (at the application level) data that is
used by tools, for example in the generation of dynamic options. Tables are
loaded and stored by names which tools use to refer to them. This allows
users to configure data tables for a local Galaxy instance without needing
to modify the tool configurations.
"""
import errno
import hashlib
import logging
import os
import os.path
import re
import string
import time
from glob import glob
from tempfile import NamedTemporaryFile
from typing import List
import refgenconf
import requests
from galaxy import util
from galaxy.util import RW_R__R__
from galaxy.util.dictifiable import Dictifiable
from galaxy.util.renamed_temporary_file import RenamedTemporaryFile
from galaxy.util.template import fill_template
log = logging.getLogger(__name__)
DEFAULT_TABLE_TYPE = 'tabular'
TOOL_DATA_TABLE_CONF_XML = """<?xml version="1.0"?>
<tables>
</tables>
"""
[docs]class ToolDataPathFiles:
[docs] def __init__(self, tool_data_path):
self.tool_data_path = os.path.abspath(tool_data_path)
self.update_time = 0
@property
def tool_data_path_files(self):
if time.time() - self.update_time > 1:
self.update_files()
return self._tool_data_path_files
[docs] def update_files(self):
try:
content = os.walk(self.tool_data_path)
self._tool_data_path_files = set(filter(os.path.exists, [os.path.join(dirpath, fn) for dirpath, _, fn_list in content for fn in fn_list if fn and fn.endswith('.loc') or fn.endswith('.loc.sample')]))
self.update_time = time.time()
except Exception:
log.exception()
self._tool_data_path_files = set()
[docs] def exists(self, path):
path = os.path.abspath(path)
if path in self.tool_data_path_files:
return True
else:
return os.path.exists(path)
[docs]class ToolDataTableManager:
"""Manages a collection of tool data tables"""
[docs] def __init__(self, tool_data_path, config_filename=None, tool_data_table_config_path_set=None, other_config_dict=None):
self.tool_data_path = tool_data_path
# This stores all defined data table entries from both the tool_data_table_conf.xml file and the shed_tool_data_table_conf.xml file
# at server startup. If tool shed repositories are installed that contain a valid file named tool_data_table_conf.xml.sample, entries
# from that file are inserted into this dict at the time of installation.
self.data_tables = {}
self.tool_data_path_files = ToolDataPathFiles(self.tool_data_path)
self.other_config_dict = other_config_dict or {}
for single_config_filename in util.listify(config_filename):
if not single_config_filename:
continue
self.load_from_config_file(single_config_filename, self.tool_data_path, from_shed_config=False)
def __getitem__(self, key):
return self.data_tables.__getitem__(key)
def __setitem__(self, key, value):
return self.data_tables.__setitem__(key, value)
def __contains__(self, key):
return self.data_tables.__contains__(key)
[docs] def load_from_config_file(self, config_filename, tool_data_path, from_shed_config=False):
"""
This method is called under 3 conditions:
1. When the ToolDataTableManager is initialized (see __init__ above).
2. Just after the ToolDataTableManager is initialized and the additional entries defined by shed_tool_data_table_conf.xml
are being loaded into the ToolDataTableManager.data_tables.
3. When a tool shed repository that includes a tool_data_table_conf.xml.sample file is being installed into a local
Galaxy instance. In this case, we have 2 entry types to handle, files whose root tag is <tables>, for example:
"""
table_elems = []
if not isinstance(config_filename, list):
config_filename = [config_filename]
for filename in config_filename:
tree = util.parse_xml(filename)
root = tree.getroot()
for table_elem in root.findall('table'):
table = ToolDataTable.from_elem(table_elem, tool_data_path, from_shed_config, filename=filename, tool_data_path_files=self.tool_data_path_files, other_config_dict=self.other_config_dict)
table_elems.append(table_elem)
if table.name not in self.data_tables:
self.data_tables[table.name] = table
log.debug("Loaded tool data table '%s' from file '%s'", table.name, filename)
else:
log.debug("Loading another instance of data table '%s' from file '%s', attempting to merge content.", table.name, filename)
self.data_tables[table.name].merge_tool_data_table(table, allow_duplicates=False) # only merge content, do not persist to disk, do not allow duplicate rows when merging
# FIXME: This does not account for an entry with the same unique build ID, but a different path.
return table_elems
[docs] def add_new_entries_from_config_file(self, config_filename, tool_data_path, shed_tool_data_table_config, persist=False):
"""
This method is called when a tool shed repository that includes a tool_data_table_conf.xml.sample file is being
installed into a local galaxy instance. We have 2 cases to handle, files whose root tag is <tables>, for example::
<tables>
<!-- Location of Tmap files -->
<table name="tmap_indexes" comment_char="#">
<columns>value, dbkey, name, path</columns>
<file path="tool-data/tmap_index.loc" />
</table>
</tables>
and files whose root tag is <table>, for example::
<!-- Location of Tmap files -->
<table name="tmap_indexes" comment_char="#">
<columns>value, dbkey, name, path</columns>
<file path="tool-data/tmap_index.loc" />
</table>
"""
error_message = ''
try:
table_elems = self.load_from_config_file(config_filename=config_filename,
tool_data_path=tool_data_path,
from_shed_config=True)
except Exception as e:
error_message = 'Error attempting to parse file {}: {}'.format(str(os.path.split(config_filename)[1]), util.unicodify(e))
log.debug(error_message, exc_info=True)
table_elems = []
if persist:
# Persist Galaxy's version of the changed tool_data_table_conf.xml file.
self.to_xml_file(shed_tool_data_table_config, table_elems)
return table_elems, error_message
[docs] def to_xml_file(self, shed_tool_data_table_config, new_elems=None, remove_elems=None):
"""
Write the current in-memory version of the shed_tool_data_table_conf.xml file to disk.
remove_elems are removed before new_elems are added.
"""
if not (new_elems or remove_elems):
log.debug('ToolDataTableManager.to_xml_file called without any elements to add or remove.')
return # no changes provided, no need to persist any changes
if not new_elems:
new_elems = []
if not remove_elems:
remove_elems = []
full_path = os.path.abspath(shed_tool_data_table_config)
# FIXME: we should lock changing this file by other threads / head nodes
try:
try:
tree = util.parse_xml(full_path)
except OSError as e:
if e.errno == errno.ENOENT:
with open(full_path, 'w') as fh:
fh.write(TOOL_DATA_TABLE_CONF_XML)
tree = util.parse_xml(full_path)
else:
raise
root = tree.getroot()
out_elems = [elem for elem in root]
except Exception as e:
out_elems = []
log.debug('Could not parse existing tool data table config, assume no existing elements: %s', e)
for elem in remove_elems:
# handle multiple occurrences of remove elem in existing elems
while elem in out_elems:
remove_elems.remove(elem)
# add new elems
out_elems.extend(new_elems)
out_path_is_new = not os.path.exists(full_path)
root = util.parse_xml_string('<?xml version="1.0"?>\n<tables></tables>')
for elem in out_elems:
root.append(elem)
with RenamedTemporaryFile(full_path, mode='w') as out:
out.write(util.xml_to_string(root, pretty=True))
os.chmod(full_path, RW_R__R__)
if out_path_is_new:
self.tool_data_path_files.update_files()
[docs] def reload_tables(self, table_names=None, path=None):
"""
Reload tool data tables. If neither table_names nor path is given, reloads all tool data tables.
"""
tables = self.get_tables()
if not table_names:
if path:
table_names = self.get_table_names_by_path(path)
else:
table_names = list(tables.keys())
elif not isinstance(table_names, list):
table_names = [table_names]
for table_name in table_names:
tables[table_name].reload_from_files()
log.debug("Reloaded tool data table '%s' from files.", table_name)
return table_names
[docs] def get_table_names_by_path(self, path):
"""Returns a list of table names given a path"""
table_names = set()
for name, data_table in self.data_tables.items():
if path in data_table.filenames:
table_names.add(name)
return list(table_names)
[docs]class ToolDataTable:
type_key: str
[docs] @classmethod
def from_elem(cls, table_elem, tool_data_path, from_shed_config, filename, tool_data_path_files, other_config_dict=None):
table_type = table_elem.get('type', 'tabular')
assert table_type in tool_data_table_types, "Unknown data table type '%s'" % table_type
return tool_data_table_types[table_type](table_elem, tool_data_path, from_shed_config=from_shed_config, filename=filename, tool_data_path_files=tool_data_path_files, other_config_dict=other_config_dict)
[docs] def __init__(self, config_element, tool_data_path, from_shed_config=False, filename=None, tool_data_path_files=None, other_config_dict=None):
self.name = config_element.get('name')
self.comment_char = config_element.get('comment_char')
self.empty_field_value = config_element.get('empty_field_value', '')
self.empty_field_values = {}
self.allow_duplicate_entries = util.asbool(config_element.get('allow_duplicate_entries', True))
self.here = filename and os.path.dirname(filename)
self.filenames = {}
self.tool_data_path = tool_data_path
self.tool_data_path_files = tool_data_path_files
self.other_config_dict = other_config_dict or {}
self.missing_index_file = None
# increment this variable any time a new entry is added, or when the table is totally reloaded
# This value has no external meaning, and does not represent an abstract version of the underlying data
self._loaded_content_version = 1
self._load_info = ([config_element, tool_data_path], {'from_shed_config': from_shed_config, 'tool_data_path_files': self.tool_data_path_files, 'other_config_dict': other_config_dict, 'filename': filename})
self._merged_load_info = []
def _update_version(self, version=None):
if version is not None:
self._loaded_content_version = version
else:
self._loaded_content_version += 1
return self._loaded_content_version
[docs] def get_empty_field_by_name(self, name):
return self.empty_field_values.get(name, self.empty_field_value)
def _add_entry(self, entry, allow_duplicates=True, persist=False, persist_on_error=False, entry_source=None, **kwd):
raise NotImplementedError("Abstract method")
[docs] def add_entry(self, entry, allow_duplicates=True, persist=False, persist_on_error=False, entry_source=None, **kwd):
self._add_entry(entry, allow_duplicates=allow_duplicates, persist=persist, persist_on_error=persist_on_error, entry_source=entry_source, **kwd)
return self._update_version()
[docs] def add_entries(self, entries, allow_duplicates=True, persist=False, persist_on_error=False, entry_source=None, **kwd):
for entry in entries:
self.add_entry(entry, allow_duplicates=allow_duplicates, persist=persist, persist_on_error=persist_on_error, entry_source=entry_source, **kwd)
return self._loaded_content_version
def _remove_entry(self, values):
raise NotImplementedError("Abstract method")
[docs] def is_current_version(self, other_version):
return self._loaded_content_version == other_version
[docs] def merge_tool_data_table(self, other_table, allow_duplicates=True, persist=False, persist_on_error=False, entry_source=None, **kwd):
raise NotImplementedError("Abstract method")
[docs] def reload_from_files(self):
new_version = self._update_version()
merged_info = self._merged_load_info
self.__init__(*self._load_info[0], **self._load_info[1])
self._update_version(version=new_version)
for (tool_data_table_class, load_info) in merged_info:
self.merge_tool_data_table(tool_data_table_class(*load_info[0], **load_info[1]), allow_duplicates=False)
return self._update_version()
[docs]class TabularToolDataTable(ToolDataTable, Dictifiable):
"""
Data stored in a tabular / separated value format on disk, allows multiple
files to be merged but all must have the same column definitions:
.. code-block:: xml
<table type="tabular" name="test">
<column name='...' index = '...' />
<file path="..." />
<file path="..." />
</table>
"""
dict_collection_visible_keys = ['name']
type_key = 'tabular'
[docs] def __init__(self, config_element, tool_data_path, from_shed_config=False, filename=None, tool_data_path_files=None, other_config_dict=None):
super().__init__(config_element, tool_data_path, from_shed_config, filename, tool_data_path_files, other_config_dict=other_config_dict)
self.config_element = config_element
self.data = []
self.configure_and_load(config_element, tool_data_path, from_shed_config)
[docs] def configure_and_load(self, config_element, tool_data_path, from_shed_config=False, url_timeout=10):
"""
Configure and load table from an XML element.
"""
self.separator = config_element.get('separator', '\t')
self.comment_char = config_element.get('comment_char', '#')
# Configure columns
self.parse_column_spec(config_element)
# store repo info if available:
repo_elem = config_element.find('tool_shed_repository')
if repo_elem is not None:
repo_info = dict(tool_shed=repo_elem.find('tool_shed').text, name=repo_elem.find('repository_name').text,
owner=repo_elem.find('repository_owner').text, installed_changeset_revision=repo_elem.find('installed_changeset_revision').text)
else:
repo_info = None
# Read every file
for file_element in config_element.findall('file'):
tmp_file = None
filename = file_element.get('path', None)
if filename is None:
# Handle URLs as files
filename = file_element.get('url', None)
if filename:
tmp_file = NamedTemporaryFile(prefix='TTDT_URL_%s-' % self.name, mode='w')
try:
tmp_file.write(requests.get(filename, timeout=url_timeout).text)
except Exception as e:
log.error('Error loading Data Table URL "%s": %s', filename, e)
continue
log.debug('Loading Data Table URL "%s" as filename "%s".', filename, tmp_file.name)
filename = tmp_file.name
tmp_file.flush()
else:
# Pull the filename from a global config
filename = file_element.get('from_config', None) or None
if filename:
filename = self.other_config_dict.get(filename, None)
filename = file_path = expand_here_template(filename, here=self.here)
found = False
if file_path is None:
log.debug("Encountered a file element (%s) that does not contain a path value when loading tool data table '%s'.", util.xml_to_string(file_element), self.name)
continue
# FIXME: splitting on and merging paths from a configuration file when loading is wonky
# Data should exist on disk in the state needed, i.e. the xml configuration should
# point directly to the desired file to load. Munging of the tool_data_tables_conf.xml.sample
# can be done during installing / testing / metadata resetting with the creation of a proper
# tool_data_tables_conf.xml file, containing correct <file path=> attributes. Allowing a
# path.join with a different root should be allowed, but splitting should not be necessary.
if tool_data_path and from_shed_config:
# Must identify with from_shed_config as well, because the
# regular galaxy app has and uses tool_data_path.
# We're loading a tool in the tool shed, so we cannot use the Galaxy tool-data
# directory which is hard-coded into the tool_data_table_conf.xml entries.
filename = os.path.split(file_path)[1]
filename = os.path.join(tool_data_path, filename)
if self.tool_data_path_files.exists(filename):
found = True
elif not os.path.isabs(filename):
# Since the path attribute can include a hard-coded path to a specific directory
# (e.g., <file path="tool-data/cg_crr_files.loc" />) which may not be the same value
# as self.tool_data_path, we'll parse the path to get the filename and see if it is
# in self.tool_data_path.
file_path, file_name = os.path.split(filename)
if file_path != self.tool_data_path:
corrected_filename = os.path.join(self.tool_data_path, file_name)
if self.tool_data_path_files.exists(corrected_filename):
filename = corrected_filename
found = True
elif not from_shed_config and self.tool_data_path_files.exists("%s.sample" % corrected_filename):
log.info("Could not find tool data %s, reading sample" % corrected_filename)
filename = "%s.sample" % corrected_filename
found = True
errors = []
if found:
self.extend_data_with(filename, errors=errors)
self._update_version()
else:
self.missing_index_file = filename
# TODO: some data tables need to exist (even if they are empty)
# for tools to load. In an installed Galaxy environment and the
# default tool_data_table_conf.xml, this will emit spurious
# warnings about missing location files that would otherwise be
# empty and we don't care about unless the admin chooses to
# populate them.
log.warning(f"Cannot find index file '{filename}' for tool data table '{self.name}'")
if filename not in self.filenames or not self.filenames[filename]['found']:
self.filenames[filename] = dict(found=found, filename=filename, from_shed_config=from_shed_config, tool_data_path=tool_data_path,
config_element=config_element, tool_shed_repository=repo_info, errors=errors)
else:
log.debug("Filename '%s' already exists in filenames (%s), not adding", filename, list(self.filenames.keys()))
# Remove URL tmp file
if tmp_file is not None:
tmp_file.close()
[docs] def merge_tool_data_table(self, other_table, allow_duplicates=True, persist=False, persist_on_error=False, entry_source=None, **kwd):
assert self.columns == other_table.columns, f"Merging tabular data tables with non matching columns is not allowed: {self.name}:{self.columns} != {other_table.name}:{other_table.columns}"
# merge filename info
for filename, info in other_table.filenames.items():
if filename not in self.filenames:
self.filenames[filename] = info
# save info about table
self._merged_load_info.append((other_table.__class__, other_table._load_info))
# If we are merging in a data table that does not allow duplicates, enforce that upon the data table
if self.allow_duplicate_entries and not other_table.allow_duplicate_entries:
log.debug('While attempting to merge tool data table "%s", the other instance of the table specified that duplicate entries are not allowed, now deduplicating all previous entries.', self.name)
self.allow_duplicate_entries = False
self._deduplicate_data()
# add data entries and return current data table version
return self.add_entries(other_table.data, allow_duplicates=allow_duplicates, persist=persist, persist_on_error=persist_on_error, entry_source=entry_source, **kwd)
[docs] def handle_found_index_file(self, filename):
self.missing_index_file = None
self.extend_data_with(filename)
[docs] def get_field(self, value):
rval = None
for i in self.get_named_fields_list():
if i['value'] == value:
rval = TabularToolDataField(i)
return rval
[docs] def get_named_fields_list(self):
rval = []
named_columns = self.get_column_name_list()
for fields in self.get_fields():
field_dict = {}
for i, field in enumerate(fields):
if i == len(named_columns):
break
field_name = named_columns[i]
if field_name is None:
field_name = i # check that this is supposed to be 0 based.
field_dict[field_name] = field
rval.append(field_dict)
return rval
[docs] def parse_column_spec(self, config_element):
"""
Parse column definitions, which can either be a set of 'column' elements
with a name and index (as in dynamic options config), or a shorthand
comma separated list of names in order as the text of a 'column_names'
element.
A column named 'value' is required.
"""
self.columns = {}
if config_element.find('columns') is not None:
column_names = util.xml_text(config_element.find('columns'))
column_names = [n.strip() for n in column_names.split(',')]
for index, name in enumerate(column_names):
self.columns[name] = index
self.largest_index = index
else:
self.largest_index = 0
for column_elem in config_element.findall('column'):
name = column_elem.get('name', None)
assert name is not None, "Required 'name' attribute missing from column def"
index = column_elem.get('index', None)
assert index is not None, "Required 'index' attribute missing from column def"
index = int(index)
self.columns[name] = index
if index > self.largest_index:
self.largest_index = index
empty_field_value = column_elem.get('empty_field_value', None)
if empty_field_value is not None:
self.empty_field_values[name] = empty_field_value
assert 'value' in self.columns, "Required 'value' column missing from column def"
if 'name' not in self.columns:
self.columns['name'] = self.columns['value']
[docs] def extend_data_with(self, filename, errors=None):
here = os.path.dirname(os.path.abspath(filename))
self.data.extend(self.parse_file_fields(filename, errors=errors, here=here))
if not self.allow_duplicate_entries:
self._deduplicate_data()
[docs] def parse_file_fields(self, filename, errors=None, here="__HERE__"):
"""
Parse separated lines from file and return a list of tuples.
TODO: Allow named access to fields using the column names.
"""
separator_char = "<TAB>" if self.separator == "\t" else self.separator
rval = []
with open(filename) as fh:
for i, line in enumerate(fh):
if line.lstrip().startswith(self.comment_char):
continue
line = line.rstrip("\n\r")
if line:
line = expand_here_template(line, here=here)
fields = line.split(self.separator)
if self.largest_index < len(fields):
rval.append(fields)
else:
line_error = "Line %i in tool data table '%s' is invalid (HINT: '%s' characters must be used to separate fields):\n%s" % ((i + 1), self.name, separator_char, line)
if errors is not None:
errors.append(line_error)
log.warning(line_error)
log.debug("Loaded %i lines from '%s' for '%s'", len(rval), filename, self.name)
return rval
[docs] def get_column_name_list(self):
rval = []
for i in range(self.largest_index + 1):
found_column = False
for name, index in self.columns.items():
if index == i:
if not found_column:
rval.append(name)
elif name == 'value':
# the column named 'value' always has priority over other named columns
rval[-1] = name
found_column = True
if not found_column:
rval.append(None)
return rval
[docs] def get_entry(self, query_attr, query_val, return_attr, default=None):
"""
Returns table entry associated with a col/val pair.
"""
rval = self.get_entries(query_attr, query_val, return_attr, default=default, limit=1)
if rval:
return rval[0]
return default
[docs] def get_entries(self, query_attr, query_val, return_attr, default=None, limit=None):
"""
Returns table entry associated with a col/val pair.
"""
query_col = self.columns.get(query_attr, None)
if query_col is None:
return default
if return_attr is not None:
return_col = self.columns.get(return_attr, None)
if return_col is None:
return default
rval = []
# Look for table entry.
for fields in self.get_fields():
if fields[query_col] == query_val:
if return_attr is None:
field_dict = {}
for i, col_name in enumerate(self.get_column_name_list()):
field_dict[col_name or i] = fields[i]
rval.append(field_dict)
else:
rval.append(fields[return_col])
if limit is not None and len(rval) == limit:
break
return rval or default
[docs] def get_filename_for_source(self, source, default=None):
if source:
# if dict, assume is compatible info dict, otherwise call method
if isinstance(source, dict):
source_repo_info = source
else:
source_repo_info = source.get_tool_shed_repository_info_dict()
else:
source_repo_info = None
filename = default
for name, value in self.filenames.items():
repo_info = value.get('tool_shed_repository', None)
if (not source_repo_info and not repo_info) or (source_repo_info and repo_info and source_repo_info == repo_info):
filename = name
break
return filename
def _add_entry(self, entry, allow_duplicates=True, persist=False, persist_on_error=False, entry_source=None, **kwd):
# accepts dict or list of columns
if isinstance(entry, dict):
fields = []
for column_name in self.get_column_name_list():
if column_name not in entry:
log.debug("Using default column value for column '%s' when adding data table entry (%s) to table '%s'.", column_name, entry, self.name)
field_value = self.get_empty_field_by_name(column_name)
else:
field_value = entry[column_name]
fields.append(field_value)
else:
fields = entry
is_error = False
if self.largest_index < len(fields):
fields = self._replace_field_separators(fields)
if (allow_duplicates and self.allow_duplicate_entries) or fields not in self.get_fields():
self.data.append(fields)
else:
log.debug("Attempted to add fields (%s) to data table '%s', but this entry already exists and allow_duplicates is False.", fields, self.name)
is_error = True
else:
log.error("Attempted to add fields (%s) to data table '%s', but there were not enough fields specified ( %i < %i ).", fields, self.name, len(fields), self.largest_index + 1)
is_error = True
filename = None
if persist and (not is_error or persist_on_error):
filename = self.get_filename_for_source(entry_source)
if filename is None:
# should we default to using any filename here instead?
log.error("Unable to determine filename for persisting data table '%s' values: '%s'.", self.name, fields)
is_error = True
else:
# FIXME: Need to lock these files for editing
log.debug("Persisting changes to file: %s", filename)
try:
data_table_fh = open(filename, 'r+b')
except OSError as e:
log.warning('Error opening data table file (%s) with r+b, assuming file does not exist and will open as wb: %s', filename, e)
data_table_fh = open(filename, 'wb')
if os.stat(filename).st_size != 0:
# ensure last existing line ends with new line
data_table_fh.seek(-1, 2) # last char in file
last_char = data_table_fh.read(1)
if last_char not in [b'\n', b'\r']:
data_table_fh.write(b'\n')
fields = "%s\n" % self.separator.join(fields)
data_table_fh.write(fields.encode('utf-8'))
return not is_error
def _remove_entry(self, values):
# update every file
for filename in self.filenames:
if os.path.exists(filename):
values = self._replace_field_separators(values)
self.filter_file_fields(filename, values)
else:
log.warning(f"Cannot find index file '{filename}' for tool data table '{self.name}'")
self.reload_from_files()
[docs] def filter_file_fields(self, loc_file, values):
"""
Reads separated lines from file and print back only the lines that pass a filter.
"""
with open(loc_file) as reader:
rval = ""
for line in reader:
if line.lstrip().startswith(self.comment_char):
rval += line
else:
line_s = line.rstrip("\n\r")
if line_s:
fields = line_s.split(self.separator)
if fields != values:
rval += line
with open(loc_file, 'w') as writer:
writer.write(rval)
return rval
def _replace_field_separators(self, fields, separator=None, replace=None, comment_char=None):
# make sure none of the fields contain separator
# make sure separator replace is different from comment_char,
# due to possible leading replace
if separator is None:
separator = self.separator
if replace is None:
if separator == " ":
if comment_char == "\t":
replace = "_"
else:
replace = "\t"
else:
if comment_char == " ":
replace = "_"
else:
replace = " "
return [x.replace(separator, replace) for x in fields]
def _deduplicate_data(self):
# Remove duplicate entries, without recreating self.data object
dup_lines = []
hash_set = set()
for i, fields in enumerate(self.data):
fields_hash = hash(self.separator.join(fields))
if fields_hash in hash_set:
dup_lines.append(i)
log.debug('Found duplicate entry in tool data table "%s", but duplicates are not allowed, removing additional entry for: "%s"', self.name, fields)
else:
hash_set.add(fields_hash)
for i in reversed(dup_lines):
self.data.pop(i)
@property
def xml_string(self):
return util.xml_to_string(self.config_element)
[docs] def to_dict(self, view='collection'):
rval = super().to_dict()
if view == 'element':
rval['columns'] = sorted(self.columns.keys(), key=lambda x: self.columns[x])
rval['fields'] = self.get_fields()
return rval
[docs]class TabularToolDataField(Dictifiable):
dict_collection_visible_keys: List[str] = []
def __getitem__(self, key):
return self.data[key]
[docs] def get_base_dir(self):
path = self.get_base_path()
if not os.path.isdir(path):
path = os.path.dirname(path)
return path
[docs] def get_filesize_map(self, rm_base_dir=False):
out = {}
for path in self.get_files():
if rm_base_dir:
out[self.clean_base_dir(path)] = os.path.getsize(path)
else:
out[path] = os.path.getsize(path)
return out
[docs] def get_fingerprint(self):
sha1 = hashlib.sha1()
fmap = self.get_filesize_map(True)
for k in sorted(fmap.keys()):
sha1.update(util.smart_str(k))
sha1.update(util.smart_str(fmap[k]))
return sha1.hexdigest()
[docs] def to_dict(self):
rval = super().to_dict()
rval['name'] = self.data['value']
rval['fields'] = self.data
rval['base_dir'] = self.get_base_dir(),
rval['files'] = self.get_filesize_map(True)
rval['fingerprint'] = self.get_fingerprint()
return rval
[docs]class RefgenieToolDataTable(TabularToolDataTable):
"""
Data stored in refgenie
.. code-block:: xml
<table name="all_fasta" type="refgenie" asset="fasta" >
<file path="refgenie.yml" />
<field name="value" template="true">${__REFGENIE_UUID__}</field>
<field name="dbkey" template="true">${__REFGENIE_GENOME__}</field>
<field name="name" template="true">${__REFGENIE_DISPLAY_NAME__}</field>
<field name="path" template="true">${__REFGENIE_ASSET__}</field>
</table>
"""
dict_collection_visible_keys = ['name']
type_key = 'refgenie'
[docs] def __init__(self, config_element, tool_data_path, from_shed_config=False, filename=None, tool_data_path_files=None, other_config_dict=None):
super().__init__(config_element, tool_data_path, from_shed_config, filename, tool_data_path_files, other_config_dict=other_config_dict)
self.config_element = config_element
self.data = []
self.configure_and_load(config_element, tool_data_path, from_shed_config)
[docs] def configure_and_load(self, config_element, tool_data_path, from_shed_config=False, url_timeout=10):
self.rg_asset = config_element.get('asset', None)
assert self.rg_asset, ValueError('You must specify an asset attribute.')
super().configure_and_load(config_element, tool_data_path, from_shed_config=from_shed_config, url_timeout=url_timeout)
[docs] def parse_column_spec(self, config_element):
self.columns = {}
self.key_map = {}
self.template_for_column = {}
self.strip_for_column = {}
self.largest_index = 0
for i, elem in enumerate(config_element.findall('field')):
name = elem.get('name', None)
assert name, ValueError('You must provide a name refgenie field element.')
value = elem.text
self.key_map[name] = value
column_index = int(elem.get('column_index', i))
empty_field_value = elem.get('empty_field_value', None)
if empty_field_value is not None:
self.empty_field_values[name] = empty_field_value
self.template_for_column[name] = util.asbool(elem.get('template', False))
self.strip_for_column[name] = util.asbool(elem.get('strip', False))
self.columns[name] = column_index
self.largest_index = max(self.largest_index, column_index)
if 'name' not in self.columns:
self.columns['name'] = self.columns['value']
[docs] def parse_file_fields(self, filename, errors=None, here="__HERE__"):
rgc = refgenconf.RefGenConf(filename)
rval = []
for genome in rgc.list_genomes_by_asset(self.rg_asset):
genome_attributes = rgc.get_genome_attributes(genome)
description = genome_attributes.get('description', None)
asset_list = rgc.list(genome, include_tags=True)[genome]
for tagged_asset in asset_list:
asset, tag = tagged_asset.rsplit(':', 1)
if asset != self.rg_asset:
continue
digest = rgc.id(genome, asset, tag=tag)
uuid = f'refgenie:{genome}/{self.rg_asset}:{tag}@{digest}'
display_name = description or f'{genome}/{tagged_asset}'
def _seek_key(key):
return rgc.seek(genome, asset, tag_name=tag, seek_key=key)
template_dict = {
'__REFGENIE_UUID__': uuid,
'__REFGENIE_GENOME__': genome,
'__REFGENIE_TAG__': tag,
'__REFGENIE_DISPLAY_NAME__': display_name,
'__REFGENIE_ASSET__': rgc.seek(genome, asset, tag_name=tag),
'__REFGENIE_ASSET_NAME__': asset,
'__REFGENIE_DIGEST__': digest,
'__REFGENIE_GENOME_ATTRIBUTES__': genome_attributes,
'__REFGENIE__': rgc,
'__REFGENIE_SEEK_KEY__': _seek_key,
}
fields = [''] * (self.largest_index + 1)
for name, index in self.columns.items():
rg_value = self.key_map[name]
# Default is hard-coded value
if self.template_for_column.get(name, False):
rg_value = fill_template(rg_value, template_dict)
if self.strip_for_column.get(name, False):
rg_value = rg_value.strip()
fields[index] = rg_value
rval.append(fields)
log.debug("Loaded %i entries from refgenie '%s' asset '%s' for '%s'", len(rval), filename, self.rg_asset, self.name)
return rval
def _remove_entry(self, values):
log.warning("Deletion from refgenie-backed '%s' data table is not supported, will only try to delete from .loc files", self.name)
# Update every non-refgenie files
super()._remove_entry(values)
[docs]def expand_here_template(content, here=None):
if here and content:
content = string.Template(content).safe_substitute({"__HERE__": here})
return content
# Registry of tool data types by type_key
tool_data_table_types = {cls.type_key: cls for cls in [TabularToolDataTable, RefgenieToolDataTable]}