#!/usr/bin/perl -w 


@dirs = qw(hs1A_50
hs1B_51
hs1C_52
hs1D_53);

@info = qw(wt+DMSO 
wt+pB
y2h2+DMSO
y2h2+pB);


@files = qw(
YKR095W-A
YPR028W
YDL191W
);

#YKR095W-A PCC1
#YPR028W YOP1
#YDL191W RPL35A
open(TG,">tmp.gnpl");
print TG "set terminal png size 1200,900 \n"; 
print TG "set output \"all.png\" \n";

print TG "set multiplot layout 2,2 \n";
#print TG "set key outside top \n" ;


for($i=0;$i<4;$i++){

$d = $dirs[$i];
$info = $info[$i];
print TG "set title '$info' \n";
print TG "set xlabel \"'Position relative to 3\' SS\" \n";
print TG "set ylabel 'Fraction splice co-transcriptionally' \n";
print TG "plot [0:600][-0.05:1.05] './smit_files/$d/YKR095W-A_splicingValues.txt' u 1:4 w l smooth sbezier title 'YKR095W PCC1' , "; 
print TG " './smit_files/$d/YPR028W_splicingValues.txt' u 1:4 w l smooth sbezier title 'YPR028W YOP1' , "; 
print TG " './smit_files/$d/YDL191W_splicingValues.txt' u 1:4 w l smooth sbezier title 'YDL191W RPL35A' \n"; 


}