annotate trinityToolWrapper.py @ 0:a9723c18d389 draft

Adding align_and_estimate_abundance tool.
author trinity_ctat
date Thu, 21 Sep 2017 15:11:18 -0400
parents
children 9e5c3f162eca
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
1 #!/usr/bin/env python
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
2
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
3
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
4 # borrowed from: http://wiki.g2.bx.psu.edu/Future/Job%20Failure%20When%20stderr and modified for use with Trinity tools.
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
5
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
6 """
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
7 Wrapper that execute a program and its arguments but reports standard error
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
8 messages only if the program exit status was not 0
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
9 Example: ./stderr_wrapper.py myprog arg1 -f arg2
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
10 """
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
11
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
12 import sys, subprocess, os
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
13
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
14 assert sys.version_info[:2] >= ( 2, 4 )
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
15
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
16 TRINITY_BASE_DIR = ""
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
17 if os.environ.has_key('TRINITY_HOME'):
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
18 TRINITY_BASE_DIR = os.environ['TRINITY_HOME'];
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
19 else:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
20 # Cicada Dennis added looking for the location of the Trinity program using the unix "which" utility.
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
21 # I tried using "command -v Trinity" but for some reason, I was getting a OS permission error with that.
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
22 try:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
23 pipe1 = subprocess.Popen(["which", "Trinity"], stdout=subprocess.PIPE)
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
24 except:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
25 t, v, tb = sys.exc_info()
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
26 sys.stderr.write("You must set the environmental variable TRINITY_HOME to the base installation directory of Trinity before running this");
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
27 raise t, v, tb
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
28 else:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
29 output1, err1 = pipe1.communicate()
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
30 # Take off the last part of the path (which is the Trinity command)
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
31 TRINITY_BASE_DIR = "/".join(output1.split("/")[0:-1])
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
32
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
33 # get bindir
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
34 bindir = sys.argv[0]
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
35 bindir = bindir.split("/")
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
36 if len(bindir) > 1:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
37 bindir.pop()
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
38 bindir = "/".join(bindir)
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
39 else:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
40 bindir = "."
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
41
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
42
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
43 ## add locations of tools to path setting.
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
44 TOOL_PATHS_FILE = bindir + "/__add_to_PATH_setting.txt";
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
45 for line in open(TOOL_PATHS_FILE):
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
46 line = line.rstrip()
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
47 os.environ['PATH'] += ":" + line
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
48
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
49
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
50 def stop_err( msg ):
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
51 sys.stderr.write( "%s\n" % msg )
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
52 sys.exit()
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
53
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
54 def __main__():
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
55 # Get command-line arguments
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
56 args = sys.argv
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
57 # Remove name of calling program, i.e. ./stderr_wrapper.py
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
58 args.pop(0)
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
59 # If there are no arguments left, we're done
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
60 if len(args) == 0:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
61 return
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
62
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
63 # If one needs to silence stdout
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
64 #args.append( ">" )
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
65 #args.append( "/dev/null" )
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
66
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
67 args[0] = "".join([TRINITY_BASE_DIR, '/', args[0]]);
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
68
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
69 cmdline = " ".join(args)
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
70
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
71
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
72
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
73 try:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
74 # Run program
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
75 err_capture = open("stderr.txt", 'w')
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
76 proc = subprocess.Popen( args=cmdline, shell=True, stderr=err_capture, stdout=sys.stdout )
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
77 returncode = proc.wait()
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
78 err_capture.close()
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
79
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
80
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
81 if returncode != 0:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
82 raise Exception
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
83
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
84 except Exception:
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
85 # Running Grinder failed: write error message to stderr
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
86 err_text = open("stderr.txt").readlines()
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
87 stop_err( "ERROR:\n" + "\n".join(err_text))
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
88
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
89
a9723c18d389 Adding align_and_estimate_abundance tool.
trinity_ctat
parents:
diff changeset
90 if __name__ == "__main__": __main__()