fix(morph-plot): handle files with single point cell types

If a file contains the definition of a point cell type, for example,
`IafTauRefCell`, we should still be able to plot the file, even if its
only a single circle/sphere. Until now, we only handled files if they
had multi-compartmental `Cell` or `Morphology` instances, and the
plotter would crash for files that had other point cell types.

Author: sanjayankur31

pyneuroml/plot/PlotMorphology.py

@@ -483,7 +483,7 @@ def plot_2D(
             radius = pop_id_vs_radii[pop_id] if pop_id in pop_id_vs_radii else 10
             color = pop_id_vs_color[pop_id] if pop_id in pop_id_vs_color else None
 
-            if cell is None:
+            if cell is None or not isinstance(cell, Cell):
                 plot_2D_point_cells(
                     offset=pos,
                     plane2d=plane2d,

pyneuroml/plot/PlotMorphologyVispy.py

@@ -607,6 +607,9 @@ def plot_interactive_3D(
         # other networks
         else:
             plottable_nml_model = nml_model
+
+        logger.debug(plottable_nml_model.info(show_contents=True))
+
     # what did we get?
     else:
         raise ValueError(f"Could not process argument: {nml_model}")
@@ -671,7 +674,7 @@ def plot_interactive_3D(
     else:
         cell = list(pop_id_vs_cell.values())[0]
 
-        if cell is not None:
+        if cell is not None and isinstance(cell, Cell):
             view_min, view_max = get_cell_bound_box(cell)
         else:
             logger.debug("Got a point cell")
@@ -790,7 +793,8 @@ def plot_interactive_3D(
             except AttributeError:
                 logging.debug(f"Plotting a point cell at {pos}")
 
-            if cell is None:
+            # a point cell component type
+            if cell is None or not isinstance(cell, Cell):
                 meshdata.append(
                     (
                         f"{radius:.1f}",

pyneuroml/utils/__init__.py

@@ -97,8 +97,12 @@ def extract_position_info(
     cell_elements = []
     popElements = []
 
-    cell_elements.extend(nml_model.cells)
-    cell_elements.extend(nml_model.cell2_ca_poolses)
+    members = nml_model.info(show_contents=True, return_format="dict")
+
+    for member, mdict in members.items():
+        if "cell" in mdict["type"].lower():
+            cell_elements.extend(mdict["members"])
+    # cell_elements.extend(nml_model.cell2_ca_poolses)
 
     # if there are no cells, look at morphologies
     if len(cell_elements) == 0:
@@ -715,7 +719,6 @@ def get_model_file_list(
             lems_def_dir = get_model_file_list(inc, filelist, rootdir, lems_def_dir)
 
     elif rootfile.endswith(".sedml"):
-
         try:
             import libsedml
         except ModuleNotFoundError:

pyneuroml/utils/plot.py

@@ -351,15 +351,14 @@ def load_minimal_morphplottable__model(
     model_members = list(vars(nml_model).keys())
     required_members = [
         "id",
-        "cells",
         "morphology",
         "cell2_ca_poolses",
         "networks",
         "populations",
         "includes",
     ]
     for m in model_members:
-        if m not in required_members:
+        if m not in required_members and "cells" not in m:
             setattr(nml_model, m, None)
             logger.debug(f"Dropped {m}")