/[PAMELA software]/rawreader/src/Read.cpp
ViewVC logotype

Diff of /rawreader/src/Read.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.3 by cafagna, Sun Dec 10 12:52:45 2006 UTC revision 1.5 by cafagna, Sun Jan 28 10:55:47 2007 UTC
# Line 6  void FindRoute(int[], int*); Line 6  void FindRoute(int[], int*);
6  short int CRC(int []);  short int CRC(int []);
7  int restrict(int [], long long int[], int);  int restrict(int [], long long int[], int);
8    
9  extern char database[80];  extern char database[120];
10  extern char MYSQL[80];  extern char MYSQL[120];
11  extern int print_flag;  extern int print_flag;
12    
13    int route_good1=0, route_good2=0, route_good3=0;
14    int pack_good1=0, pack_good2=0, pack_good3=0;
15    int number_good1=0, number_good2=0, number_good3=0;
16    
17  void read(char fni[], char fno[], char INF[],  void read(char fni[], char fno[], char INF[],
18            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)
19  {  {
# Line 20  void read(char fni[], char fno[], char I Line 24  void read(char fni[], char fno[], char I
24   char buf[1024];   char buf[1024];
25   //  char stop1;   //  char stop1;
26   char ch1,ch2,ch3,ch, ch4[80], ch5[5];   char ch1,ch2,ch3,ch, ch4[80], ch5[5];
27   char file_route[80];  /*protocol binary filesm nnnnn_mmm_xxx_cln1Er.log*/   char file_route[100];  /*protocol binary filesm nnnnn_mmm_xxx_cln1Er.log*/
28   char miss_file[80];   /*protocol files nnnnn_mmm_xxx_mis_cln1.dat*/   char miss_file[100];   /*protocol files nnnnn_mmm_xxx_mis_cln1.dat*/
29   char cln1_file[80];   /*route binary data files nnnnn_mmm_xxx_cln1.pam*/   char cln1_file[100];   /*route binary data files nnnnn_mmm_xxx_cln1.pam*/
30   char cln1_kadr[80];   /*route binary data files nnnnn_mmm_xxx_cln1.bin*/   char cln1_kadr[100];   /*route binary data files nnnnn_mmm_xxx_cln1.bin*/
31   char log_kadr[80];    /*log binary data files nnnnn_mmm_xxx_cln1.log*/   char log_kadr[100];    /*log binary data files nnnnn_mmm_xxx_cln1.log*/
32   char log_route[80];   /*log binary file for number of routes*/   char log_route[100];   /*log binary file for number of routes*/
33   char Short_Buf[8], Name_Buf[4];   char Short_Buf[8], Name_Buf[4];
34   char temp[10], temp1[80], connection[80], command[80], name11[80];   char temp[10], temp1[80], connection[120], command[80], name11[80];
35   short int CalcCheckSum;   short int CalcCheckSum;
36    
37   int LenthBuf=1015, flag_help=0;   int LenthBuf=1015, flag_help=0;
# Line 49  void read(char fni[], char fno[], char I Line 53  void read(char fni[], char fno[], char I
53   long long int icount_total=0;   /* Total number of bytes*/   long long int icount_total=0;   /* Total number of bytes*/
54   long long int packet_first[100];/* First cadre number for each route*/   long long int packet_first[100];/* First cadre number for each route*/
55   long long int packet_last[100]; /* Second cadre number for each route*/   long long int packet_last[100]; /* Second cadre number for each route*/
56   long int count_missed[100];     /* Number of missed cadres in each route*/   long long int count_missed[100];/* Number of missed cadres in each route*/
57   long int count_sec_missed[100];   long long int count_sec_missed[100];
58   unsigned long long int count_ident[100]; /* Number of good bytes in each route*/   unsigned long long int count_ident[100]; /* Number of good bytes in each route*/
59   unsigned long long int count_tot[100];   /* Total number of bytes in each route*/   unsigned long long int count_tot[100];   /* Total number of bytes in each route*/
60   int NumRout[1000];              /* Number of route*/   int NumRout[1000];              /* Number of route*/
# Line 63  void read(char fni[], char fno[], char I Line 67  void read(char fni[], char fno[], char I
67   int adr2=0;                     /* Last byte of cadre in cln2.bin*/   int adr2=0;                     /* Last byte of cadre in cln2.bin*/
68   int tmp=0;   int tmp=0;
69   int numPack[4];                 /* Number of cadre in given route*/   int numPack[4];                 /* Number of cadre in given route*/
70   int good_pack=0;  
  long long int good_number=0;  
  int good_route=0;  
71   int adress[4];   int adress[4];
72   int R_er[100];           /* Route errors number in given route*/   int R_er[100];           /* Route errors number in given route*/
73   int R_er_interf1[100];   /* Interface Route error at the begining of the route*/   int R_er_interf1[100];   /* Interface Route error at the begining of the route*/
# Line 84  void read(char fni[], char fno[], char I Line 86  void read(char fni[], char fno[], char I
86   int flag_CRC_cor[3];     /* Flag for CRC error correction*/   int flag_CRC_cor[3];     /* Flag for CRC error correction*/
87   int flag_tr[3];          /* Truncate flag*/   int flag_tr[3];          /* Truncate flag*/
88   short int flags;         /* Number for all errors*/   short int flags;         /* Number for all errors*/
89     int num_rows=0;
90     int control[1];
91     control[0]=0;
92     control[1]=0;
93   for (j=0; j<4; j++) {adress[j]=0; numPack[j]=0;}   for (j=0; j<4; j++) {adress[j]=0; numPack[j]=0;}
94   for (j=0; j<1024; j++)   for (j=0; j<1024; j++)
95   {buffer1[j]=0; buffer2[j]=0; buffer3[j]=0;}   {buffer1[j]=0; buffer2[j]=0; buffer3[j]=0;}
# Line 114  void read(char fni[], char fno[], char I Line 119  void read(char fni[], char fno[], char I
119    
120   strcpy(connection,MYSQL);   strcpy(connection,MYSQL);
121   strcat(connection,database);   strcat(connection,database);
122   TSQLServer *sqlServer = TSQLServer::Connect(connection, "pamelaprod", "rawprod");   TSQLServer *sqlServer = TSQLServer::Connect(connection,"pamelaprod","rawprod");
123     if (sqlServer==0) {cerr<<"can not connect with MYSQL sever\n"; exit(1);}
124   // open input file "nnnnnmmm.pam"   // open input file "nnnnnmmm.pam"
125   fin.open(fni, ifstream::binary);   fin.open(fni, ifstream::binary);
126   if(fin.good()==0) {cerr<<"can not open input file"; exit(0);}   if(fin.good()==0) {cerr<<"can not open input file"; exit(1);}
127   //open output file "nnnnnmmm_work.dat"   //open output file "nnnnnmmm_work.dat"
128   fout=fopen(fno, "wb");   fout=fopen(fno, "wb");
129   if(fout==0) {cerr<<"can not open work file"; exit(0);}   if(fout==0) {cerr<<"can not open work file"; exit(1);}
130    
131   strcpy(ch5,"info");   strcpy(ch5,"info");
132   strcpy(ch4,fno);   strcpy(ch4,fno);
133   for(i=0;i<4;i++){ch4[strlen(ch4)-8+i]=ch5[i];}   for(i=0;i<4;i++){ch4[strlen(ch4)-8+i]=ch5[i];}
134   // open input file "nnnnnmmm_info.dat"   // open input file "nnnnnmmm_info.dat"
135   fo_info.open(ch4, ifstream::binary);   fo_info.open(ch4, ifstream::binary);
136   if(fo_info.good()==0) {cerr<<"can not open info file"; exit(0);}   if(fo_info.good()==0) {cerr<<"can not open info file"; exit(1);}
137    
138   //open output file "nnnnnmmm_protocol.txt"   //open output file "nnnnnmmm_protocol.txt"
139   fo=fopen(protocol, "wb");   fo=fopen(protocol, "wb");
140   if(fo==0) {cerr<<"can not open protocol file"; exit(0);}   if(fo==0) {cerr<<"can not open protocol file"; exit(1);}
141   //open output file "nnnnnmmm_route.log"   //open output file "nnnnnmmm_route.log"
142   strncpy(log_route,INF,strlen(INF)-9);   strncpy(log_route,INF,strlen(INF)-9);
143   log_route[strlen(INF)-9]='\0';   log_route[strlen(INF)-9]='\0';
144   strcat(log_route,"_route.log");   strcat(log_route,"_route.log");
145   fo_route=fopen(log_route, "wb");   fo_route=fopen(log_route, "wb");
146   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);}
147   /////////////////////////////////////////////////////////////////   /////////////////////////////////////////////////////////////////
148   // cin>>stop1;   // cin>>stop1;
149   k=0;   k=0;
150   while(fin.eof()==0)   while(fin.eof()==0)
151   {   {
152    counter=0;          //pounter for position of cadre truncation    counter=0;          //pounter for position of cadre truncation
# Line 163  void read(char fni[], char fno[], char I Line 169  void read(char fni[], char fno[], char I
169    count_ident[Rcount-1]+=8; count_tot[Rcount-1]+=8;    count_ident[Rcount-1]+=8; count_tot[Rcount-1]+=8;
170    
171    while((ch1!=(char)(CODE1)&&ch2!=(char)(CODE2)&&    while((ch1!=(char)(CODE1)&&ch2!=(char)(CODE2)&&
172    (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]))))
173               ||               ||
174    (ch2!=(char)(CODE2)&&ch3!=(char)(CODE3)&&(buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])||(numP[2]+1)!=tmpP))    (ch2!=(char)(CODE2)&&ch3!=(char)(CODE3)&&(buffer3[1023]!=(int)(unsigned char)(Short_Buf[4])||(numP[2]+1)!=tmpP))
175               ||               ||
176    (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])||
177       (buffer3[5]!=(int)(unsigned char)(Short_Buf[2])&& buffer3[6]!=(int)(unsigned char)(Short_Buf[3]))))
178               ||               ||
179    (ch1!=(char)(CODE1)&&    (ch1!=(char)(CODE1)&&
180     (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])&&
181        buffer3[5]!=(int)(unsigned char)(Short_Buf[2])&& buffer3[6]!=(int)(unsigned char)(Short_Buf[3])))
182               ||               ||
183    (ch2!=(char)(CODE2)&&    (ch2!=(char)(CODE2)&&
184     (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])&&
185        buffer3[5]!=(int)(unsigned char)(Short_Buf[2])&& buffer3[6]!=(int)(unsigned char)(Short_Buf[3])))
186               ||               ||
187    (ch3!=(char)(CODE3)&&    (ch3!=(char)(CODE3)&&
188     (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])&&
189        buffer3[5]!=(int)(unsigned char)(Short_Buf[2])&& buffer3[6]!=(int)(unsigned char)(Short_Buf[3])))  
190               ||               ||
191    (ch1!=(char)(CODE1)&&ch2!=(char)(CODE2)&&ch3!=(char)(CODE3)))    (ch1!=(char)(CODE1)&&ch2!=(char)(CODE2)&&ch3!=(char)(CODE3)))
192    {fin.get(ch);    {fin.get(ch);
# Line 196  void read(char fni[], char fno[], char I Line 206  void read(char fni[], char fno[], char I
206     fclose(fo);     fclose(fo);
207     remove(protocol);     remove(protocol);
208     remove(INF);     remove(INF);
209     exit(0);     exit(1);
210    }    }
   
