rdown pe

code_new/Sumit/rdown_pe.py

@@ -65,7 +65,7 @@ class Ringdown_PE:
 
         sigma2 = self.errors2 + self.l_int*(self.datare** 2+self.dataim**2) * np.exp(2 * theta[-1])
 
-        result = -0.5*np.sum(((self.datare - modelevre)**2+(self.dataim - modelevim)**2)/sigma2+self.l_int*(2*np.log(sigma2)))-self.l_int*self.norm_factor
+        result = -0.5*np.sum(((self.datare - modelevre)**2+(self.dataim - modelevim)**2)/sigma2+self.l_int*(np.log(2*np.pi*sigma2)))-self.l_int*self.norm_factor
 
         if np.isnan(result):
             return -np.inf