--- rawreader/src/Read.cpp	2006/12/10 12:21:23	1.2
+++ rawreader/src/Read.cpp	2007/01/28 10:55:47	1.5
@@ -6,79 +6,90 @@
 short int CRC(int []);
 int restrict(int [], long long int[], int);
 
+extern char database[120];
+extern char MYSQL[120];
+extern int print_flag;
+
+int route_good1=0, route_good2=0, route_good3=0;
+int pack_good1=0, pack_good2=0, pack_good3=0;
+int number_good1=0, number_good2=0, number_good3=0;
+
 void read(char fni[], char fno[], char INF[],
-		  char protocol[], char path[], char name[], char nnnnn_TMP[], float quality[], int real_route[], int &Rcount)
+          char protocol[], char path[], char name[], char nnnnn_TMP[], float quality[], int real_route[], int &Rcount)
 {
- FILE *fo, *fout, *foRout, *foMiss, *fo_log_kadr, *fo_route;
+ FILE *fo, *fout, *foRout, *foMiss, *fo_log_kadr, *fo_route, *fo_tmp;
  ifstream fin;
- ofstream fo_cln1, fokadr_cln1;
+ ofstream fo_cln1, fokadr_cln1, fo_info;
  
  char buf[1024];
  //  char stop1;
- char ch1,ch2,ch3, ch;
- char file_route[80];  /*protocol binary filesm mmm_counter_nnnnn_cln1Er.log*/
- char miss_file[80];   /*protocol files mmm_counter_nnnnn_mis_cln1.dat*/
- char cln1_file[80];   /*route binary data files mmm_counter_nnnnn_cln1.pam*/
- char cln1_kadr[80];   /*route binary data files mmm_counter_nnnnn_cln1.bin*/
- char log_kadr[80];    /*log binary data files mmm_counter_nnnnn_cln1.log*/
- char log_route[80];   /*log binary file for number of routes*/
+ char ch1,ch2,ch3,ch, ch4[80], ch5[5];
+ char file_route[100];  /*protocol binary filesm nnnnn_mmm_xxx_cln1Er.log*/
+ char miss_file[100];   /*protocol files nnnnn_mmm_xxx_mis_cln1.dat*/
+ char cln1_file[100];   /*route binary data files nnnnn_mmm_xxx_cln1.pam*/
+ char cln1_kadr[100];   /*route binary data files nnnnn_mmm_xxx_cln1.bin*/
+ char log_kadr[100];    /*log binary data files nnnnn_mmm_xxx_cln1.log*/
+ char log_route[100];   /*log binary file for number of routes*/
  char Short_Buf[8], Name_Buf[4]; 
+ char temp[10], temp1[80], connection[120], command[80], name11[80];
  short int CalcCheckSum;
 
  int LenthBuf=1015, flag_help=0;
  int i, j, k, file_length=0, Hex[3];
- int N_er=0;              /*Number of kadr with name error*/
- int CRC0_er=0;           /*Number of kadr with CRC0 error*/
- int P_er=0;              /*Number of kadr with packet error number*/
- int M_er=0;              /*Number of kadr with rout error number*/
+ int N_er=0;              /*Number of cadre with name error*/
+ int CRC0_er=0;           /*Number of cadre with CRC0 error*/
+ int P_er=0;              /*Number of cadre with packet error number*/
  int end=0;               /*End of last pamela packet*/
- // Buffers for Kadr analysis
+ // Buffers for Cadre analysis
  int buffer1[1024], buffer2[1024], buffer3[1024]; 
  int rr=0;
-
+ int cirquit_no=0, route=0;
+ 
  int header1, header2, header3, CRC_0[3], CRC_1[3];
  long long int icount_int;
  long long int icount_frac, tmpP, numP[3];	      
- long int N1=0;                  /* Number of identified kadrs*/
- long long int icount=0;         /* Number of bytes in identified kadrs*/ 
+ long int N1=0;                  /* Number of identified cadres*/
+ long long int icount=0;         /* Number of bytes in identified cadres*/ 
  long long int icount_total=0;   /* Total number of bytes*/
- long long int packet_first[100];/*First kadr number for each route*/
- long long int packet_last[100]; /*Second kadr number for each route*/
- long int count_missed[100];     /* Number of missed kadrs in each route*/
- long int count_sec_missed[100];
- unsigned long long int count_ident[100]; /* Number of good bytes in each rout*/
- unsigned long long int count_tot[100];   /* Total number of bytes in each rout*/
- int NumRout[1000];              /* Number of Route*/
- long long int Num_kadrs[1000];  /* Number of identified kadrs in each rout*/
- long int Num_good_kadrs[1000];  /* Number of normal kadrs in each rout*/
- long long int Num_bad_kadrs[1000]; /* Bad kadrs number in each rout*/
-// float quality[100];           /* File quality pointer*/
- int counter=0;                  /* Last byte position in truncated kadr*/
- int adr1=0;                     /* Initial byte of kadr in cln1.bin*/
- int adr2=0;                     /* Last byte of kadr in cln2.bin*/
+ long long int packet_first[100];/* First cadre number for each route*/
+ long long int packet_last[100]; /* Second cadre number for each route*/
+ long long int count_missed[100];/* Number of missed cadres in each route*/
+ long long int count_sec_missed[100];
+ unsigned long long int count_ident[100]; /* Number of good bytes in each route*/
+ unsigned long long int count_tot[100];   /* Total number of bytes in each route*/
+ int NumRout[1000];              /* Number of route*/
+ long long int Num_kadrs[1000];  /* Number of identified cadres in each route*/
+ long int Num_good_kadrs[1000];  /* Number of normal cadres in each route*/
+ long long int Num_bad_kadrs[1000]; /* Bad cadres number in each route*/
+
+ int counter=0;                  /* Last byte position in truncated cadre*/
+ int adr1=0;                     /* Initial byte of cadre in cln1.bin*/
+ int adr2=0;                     /* Last byte of cadre in cln2.bin*/
  int tmp=0;
- int numPack[4];                 /* Number of kadr in given route*/
- int good_pack=0;
- int good_number=0;
- int good_route=0;
+ int numPack[4];                 /* Number of cadre in given route*/
+
  int adress[4];
- int R_er[100];           /* Rout errors number in given rout*/
- int R_er_interf1[100];   /* Interface Rout error at begin of rout*/
- int R_er_interf2[100];   /* Interface Rout error at end of rout*/
- int Mis_group[1000];     /* Number of missed kadrs group
+ int R_er[100];           /* Route errors number in given route*/
+ int R_er_interf1[100];   /* Interface Route error at the begining of the route*/
+ int R_er_interf2[100];   /* Interface Route error at end of route*/
+ int Mis_group[1000];     /* Number of missed cadres group
                              for each route*/
+ float first_missed[100];			     
  int Pack_Er[1000];       /* Number of packets with packet number
                              error for each route*/
- int CRC_bad_kadr[3];     /* Current CRC bad kadr number*/
+ int CRC_bad_kadr[3];     /* Current CRC bad cadre number*/
  int flag_pack[3];        /* Flag for error in packet number*/
- int flag_rout[3];        /* Flag for error in rout number*/
+ int flag_rout[3];        /* Flag for error in route number*/
  int flag_name[3];        /* Flag for name error*/
  int flag_CRC0[3];        /* Flag for CRC0 error*/
  int flag_CRC[3];         /* Flag for CRC error*/
  int flag_CRC_cor[3];     /* Flag for CRC error correction*/
  int flag_tr[3];          /* Truncate flag*/
  short int flags;         /* Number for all errors*/
-
+ int num_rows=0;
+ int control[1];
+ control[0]=0;
+ control[1]=0;
  for (j=0; j<4; j++) {adress[j]=0; numPack[j]=0;}
  for (j=0; j<1024; j++) 
  {buffer1[j]=0; buffer2[j]=0; buffer3[j]=0;}
@@ -89,6 +100,7 @@
   packet_last[i]=0;
   count_missed[i]=0;
   count_sec_missed[i]=0;
+  first_missed[i]=0;
   count_ident[i]=0;
   count_tot[i]=0;
   R_er[i]=0;
@@ -102,29 +114,43 @@
   flag_pack[i]=0; flag_rout[i]=0;
   flag_name[i]=0; flag_CRC0[i]=0; flag_tr[i]=0;
   flag_CRC[i]=0; flag_CRC_cor[i]=0;} 
