<tool id="deeptools_plot_enrichment" name="plotEnrichment" version="@WRAPPER_VERSION@.0" profile="18.01">
    <description>plots read/fragment coverage over sets of regions</description>
    <macros>
        <token name="@BINARY@">plotEnrichment</token>
        <import>deepTools_macros.xml</import>
    </macros>
    <expand macro="requirements" />
    <command>
<![CDATA[
        @multiple_input_bams@
        @BINARY@
            @THREADS@

            --plotFile '$outFileName'

            --bamfiles #echo " ".join($files)#
            --labels #echo " ".join($labels)#

            @multiple_bed@

            #if $outRawCounts:
                --outRawCounts '$outFileRawCounts'
            #end if

            --plotFileFormat '$outFileFormat'

            #if str($region).strip() != "":
                --region '$region'
            #end if

            #if $advancedOpt.showAdvancedOpt == "yes"
                #if $advancedOpt.attributeKey:
                    --attributeKey '$advancedOpt.attributeKey'
                #end if

                #if $advancedOpt.labels and str($advancedOpt.labels).strip() != "":
                    --labels $advancedOpt.labels
                #end if

                #if $advancedOpt.regionLabels and str($advancedOpt.regionLabels).strip() != "":
                    --regionLabels $advancedOpt.regionLabels
                #end if

                #if $advancedOpt.plotTitle and str($advancedOpt.plotTitle.value) != "":
                    --plotTitle '$advancedOpt.plotTitle'
                #end if

                $advancedOpt.perSample

                $advancedOpt.variableScales
                
                $advancedOpt.perSample

                $advancedOpt.variableScales

                --plotWidth '$advancedOpt.plotWidth'

                --plotHeight '$advancedOpt.plotHeight'

                #if str($advancedOpt.colors).strip() != "":
                    --colors #echo ' '.join( ["'%s'" % $color for $color in $advancedOpt.colors.split()] )#
                #end if

                --numPlotsPerRow '$advancedOpt.numPlotsPerRow'

                --alpha '$advancedOpt.alpha'

                @ADVANCED_OPTS_READ_PROCESSING@

                #if $advancedOpt.Offset:
                    --Offset $advancedOpt.Offset
                #end if

                @blacklist@

                $advancedOpt.keepExons

            #end if
]]>
    </command>

    <inputs>
        <expand macro="multiple_input_bams" MIN="1"/>
        <expand macro="custom_sample_labels" />
        <expand macro="multiple_bed" />

        <expand macro="region_limit_operation" />
        <expand macro="input_image_file_format" />

        <expand macro="advancedOpt_scaffold">
            <param argument="--attributeKey" type="text" optional="true" label="Optional attribute key"
                help="Instead of deriving the feature label from the feature column, use the value of the given
                attribute key. For example, the gene_biotype. Note that 'None' is used for BED files
                or entries where the attributeKey is not found." />

            <param argument="--plotHeight" type="integer" value="20" min="3"
                label="Plot height"
                help="Height in cm. The default for the plot height is 20 centimeters. The minimum value is 3 cm." />

            <param argument="--plotWidth" type="integer" value="20" min="1"
                label="Plot width"
                help="Width in cm. The default value is 20 centimeters. The minimum value is 1 cm." />

            <param argument="--labels" type="text" size="30"
                label="Labels for the samples (each BAM file) plotted"
                help="The default is to use the file name of the sample. The sample labels should be separated by
                spaces and quoted if a label itself contains a space, e.g., label-1 &quot;label 2&quot;">
                <sanitizer>
                    <valid initial="string.printable" />
                </sanitizer>
            </param>

            <param argument="--regionLabels" type="text" size="30"
                label="Labels for the features"
                help="For BED files, this sets the label for all entries in each file.
                      For GTF files this is ignored." />

            <expand macro="plotTitle" />

            <param argument="--alpha" type="float" value="0.9" min="0.0" max="1.0"
                label="Alpha channel"
                help="The alpha channel value (between 0 and 1). A value of 0 is transparent. Default: 0.9" />

            <param argument="--colors" type="text" value="" size="40"
                label="List of colors to use for the plotted bars"
                help="Color names and html hex strings (e.g. #eeff22) are accepted.
                The color names should be separated by spaces. (--colors red blue green)">
                <validator type="expression"
                    message="Only numbers, digits, '#' and spaces are allowed.">all(c in ' #abcdefghijklmnopqrstuvwxyz0123456789' for c in value)</validator>
            </param>

            <param argument="--perSample" type="boolean" truevalue="--perSample" falsevalue=""
                label="Create one plot per-sample"
                help="Plot per-sample, rather than per-feature." />

            <param argument="--variableScales" type="boolean" truevalue="--variableScales" falsevalue=""
                label="Allow variable scales"
                help="By default, the y-axis always goes from 0 to 100. If this option is selected, then the maximum value will vary with the dataset." />

            <param argument="--numPlotsPerRow" type="integer" value="4" min="1"
                label="Number of plots/row"
                help="" />

            <expand macro="read_processing_options" />

            <param argument="--keepExons" type="boolean" truevalue="--keepExons" falsevalue=""
                label="Include BED12 exons"
                help="Only for BED12 files, include columns 10-12 rather than just 2 and 3." />

            <param argument="--Offset" type="text"  value="" optional="True"
                label="Offset inside each alignment to use for the signal location."
                help="Uses this offset inside of each read as the signal. This is useful in
                    cases like RiboSeq or GROseq, where only the 12th, 15th or 1st base aligned
                    should be used to denote where the signal is (rather than the span of the
                    whole alignment). This can be paired with the --filterRNAstrand option. Note
                    that negative values indicate offsets from the end of each read. A value of
                    1 indicates the first base of the alignment (taking alignment orientation
                    into account). Likewise, a value of -1 is the last base of the alignment. An
                    offset of 0 is not permitted. If two values (separated by spaces) are specified, then they will be
                    used to specify a range of positions. Note that specifying something like
                    --Offset 5 -1 will result in the 5th through last position being used, which
                    is equivalent to trimming 4 bases from the 5-prime end of alignments." />

            <expand macro="blacklist" />

        </expand>

        <param argument="--outRawCounts" type="boolean" label="Save percentages to a file" help=""/>

    </inputs>
    <outputs>
        <expand macro="output_image_file_format_not_nested" />
        <data format="tabular" name="outFileRawCounts" label="${tool.name} on ${on_string}: percentages">
            <filter>outRawCounts is True</filter>
        </data>
    </outputs>
    <tests>
        <test>
            <param name="bamfiles" value="bowtie2 test1.bam,bowtie2 test1.bam" ftype="bam" />
            <param name="BED" value="multiBamSummary_regions.bed,multiBamSummary_regions.bed" ftype="bed" />
            <param name="outRawCounts" value="true" />
            <param name="showAdvancedOpt" value="yes" />
            <param name="minMappingQuality" value="0" />
            <param name="regionLabels" value="up down" />
            <output name="outFileName" file="plotEnrichment_output.png" ftype="png" compare="sim_size" delta="1000" />
            <output name="outFileRawCounts" file="plotEnrichment_output.txt" ftype="tabular" />
        </test>
    </tests>

    <help>
<![CDATA[

What it does
-------------

This tool determines read/fragment coverage of regions. The regions (e.g., exons, genes or peaks) can be specified in one or more BED or GTF files. For GTF files, the feature type is taken from column 3. For BED files, the file name is used. For BED files, the feature labels can be changed.

Output
--------

The output file is a plot in the format specified. A table of the percentages and raw counts can also be created.

.. image:: $PATH_TO_IMAGES/plotEnrichment_output.png
   :width: 600
   :height: 600

-----

@REFERENCES@
]]>
    </help>
    <expand macro="citations" />
</tool>