def trunk_vectors(nrn, neurite_type=NeuriteType.all):
neurite_filter = is_type(neurite_type)<EOL>nrns = neuron_population(nrn)<EOL>return np.array([morphmath.vector(s.root_node.points[<NUM_LIT:0>], n.soma.center)<EOL>for n in nrns<EOL>for s in n.neurites if neurite_filter(s)])<EOL>
Calculates the vectors between all the trunks of the neuron
    and the soma center.