code: def chi_squared(self, p=None):
          
          chi2s = self.chi_squareds(p)
          if chi2s == None: return None
          
          return sum(self.chi_squareds(p))
docs: Returns the total chi squared (summed over all massaged data sets).
          
          p=None means use the fit results.