211    /*****************************************************************/    /*****************************************************************/
212    /*********************data reading********************************/    /*********************data reading********************************/
213    fin.read(buf,LenthBuf);    fin.read(buf,LenthBuf);
# Line 257  void read(char fni[], char fno[], char I Line 266  void read(char fni[], char fno[], char I
266                     for (i=0; i<2; i++)                     for (i=0; i<2; i++)
267                     {CRC_1[i]=CRC_1[2];}                     {CRC_1[i]=CRC_1[2];}
268             packet_first[Rcount-1]=numP[2];             packet_first[Rcount-1]=numP[2];
269               control[0]=0; control[1]=0;
270                     break;                     break;
271            }            }
272     default:{     default:{
# Line 294  void read(char fni[], char fno[], char I Line 304  void read(char fni[], char fno[], char I
304              CRC_0[2]=(int)(unsigned char)(Short_Buf[4]);              CRC_0[2]=(int)(unsigned char)(Short_Buf[4]);
305                      for (i=0; i<2; i++) {CRC_1[i]=CRC_1[i+1];}                      for (i=0; i<2; i++) {CRC_1[i]=CRC_1[i+1];}
306              CRC_1[2]=(int)(unsigned char)(ch);                  CRC_1[2]=(int)(unsigned char)(ch);    
307                control[0]=control[1];
308                control[1]=0;
309             }             }
310    }      }  
311    if (print_flag>0 && fmod(N1,10000)==0.) cout<<N1<<"\n";            //  if (print_flag>0 && fmod(N1,10000)==0) cout<<N1<<"\n";      
312      if (print_flag>0 && !(N1%10000)) cout<<N1<<"\n";
313    
314  /******************CRC ERROR ANALYSIS***************************/  /******************CRC ERROR ANALYSIS***************************/
315    CalcCheckSum=CRC(buffer3);    CalcCheckSum=CRC(buffer3);
# Line 305  void read(char fni[], char fno[], char I Line 318  void read(char fni[], char fno[], char I
318  //write in file nnnnnmmm_info.dat      //write in file nnnnnmmm_info.dat    
319     fo_info<<"real CRC = "<<hex<<buffer3[1023]<<" (calculated CRC ="<<CalcCheckSum<<") in cadre "<<dec<<numP[2]<<"\n";     fo_info<<"real CRC = "<<hex<<buffer3[1023]<<" (calculated CRC ="<<CalcCheckSum<<") in cadre "<<dec<<numP[2]<<"\n";
320  //write in file nnnnnmmm_work.dat        //write in file nnnnnmmm_work.dat      
321     fprintf(fout,"%i %x (%x)", numP[2], CalcCheckSum, buffer3[1023]);     fprintf(fout,"%10i %4i %4i\n", numP[2], CalcCheckSum, buffer3[1023]);
322    
323     flag_CRC[2]=1;     flag_CRC[2]=1;
324     Num_bad_kadrs[Rcount-1]+=1;     Num_bad_kadrs[Rcount-1]+=1;
# Line 314  void read(char fni[], char fno[], char I Line 327  void read(char fni[], char fno[], char I
327     if (counter!=0)     if (counter!=0)
328     {fo_info<<" truncate byte="<<counter<<" in cadre "<<numP[2]<<"\n";     {fo_info<<" truncate byte="<<counter<<" in cadre "<<numP[2]<<"\n";
329      if (print_flag>0)      if (print_flag>0)
330      {cout<<" truncate counter="<<counter<<" route="<<buffer2[3]<<" "<<buffer3[3]<<" kadr="<<numP[2]<<" "<<numP[1]<<" "<<numP[0]<<"\n";}      {cout<<" truncate counter="<<counter<<" route="<<buffer2[3]<<" "<<buffer3[3]<<" cadre="<<numP[2]<<" "<<numP[1]<<" "<<numP[0]<<"\n";}
331      fin.seekg(-(1024-counter),ios::cur);      fin.seekg(-(1024-counter),ios::cur);
332      icount-=1024;      icount-=1024;
333      count_ident[Rcount-1]-=1024;      count_ident[Rcount-1]-=1024;
# Line 338  void read(char fni[], char fno[], char I Line 351  void read(char fni[], char fno[], char I
351    {flag_CRC0[1]=1; CRC0_er=N1-1;}    {flag_CRC0[1]=1; CRC0_er=N1-1;}
352    /****************************************************************/    /****************************************************************/
353    
354    /******************KADR NUMBER ANALYSIS**************************/    /******************CADRE NUMBER ANALYSIS*************************/
355    if ((numP[1]!=(numP[0]+1))&&(numP[1]!=(numP[2]-1))&&(N1!=1)&&(CRC_bad_kadr[1]==numP[1])&&(numP[1]!=0))    if ((numP[1]!=(numP[0]+1))&&(numP[1]!=(numP[2]-1))&&(N1!=1)&&(CRC_bad_kadr[1]==numP[1])&&(numP[1]!=0))
356    {P_er=N1-1; flag_pack[1]=1;    {P_er=N1-1; flag_pack[1]=1;
357    //Correct error packet number//    //Correct error packet number//
# Line 348  void read(char fni[], char fno[], char I Line 361  void read(char fni[], char fno[], char I
361     {if (tmpP/(long int)(pow(256.,2-i))>=1)     {if (tmpP/(long int)(pow(256.,2-i))>=1)
362      {Hex[2-i]=(int)(tmpP/(long int)(pow(256.,2-i)));      {Hex[2-i]=(int)(tmpP/(long int)(pow(256.,2-i)));
363       tmpP=tmpP-(long int)pow(256.,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)));}}
364     for (i=4;i<7;i++) (buffer2[i]=Hex[i-4]);  //   for (i=4;i<7;i++) (buffer2[i]=Hex[i-4]);
365     if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;     if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
366    
367     tmpP=numP[2]-1;     tmpP=numP[2]-1;
# Line 357  void read(char fni[], char fno[], char I Line 370  void read(char fni[], char fno[], char I
370     {if (tmpP/(long int)(pow(256.,2-i))>=1)     {if (tmpP/(long int)(pow(256.,2-i))>=1)
371      {Hex[2-i]=(int)(tmpP/(long int)(pow(256.,2-i)));      {Hex[2-i]=(int)(tmpP/(long int)(pow(256.,2-i)));
372           tmpP=tmpP-(long int)pow(256.,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)));}}
373     for (i=4;i<7;i++) (buffer2[i]=Hex[i-4]);  //   for (i=4;i<7;i++) (buffer2[i]=Hex[i-4]);
374     if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;     if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
375    }    }
376    /////////////////////////////////////////////////////////////////////    /////////////////////////////////////////////////////////////////////
# Line 368  void read(char fni[], char fno[], char I Line 381  void read(char fni[], char fno[], char I
381    /**cadre number:  21 (22 17 18) 19 ******************************/    /**cadre number:  21 (22 17 18) 19 ******************************/
382    /**missed cadres: 21  20 19 18  17 ******************************/      /**missed cadres: 21  20 19 18  17 ******************************/  
383    if ((numP[1]<numP[0])&&(numP[1]==numP[2]-1)&&(buffer1[3]==buffer3[3])&&(flag_tr[1]==0)&&    if ((numP[1]<numP[0])&&(numP[1]==numP[2]-1)&&(buffer1[3]==buffer3[3])&&(flag_tr[1]==0)&&
384            (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))
385    {Mis_group[Rcount-1]+=1;    {Mis_group[Rcount-1]+=1;
386     count_sec_missed[Rcount-1]+=numP[0]-1-numP[1];     count_sec_missed[Rcount-1]+=numP[0]-1-numP[1];
387     fprintf(foMiss,"secondary transmitted cadres may be %lli - %lli number of cadres=%lli\n",numP[1],numP[0],numP[0]-1-numP[1]);     fprintf(foMiss,"secondary transmitted cadres may be %lli - %lli number of cadres=%lli\n",
388               numP[1],numP[0],numP[0]-1-numP[1]);
389    }      }  
     
390    /**route number: 10 (10 19 10) 10 ******************************/    /**route number: 10 (10 19 10) 10 ******************************/
391    /**cadre number: 11 (12 15 20) 21 ******************************/    /**cadre number: 11 (12 15 20) 21 ******************************/
392    /**missed cadres:   16 17 18 19   ******************************/      /**missed cadres:    16 17 18 19   ******************************/  
393    if ((numP[1]>numP[0]+1)&&(numP[1]<numP[2]-1)&&(buffer1[3]==buffer3[3])&&(flag_tr[1]==0)&&    if ((numP[1]>numP[0]+1)&&(numP[1]<numP[2]-1)&&(buffer1[3]==buffer3[3])&&(flag_tr[1]==0)&&
394            (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0))        (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0))
395    {Mis_group[Rcount-1]+=1;    {Mis_group[Rcount-1]+=1;
396     fprintf(foMiss,"1) %9lli%9lli\n",numP[1]+1,numP[2]-1);     fprintf(foMiss,"1) %9lli%9lli\n",numP[1]+1,numP[2]-1);
397       if (print_flag>0) cout<<"1) "<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";      
398     count_missed[Rcount-1]+=numP[2]-numP[1]-1;}     count_missed[Rcount-1]+=numP[2]-numP[1]-1;}
399    /**route number: (10 10 10) 10 10 ******************************/    /**route number: (10 10 10) 10 10 ******************************/
400    /**cadre number: (11 12 15) 20 21 ******************************/    /**cadre number: (11 29 15) 14 15 ******************************/
   /**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,"2) %9lli%9lli\n",numP[0]+2,numP[2]-1);  
    count_missed[Rcount-1]+=numP[2]-numP[0]-1-1;}  
   /**route number: (10 10 10) 10 10 ******************************/  
   /**cadre number:  (11 29 15) 14 15 ******************************/  
401    /**number of missed cadres:2       ******************************/      /**number of missed cadres:2       ******************************/  
402    else if (((numP[1]<numP[0]+1)||(numP[1]>numP[2]-1))&&(numP[2]-numP[0]>2)&&    else if (((numP[1]<numP[0]+1)||(numP[1]>numP[2]-1))&&(numP[2]-numP[0]>2)&&
403             (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&(flag_tr[1]==0)&&             (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&(flag_tr[1]==0)&&
404             (flag_pack[2]==0)&&(flag_pack[1]==1)&&(flag_pack[0]==0))             (flag_pack[2]==0)&&(flag_pack[1]==1)&&(flag_pack[0]==0))
405    {Mis_group[Rcount-1]+=1;    {Mis_group[Rcount-1]+=1;
406     fprintf(foMiss,"3) %9lli%9lli -1\n",numP[0]+1,numP[2]-1);     fprintf(foMiss,"2) %9lli%9lli\n",numP[0]+1,numP[2]-1);
407       if (print_flag>0) cout<<"2) "<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";  
408     count_missed[Rcount-1]+=numP[2]-numP[0]-1-1;}     count_missed[Rcount-1]+=numP[2]-numP[0]-1-1;}
409      else if (((numP[1]<numP[0]+1)||(numP[1]>numP[2]-1))&&((numP[2]-numP[0])==2)&&
410               (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&(flag_tr[1]==0))
411      {control[1]=1;}  
412      /**route number: (10 10 10) 10 10 ******************************/
413      /**cadre number: (11 15 16) 17 18 ******************************/
414      /**missed kadrs:  12 13 14        ******************************/  
415      else if ((numP[1]>numP[0]+1)&&(numP[1]<=numP[2]-1)&&(buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&
416               (flag_tr[1]==0)&&(flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0)&&(control[0]==0))
417      {Mis_group[Rcount-1]+=1;
418       fprintf(foMiss,"3) %9lli%9lli\n",numP[0]+1,numP[1]-1);
419       if (print_flag>0) cout<<"3) "<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";      
420       count_missed[Rcount-1]+=numP[1]-numP[0]-1;}
421    /**route number:(10 10 )10 10 10 10 ******************************/    /**route number:(10 10 )10 10 10 10 ******************************/
422    /**cadre number:( 1 11 )12 13 14 15 ******************************/    /**cadre number:( 1 11 )12 13 14 15 ******************************/
423    /**number of missed cadres: 10       ******************************/      /**number of missed cadres: 10      ******************************/  
424    else if ((N1>1)&&(abs(numP[2]-good_pack-2)>abs(Num_kadrs[Rcount-1]-good_number))&&    else if ((N1>2)&&(abs(numP[1]-pack_good3)>abs(Num_kadrs[Rcount-1]-number_good3-1))&&
425             (buffer2[3]==good_route)&&(flag_tr[2]==0)&&(flag_pack[2]==0))             (buffer2[3]==route_good3)&&(flag_tr[2]==0)&&(flag_pack[2]==0)&&(numP[2]>numP[0]+2)&&(control[0]==0))
426    {Mis_group[Rcount-1]+=1;      {Mis_group[Rcount-1]+=1;
427     count_missed[Rcount-1]+=abs(numP[2]-good_pack)-1;     if (pack_good3<numP[2]) count_missed[Rcount-1]+=abs(numP[2]-pack_good3)-1;
428     fprintf(foMiss,"4) %9i%9lli number of missed cadres=%i\n", good_pack, numP[2], abs(numP[2]-good_pack)-1);     fprintf(foMiss,"4) %9lli%9lli number of missed cadres=%i\n", pack_good3, numP[1], abs(numP[1]-pack_good3-1));
429    }     if (print_flag>0) cout<<"P1="<<numP[0]<<" P2="<<numP[1]<<" P3="<<numP[2]<<"\n";
430    //   cout<<"pack_good1="<<pack_good1<<" pack_good2="<<pack_good2<<" pack_good3="<<pack_good3<<"\n";
431    //   cout<<"number_good1="<<number_good1<<" number_good2="<<number_good2<<" number_good3="<<number_good3<<"\n";
432      }  
433    //truncate cadres    //truncate cadres
434    /**route number: 10 (10 10 10) 10 ******************************/    /**route number: 10 (10 10 10) 10  ******************************/
435    /**cadre number:  11 (12 13 14) 15 ******************************/    /**cadre number:  11 (12 13 14) 15 ******************************/
436    /**missed truncate cadre:13        ******************************/      /**missed truncate cadre:13        ******************************/  
437    else if ((numP[0]<numP[2]-1)&&(buffer1[3]==buffer3[3])&&    else if ((numP[0]<numP[2]-1)&&(buffer1[3]==buffer3[3])&&
# Line 421  void read(char fni[], char fno[], char I Line 442  void read(char fni[], char fno[], char I
442     count_missed[Rcount-1]+=numP[2]-numP[0]-1;}     count_missed[Rcount-1]+=numP[2]-numP[0]-1;}
443    /***************************************************************/    /***************************************************************/
444    
445    if (flag_CRC[2]!=1) {good_pack=numP[2]; good_route=buffer3[3]; good_number=Num_kadrs[Rcount-1];}    if (flag_CRC[2]!=1)
446      {
447       pack_good1=pack_good2;
448       pack_good2=pack_good3;
449       route_good1=route_good2;
450       route_good2=route_good3;
451       number_good1=number_good2;
452       number_good2=number_good3;    
453    
454       pack_good3=numP[2];      
455       route_good3=buffer3[3];
456       number_good3=Num_kadrs[Rcount-1];  
457      }
458    /******************ROUTE NUMBER ANALYSIS**************************/    /******************ROUTE NUMBER ANALYSIS**************************/
459    if ((buffer2[3]!=buffer1[3])&&(buffer2[3]!=buffer3[3])&&(buffer1[3]!=buffer3[3])&&    if ((buffer2[3]!=buffer1[3])&&(buffer2[3]!=buffer3[3])&&(buffer1[3]!=buffer3[3])&&
460        (flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))        (flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))
# Line 437  void read(char fni[], char fno[], char I Line 469  void read(char fni[], char fno[], char I
469     packet_last[Rcount-1]=numP[2];     packet_last[Rcount-1]=numP[2];
470     NumRout[Rcount-1]=1;     NumRout[Rcount-1]=1;
471     // CREATE FILE NAMES     // CREATE FILE NAMES
472     filename(file_route, miss_file, cln1_file, path, NumRout[Rcount-1], name);                         filename(file_route, miss_file, cln1_file, path, NumRout[Rcount-1], name);
473     strncpy(cln1_kadr,cln1_file,strlen(cln1_file)-3);     strncpy(cln1_kadr,cln1_file,strlen(cln1_file)-3);
474     cln1_kadr[strlen(cln1_file)-3]='\0';     cln1_kadr[strlen(cln1_file)-3]='\0';
475     strcat(cln1_kadr,"bin");     strcat(cln1_kadr,"bin");
# Line 473  void read(char fni[], char fno[], char I Line 505  void read(char fni[], char fno[], char I
505     fokadr_cln1.open(cln1_kadr,ios::binary); //nnnnn_mmm_xxx_cln1.bin     fokadr_cln1.open(cln1_kadr,ios::binary); //nnnnn_mmm_xxx_cln1.bin
506     if (!fokadr_cln1) cout<<"can not open "<<cln1_kadr<<"\n";     if (!fokadr_cln1) cout<<"can not open "<<cln1_kadr<<"\n";
507     fo_log_kadr=fopen(log_kadr,"wb");        //nnnnn_mmm_xxx_cln1_bin.log     fo_log_kadr=fopen(log_kadr,"wb");        //nnnnn_mmm_xxx_cln1_bin.log
508    }    }
509    /***************New route****************************************/    /***************New route****************************************/
510    /**route number:   7 ( 7  0 8) 8*********************************/    /**route number:   7 ( 7  0 8) 8*********************************/
511    /**cadre number: 199 (200 0 1) 2*********************************/    /**cadre number: 199 (200 0 1) 2*********************************/
512    /**route number:   7 ( 7  8 8) 8*********************************/    /**route number:   7 ( 7  8 8) 8*********************************/
513    /**cadre number: 199 (200 0 1) 2*********************************/    /**cadre number: 199 (200 0 1) 2*********************************/
514    if ( ((buffer1[3]!=buffer3[3])&&(buffer2[3]==0)&&(N1>3)&&(numP[1]!=numP[0]+1)&&(numP[1]<10))||    if ( ((buffer1[3]!=buffer3[3])&&(buffer2[3]==0)&&(N1>3)&&(numP[1]!=numP[0]+1)&&(numP[1]<10))||
515         ((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)&&
516            (flag_pack[0]==0)&&(flag_pack[1]==0))||
517         ((numP[2]==numP[1]+1)&&(numP[1]==0)&&(Num_kadrs[Rcount-1]>1000000)) )         ((numP[2]==numP[1]+1)&&(numP[1]==0)&&(Num_kadrs[Rcount-1]>1000000)) )
518    {Rcount+=1;    {Rcount+=1;
519     packet_last[Rcount-2]=numP[0];     packet_last[Rcount-2]=numP[0];
# Line 573  void read(char fni[], char fno[], char I Line 606  void read(char fni[], char fno[], char I
606                                   pow((float)(2*flag_CRC_cor[1]),6));                                   pow((float)(2*flag_CRC_cor[1]),6));
607     fprintf(foRout,"%1c%1c%1c%1c",(char)(buffer2[4]),     fprintf(foRout,"%1c%1c%1c%1c",(char)(buffer2[4]),
608             (char)(buffer2[5]),(char)(buffer2[6]),(char)(flags));             (char)(buffer2[5]),(char)(buffer2[6]),(char)(flags));
609     fo_info<<"CRC error in cadre "<<CRC_bad_kadr[1]<<" error flag="<<flags<<"\n";}     fo_info<<" error flag="<<flags<<" in cadre "<<numP[1]<<"\n";}
610    
611    /***************************************************************/    /***************************************************************/
612    //*Packet analysis**********************************************    //*Packet analysis**********************************************
# Line 626  void read(char fni[], char fno[], char I Line 659  void read(char fni[], char fno[], char I
659   }   }
660  /***************************************************************/  /***************************************************************/
661  /***************************************************************/  /***************************************************************/
662   if (flag_CRC[2]==0||(flag_rout[2]==0&&flag_pack[2]==0)||   if (flag_CRC[2]==0||(flag_rout[2]==0&&flag_pack[2]==0)||(flag_CRC[2]!=0&&flag_CRC_cor[2]==1))
      (flag_CRC[2]!=0&&flag_CRC_cor[2]==1))  
663   {end=0;   {end=0;
664    for (i=0;i<1024;i++) fokadr_cln1<<(unsigned char)buffer3[i];    for (i=0;i<1024;i++) fokadr_cln1<<(unsigned char)buffer3[i];
665    //Write in file nnnnn_mmm_xxx_cln1_bin.log    //Write in file nnnnn_mmm_xxx_cln1_bin.log
# Line 690  void read(char fni[], char fno[], char I Line 722  void read(char fni[], char fno[], char I
722    oss << "INSERT INTO RAW_FILES (cirquit_n, route_n, counter, NumCadr, CRCerr, cln1_name, cln1_size)"    oss << "INSERT INTO RAW_FILES (cirquit_n, route_n, counter, NumCadr, CRCerr, cln1_name, cln1_size)"
723    << " VALUES ("<< cirquit_no << ","<< route << ","<< i+1 <<","<< Num_kadrs[i]<< "," << Num_bad_kadrs[i] << ",'"    << " VALUES ("<< cirquit_no << ","<< route << ","<< i+1 <<","<< Num_kadrs[i]<< "," << Num_bad_kadrs[i] << ",'"
724    << temp << "'," << count_ident[i]-Num_kadrs[i]*9 << ")";    << temp << "'," << count_ident[i]-Num_kadrs[i]*9 << ")";
725    if (print_flag>0) cout<<oss.str().c_str()<<" table RAW_FILES\n";          if (print_flag>0) cout<<oss.str().c_str()<<" table RAW_FILES\n";
726    TSQLResult *res=sqlServer->Query(oss.str().c_str());    TSQLResult *res=sqlServer->Query(oss.str().c_str());
727    
728    Num_good_kadrs[i]=Num_kadrs[i]-Num_bad_kadrs[i];    Num_good_kadrs[i]=Num_kadrs[i]-Num_bad_kadrs[i];
# Line 706  void read(char fni[], char fno[], char I Line 738  void read(char fni[], char fno[], char I
738    fprintf(fo," total size of route=%lli bytes;\n",count_tot[i]);    fprintf(fo," total size of route=%lli bytes;\n",count_tot[i]);
739    fprintf(fo,"first cadre number= %lli:", packet_first[i]);    fprintf(fo,"first cadre number= %lli:", packet_first[i]);
740    fprintf(fo," last cadre number= %lli:", packet_last[i]);    fprintf(fo," last cadre number= %lli:", packet_last[i]);
741    fprintf(fo," Number of missed cadres= %lli:\n", count_missed[i]);      fprintf(fo," Number of missed cadres= %lli: Number of double transmitted cadres= %lli\n",
742                 count_missed[i], count_sec_missed[i]);  
743    fprintf(fo,"number of identified cadres=%lli: quality index=%f: completness=%f \n", Num_kadrs[i], quality[i],    fprintf(fo,"number of identified cadres=%lli: quality index=%f: completness=%f \n", Num_kadrs[i], quality[i],
744            first_missed[i]);              first_missed[i]);  
745   }   }
746    
747   //***route's number log file**************************************   //***route number log file**************************************
748   for (i=0;i<Rcount;i++) {fprintf(fo_route,"%1c",(char)(NumRout[i]));}   for (i=0;i<Rcount;i++) {fprintf(fo_route,"%1c",(char)(NumRout[i]));}
749    
750                    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.23