Mercurial > repos > galaxyp > msconvert

annotate msconvert_win/msconvert_wrapper.py @ 6:6b6bba73eadb draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
author galaxyp
date Tue, 14 Mar 2017 16:52:39 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
6
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
1 #!/usr/bin/env python
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
2 import optparse
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
3 import os
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
4 import sys
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
5 import tempfile
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
6 import shutil
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
7 import subprocess
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
8 import re
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
9 import logging
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
10 import shlex
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
11
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
12 assert sys.version_info[:2] >= (2, 6)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
13
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
14 log = logging.getLogger(__name__)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
15 working_directory = os.getcwd()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
16 tmp_stderr_name = tempfile.NamedTemporaryFile(dir=working_directory, suffix='.stderr').name
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
17 tmp_stdout_name = tempfile.NamedTemporaryFile(dir=working_directory, suffix='.stdout').name
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
18
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
19
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
20 def stop_err(msg):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
21 sys.stderr.write("%s\n" % msg)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
22 sys.exit()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
23
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
24
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
25 def read_stderr():
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
26 stderr = ''
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
27 if(os.path.exists(tmp_stderr_name)):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
28 with open(tmp_stderr_name, 'rb') as tmp_stderr:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
29 buffsize = 1048576
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
30 try:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
31 while True:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
32 stderr += tmp_stderr.read(buffsize)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
33 if not stderr or len(stderr) % buffsize != 0:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
34 break
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
35 except OverflowError:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
36 pass
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
37 return stderr
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
38
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
39
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
40 def execute(command, stdin=None):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
41 try:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
42 with open(tmp_stderr_name, 'wb') as tmp_stderr:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
43 with open(tmp_stdout_name, 'wb') as tmp_stdout:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
44 args = shlex.split(command) # handle proper splitting of quoted args
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
45 proc = subprocess.Popen(args=args, shell=False, stderr=tmp_stderr.fileno(), stdout=tmp_stdout.fileno(), stdin=stdin, env=os.environ)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
46 returncode = proc.wait()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
47 if returncode != 0:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
48 raise Exception("Program returned with non-zero exit code %d. stderr: %s" % (returncode, read_stderr()))
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
49 finally:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
50 print(( open(tmp_stderr_name, "r").read() ))
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
51 print(( open(tmp_stdout_name, "r").read() ))
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
52
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
53
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
54 def delete_file(path):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
55 if os.path.exists(path):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
56 try:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
57 os.remove(path)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
58 except:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
59 pass
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
60
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
61
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
62 def delete_directory(directory):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
63 if os.path.exists(directory):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
64 try:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
65 shutil.rmtree(directory)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
66 except:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
67 pass
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
68
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
69
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
70 def symlink(source, link_name):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
71 import platform
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
72 if platform.system() == 'Windows':
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
73 try:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
74 import win32file
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
75 win32file.CreateSymbolicLink(source, link_name, 1)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
76 except:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
77 shutil.copy(source, link_name)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
78 else:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
79 os.symlink(source, link_name)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
80
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
81
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
82 def copy_to_working_directory(data_file, relative_path):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
83 if os.path.abspath(data_file) != os.path.abspath(relative_path):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
84 symlink(data_file, relative_path)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
85 return relative_path
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
86
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
87
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
88 def __main__():
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
89 run_script()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
90
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
91 #ENDTEMPLATE
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
92
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
93 to_extensions = ['mzML', 'mzXML', 'unindexed_mzML', 'unindexed_mzXML', 'mgf', 'mz5', 'txt', 'ms2', 'cms2']
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
94
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
95
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
96 def str_to_bool(v):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
97 """ From http://stackoverflow.com/questions/715417/converting-from-a-string-to-boolean-in-python """
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
98 return v.lower() in ["yes", "true", "t", "1"]
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
99
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
100
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
101 def _add_filter(filters_file, contents):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
102 filters_file.write("filter=%s\n" % contents)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
103
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
104
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
105 def _skip_line(options, file_num, line_parts):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
106 file_num_column = options.filter_table_file_column
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
107 if not file_num_column:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
108 return False
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
109 else:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
110 target_file_num_val = str(file_num).strip()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
111 query_file_num_val = line_parts[int(file_num_column) - 1].strip()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
112 #print "target %s, query %s" % (target_file_num_val, query_file_num_val)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
113 return target_file_num_val != query_file_num_val
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
114
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
115
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
116 def _read_table_numbers(path, options, file_num=None):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
117 unique_numbers = set([])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
118 column_num = options.filter_table_column
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
119 input = open(path, "r")
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
120 first_line = True
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
121 for line in input:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
122 if not line:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
123 continue
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
124 line = line.strip()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
125 if line.startswith("#"):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
126 first_line = False
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
127 continue
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
128 if column_num == None:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
129 column = line
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
130 else:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
131 line_parts = line.split("\t")
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
132 if _skip_line(options, file_num, line_parts):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
133 continue
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
134 column = line_parts[int(column_num) - 1]
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
135 match = re.match("\d+", column)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
136 if match:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
137 unique_numbers.add(int(match.group()))
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
138 first_line = False
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
139 return unique_numbers
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
140
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
141
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
142 def shellquote(s):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
143 return '"' + s.replace('"', '\\"') + '"'
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
144
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
145
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
146 def _add_filter_line_from_file(filter_file, options, file_num=None):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
147 file = options.filter_table
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
148 if not file:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
149 return
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
150 numbers = _read_table_numbers(file, options, file_num)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
151 msconvert_int_set = " ".join([str(number) for number in numbers])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
152 filter_type = options.filter_table_type
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
153 if filter_type == 'number':
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
154 filter_prefix = 'scanNumber'
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
155 else:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
156 filter_prefix = 'index'
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
157 _add_filter(filter_file, "%s %s" % (filter_prefix, msconvert_int_set))
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
158
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
159
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
160 def _create_filters_file(options, file_num=None, debug=False):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
161 suffix = "" if not file_num else str(file_num)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
162 filters_file_path = "filters%s" % suffix
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
163 filters_file = open(filters_file_path, "w")
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
164 if options.filters_file:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
165 filters_file.write(open(options.filters_file, "r").read())
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
166 for filter in options.filter:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
167 _add_filter(filters_file, filter)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
168 _add_filter_line_from_file(filters_file, options, file_num=file_num)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
169
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
170 filters_file.close()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
171 if debug:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
172 print(( open(filters_file_path, "r").read() ))
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
173 return filters_file_path
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
174
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
175
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
176 def _build_base_cmd(options,args=None):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
177 to_extension = options.toextension
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
178 if to_extension.startswith("unindexed_"):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
179 to_extension = to_extension[len("unindexed_"):]
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
180 to_params = "--noindex"
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
181 else:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
182 to_params = ""
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
183 cmd = "msconvert --%s %s" % (to_extension, to_params)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
184 if args:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
185 cmd = "%s %s" % (cmd, ' '.join(args))
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
186 if str_to_bool(options.zlib):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
187 cmd = "%s %s" % (cmd, "--zlib")
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
188 if options.binaryencoding:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
189 cmd = "%s --%s" % (cmd, options.binaryencoding)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
190 if options.mzencoding:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
191 cmd = "%s --mz%s" % (cmd, options.mzencoding)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
192 if options.intensityencoding:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
193 cmd = "%s --inten%s" % (cmd, options.intensityencoding)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
194 return cmd
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
195
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
196
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
197 def _run(base_cmd, output_dir='output', inputs=[], debug=False):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
198 inputs_as_str = " ".join(['%s' % shellquote(input) for input in inputs])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
199 os.mkdir(output_dir)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
200 cmd = "%s -o %s %s" % (base_cmd, shellquote(output_dir), inputs_as_str)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
201 if debug:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
202 print(cmd)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
203 execute(cmd)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
204 output_files = os.listdir(output_dir)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
205 assert len(output_files) == 1
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
206 output_file = output_files[0]
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
207 return os.path.join(output_dir, output_file)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
208
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
209
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
210 def run_script():
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
211 parser = optparse.OptionParser()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
212 parser.add_option('--input', dest='inputs', action='append', default=[])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
213 parser.add_option('--input_name', dest='input_names', action='append', default=[])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
214 parser.add_option('--implicit', dest='implicits', action='append', default=[], help='input files that should NOT be on the msconvert command line.')
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
215 parser.add_option('--ident', dest='idents', action='append', default=[])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
216 parser.add_option('--ident_name', dest='ident_names', action='append', default=[])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
217 parser.add_option('--output', dest='output')
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
218 parser.add_option('--refinement', dest='refinement')
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
219 parser.add_option('--fromextension', dest='fromextension')
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
220 parser.add_option('--toextension', dest='toextension', default='mzML', choices=to_extensions)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
221 parser.add_option('--binaryencoding', dest='binaryencoding', choices=['32', '64'])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
222 parser.add_option('--mzencoding', dest='mzencoding', choices=['32', '64'])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
223 parser.add_option('--intensityencoding', dest='intensityencoding', choices=['32', '64'])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
224 parser.add_option('--zlib', dest='zlib', default="false")
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
225 parser.add_option('--filter', dest='filter', action='append', default=[])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
226 parser.add_option('--filters_file', dest='filters_file', default=None)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
227 parser.add_option('--filter_table', default=None)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
228 parser.add_option('--filter_table_type', default='index', choices=['index', 'number'])
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
229 parser.add_option('--filter_table_column', default=None)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
230 parser.add_option('--filter_table_file_column', default=None)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
231 parser.add_option('--debug', dest='debug', action='store_true', default=False)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
232
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
233 (options, args) = parser.parse_args()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
234 if len(options.inputs) < 1:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
235 stop_err("No input files to msconvert specified")
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
236 if len(options.input_names) > 0 and len(options.input_names) != len(options.inputs):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
237 stop_err("Number(s) of supplied input names and input files do not match")
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
238 if not options.output:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
239 stop_err("Must specify output location")
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
240 input_files = []
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
241 for i, input in enumerate(options.inputs):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
242 input_base = None
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
243 if len(options.input_names) > i:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
244 input_base = options.input_names[i]
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
245 input_base = input_base.replace("'", "").replace("\"", "")
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
246 print("1- input_base: %s" % input_base)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
247 if not input_base:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
248 input_base = 'input%s' % i
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
249 print("2- input_base: %s" % input_base)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
250 if not input_base.lower().endswith('.%s' % options.fromextension.lower()) and input not in options.implicits:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
251 input_file = '%s.%s' % (input_base, options.fromextension)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
252 print("3- input_base: %s" % input_base)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
253 print("3- input_file: %s" % input_file)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
254 else:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
255 input_file = input_base
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
256 print("4- input_base: %s" % input_base)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
257 print("4- input_file: %s" % input_file)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
258 input_file = input_file
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
259 copy_to_working_directory(input, input_file)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
260 if input in options.implicits:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
261 continue
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
262 input_files.append(input_file)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
263 for i, ident in enumerate(options.idents):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
264 ident_file = options.ident_names[i]
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
265 copy_to_working_directory(ident, ident_file)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
266
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
267 cmd = _build_base_cmd(options,args=args)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
268 file_column = options.filter_table_file_column
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
269 if not file_column:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
270 # Apply same filters to all files, just create a unviersal filter files
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
271 # and run msconvert once.
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
272 filters_file_path = _create_filters_file(options, debug=options.debug)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
273 cmd = "%s -c %s" % (cmd, filters_file_path)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
274 else:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
275 # Dispatching on a column to filter different files differently, need to filter
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
276 # each input once with msconvert and then merge once.
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
277 filtered_files = []
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
278 for index, input_file in enumerate(input_files):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
279 filters_file_path = _create_filters_file(options, index + 1, debug=options.debug)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
280 filter_cmd = "%s -c %s" % (cmd, filters_file_path)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
281 filtered_output_file = _run(filter_cmd, output_dir='output%d' % index, inputs=[input_file], debug=options.debug)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
282 filtered_files.append(filtered_output_file)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
283 input_files = filtered_files
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
284 if len(input_files) > 1:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
285 cmd = "%s --merge" % cmd
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
286 output_file = _run(cmd, output_dir='output', inputs=input_files, debug=options.debug)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
287 shutil.copy(output_file, options.output)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
288 if options.refinement:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
289 # .mzRefinement.tsv
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
290 files = os.listdir(os.getcwd())
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
291 for fname in files:
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
292 if fname.endswith('.mzRefinement.tsv'):
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
293 shutil.copy(fname, options.refinement)
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
294 break
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
295
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
296 def __main__():
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
297 run_script()
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
298
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
299 if __name__ == '__main__':
6b6bba73eadb planemo upload commit d56659dd48f8c554a832787e71aca6ae65c90848
galaxyp
parents:
diff changeset
300 __main__()