#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<4; type++) {
//---------------------------------------
switch(type){
 case 1: sprintf(cid,"meanHVchVolt_E_%02d-%02d-%d.txt",DAY,MONTH,YEAR); break;
 case 2: sprintf(cid,"meanHVchCurr_E_%02d-%02d-%d.txt",DAY,MONTH,YEAR); break;
 case 3: sprintf(cid,"meanMOch_Tm0_E_%02d-%02d-%d.txt",DAY,MONTH,YEAR); break;
}
ofstream fout(cid);
cout<<cid<<" created"<<endl;
fout<<"bec,eta,disk,phi,ndata, avr, sig, delta "<<endl;
//-------------------
char name[4][6]={"CIM","CO","AIM","AO"};
int NquadModule[4] = {10,13,10,13};
int bec,Q,disk,B302,Loop,module,eta,day,month,year,hour,min,sec,subsec;
int bec0,Q0,disk0,B3020,Loop0,module0,eta0, phi0;
//----------------------
for (int Layer=0; Layer<4 ; Layer++) {
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_%s_%02d-%02d-%d.txt",name[Layer],DAY,MONTH,YEAR);
         break;
 case 2: sprintf(cid,"dcsHVchCurr_%s_%02d-%02d-%d.txt",name[Layer],DAY,MONTH,YEAR);
         break;
 case 3: sprintf(cid,"dcsMOch_Tm0_%s_%02d-%02d-%d.txt",name[Layer],DAY,MONTH,YEAR);
         break;
}
ifstream fin(cid);
fin.getline(line,sizeof(line));   // skip 1st line
while(fin>>bec>>Q>>disk>>B302>>Loop>>eta>>module>>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 || disk0 != disk || B3020!=B302 
     || Loop0!=Loop || module0!=module || 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.));
         phi0 = NquadModule[Layer] * (Q0-1) + module0 ;
         fout<<setw(3)<<bec0<<setw(3)<<eta0<<setw(4)<<disk0
           <<setw(5)<<phi0<<setw(5)<< n<<"  "<<setw(7)<< avr<<"  "
           << setw(4)<<sig<<"  "<<setw(4)<<delta<<endl;
         if((type==1 && delta>2000) || (type==3 && delta>20.) )
           cout<<"abnormal: "<<YEAR<<setw(3)<<bec0<<setw(3)<<eta0<<setw(4)
           <<disk0 <<setw(5)<<phi0<<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 ; disk0 = disk ; B3020=B302;
     Loop0=Loop ; module0=module ; 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.));
           phi0 = NquadModule[Layer] * (Q0-1) + module0 ;
           fout<<setw(3)<<bec0<<setw(3)<<eta0<<setw(4)<<disk0
           <<setw(5)<<phi0<<setw(5)<< n<<"  "<<setw(7)<< avr<<"  "
           << setw(4)<<sig<<"  "<<setw(4)<<delta<<endl;
         if((type==1 && delta>2000) || (type==3 && delta>20.) )
           cout<<"abnormal: "<<YEAR<<setw(3)<<bec0<<setw(3)<<eta0<<setw(4)
           <<disk0 <<setw(5)<<phi0<<setw(5)<< n<<"  "<<setw(7)<< avr<<"  "
           << setw(4)<<sig<<"  "<<setw(4)<<delta<<endl;
}
fin.close();
}
fout.close();
}
return(0);
}