[Bio] / MGRASTBackend / MGRASTPipeline.pm Repository:
ViewVC logotype

Annotation of /MGRASTBackend/MGRASTPipeline.pm

Parent Directory Parent Directory | Revision Log Revision Log


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

1 : jared 1.1 package MGRASTBackend::MGRASTPipeline;
2 :    
3 :     use ClusterStage;
4 :     use Filesys::DfPortable;
5 :     use Mail::Mailer;
6 :     use Job48;
7 :     use FIG;
8 :    
9 :     use strict;
10 :     use warnings;
11 :    
12 :     use base qw( Exporter );
13 :     our @EXPORT = qw ( check_free_space get_jobs get_pipeline_for_job process_upload mark_job_done);
14 :    
15 :     1;
16 :    
17 :    
18 :     ########################################
19 :     #
20 :     # Basic methods
21 :     #
22 :     #########################################
23 :    
24 :     sub check_free_space {
25 :     my ($job_spool_dir) = @_;
26 :     my $df = dfportable($job_spool_dir, 1024*1024*1024);
27 :     if (!defined($df)){
28 :     die "dfportable call failed on $job_spool_dir: $!";
29 :     }
30 :     if ($df->{bavail} < 10){
31 :     die sprintf "Not enough free space available (%.1f GB) in $job_spool_dir", $df->{bavail};
32 :     }
33 :     return;
34 :     }
35 :    
36 :     sub get_jobs {
37 :     my ($job_spool_dir) = @_;
38 :     opendir(D, $job_spool_dir) or die "Cannot open job directory $job_spool_dir: $!\n";
39 :     my @jobs = sort { $a <=> $b } grep { /^\d+$/ and -d "$job_spool_dir/$_" } readdir(D);
40 :     return \@jobs;
41 :     }
42 :    
43 :     sub get_pipeline_for_job {
44 :     my ($job_dir) = @_;
45 : jared 1.4 if(! -f "$job_dir/PIPELINE"){
46 : jared 1.1 return &{&pipeline("default")};
47 :     } else {
48 : jared 1.4 my $pipe = &FIG::file_head("$job_dir/PIPELINE", 1);
49 : jared 1.1 chomp $pipe;
50 :     if (ref(&pipeline($pipe)) eq 'CODE'){
51 :     return &{&pipeline($pipe)};
52 :     } else {
53 :     return undef;
54 :     }
55 :     }
56 :     }
57 :    
58 :     ########################################
59 :     #
60 :     # Pipelines
61 :     #
62 :     #########################################
63 :    
64 :     sub pipeline {
65 :     my ($name) = @_;
66 :     my %pipelines = (
67 :     default => \&default_process,
68 :     );
69 :    
70 :     if(defined $pipelines{$name}){
71 :     return $pipelines{$name};
72 :     } else {
73 : jared 1.7 return $pipelines{default};
74 : jared 1.1 }
75 :     }
76 :    
77 :     sub default_process {
78 :     return [[uploaded => \&process_upload],
79 :     [preprocess => ClusterStage->new('mg_preprocess',
80 :     sge_flag => "-l mg_preprocess",
81 :     )],
82 :     [sims => ClusterStage->new('mg_sims',
83 :     start_locally => 1,
84 :     )],
85 :     [check_sims => ClusterStage->new('mg_check_sims',
86 :     start_locally => 1,
87 :     )],
88 :     [create_seed_org => ClusterStage->new('mg_create_seed_org',
89 :     sge_flag => "-l mg_postproc_taxa_sims",
90 :     )],
91 : jared 1.6 [export_to_genbank => ClusterStage->new('mg_sims_to_gff_and_gbk',
92 : jared 1.3 start_locally => 1,
93 : jared 1.1 )],
94 :     ];
95 :     }
96 :    
97 :     ########################################
98 :     #
99 :     # Pipeline Stages
100 :     #
101 :     #########################################
102 :    
103 :     sub process_upload
104 :     {
105 :     return;
106 :     }
107 :    
108 :     sub mark_job_done
109 :     {
110 :    
111 :     print "Attempting to Email user\n";
112 : jared 1.2
113 : jared 1.1 my($job_id, $job_dir, $meta, $req) = @_;
114 :    
115 :     if (open(D, ">$job_dir/DONE"))
116 :     {
117 :     print D time . "\n";
118 :     close(D);
119 :     }
120 :     else
121 :     {
122 :     warn "Error opening $job_dir/DONE: $!\n";
123 :     }
124 :    
125 :     my $job = new Job48($job_id);
126 :    
127 :     my $userobj = $job->getUserObject();
128 :    
129 :     print "setting meta $meta\n";
130 :     $meta->set_metadata("status.final","complete");
131 :     print "setting meta $meta .. done\n";
132 :    
133 :     if ($userobj)
134 :     {
135 :     my $email = $userobj->email();
136 :     my $name = join(" " , $userobj->firstname(), $userobj->lastname());
137 :    
138 :     my $full = $name ? "$name <$email>" : $email;
139 :     print "send email to $full\n";
140 :    
141 :     my $mail = Mail::Mailer->new();
142 :     $mail->open({
143 :     To => $full,
144 :     From => 'Metagenome RAST server <mg-rast@mcs.anl.gov>',
145 :     Subject => "MG-RAST job completed"
146 :     });
147 :    
148 :     my $gname = $job->genome_name;
149 :     my $entry = $FIG_Config::fortyeight_home;
150 :     $entry = "http://metagenomics.nmpdr.org/" if $entry eq '';
151 :     print $mail "The annotation job that you submitted for $gname has completed.\n";
152 :     print $mail "It is available for browsing at $entry as job number $job_id.\n";
153 :     $mail->close();
154 :     }
155 :     }

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3