annotate protein_prophet_wrapper.rb @ 0:adf8ac7aa6c8

Uploaded
author iracooke
date Thu, 06 Jun 2013 01:33:00 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
1 require 'pathname'
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
2
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
3 $VERBOSE=nil
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
4
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
5 actual_output_path_string=ARGV.shift
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
6
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
7 # Second argument is the original input file name ... we'll change this below
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
8 original_input_file=ARGV[0]
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
9
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
10 # Before doing anything we append create a link to the input file in our working dir with ".pep.xml" appended to the input
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
11 # name because peptide prophet can't handle anything else
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
12
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
13 wd= Dir.pwd
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
14
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
15 original_input_path=Pathname.new("#{original_input_file}")
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
16 actual_input_path_string="#{wd}/#{original_input_path.basename}.pep.xml"
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
17
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
18 cmd = "ln -s #{original_input_file} #{actual_input_path_string};"
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
19
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
20 cmd << "rvm 1.9.3@protk-1.2.0 do protein_prophet.rb"
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
21
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
22
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
23 ARGV[0]="#{actual_input_path_string}"
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
24
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
25 ARGV.each { |a|
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
26
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
27 cmd << " #{a}"
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
28 }
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
29
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
30 cmd << " -o protein_prophet_results.prot.xml"
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
31
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
32 cmd << ";ruby -pi -e \"gsub('#{actual_input_path_string}', '#{original_input_file}.pep.xml')\" protein_prophet_results.prot.xml"
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
33
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
34 %x[#{cmd}]
adf8ac7aa6c8 Uploaded
iracooke
parents:
diff changeset
35