Mercurial > repos > modencode-dcc > spp_package

changeset 10:c57381059179 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author modencode-dcc
date Thu, 07 Feb 2013 23:06:21 -0500
parents dfdaa943f855
children 495a6d033ca1
files spp_wrapper.py
diffstat 1 files changed, 9 insertions(+), 7 deletions(-) [+]
line wrap: on
line diff
--- a/spp_wrapper.py	Wed Feb 06 21:23:05 2013 -0500
+++ b/spp_wrapper.py	Thu Feb 07 23:06:21 2013 -0500
@@ -47,7 +47,8 @@
 
     #run cmdline
     #========================================================================================
-    tmp_dir = tempfile.mkdtemp()
+    #tmp_dir = tempfile.mkdtemp()
+    tmp_dir = os.path.dirname(options['chip_file'])
     stderr_name = tempfile.NamedTemporaryFile().name
     proc = subprocess.Popen( args=cmdline, shell=True, cwd=tmp_dir, stderr=open( stderr_name, 'wb' ) )
     proc.wait()
@@ -66,17 +67,18 @@
 
     #determine if the outputs are there, copy them to the appropriate dir and filename
     #========================================================================================
+    chip_name = os.path.basename(options['chip_file'])
+    input_name = os.path.basename(options['input_file'])
+
     created_default_file =  os.path.join( tmp_dir, "default_output.txt" )
     if os.path.exists( created_default_file ):
         shutil.move( created_default_file, output_default_file )
 
-    chip_name = os.path.basename(options['chip_file'])
-    input_name = os.path.basename(options['input_file'])
-    created_narrow_peak =  os.path.join( tmp_dir, "%s_VS_%s.narrowPeak.gz" % (chip_name, input_name) )
+    created_narrow_peak =  os.path.join( tmp_dir, "%s_VS_%s.narrowPeak" % (chip_name, input_name) )
     if os.path.exists( created_narrow_peak ):
         shutil.move( created_narrow_peak, output_narrow_peak )
  
-    created_region_peak =  os.path.join( tmp_dir, "%s_VS_%s.regionPeak.gz" % (chip_name, input_name) )
+    created_region_peak =  os.path.join( tmp_dir, "%s_VS_%s.regionPeak" % (chip_name, input_name) )
     if os.path.exists( created_region_peak ):
         shutil.move( created_region_peak, output_region_peak )
 
@@ -84,11 +86,11 @@
     if os.path.exists( created_peakshift_file ):
         shutil.move( created_peakshift_file, output_peakshift_file )
 
-    created_rdata_file =  os.path.join( tmp_dir, "%s.Rdata" % options['chip_file'] )
+    created_rdata_file =  os.path.join( tmp_dir, "%s.Rdata" % chip_name )
     if os.path.exists( created_rdata_file ):
         shutil.move( created_rdata_file, output_rdata_file )
 
-    created_plot_file =  os.path.join( tmp_dir, "%s.pdf" % options['chip_file'] )
+    created_plot_file =  os.path.join( tmp_dir, "%s.pdf" % chip_name )
     if os.path.exists( created_plot_file ):
         shutil.move( created_plot_file, output_plot_file )