+
+ stringstream oss;
+
+ strcpy(connection,MYSQL);
+ strcat(connection,database);
+ TSQLServer *sqlServer = TSQLServer::Connect(connection,"pamelaprod","rawprod");
+ if (sqlServer==0) {cerr<<"can not connect with MYSQL sever\n"; exit(1);}
  // open input file "nnnnnmmm.pam"
- fin.open(fni, ifstream::binary);
- if(fin.good()==0) {cerr<<"can not open input file"; exit(0);}
+ fin.open(fni, ifstream::binary); 
+ if(fin.good()==0) {cerr<<"can not open input file"; exit(1);}
  //open output file "nnnnnmmm_work.dat"
  fout=fopen(fno, "wb");
- if(fout==0) {cerr<<"can not open outputfile"; exit(0);}
+ if(fout==0) {cerr<<"can not open work file"; exit(1);}
+
+ strcpy(ch5,"info");
+ strcpy(ch4,fno);
+ for(i=0;i<4;i++){ch4[strlen(ch4)-8+i]=ch5[i];} 
+ // open input file "nnnnnmmm_info.dat"
+ fo_info.open(ch4, ifstream::binary); 
+ if(fo_info.good()==0) {cerr<<"can not open info file"; exit(1);} 
+
  //open output file "nnnnnmmm_protocol.txt"
  fo=fopen(protocol, "wb");
- if(fo==0) {cerr<<"can not open protocol file"; exit(0);}
+ if(fo==0) {cerr<<"can not open protocol file"; exit(1);}
  //open output file "nnnnnmmm_route.log"
  strncpy(log_route,INF,strlen(INF)-9);
  log_route[strlen(INF)-9]='\0';
  strcat(log_route,"_route.log");
  fo_route=fopen(log_route, "wb");
- if(fo_route==0) {cerr<<"can not open file for routes number"; exit(0);}
+ if(fo_route==0) {cerr<<"can not open file for routes number"; exit(1);}
  /////////////////////////////////////////////////////////////////
  // cin>>stop1;
- k=0;
+ k=0; 
  while(fin.eof()==0)
  {
-  counter=0;          //pounter for position of kadr truncation
+  counter=0;          //pounter for position of cadre truncation
   N_er=0;           
-  M_er=0;
   P_er=0;
   CRC0_er=0;
   if (N1==0) Rcount=1;
@@ -143,20 +169,24 @@
   count_ident[Rcount-1]+=8; count_tot[Rcount-1]+=8;
 
   while((ch1!=(char)(CODE1)&&ch2!=(char)(CODE2)&&
-  (N1==0||buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])||buffer3[3]!=(int)(unsigned char)(Short_Buf[0])))
+  (N1==0||(buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])||buffer3[3]!=(int)(unsigned char)(Short_Buf[0]))))
 	     ||
   (ch2!=(char)(CODE2)&&ch3!=(char)(CODE3)&&(buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])||(numP[2]+1)!=tmpP))
 	     ||
