post:getlineageindices [FAST]

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
====== getLineageIndices.m ======

**Class:** Function

**Description:** Allows you to extract the indices and generations of all descendants of a target cell from the procTracks structure. 

**Author:** Oliver J. Meacock

<code matlab getLineageIndices.m>
function [linInds,linLens] = getLineageIndices(procTracks,motherID,callDepth)
%GETLINEAGEINDICES is a recursive function that allows you to extract all of
%the indicies from a given tree-structured lineage. Now also returns the 
%distance of each lineage from the founder cell.
%
%   INPUTS:
%       -procTracks: A track structure which has been output by the
%       Divisions module of the FAST GUI. Includes fields D1, D2 and M,
%       indicating the two daughters of this cell and the mother of this
%       cell.
%       -motherID: The ID of the track that is the founder of the current
%       lineage.
%       -callDepth: The current depth into the recursion you are. Starts at
%       0 upon first call.
%
%   OUTPUTS:
%       -linInds: A vector of all the descendants that arise from the initial mother ID
%       -linLens: The distance from the foundational track each track is.
%
%   Author: Oliver J. Meacock, (c) 2019

D1ID = procTracks(motherID).D1;
D2ID = procTracks(motherID).D2;

if ~isempty(D1ID)
    [D1lins,D1lens] = getLineageIndices(procTracks,D1ID,callDepth + 1);
else
    D1lins = [];
    D1lens = [];
end
if ~isempty(D2ID)
    [D2lins,D2lens] = getLineageIndices(procTracks,D2ID,callDepth + 1);
else
    D2lins = [];
    D2lens = [];
end

linInds = [motherID;D1lins;D2lins];
linLens = [callDepth;D1lens;D2lens];
</code>