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