<tool id="mdanalysis_angle" name="Angle Analysis" version="@TOOL_VERSION@+galaxy@GALAXY_VERSION@">
    <description>- time series of Angles</description>
    <macros>
        <import>macros.xml</import>
        <token name="@GALAXY_VERSION@">0</token>
    </macros>   
    <expand macro="requirements" />
    <command detect_errors="exit_code">
<![CDATA[
    python '$__tool_directory__/angle.py' 
        --itraj '$trajin'
        --istr '$strin'
        --itrajext '$trajin.ext'
        --istrext '$strin.ext'
        --isegid1 '$segid1'
        --iresid1 '$resid1'
        --iname1 '$name1'
        --isegid2 '$segid2'
        --iresid2 '$resid2'
        --iname2 '$name2'
        --isegid3 '$segid3'
        --iresid3 '$resid3'
        --iname3 '$name3'
        --output '$output'
        --oangle_plot '$angle_plot'
    2>&1
]]></command>
    <inputs>
        <expand macro="analysis_inputs"/>
        <param name="segid1"  type="text" value="PRO" label="Segment ID of atom 1"/>
        <param name="resid1"  type="text" value="212" label="Residue ID of atom 1"/>
        <param name="name1"  type="text" value="OE2" label="Atom name of atom 1"/>
        <param name="segid2"  type="text" value="HET" label="Segment ID of atom 2"/>
        <param name="resid2"  type="text" value="3" label="Residue ID of atom 2"/>
        <param name="name2"  type="text" value="C1" label="Atom name of atom 2"/>
        <param name="segid3"  type="text" value="HET" label="Segment ID of atom 3"/>
        <param name="resid3"  type="text" value="3" label="Residue ID of atom 3"/>
        <param name="name3"  type="text" value="C2" label="Atom name of atom 3"/>
    </inputs>
    <outputs>
        <data format="tabular" name="output" label="Angle Analysis raw data"/>
        <data format="png" name="angle_plot" label="Angle Analysis Plot"/>
    </outputs>
    <tests>
        <test>
            <expand macro="tests_inputs"/>
            <param name="segid1" value="PRO"/>
            <param name="resid1" value="212"/>
            <param name="name1" value="OE2"/>
            <param name="segid2" value="HET"/>
            <param name="resid2" value="3"/>
            <param name="name2" value="C1"/>
            <param name="segid3" value="HET"/>
            <param name="resid3" value="3"/>
            <param name="name3" value="C2"/>
            <output name="output" file="Angle_Analysis_raw_data.tabular" />
        </test>
        <test>
            <expand macro="tests_inputs_gmx"/>
            <param name="segid1" value="SYSTEM"/>
            <param name="resid1" value="212"/>
            <param name="name1" value="OE2"/>
            <param name="segid2" value="SYSTEM"/>
            <param name="resid2" value="3"/>
            <param name="name2" value="C1"/>
            <param name="segid3" value="SYSTEM"/>
            <param name="resid3" value="3"/>
            <param name="name3" value="C2"/>
            <output name="output">
                <assert_contents>
                    <has_n_columns n="2" />
                    <has_line_matching expression="0.0\t70.*" />
                    <has_line_matching expression="10.0\t71.*" />
                    <has_line_matching expression="11.0\t81.*" />
                </assert_contents>
            </output>
        </test>
    </tests>
    <help><![CDATA[
.. class:: infomark

**What it does**
        
This tool calculates and plots the angle between three atoms.

_____


.. class:: infomark

**Input**

       - Trajectory file  (DCD).
       - PDB file.
       - Segment IDs, residue IDs and names of the three atoms for calculating angles.

Note that a MDAnalysis 'segment' is a larger organizational unit,  for example one protein or all the solvent molecules or simply the whole system.

_____

        
.. class:: infomark

**Output**

       - Tab-separated file of raw data of the angle between three atoms calculated for each frame.
       - Image (as png) of the time series graph.

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