<tool id="compare_humann2_output" name="Compare outputs of HUMAnN2 for several samples" version="0.1.0">
    <description>and extract information</description>

    <requirements>
    </requirements>

    <stdio>
        <exit_code range="1:" />
        <exit_code range=":-1" />
    </stdio>

    <version_command></version_command>

    <command><![CDATA[
        mkdir specifics 
        &&

        python $__tool_directory__/compare_humann2_output.py
            #for $sample in $samples:
                --sample_name "${sample.sample_name}"
                --charact_input_file "${sample.input}"
                --specific_output_file "specifics/specific_to_${sample.sample_name}.txt"
            #end for

            --most_abundant_characteristics_to_extract $charact_nb
            --more_abundant_output_file $more_abundant_output_file
            --similar_output_file $similar_output_file
            > $log
                
    ]]></command>

    <inputs>
        <repeat name="samples" title="Add sample and input file (HUMAnN2 output after normalization)" >
            <param name="sample_name" type="text" label="Name of the sample" help="(--sample_name)"/>
            <param name="input" format="txt,tabular" type="data" label="Input file corresponding to HUMAnN2 output" help="The HUMAnN2 output file contains relative abundance of gene families or pathways (after normalization, --charact_input_file)"/>
        </repeat>

        <param name="charact_nb" type="integer" value="10" label="Number of most abundant characteristics to extract for each sample" help="(--most_abundant_characteristics_to_extract)"/>
    </inputs>

    <outputs>
        <data name="more_abundant_output_file" format="tabular"
            label="${tool.name} on ${on_string}: More abundant characteristics for each sample" />
        <data name="similar_output_file" format="tabular"
            label="${tool.name} on ${on_string}: Similar characteristics and the relative abundances for all samples" />
        <data name="log" format="txt"
            label="${tool.name} on ${on_string}: Log" />
        <collection name="specific_files" type="list">
            <discover_datasets pattern="__designation_and_ext__" directory="specifics"/>
        </collection>

    </outputs>

    <tests>
    </tests>

    <help><![CDATA[
        **What it does**

        This tool compare HUMANnN2 outputs with gene families or pathways relative abundances of several samples and extract several files:

        - Similar gene families or pathways between the samples and the relative abundances and these similar characteristics
        - Most abundant gene families or pathways for each sample and the corresponding relative abundance in all samples

    ]]></help>

    <citations>
    </citations>
</tool>