def norm(self): x = centered1d(self.x) y = centered1d(self.y) return scipy.sqrt(trapz2(self.intensity(), x=x, y=y))