From: Alex Schmidt Date: Tue, 29 May 2012 10:09:56 +0000 (+0200) Subject: wrong variances\! X-Git-Url: http://git.treefish.org/~alex/phys/latlib.git/commitdiff_plain/759171112bbb5c0121f2ba762a2b1d9ad96a271d wrong variances\! --- diff --git a/obs.hpp b/obs.hpp index 99d97fa..dd9b23d 100644 --- a/obs.hpp +++ b/obs.hpp @@ -74,7 +74,7 @@ void obs::mean(const string& compid, vector< vector > *meas, co for(vector< vector >::iterator measIt = meas->begin(); measIt != meas->end(); ++measIt) computations[compid].err += pow( (*measIt)[ival] - computations[compid].val, 2 ); - computations[compid].err /= nmeas*(nmeas-1); + computations[compid].err /= nmeas-1; computations[compid].err = sqrt( computations[compid].err ); } @@ -92,8 +92,8 @@ void obs::mean(const string& compid, vector< vector > *meas, const for(vector< vector >::iterator measIt = meas->begin(); measIt != meas->end(); ++measIt) computations[compid].err += pow( (*measIt)[ival] - computations[compid].val, 2 ); - computations[compid].err /= nmeas*(nmeas-1); - computations[compid].err = sqrt( computations[compid].err ); + computations[compid].err /= nmeas-1; + computations[compid].err = sqrt( computations[compid].err ); } template