#!/usr/bin/perl
# Description: Extract pre-calculated SIDD values of region within a prokaryotic genome sequence
# Author: Peter Fischer Hallin
# Email: pfh@cbs.dtu.dk
# Version: 1.0a ws0
# Date: 2009-05-04

use Getopt::Long;
use Data::Dumper;

my ($string,$bool);

my ($w,$e,$s,$g,$a,$f,$t);

&GetOptions (
'w:s'  , \$w,
'e:s'  , \$e,
's:s'  , \$s,
'g:s'  , \$g,
'a:s'  , \$a,
'f:s'  , \$f,
't:s'  , \$t,
);


require "xml-compile.pl";

my $proxy = WSDLclient ( 'http://www.cbs.dtu.dk/ws/SIDDbase/SIDDbase_1_0a_ws0.wsdl' );

# input parameters for retrieving the sequence information
my $response = $proxy->{getSIDD}->( 
 parameters => { 
  parameters => {
   accession => $a,
   energetics => $e , 
   sigma => $s , 
   from => $f,
   to => $t,
   weight => $w ,
   Gformat => '%0.2f' , 
   Pformat => '%0.6f' , 
   format => 'string' , 
  }
 }
);
my @POS;
die "unable to obtain response inless $response\n" unless defined $response;
printf STDERR "version: %s\n"  , $response->{parameters}->{output}->{version};
printf STDERR "organism_name: %s\n"  , $response->{parameters}->{output}->{organism_name};
printf STDERR "segment: %s\n"  , $response->{parameters}->{output}->{segment};
printf STDERR "entrez_pid: %s\n"  , $response->{parameters}->{output}->{entrez_pid};
printf STDERR "method: %s\n"  , $response->{parameters}->{output}->{method};
printf STDERR "digest: %s\n"  , $response->{parameters}->{output}->{digest};
printf STDERR "geometry: %s\n"  , $response->{parameters}->{output}->{geometry};
foreach ( split ",",$response->{parameters}->{output}->{string}->{G}){
print "$_\n";
}