From b6fd50ed91cb2001e2b564c5b2da1da576c3e54b Mon Sep 17 00:00:00 2001
From: martinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Date: Fri, 5 Dec 2008 02:52:08 +0000
Subject: [PATCH] added bed2fixedstep.c as backend to calc_fixedstep

git-svn-id: http://biopieces.googlecode.com/svn/trunk@330 74ccb610-7750-0410-82ae-013aeee3265d
---
 code_perl/Maasha/Biopieces.pm | 14 +++++---------
 1 file changed, 5 insertions(+), 9 deletions(-)

diff --git a/code_perl/Maasha/Biopieces.pm b/code_perl/Maasha/Biopieces.pm
index 5f6dd07..59b7df7 100644
--- a/code_perl/Maasha/Biopieces.pm
+++ b/code_perl/Maasha/Biopieces.pm
@@ -457,13 +457,6 @@ sub get_options
             percent|p
         );
     }
-    elsif ( $script eq "calc_fixedstep" )
-    {
-        @options = qw(
-            score|S
-            log10|L
-        );
-    }
     elsif ( $script eq "transliterate_seq" )
     {
         @options = qw(
@@ -2680,7 +2673,10 @@ sub script_calc_fixedstep
 
     foreach $chr ( sort keys %{ $file_hash } )
     {
-        $fixedstep_file = Maasha::UCSC::Wiggle::fixedstep_calc( $file_hash->{ $chr }, $chr, $options->{ 'score' }, $options->{ 'log10' } );
+        $bed_file       = $file_hash->{ $chr };
+        $fixedstep_file = "$bed_file.fixedstep";
+        
+        Maasha::Common::run( "bed2fixedstep", "< $bed_file > $fixedstep_file" );
 
         $fh_in = Maasha::Filesys::file_read_open( $fixedstep_file );
 
@@ -2693,8 +2689,8 @@ sub script_calc_fixedstep
 
         close $fh_in;
 
-        unlink $file_hash->{ $chr };
         unlink $bed_file;
+        unlink $fixedstep_file;
     }
 }
 
-- 
2.39.5