#include <iostream>
#include <cmath>
#include <fstream>
using namespace std;

int main()
{
char cid[30];
int Layer=3;
int bec,Q,layer,B302,Loop,pipe,eta,day,month,year,hour,min,sec,subsec;
int bec0,Q0,layer0,B3020,Loop0,pipe0,eta0;
int n=0;
double dcsvalue,sum=0., sum2=0., avr,sig;
char line[255]; 
//----------------------------------------
sprintf(cid,"meanHVchCurr_B%d.txt",Layer);
ofstream fout(cid);
fout<<"bec,layer,pipe, eta, n, avr, sig "<<endl;
//------------------readin data ---------------------------
sprintf(cid,"dcsHVchCurr_B%d.txt",Layer);
ifstream fin(cid);
fin.getline(line,sizeof(line));   // skip 1st line
while(fin>>bec>>Q>>layer>>B302>>Loop>>pipe>>eta>>dcsvalue>>day>>month>>year>>hour>>min>>sec>>subsec)
{
if(bec0 != bec && Q0!=Q || layer0 != layer || B3020!=B302 
     || Loop0!=Loop || pipe0!=pipe || eta0 != eta) {
     if(n != 0) {
         avr = sum/n; sig = 0.;
         if (n > 1) sig = sqrt((sum2 - n*avr*avr)/(n-1.));
         fout<<bec0<<" "<<layer0<<" "<<pipe0<<" "<< eta0<<"    "<< n<<"   "<< avr<<"   "<< sig<<endl;
     n=0; sum=0.; sum2=0.;
     }
     bec0 = bec ; Q0=Q ; layer0 = layer ; B3020=B302; 
     Loop0=Loop ; pipe0=pipe ; eta0 = eta;n++;
     sum += dcsvalue; sum2 += dcsvalue* dcsvalue;
} else {
    sum += dcsvalue;
    sum2 += dcsvalue* dcsvalue;
   n++;
}
}
if(n != 0) {
         avr = sum/n; sig = 0.;
         if (n > 1) sig = sqrt((sum2 - n*avr*avr)/(n-1.));
         fout<<bec0<<" "<<layer0<<" "<<pipe0<<" "<< eta0<<"    "<< n<<"   "<< avr<<"   "<< sig<<endl; }
fin.close();
fout.close();
cout<<cid<<" created"<<endl;
return(0);
}