-  (ch3!=(char)(CODE3)&&ch1!=(char)(CODE1)&&((numP[2]+1)!=tmpP||buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])))
+  (ch3!=(char)(CODE3)&&ch1!=(char)(CODE1)&&(buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])||
+   (buffer3[5]!=(int)(unsigned char)(Short_Buf[2])&& buffer3[6]!=(int)(unsigned char)(Short_Buf[3]))))
 	     ||
   (ch1!=(char)(CODE1)&&
-   (buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])||buffer3[3]!=(int)(unsigned char)(Short_Buf[0])||(numP[2]+1)!=tmpP))
+   (buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])&& buffer3[3]!=(int)(unsigned char)(Short_Buf[0])&&
+    buffer3[5]!=(int)(unsigned char)(Short_Buf[2])&& buffer3[6]!=(int)(unsigned char)(Short_Buf[3])))
   	     ||
   (ch2!=(char)(CODE2)&&
-   (buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])||buffer3[3]!=(int)(unsigned char)(Short_Buf[0])||(numP[2]+1)!=tmpP))
+   (buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])&& buffer3[3]!=(int)(unsigned char)(Short_Buf[0])&&
+    buffer3[5]!=(int)(unsigned char)(Short_Buf[2])&& buffer3[6]!=(int)(unsigned char)(Short_Buf[3]))) 
 	     ||
   (ch3!=(char)(CODE3)&&
-   (buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])||buffer3[3]!=(int)(unsigned char)(Short_Buf[0])||(numP[2]+1)!=tmpP))
+   (buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])&& buffer3[3]!=(int)(unsigned char)(Short_Buf[0])&&
+    buffer3[5]!=(int)(unsigned char)(Short_Buf[2])&& buffer3[6]!=(int)(unsigned char)(Short_Buf[3])))  
 	     ||
   (ch1!=(char)(CODE1)&&ch2!=(char)(CODE2)&&ch3!=(char)(CODE3)))
   {fin.get(ch);
@@ -171,13 +201,13 @@
   if ((fin.eof())&&(N1==0)) 
   {
    fclose(fout);
+   fo_info.close();   
    remove(fno);
    fclose(fo);
    remove(protocol);
    remove(INF);
-   exit(0);
+   exit(1);
   }
-
   /*****************************************************************/
   /*********************data reading********************************/
   fin.read(buf,LenthBuf); 
@@ -222,7 +252,7 @@
 		   CRC_bad_kadr[2]=0;
 		   flag_pack[2]=0; flag_rout[2]=0;
 		   flag_name[2]=0; flag_CRC0[2]=0; flag_tr[2]=0;
-		   flag_CRC[2]=0; flag_CRC_cor[2]=0;
+		   flag_CRC[2]=0; flag_CRC_cor[2]=0;	   
            for (i=0; i<2; i++)
 		   {numP[i]=numP[2];CRC_bad_kadr[i]=CRC_bad_kadr[2];
 		    flag_pack[i]=flag_pack[2]; flag_rout[i]=flag_rout[2];
@@ -236,6 +266,7 @@
 		   for (i=0; i<2; i++)
 		   {CRC_1[i]=CRC_1[2];}
            packet_first[Rcount-1]=numP[2];
+	   control[0]=0; control[1]=0;
 		   break;
           }
    default:{
@@ -273,29 +304,30 @@
             CRC_0[2]=(int)(unsigned char)(Short_Buf[4]);
 		    for (i=0; i<2; i++) {CRC_1[i]=CRC_1[i+1];}
             CRC_1[2]=(int)(unsigned char)(ch);    
+	    control[0]=control[1];
+	    control[1]=0;
            }
   }  
-//  if (fmod(N1,10000)==0.) cout<<N1<<"\n";	  
+  //  if (print_flag>0 && fmod(N1,10000)==0) cout<<N1<<"\n";	  
+  if (print_flag>0 && !(N1%10000)) cout<<N1<<"\n";
 
 /******************CRC ERROR ANALYSIS***************************/
   CalcCheckSum=CRC(buffer3);
   if (CalcCheckSum!=(short int)(buffer3[1023]))
