docs: def com(self):
          """Center of mass, calculated from the mean of the PMT positions"""
          if self._com is None:
              self._com = np.mean(self.pmts.pos, axis=0)
          return self._com