X-Git-Url: http://git.treefish.org/~alex/phys/latlib.git/blobdiff_plain/f88ac5417385884367cf8cd5fc105ec4ffda8904..027b93b16fe0d77104e0e711377c5728b0934035:/obs.hpp?ds=inline

diff --git a/obs.hpp b/obs.hpp
index 2501f14..99d97fa 100644
--- a/obs.hpp
+++ b/obs.hpp
@@ -35,6 +35,7 @@ private:
   map<string,result> computations;
   
   void mean(const string& compid, vector< vector <double> > *meas, const int& ival);
+  void mean(const string& compid, vector< vector <int> > *meas, const int& ival);
 };
 
 
@@ -75,6 +76,24 @@ void obs<obstype>::mean(const string& compid, vector< vector<double> > *meas, co
     computations[compid].err += pow( (*measIt)[ival] - computations[compid].val, 2 );
   computations[compid].err /= nmeas*(nmeas-1);
   computations[compid].err = sqrt( computations[compid].err );	 
+}
+
+template <typename obstype>
+void obs<obstype>::mean(const string& compid, vector< vector<int> > *meas, const int& ival)
+{
+  computations[compid].val = 0;
+  computations[compid].err = 0;
+  
+  int nmeas = meas->size();
+
+  for(vector< vector<int> >::iterator measIt = meas->begin(); measIt != meas->end(); ++measIt)
+    computations[compid].val += (*measIt)[ival];
+  computations[compid].val /= nmeas;
+  
+  for(vector< vector<int> >::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 );	 
 }				 
     
 template <typename obstype>