-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathdock_queue_1.pl
40 lines (29 loc) · 965 Bytes
/
dock_queue_1.pl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
#! /usr/bin/perl
$projectname = $ARGV[0];
$max_run = $ARGV[1];
$thorough = $ARGV[2];
$script_home = $ENV{'PWD'};
$results =$ENV{'PWD'}.'/'.$projectname.'/';
%ligand = (1 => "d_leucine", 2 => "l_leucine", 3 => "d_isoleucine", 4 => "l_isoleucine", 5 => "d_valine");
$NoLigand = scalar(keys %ligand);
for($j=1;$j<=$NoLigand;$j++)
{
chdir $script_home;
$command= "perl Multiple_ICM_Dockings.pl $projectname $ligand{$j} $max_run $thorough";
system($command);
}
open(W,'>',$results."Total_Results.log") || die "Please give me output filename $!";
for($f=1;$f<=$NoLigand;$f++)
{
open(Readlog,'<', "$results$ligand{$f}".'/'.$projectname.'_'.$ligand{$f}.'.log') || die $!;
while($line=<Readlog>)
{
chomp($line);
foreach($line) { s/^\s+//;s/\s+$//; s/\s+/ /g; }
my @temp=split(' ',$line);
$ICMscore = $temp[1];
close(Readlog)||die $!;
}
printf W "%-15s\t%15s\n",$ligand{$f},$ICMscore;
}
close(W)||die $!;