[Bio] / FortyEight / rp_CELLO_attribute_generation.pl Repository:
ViewVC logotype

Annotation of /FortyEight/rp_CELLO_attribute_generation.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (view) (download) (as text)

1 : mkubal 1.1
2 :     #
3 :     # Compute cell localization of proteins using CELLO server and create attribute data.
4 :     #
5 :    
6 :     use Data::Dumper;
7 :     use Carp;
8 :     use strict;
9 :     use FIG;
10 :     use FIG_Config;
11 :     use File::Basename;
12 :     use GenomeMeta;
13 :     use Sim;
14 :    
15 :     @ARGV == 1 or die "Usage: $0 job-dir\n";
16 :    
17 :     my $jobdir = shift;
18 :    
19 :     -d $jobdir or die "$0: job dir $jobdir does not exist\n";
20 :    
21 :     my $hostname = `hostname`;
22 :     chomp $hostname;
23 :    
24 :     my $genome = &FIG::file_head("$jobdir/GENOME_ID");
25 :     chomp $genome;
26 :     $genome =~ /^\d+\.\d+/ or die "$0: Cannnot find genome ID for jobdir $jobdir\n";
27 :    
28 :     my $meta = new GenomeMeta($genome, "$jobdir/meta.xml");
29 :    
30 :     my $genome_dir = "$jobdir/rp/$genome";
31 :    
32 :     my @apps = qw(rp_CELLO_search
33 :     rp_make_CELLO_attributes
34 :     );
35 :    
36 :     $meta->set_metadata("cello.hostname", $hostname);
37 :    
38 :     for my $app (@apps)
39 :     {
40 :     $meta->add_log_entry($0, "start $app on $hostname in $jobdir");
41 :    
42 :     my $cmd = "$FIG_Config::bin/$app $jobdir $genome > $jobdir/rp.errors/$app.stderr 2>&1";
43 :     warn "Compute: $cmd\n";
44 :     my $rc = system($cmd);
45 :     if ($rc != 0)
46 :     {
47 :     &fatal("$app computation failed with rc=$rc");
48 :     }
49 :     }
50 :    
51 :     $meta->add_log_entry($0, "finish cello computation on $jobdir");
52 :     $meta->set_metadata("status.cello", "complete");
53 :     exit(0);
54 :    
55 :     sub fatal
56 :     {
57 :     my($msg) = @_;
58 :    
59 :     if ($meta)
60 :     {
61 :     $meta->add_log_entry($0, ['fatal error', $msg]);
62 :     $meta->set_metadata("status.cello", "error");
63 :     }
64 :    
65 :     croak "$0: $msg";
66 :     }
67 :    

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3