dNdEta_Run2023/analysis_INTT_CW/vector_stddev.h

0001 #ifndef VECTOR_STDDEV_H
0002 #define VECTOR_STDDEV_H
0003
0004 double  vector_average (vector <double> input_vector) {
0005     return accumulate( input_vector.begin(), input_vector.end(), 0.0 ) / double(input_vector.size());
0006 }
0007
0008 double vector_stddev (vector <double> input_vector){
0009
0010     double sum_subt = 0;
0011     double average  = accumulate( input_vector.begin(), input_vector.end(), 0.0 ) / double(input_vector.size());
0012
0013     // cout<<"average is : "<<average<<endl;
0014
0015     for (int i=0; i<input_vector.size(); i++){ sum_subt += pow((input_vector[i] - average),2); }
0016
0017     //cout<<"sum_subt : "<<sum_subt<<endl;
0018     // cout<<"print from the function, average : "<<average<<" std : "<<stddev<<endl;
0019
0020     return sqrt( sum_subt / double(input_vector.size()-1) );
0021 }
0022
0023 #endif