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

int main()
{
int YEAR=2014, MONTH=11, DAY=19;
char cid[30];
for (int type=1; type<5; type++) {
//---------------------------------------
switch(type){
 case 1: sprintf(cid,"meanHVchVolt_B_%02d-%02d-%d.txt",DAY,MONTH,YEAR); break;
 case 2: sprintf(cid,"meanHVchCurr_B_%02d-%02d-%d.txt",DAY,MONTH,YEAR); break;
 case 3: sprintf(cid,"meanMOch_Tm0_B_%02d-%02d-%d.txt",DAY,MONTH,YEAR); break;
 case 4: sprintf(cid,"meanMOch_Tm1_B_%02d-%02d-%d.txt",DAY,MONTH,YEAR); break;
}
ofstream fout(cid);
cout<<cid<<" created"<<endl;
fout<<"bec,layer,phi, eta, n, avr, sig, delta "<<endl;
//-------------------
for (int Layer=3; Layer<7 ; Layer++) {
int bec,Q,layer,B302,Loop,phi,eta,day,month,year,hour,min,sec,subsec;
int bec0,Q0,layer0,B3020,Loop0,phi0,eta0;
int n=0;
double dcsvalue,sum=0., sum2=0., max=-999., mini=99999999.,avr, sig,delta;
char line[255]; 
//----------------------------------------
//------------------readin data ---------------------------
switch(type){
 case 1: sprintf(cid,"dcsHVchVolt_B%d_%02d-%02d-%d.txt",Layer,DAY,MONTH,YEAR);
         break;
 case 2: sprintf(cid,"dcsHVchCurr_B%d_%02d-%02d-%d.txt",Layer,DAY,MONTH,YEAR);
         break;
 case 3: sprintf(cid,"dcsMOch_Tm0_B%d_%02d-%02d-%d.txt",Layer,DAY,MONTH,YEAR);
         break;
 case 4: sprintf(cid,"dcsMOch_Tm1_B%d_%02d-%02d-%d.txt",Layer,DAY,MONTH,YEAR);
         break;
}
ifstream fin(cid);
fin.getline(line,sizeof(line));   // skip 1st line
while(fin>>bec>>Q>>layer>>B302>>Loop>>phi>>eta>>dcsvalue>>day>>month>>year>>hour>>min>>sec>>subsec)
{
//remove bad dcsdata (Volt and Curr only)
if((type==1 || type==2) && dcsvalue==0 )continue;  
if(bec0 != bec && Q0!=Q || layer0 != layer || B3020!=B302 
     || Loop0!=Loop || phi0!=phi || eta0 != eta) {
     if(n != 0) {
         avr = sum/n; sig = 0.; delta=max-mini;
         if (n > 1 && delta > 0. ) sig = sqrt((sum2 - n*avr*avr)/(n-1.));
           fout<<setw(3)<<bec0<<setw(3)<<layer0<<setw(4)<<phi0
           <<setw(5)<< eta0<<setw(5)<< n<<"  "<<setw(7)<< avr<<"  "
           << setw(4)<<sig<<"  "<<setw(4)<<delta<<endl; 
         if(type==1 && delta>2000)
           cout<<"abnormal: "<<YEAR<<setw(3)<<bec0<<setw(3)<<layer0
           <<setw(4)<<phi0 <<setw(5)<< eta0<<setw(5)<< n<<"  "<<setw(7)
           << avr<<"  " << setw(4)<<sig<<"  "<<setw(4)<<delta<<endl;
         n=0; sum=0.; sum2=0.; max=-999.; mini=99999999.;
     }
     bec0 = bec ; Q0=Q ; layer0 = layer ; B3020=B302; 
     Loop0=Loop ; phi0=phi ; eta0 = eta;
     n++; sum += dcsvalue; sum2 += dcsvalue* dcsvalue;
     if(dcsvalue>max) max = dcsvalue;
     if(dcsvalue<mini)mini= dcsvalue;
} else {
     sum += dcsvalue;
     sum2 += dcsvalue* dcsvalue;
     if(dcsvalue>max) max = dcsvalue;
     if(dcsvalue<mini)mini= dcsvalue;
     n++;
}
}
if(n != 0) {
    avr = sum/n; sig = 0.; delta=max-mini;
    if (n > 1 && delta > 0. ) sig = sqrt((sum2 - n*avr*avr)/(n-1.));
           fout<<setw(3)<<bec0<<setw(3)<<layer0<<setw(4)<<phi0
           <<setw(5)<< eta0<<setw(5)<< n<<"  "<<setw(7)<< avr<<"  "
           << setw(4)<<sig<<"  "<<setw(4)<<delta<<endl; 
         if(type==1 && delta>2000)
           cout<<"abnormal: "<<YEAR<<setw(3)<<bec0<<setw(3)<<layer0
           <<setw(4)<<phi0 <<setw(5)<< eta0<<setw(5)<< n<<"  "<<setw(7)
           << avr<<"  " << setw(4)<<sig<<"  "<<setw(4)<<delta<<endl;
}
fin.close();
}
fout.close();
}
return(0);
}