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