Bug in new prejack observables.

[phys/u1casc.git] / u1casc-ordinary / obs_phi2.hpp

diff --git a/u1casc-ordinary/obs_phi2.hpp b/u1casc-ordinary/obs_phi2.hpp
index 3e4fe4878b87afb8df0526776ef9a184c2911f11..0cad20eb9e15d88c390f73bfebd9fee16edff738 100644 (file)
--- a/u1casc-ordinary/obs_phi2.hpp
+++ b/u1casc-ordinary/obs_phi2.hpp
@@ -97,10 +97,9 @@ void obs_phi2::prePhi2Sus(vector< vector <double> > *allVals, vector <double> *p
 }
 
 double obs_phi2::phi2Sus(vector <double> *preCalculated, vector<double> *excludedMeas, int nmeas, void *para) {
 }
 
 double obs_phi2::phi2Sus(vector <double> *preCalculated, vector<double> *excludedMeas, int nmeas, void *para) {

-  return (( (*preCalculated)[1] - pow((*excludedMeas)[0], 2) )
-         - pow( (*preCalculated)[0] - (*excludedMeas)[0], 2 ))
-    * *(int*)para
-    / (nmeas-1);
+  return (( (*preCalculated)[1] - pow((*excludedMeas)[0], 2) ) / (nmeas-1)
+         - pow( ( (*preCalculated)[0] - (*excludedMeas)[0] ) / (nmeas-1), 2 ))
+    * *(int*)para;

 }
 
 #endif
 }
 
 #endif