#!/usr/bin/env perl

$conf = $ARGV[0];
$n1 = $ARGV[1];
$n2 = $ARGV[2];

$esum0 = 0;
$esum1 = 0;
for($i=$n1;$i<=$n2;$i++){

    $i = sprintf("%0.3d",$i);


    system("analyze $conf.$i E | grep \"ACE\"  -A 1 > log ");


    open(IN,"log");
    @temp = split(" ",<IN>);
    $e0 = $temp[2];
    @temp = split(" ",<IN>);
    $e1 = $temp[1];
    close IN;



   print "$i  $e0 $e1 
";


#---average-calculation---

   $esum0 = $esum0 + $e0;
   $esum1 = $esum1 + $e1;

}

{
   $esum0=$esum0/$n2;
   $esum1=$esum1/$n2;
    $esum0= sprintf("%0.4f",$esum0);
   $esum1= sprintf("%0.4f",$esum1);

}

#averages
my@coolbeans = ('averages',  " $esum0", "$esum1");

print "
@coolbeans\n";




#---standard-deviation---


for($i=$n1;$i<=$n2;$i++){

        $estd0=($esum0-$e0) ** 2;
        $estd0=($estd0/$n2) **0.5;
        $estd0= sprintf ("%0.4f", $estd0);

        $estd1=($esum1-$e1) ** 2;
        $estd1=($estd1/$n2) **0.5;
        $estd1= sprintf ("%0.4f", $estd1);
}

#stdev
my@arroba = ('stdev', "$estd0", "$estd1");
print "
@arroba\n";