Merge pull request #67 from biowdl/controlTreat

fix CaseControl and macs2 tasks

biopet/sampleconfig.wdl

@@ -85,12 +85,13 @@ task CaseControl {
         File? toolJar
         String? preCommand
         Array[File]+ inputFiles
+        Array[File]+ inputIndexFiles
         Array[File]+ sampleConfigs
         String outputPath
         String controlTag = "control"
 
         Int memory = 4
-        Float memoryMultiplier = 1.5
+        Float memoryMultiplier = 2.0
     }
 
     String toolCommand = if defined(toolJar)
@@ -101,9 +102,9 @@ task CaseControl {
         set -e -o pipefail
         ~{preCommand}
         mkdir -p $(dirname ~{outputPath})
-        ~{toolCommand} CromwellArrays \
-        -i ~{sep="-i " inputFiles} \
-        -s ~{sep="-s " sampleConfigs} \
+        ~{toolCommand} CaseControl \
+        -i ~{sep=" -i " inputFiles} \
+        -s ~{sep=" -s " sampleConfigs} \
         ~{"-o " + outputPath} \
         ~{"--controlTag " + controlTag}
     }

macs2.wdl

@@ -7,6 +7,8 @@ task PeakCalling {
         String? preCommand
         Array[File]+ inputBams
         Array[File]+ inputBamsIndex
+        Array[File]+? controlBams
+        Array[File]+? controlBamsIndex
         String outDir
         String sampleName
         Int threads = 1
@@ -19,6 +21,7 @@ task PeakCalling {
         ~{preCommand}
         macs2 callpeak \
         --treatment ~{sep = ' ' inputBams} \
+        ~{true="--control" false="" defined(controlBams)} ~{sep = ' ' controlBams} \
         --outdir ~{outDir} \
         --name ~{sampleName} \
         ~{true='--nomodel' false='' nomodel}