-  {
-   /*if (numP[2]==1026) 
-   {for (i=0;i<8;i++) printf(" %x ",buffer3[i]); printf("\n");
-    for (i=8;i<16;i++) printf(" %x ",buffer3[i]); printf("\n");
-    for (i=1000;i<1008;i++) printf(" %x ",buffer3[i]); printf("\n");
-	printf("%x %x %x %x",buffer3[1020],buffer3[1021],buffer3[1022],buffer3[1023]);
-    cin>>i;}*/
-   fprintf(fout,"CRC = %x (%x) error in kadr %i\n", buffer3[1023],CalcCheckSum,numP[2]);     //write in file nnnnn_work.dat
+  {    
+//write in file nnnnnmmm_info.dat    
+   fo_info<<"real CRC = "<<hex<<buffer3[1023]<<" (calculated CRC ="<<CalcCheckSum<<") in cadre "<<dec<<numP[2]<<"\n";
+//write in file nnnnnmmm_work.dat       
+   fprintf(fout,"%10i %4i %4i\n", numP[2], CalcCheckSum, buffer3[1023]);
 
    flag_CRC[2]=1;
    Num_bad_kadrs[Rcount-1]+=1;
    CRC_bad_kadr[2]=numP[2];
    counter=restrict(buffer3,numP,buffer2[3]);
    if (counter!=0)
-   {fprintf(fout," truncate byte=%i",counter);
-//    cout<<" truncate counter="<<counter<<" route="<<buffer2[3]<<" "<<buffer3[3]<<" kadr="<<numP[2]<<" "<<numP[1]<<" "<<numP[0]<<"\n";
+   {fo_info<<" truncate byte="<<counter<<" in cadre "<<numP[2]<<"\n";
+    if (print_flag>0) 
+    {cout<<" truncate counter="<<counter<<" route="<<buffer2[3]<<" "<<buffer3[3]<<" cadre="<<numP[2]<<" "<<numP[1]<<" "<<numP[0]<<"\n";}
     fin.seekg(-(1024-counter),ios::cur);
     icount-=1024;
     count_ident[Rcount-1]-=1024;
@@ -319,7 +351,7 @@
   {flag_CRC0[1]=1; CRC0_er=N1-1;}
   /****************************************************************/
 
-  /******************KADR NUMBER ANALYSIS**************************/
+  /******************CADRE NUMBER ANALYSIS*************************/
   if ((numP[1]!=(numP[0]+1))&&(numP[1]!=(numP[2]-1))&&(N1!=1)&&(CRC_bad_kadr[1]==numP[1])&&(numP[1]!=0)) 
   {P_er=N1-1; flag_pack[1]=1;
   //Correct error packet number//
@@ -329,7 +361,7 @@
    {if (tmpP/(long int)(pow(256.,2-i))>=1) 
     {Hex[2-i]=(int)(tmpP/(long int)(pow(256.,2-i)));
      tmpP=tmpP-(long int)pow(256.,2-i)*(int)(tmpP/(long int)(pow(256.,2-i)));}}
-   for (i=4;i<7;i++) (buffer2[i]=Hex[i-4]);
+//   for (i=4;i<7;i++) (buffer2[i]=Hex[i-4]);
    if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
 
    tmpP=numP[2]-1;
@@ -338,70 +370,91 @@
    {if (tmpP/(long int)(pow(256.,2-i))>=1) 
     {Hex[2-i]=(int)(tmpP/(long int)(pow(256.,2-i)));
 	 tmpP=tmpP-(long int)pow(256.,2-i)*(int)(tmpP/(long int)(pow(256.,2-i)));}}
-   for (i=4;i<7;i++) (buffer2[i]=Hex[i-4]);
+//   for (i=4;i<7;i++) (buffer2[i]=Hex[i-4]);
    if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
   }
   /////////////////////////////////////////////////////////////////////
   if ((numP[2]!=(numP[1]+1))&&(N1!=1)&&(buffer3[3]==buffer2[3])&&(CRC_bad_kadr[2]==numP[2])&&(numP[2]!=0)) flag_pack[2]=1;
   if (flag_pack[2]!=1) packet_last[Rcount-1]=numP[2];
-  ////////Missing kadr analysis////////////////////////////////////
-  /**route number: 10 (10 19 10) 10 ******************************/
-  /**kadr number:  21 (22 17 18) 19 ******************************/
-  /**missed kadrs: 21 20 19 18 17   ******************************/  
+  ////////Missing cadre analysis////////////////////////////////////
+  /**route number:  10 (10 19 10) 10 ******************************/
+  /**cadre number:  21 (22 17 18) 19 ******************************/
+  /**missed cadres: 21  20 19 18  17 ******************************/  
   if ((numP[1]<numP[0])&&(numP[1]==numP[2]-1)&&(buffer1[3]==buffer3[3])&&(flag_tr[1]==0)&&
- 	  (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0))
+ 	  (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0)&&(control[0]==0))
   {Mis_group[Rcount-1]+=1;
-   fprintf(foMiss,"secondary transmitted kadres may be %9i%9i\n",numP[1]+1,numP[0]-1);
-   count_sec_missed[Rcount-1]+=numP[0]-1-numP[1]-1;}  
+   count_sec_missed[Rcount-1]+=numP[0]-1-numP[1];
+   fprintf(foMiss,"secondary transmitted cadres may be %lli - %lli number of cadres=%lli\n",
+           numP[1],numP[0],numP[0]-1-numP[1]);
+  }  
   /**route number: 10 (10 19 10) 10 ******************************/
-  /**kadr number:  11 (12 15 20) 21 ******************************/
-  /**missed kadrs:    16 17 18 19   ******************************/  
+  /**cadre number: 11 (12 15 20) 21 ******************************/
+  /**missed cadres:    16 17 18 19   ******************************/  
   if ((numP[1]>numP[0]+1)&&(numP[1]<numP[2]-1)&&(buffer1[3]==buffer3[3])&&(flag_tr[1]==0)&&
- 	  (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0))
+      (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0))
   {Mis_group[Rcount-1]+=1;
-   fprintf(foMiss,"%9i%9i\n",numP[1]+1,numP[2]-1);
+   fprintf(foMiss,"1) %9lli%9lli\n",numP[1]+1,numP[2]-1);
+   if (print_flag>0) cout<<"1) "<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";      
    count_missed[Rcount-1]+=numP[2]-numP[1]-1;}
   /**route number: (10 10 10) 10 10 ******************************/
-  /**kadr number:  (11 12 15) 20 21 ******************************/
-  /**missed kadrs:     13 14        ******************************/  
-  else if ((numP[1]>=numP[0]+1)&&(numP[1]<numP[2]-1)&&
-	       (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&
-		   (flag_tr[1]==0)&&(flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))
-  {Mis_group[Rcount-1]+=1;
-   fprintf(foMiss,"%9i%9i\n",numP[0]+2,numP[2]-1);
-   count_missed[Rcount-1]+=numP[2]-numP[0]-1-1;}
-  /**route number: (10 10 10) 10 10 ******************************/
-  /**kadr number:  (11 29 15) 14 15 ******************************/
-  /**number of missed kadrs:2       ******************************/  
+  /**cadre number: (11 29 15) 14 15 ******************************/
+  /**number of missed cadres:2       ******************************/  
   else if (((numP[1]<numP[0]+1)||(numP[1]>numP[2]-1))&&(numP[2]-numP[0]>2)&&
            (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&(flag_tr[1]==0)&&
      	   (flag_pack[2]==0)&&(flag_pack[1]==1)&&(flag_pack[0]==0))
   {Mis_group[Rcount-1]+=1;
-   fprintf(foMiss,"%9i%9i -1\n",numP[0]+1,numP[2]-1);
+   fprintf(foMiss,"2) %9lli%9lli\n",numP[0]+1,numP[2]-1);
+   if (print_flag>0) cout<<"2) "<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";   
    count_missed[Rcount-1]+=numP[2]-numP[0]-1-1;}
-  /**route number:10 (10 10 10) 10 10 ******************************/
-  /**kadr number:  1 (11 12 13) 14 15 ******************************/
-  /**number of missed kadrs: 10       ******************************/   
-  else if ((N1>1)&&(abs(numP[2]-good_pack-2)>abs(Num_kadrs[Rcount-1]-good_number))&&
-           (buffer2[3]==good_route)&&(flag_tr[2]==0)&&(flag_pack[2]==0))
+  else if (((numP[1]<numP[0]+1)||(numP[1]>numP[2]-1))&&((numP[2]-numP[0])==2)&&
+           (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&(flag_tr[1]==0))
+  {control[1]=1;}   
+  /**route number: (10 10 10) 10 10 ******************************/
+  /**cadre number: (11 15 16) 17 18 ******************************/
+  /**missed kadrs:  12 13 14        ******************************/  
+  else if ((numP[1]>numP[0]+1)&&(numP[1]<=numP[2]-1)&&(buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&
+	   (flag_tr[1]==0)&&(flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0)&&(control[0]==0))
   {Mis_group[Rcount-1]+=1;
-   fprintf(foMiss,"%9i%9i \n",good_pack+1,numP[2]-1);
-   count_missed[Rcount-1]+=numP[2]-good_pack-1;
-  }
-  //truncate kadr
-  /**route number: 10 (10 10 10) 10 ******************************/
-  /**kadr number:  11 (12 13 14) 15 ******************************/
-  /**missed truncate kadr:13        ******************************/  
+   fprintf(foMiss,"3) %9lli%9lli\n",numP[0]+1,numP[1]-1);
+   if (print_flag>0) cout<<"3) "<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";      
+   count_missed[Rcount-1]+=numP[1]-numP[0]-1;}
+  /**route number:(10 10 )10 10 10 10 ******************************/
+  /**cadre number:( 1 11 )12 13 14 15 ******************************/
+  /**number of missed cadres: 10      ******************************/  
+  else if ((N1>2)&&(abs(numP[1]-pack_good3)>abs(Num_kadrs[Rcount-1]-number_good3-1))&&
+           (buffer2[3]==route_good3)&&(flag_tr[2]==0)&&(flag_pack[2]==0)&&(numP[2]>numP[0]+2)&&(control[0]==0))
+  {Mis_group[Rcount-1]+=1;
+   if (pack_good3<numP[2]) count_missed[Rcount-1]+=abs(numP[2]-pack_good3)-1;
+   fprintf(foMiss,"4) %9lli%9lli number of missed cadres=%i\n", pack_good3, numP[1], abs(numP[1]-pack_good3-1));
+   if (print_flag>0) cout<<"P1="<<numP[0]<<" P2="<<numP[1]<<" P3="<<numP[2]<<"\n";
+//   cout<<"pack_good1="<<pack_good1<<" pack_good2="<<pack_good2<<" pack_good3="<<pack_good3<<"\n";
+//   cout<<"number_good1="<<number_good1<<" number_good2="<<number_good2<<" number_good3="<<number_good3<<"\n";
+  }  
+  //truncate cadres
+  /**route number: 10 (10 10 10) 10  ******************************/
+  /**cadre number:  11 (12 13 14) 15 ******************************/
+  /**missed truncate cadre:13        ******************************/  
   else if ((numP[0]<numP[2]-1)&&(buffer1[3]==buffer3[3])&&
 	       (flag_tr[0]==0)&&(flag_tr[1]==1)&&(flag_tr[2]==0)&&
 		   (flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0)) 
   {Mis_group[Rcount-1]+=1;
-   fprintf(foMiss,"%9i%9i\n",numP[0]+1,numP[2]-1);
+   fprintf(foMiss,"5) %9lli%9lli\n",numP[0]+1,numP[2]-1);
    count_missed[Rcount-1]+=numP[2]-numP[0]-1;}
   /***************************************************************/
 
-  if (flag_CRC[2]!=1) {good_pack=numP[2]; good_route=buffer3[3]; good_number=Num_kadrs[Rcount-1];}
-
+  if (flag_CRC[2]!=1)
+  {
+   pack_good1=pack_good2;
+   pack_good2=pack_good3;
+   route_good1=route_good2;
+   route_good2=route_good3;
+   number_good1=number_good2;
+   number_good2=number_good3;    
+
+   pack_good3=numP[2];      
+   route_good3=buffer3[3];
+   number_good3=Num_kadrs[Rcount-1];  
+  }
   /******************ROUTE NUMBER ANALYSIS**************************/
   if ((buffer2[3]!=buffer1[3])&&(buffer2[3]!=buffer3[3])&&(buffer1[3]!=buffer3[3])&&
       (flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))
@@ -410,14 +463,13 @@
   }
   /**************First route***************************************/
   /**route number: (0 4) 4*****************************************/
-  /**kadr number:  (0 1) 2*****************************************/
+  /**cadre number: (0 1) 2*****************************************/
   if (N1==2)
   {packet_first[Rcount-1]=numP[1];
    packet_last[Rcount-1]=numP[2];
    NumRout[Rcount-1]=1;
    // CREATE FILE NAMES
-   filename(file_route, miss_file, cln1_file, path,
-            NumRout[Rcount-1], name);			 
+   filename(file_route, miss_file, cln1_file, path, NumRout[Rcount-1], name);
    strncpy(cln1_kadr,cln1_file,strlen(cln1_file)-3);
    cln1_kadr[strlen(cln1_file)-3]='\0';
    strcat(cln1_kadr,"bin");
@@ -425,28 +477,51 @@
    log_kadr[strlen(file_route)-6]='\0';
    strcat(log_kadr,"_bin.log");
    //Open files for new route
-   foRout=fopen(file_route, "wb");          //xxx_mmm_nnnnn_cln1Er.log
-   foMiss=fopen(miss_file, "wb");           //xxx_mmm_nnnnn_mis_cln1.dat 
-   fo_cln1.open(cln1_file, ios::binary);    //xxx_mmm_nnnnn_cln1.pam
+   foRout=fopen(file_route, "wb");          //nnnnn_mmm_xxx_cln1Er.log
+   foMiss=fopen(miss_file, "wb");           //nnnnn_mmm_xxx_mis_cln1.dat 
+   if (CLN1_MEPHI>0)
+   {strcpy(ch4,"");
+    strncat(ch4,path,strlen(path)-6);
+    strcat(strcpy(temp1,ch4),"Cln1/tmp.dat");
+    if ((fo_tmp=fopen(temp1,"w"))==NULL)
+    {strcat(strcat(strcpy(command,"mkdir "),ch4),"Cln1");
+     //Creata special directory for cln1 files
+     system(command);
+    }
+    else
+    {
+     fclose(fo_tmp);
+     remove(temp1);
+    }
+    for(i=0;i<22;i++)
+    {name11[i]=cln1_file[i+strlen(cln1_file)-22];}
+    name11[22]='\0';
+    strcat(strcat(ch4,"Cln1/"),name11);
+    strcpy(cln1_file,ch4);
+   }
+   fo_cln1.open(cln1_file, ios::binary);   //nnnnn_mmm_xxx_cln1.pam
+   
    if (!fo_cln1) cout<<"can not open "<<cln1_file<<"\n";
-   fokadr_cln1.open(cln1_kadr,ios::binary); //xxx_mmm_nnnnn_cln1.bin
+   fokadr_cln1.open(cln1_kadr,ios::binary); //nnnnn_mmm_xxx_cln1.bin
    if (!fokadr_cln1) cout<<"can not open "<<cln1_kadr<<"\n";
-   fo_log_kadr=fopen(log_kadr,"wb");        //xxx_mmm_nnnnn_cln1_bin.log
-  } 
-  /***************New route***************************************/
-  /**route number: 7  ( 7  0 8) 8*********************************/
-  /**kadr number: 199 (200 0 1) 2*********************************/
-  /**route number: 7  ( 7  8 8) 8*********************************/
-  /**kadr number: 199 (200 0 1) 2*********************************/
+   fo_log_kadr=fopen(log_kadr,"wb");        //nnnnn_mmm_xxx_cln1_bin.log
+  }
+  /***************New route****************************************/
+  /**route number:   7 ( 7  0 8) 8*********************************/
+  /**cadre number: 199 (200 0 1) 2*********************************/
+  /**route number:   7 ( 7  8 8) 8*********************************/
+  /**cadre number: 199 (200 0 1) 2*********************************/
   if ( ((buffer1[3]!=buffer3[3])&&(buffer2[3]==0)&&(N1>3)&&(numP[1]!=numP[0]+1)&&(numP[1]<10))||
-       ((buffer1[3]!=buffer2[3])&&(buffer2[3]==buffer3[3])&&(buffer1[3]!=0)&&(numP[1]!=numP[0]+1)&&(numP[1]<10))||
+       ((buffer1[3]!=buffer2[3])&&(buffer2[3]==buffer3[3])&&(buffer1[3]!=0)&&(numP[1]!=numP[0]+1)&&(numP[1]<10)&&
+        (flag_pack[0]==0)&&(flag_pack[1]==0))||
        ((numP[2]==numP[1]+1)&&(numP[1]==0)&&(Num_kadrs[Rcount-1]>1000000)) )
   {Rcount+=1;
    packet_last[Rcount-2]=numP[0]; 
    packet_first[Rcount-1]=numP[1]; 
    NumRout[Rcount-1]=Rcount;
-   cout<<"New route="<<Rcount<<" "<<buffer1[3]<<" "<<buffer2[3]<<" "<<buffer3[3]<<"\n";
-   cout<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";
+   if (print_flag>0) 
+   {cout<<"New route="<<Rcount<<" "<<buffer1[3]<<" "<<buffer2[3]<<" "<<buffer3[3]<<"\n";
+    cout<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";}
    if (flag_tr[1]==0)
    {Num_kadrs[Rcount-2]-=2;
     Num_kadrs[Rcount-1]+=2;
@@ -458,8 +533,7 @@
     count_ident[Rcount-1]+=1024; count_ident[Rcount-2]-=1024;
     count_tot[Rcount-1]+=1024; count_tot[Rcount-2]-=1024;}
    // CREATE FILE NAMES
-   filename(file_route, miss_file, cln1_file, path,
-            NumRout[Rcount-1], name);			 
+   filename(file_route, miss_file, cln1_file, path, NumRout[Rcount-1], name);			 
    strncpy(cln1_kadr,cln1_file,strlen(cln1_file)-3);
    cln1_kadr[strlen(cln1_file)-3]='\0';
    strcat(cln1_kadr,"bin");
@@ -472,11 +546,11 @@
    fokadr_cln1.close();
    fclose(fo_log_kadr);
    //Open files for new route
-   foRout=fopen(file_route, "wb");          //xxx_mmm_nnnnn_cln1Er.log
-   foMiss=fopen(miss_file, "wb");           //xxx_mmm_nnnnn_mis_cln1.dat    
-   fo_cln1.open(cln1_file, ios::binary);    //xxx_mmm_nnnnn_cln1.pam
-   fokadr_cln1.open(cln1_kadr,ios::binary); //xxx_mmm_nnnnn_cln1.bin
-   fo_log_kadr=fopen(log_kadr,"wb");        //xxx_mmm_nnnnn_cln1_bin.log
+   foRout=fopen(file_route, "wb");          //nnnnn_mmm_xxx_cln1Er.log
+   foMiss=fopen(miss_file, "wb");           //nnnnn_mmm_xxx_mis_cln1.dat    
+   fo_cln1.open(cln1_file, ios::binary);    //nnnnn_mmm_xxx_cln1.pam
+   fokadr_cln1.open(cln1_kadr,ios::binary); //nnnnn_mmm_xxx_cln1.bin
+   fo_log_kadr=fopen(log_kadr,"wb");        //nnnnn_mmm_xxx_cln1_bin.log
    if ((CRC_bad_kadr[2]==numP[2])&&(numP[2]!=0)) {Num_bad_kadrs[Rcount-1]+=1; Num_bad_kadrs[Rcount-2]-=1;}
    if ((CRC_bad_kadr[1]==numP[1])&&(numP[1]!=0))
    {Num_bad_kadrs[Rcount-1]+=1; Num_bad_kadrs[Rcount-2]-=1;
@@ -487,8 +561,8 @@
 	       pow((float)(2*flag_CRC[1]),5)+
 	       pow((float)(2*flag_CRC_cor[1]),6));
     fprintf(foRout,"%1c%1c%1c%1c",(char)(buffer2[4]),(char)(buffer2[5]),(char)(buffer2[6]),(char)(flags));
-    fprintf(fout,"%i New route %i\n",CRC_bad_kadr[1],flags);
-    CRC_bad_kadr[1]=-1;} 
+    fo_info<<"New route: cadre="<<CRC_bad_kadr[1]<<" error flag="<<flags<<"\n";
+    CRC_bad_kadr[1]-=1;}
   }
   
   if (((numP[1]==0)&&(N1>1))||(N1==2))
@@ -507,13 +581,18 @@
   if (CRC_bad_kadr[1]==numP[1])
   {
    //Correct CRC0 error/////////////////////////////////////////////////
-   buffer2[7]=buffer1[1023];
-   if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
+   if (flag_CRC0[1]!=0)
+   {
+    buffer2[7]=buffer1[1023];
+    if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
+   }
    //Correct name error/////////////////////////////////////////////////
-   buffer2[0]=(int)(CODE1);
-   buffer2[1]=(int)(CODE2);
-   buffer2[2]=(int)(CODE3);
-   if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
+   if (flag_name[1]!=0)
+   {buffer2[0]=(int)(CODE1);
+    buffer2[1]=(int)(CODE2);
+    buffer2[2]=(int)(CODE3);
+    if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
+   }
   }
   ////////////////////////////////////////////////////////////////////
   /*******************************************************************/
@@ -527,14 +606,13 @@
 			         pow((float)(2*flag_CRC_cor[1]),6));
    fprintf(foRout,"%1c%1c%1c%1c",(char)(buffer2[4]),
            (char)(buffer2[5]),(char)(buffer2[6]),(char)(flags));
