annotate link_scan_datasets.py @ 3:7f02fc51bddf draft

planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
author jjohnson
date Wed, 25 Jul 2018 10:58:17 -0400
parents
children 274fdc50169b
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
1 #!/usr/bin/env python
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
2
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
3 from __future__ import print_function
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
4
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
5 import argparse
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
6 import difflib
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
7 from difflib import SequenceMatcher
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
8 import os
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
9 import os.path
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
10 import sys
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
11 import xml.sax
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
12
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
13
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
14
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
15 def __main__():
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
16 parser = argparse.ArgumentParser(
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
17 description='link spectrum datasets to the name used' +
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
18 ' in the identification dataset')
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
19 parser.add_argument(
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
20 'ident_files', nargs='+',
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
21 help='Pepxml or mzIdentML')
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
22 parser.add_argument(
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
23 '-n', '--scan_name', default=[], action='append',
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
24 help='Name for scan file')
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
25 parser.add_argument(
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
26 '-f', '--scan_file', default=[], action='append',
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
27 help='Path for scan file')
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
28 args = parser.parse_args()
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
29
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
30 class MzidHandler( xml.sax.ContentHandler):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
31 def __init__(self):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
32 xml.sax.ContentHandler.__init__(self)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
33 self.spectraDataFiles = []
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
34 self.spectraDataNames = []
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
35 self.searchDatabaseFiles = []
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
36 self.searchDatabaseNames = []
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
37 def startElement(self, tag, attrs):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
38 if tag == 'SpectraData':
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
39 id = attrs['id']
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
40 path = attrs['location']
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
41 filename = os.path.basename(path)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
42 name = attrs['name'] if 'name' in attrs else None
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
43 self.spectraDataFiles.append(filename)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
44 self.spectraDataNames.append(name if name else id)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
45 print ("SpectraData: %s %s" % (name if name else id, path))
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
46 if tag == 'SearchDatabase':
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
47 id = attrs['id']
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
48 path = attrs['location']
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
49 filename = os.path.basename(path)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
50 name = attrs['name'] if 'name' in attrs else None
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
51 self.searchDatabaseFiles.append(filename)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
52 self.searchDatabaseNames.append(name if name else id)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
53 print ("SearchDatabase: %s %s" % (name if name else id, path))
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
54 def endElement( self, name):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
55 pass
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
56 def characters( self, data):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
57 pass
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
58
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
59 class PepXmlHandler( xml.sax.ContentHandler):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
60 def __init__(self):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
61 xml.sax.ContentHandler.__init__(self)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
62 self.spectraDataFiles = []
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
63 self.spectraDataNames = []
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
64 def startElement(self, tag, attrs):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
65 if tag == 'msms_run_summary':
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
66 basename = attrs['base_name']
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
67 name = os.path.basename(basename)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
68 ext = attrs['raw_data']
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
69 path = '%s%s' % (basename,ext)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
70 filename = os.path.basename(path)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
71 self.spectraDataFiles.append(filename)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
72 self.spectraDataNames.append(name)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
73 print ("SpectraData: %s %s" % (name, path))
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
74 def endElement( self, name):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
75 pass
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
76 def characters( self, data):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
77 pass
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
78
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
79 parser = xml.sax.make_parser()
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
80 parser.setFeature(xml.sax.handler.feature_namespaces, 0)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
81 handler = PepXmlHandler()
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
82 parser.setContentHandler( handler )
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
83 for ident in args.ident_files:
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
84 parser.parse(ident)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
85
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
86 spectra_names = handler.spectraDataFiles
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
87
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
88 def best_match(name):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
89 if name in spectra_names:
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
90 return name
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
91 try:
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
92 r = [SequenceMatcher(None, name, spectra_names[x]).ratio() for x in range(len(spectra_names))]
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
93 return spectra_names[r.index(max(r))]
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
94 except Exception, e:
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
95 print ("best_match: %s %s" % (name, e))
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
96
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
97 for i,name in enumerate(args.scan_name):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
98 path = args.scan_file[i] if len(args.scan_file) > i else ''
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
99 (root, ext) = os.path.splitext(name)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
100 print ("SpectraFile: %s %s" % (name, path))
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
101 iname = best_match(name)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
102 print ("IdentName: %s %s" % (name, iname))
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
103 if not os.path.exists(iname) and os.path.exists(path):
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
104 os.symlink(path, iname)
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
105 print ("%s -> %s" % (iname, path))
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
106
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
107
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
108 if __name__ == "__main__":
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
109 __main__()
7f02fc51bddf planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/spectrast commit 379705f578f9a0465f497894c7d2b5f68b6a55e6-dirty
jjohnson
parents:
diff changeset
110