-   fprintf(fout,"%i flags=%i\n",CRC_bad_kadr[1],flags);}
-   
+   fo_info<<" error flag="<<flags<<" in cadre "<<numP[1]<<"\n";}
+ 
   /***************************************************************/
   //*Packet analysis**********************************************
-  if (N_er!=0) fprintf(fout,"Name error %i in kadr %i\n",numP[2],N_er);
-  if (M_er!=0) fprintf(fout,"Root number error %i in kadr %i\n",numP[1],M_er);
-  if (P_er!=0) fprintf(fout,"Number packet error %i in kadr %i\n",numP[1], P_er);
-  if (CRC0_er!=0) fprintf(fout,"CRC0 error %i in kadr %i\n",numP[1], CRC0_er);
+  if (N_er!=0) fo_info<<"Name error "<<N_er<<" in cadre "<<numP[2]<<"\n";
+  if (P_er!=0) fo_info<<"Number packet error "<<numP[1]<<" in cadre "<<P_er<<"\n";
+  if (CRC0_er!=0) fo_info<<"CRC0 error "<<numP[1]<<" in cadre "<<CRC0_er<<"\n";
 
   /***************PAMELA data extracting for each route***********/
 //  if (N1!=1&&(flag_CRC[1]==0||(flag_rout[1]==0&&flag_pack[1]==0)||
@@ -552,7 +630,7 @@
    }
    else {for(i=8;i<1023;i++) fo_cln1<<(unsigned char)buffer2[i];}
    for (i=0;i<1024;i++) {fokadr_cln1<<(unsigned char)buffer2[i];}
-   //Write in file mmm_counter_nnnnn_cln1_bin.log
+   //Write in file nnnnn_mmm_xxx_cln1_bin.log
    tmp=numP[1]; 
    for (k=0;k<4;k++) numPack[k]=0;
    for (k=0;k<3;k++)
@@ -581,11 +659,10 @@
  }
 /***************************************************************/
 /***************************************************************/
- if (flag_CRC[2]==0||(flag_rout[2]==0&&flag_pack[2]==0)||
-     (flag_CRC[2]!=0&&flag_CRC_cor[2]==1)) 
+ if (flag_CRC[2]==0||(flag_rout[2]==0&&flag_pack[2]==0)||(flag_CRC[2]!=0&&flag_CRC_cor[2]==1)) 
  {end=0;
   for (i=0;i<1024;i++) fokadr_cln1<<(unsigned char)buffer3[i];
-  //Write in file mmm_counter_nnnnn_cln1_bin.log
+  //Write in file nnnnn_mmm_xxx_cln1_bin.log
   tmp=numP[2];
   for (k=0;k<4;k++) numPack[k]=0;
   for (k=0;k<3;k++)
@@ -631,26 +708,50 @@
  fprintf(fo,"file size= %i,%i kByte \n", icount_int, icount_frac);
  fprintf(fo,"number identified cadres= %i\n", N1);
  fprintf(fo,"number of routes= %i\n", Rcount); 
- 
+
+ for(i=0;i<5;i++) {temp[i]=name[i]; temp[5]='\0';}
+ sscanf(temp,"%5d", &cirquit_no);
+ for(i=0;i<3;i++) {temp[i]=name[i+5]; temp[3]='\0';}
+ sscanf(temp,"%3i", &route);
+  
  for (i=0;i<Rcount;i++) 
  {
+  for(j=0;j<22;j++) temp[j]=cln1_file[j+strlen(cln1_file)-22]; 
+  temp[22]='\0';
+  oss.str("");
+  oss << "INSERT INTO RAW_FILES (cirquit_n, route_n, counter, NumCadr, CRCerr, cln1_name, cln1_size)"
+  << " VALUES ("<< cirquit_no << ","<< route << ","<< i+1 <<","<< Num_kadrs[i]<< "," << Num_bad_kadrs[i] << ",'"
+  << temp << "'," << count_ident[i]-Num_kadrs[i]*9 << ")"; 
+  if (print_flag>0) cout<<oss.str().c_str()<<" table RAW_FILES\n";
+  TSQLResult *res=sqlServer->Query(oss.str().c_str());
+ 
   Num_good_kadrs[i]=Num_kadrs[i]-Num_bad_kadrs[i];
+  if (packet_first[i]>1) 
+  {first_missed[i]=(float)(1.)-float(packet_first[i]-1.)/(float)(packet_last[i]);}
+  else if (packet_first[i]==0)
+  {first_missed[i]=1.;}
+  
+
   quality[i]=(float)(1.)-(float)(Num_bad_kadrs[i])/(float)(Num_kadrs[i]);
   fprintf(fo,"\nRoute number=%i: number of normal cadres= %i: number of errors=%i; \n",NumRout[i],Num_good_kadrs[i],Num_bad_kadrs[i]); 
-  fprintf(fo," size of identified cadres =%lli bytes;",count_ident[i]);
+  fprintf(fo,"size of identified cadres =%lli bytes;",count_ident[i]);
   fprintf(fo," total size of route=%lli bytes;\n",count_tot[i]);
-  fprintf(fo," first cadre number= %lli:", packet_first[i]);
+  fprintf(fo,"first cadre number= %lli:", packet_first[i]);
   fprintf(fo," last cadre number= %lli:", packet_last[i]);
-  fprintf(fo," Number of missed cadres= %i:\n", count_missed[i]);  
-  fprintf(fo,"number of identified cadres=%lli: quality index=%f\n", Num_kadrs[i], quality[i]);  
+  fprintf(fo," Number of missed cadres= %lli: Number of double transmitted cadres= %lli\n",
+             count_missed[i], count_sec_missed[i]);  
+  fprintf(fo,"number of identified cadres=%lli: quality index=%f: completness=%f \n", Num_kadrs[i], quality[i], 
+          first_missed[i]);  
  }
 
- //***route's number log file**************************************
+ //***route number log file**************************************
  for (i=0;i<Rcount;i++) {fprintf(fo_route,"%1c",(char)(NumRout[i]));}
+
 	 
 //***************************************************************
   fin.close();
-  fclose(fout);
+  fclose(fout);  
+  fo_info.close();
   fclose(fo);
   fclose(foRout);
   fclose(foMiss);