/[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.1 by cafagna, Fri Sep 1 10:56:22 2006 UTC revision 1.3 by cafagna, Sun Dec 10 12:52:45 2006 UTC
# Line 1  Line 1 
1  #include "descript.h"  #include "descript.h"
2    #include "flag.h"
3    
4  void filename(char [], char [], char [], char[], int, char[]);  void filename(char [], char [], char [], char[], int, char[]);
5  void FindRoute(int[], int*);  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];
10    extern char MYSQL[80];
11    extern int print_flag;
12    
13  void read(char fni[], char fno[], char INF[],  void read(char fni[], char fno[], char INF[],
14                    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)
15  {  {
16   FILE *fo, *fo1, *fout, *foRout, *foMiss, *fo_log_kadr, *fo_route;   FILE *fo, *fout, *foRout, *foMiss, *fo_log_kadr, *fo_route, *fo_tmp;
17   ifstream fin;   ifstream fin;
18   ofstream fo_cln1, fokadr_cln1;   ofstream fo_cln1, fokadr_cln1, fo_info;
19    
20   char buf[1024];   char buf[1024];
21   //  char stop1;   //  char stop1;
22   char ch1,ch2,ch3, ch, file_name_prot[80];   char ch1,ch2,ch3,ch, ch4[80], ch5[5];
23   char file_route[80];  /*protocol binary filesm mmm_counter_nnnnn_cln1Er.log*/   char file_route[80];  /*protocol binary filesm nnnnn_mmm_xxx_cln1Er.log*/
24   char miss_file[80];   /*protocol files mmm_counter_nnnnn_mis_cln1.dat*/   char miss_file[80];   /*protocol files nnnnn_mmm_xxx_mis_cln1.dat*/
25   char cln1_file[80];   /*route binary data files mmm_counter_nnnnn_cln1.pam*/   char cln1_file[80];   /*route binary data files nnnnn_mmm_xxx_cln1.pam*/
26   char cln1_kadr[80];   /*route binary data files mmm_counter_nnnnn_cln1.bin*/   char cln1_kadr[80];   /*route binary data files nnnnn_mmm_xxx_cln1.bin*/
27   char log_kadr[80];    /*log binary data files mmm_counter_nnnnn_cln1.log*/   char log_kadr[80];    /*log binary data files nnnnn_mmm_xxx_cln1.log*/
28   char log_route[80];   /*log binary file for number of routes*/   char log_route[80];   /*log binary file for number of routes*/
29   char Short_Buf[8], Name_Buf[4];   char Short_Buf[8], Name_Buf[4];
30     char temp[10], temp1[80], connection[80], command[80], name11[80];
31   short int CalcCheckSum;   short int CalcCheckSum;
32    
33   int LenthBuf=1015, flag_help=0;   int LenthBuf=1015, flag_help=0;
34   int i, j, k, file_length=0, Hex[3];   int i, j, k, file_length=0, Hex[3];
35   int N_er=0;              /*Number of kadr with name error*/   int N_er=0;              /*Number of cadre with name error*/
36   int CRC0_er=0;           /*Number of kadr with CRC0 error*/   int CRC0_er=0;           /*Number of cadre with CRC0 error*/
37   int P_er=0;              /*Number of kadr with packet error number*/   int P_er=0;              /*Number of cadre with packet error number*/
  int M_er=0;              /*Number of kadr with rout error number*/  
38   int end=0;               /*End of last pamela packet*/   int end=0;               /*End of last pamela packet*/
39   // Buffers for Kadr analysis   // Buffers for Cadre analysis
40   int buffer1[1024], buffer2[1024], buffer3[1024];   int buffer1[1024], buffer2[1024], buffer3[1024];
41   int rr=0;   int rr=0;
42     int cirquit_no=0, route=0;
43    
44   int header1, header2, header3, CRC_0[3], CRC_1[3];   int header1, header2, header3, CRC_0[3], CRC_1[3];
45   long long int icount_int;   long long int icount_int;
46   long long int icount_frac, tmpP, numP[3];               long long int icount_frac, tmpP, numP[3];            
47   long int N1=0;                  /* Number of identified kadrs*/   long int N1=0;                  /* Number of identified cadres*/
48   long long int icount=0;         /* Number of bytes in identified kadrs*/   long long int icount=0;         /* Number of bytes in identified cadres*/
49   long long int icount_total=0;   /* Total number of bytes*/   long long int icount_total=0;   /* Total number of bytes*/
50   long long int packet_first[100];/*First kadr number for each route*/   long long int packet_first[100];/* First cadre number for each route*/
51   long long int packet_last[100]; /*Second kadr number for each route*/   long long int packet_last[100]; /* Second cadre number for each route*/
52   long int count_missed[100];     /* Number of missed kadrs in each route*/   long int count_missed[100];     /* Number of missed cadres in each route*/
53   long int count_sec_missed[100];   long int count_sec_missed[100];
54   unsigned long long int count_ident[100]; /* Number of good bytes in each rout*/   unsigned long long int count_ident[100]; /* Number of good bytes in each route*/
55   unsigned long long int count_tot[100];   /* Total number of bytes in each rout*/   unsigned long long int count_tot[100];   /* Total number of bytes in each route*/
56   int NumRout[1000];              /* Number of Route*/   int NumRout[1000];              /* Number of route*/
57   long long int Num_kadrs[1000];  /* Number of identified kadrs in each rout*/   long long int Num_kadrs[1000];  /* Number of identified cadres in each route*/
58   long int Num_good_kadrs[1000];  /* Number of normal kadrs in each rout*/   long int Num_good_kadrs[1000];  /* Number of normal cadres in each route*/
59   long long int Num_bad_kadrs[1000]; /* Bad kadrs number in each rout*/   long long int Num_bad_kadrs[1000]; /* Bad cadres number in each route*/
60  // float quality[100];           /* File quality pointer*/  
61   int counter=0;                  /* Last byte position in truncated kadr*/   int counter=0;                  /* Last byte position in truncated cadre*/
62   int adr1=0;                     /* Initial byte of kadr in cln1.bin*/   int adr1=0;                     /* Initial byte of cadre in cln1.bin*/
63   int adr2=0;                     /* Last byte of kadr in cln2.bin*/   int adr2=0;                     /* Last byte of cadre in cln2.bin*/
64   int tmp=0;   int tmp=0;
65   int numPack[4];                 /* Number of kadr in given route*/   int numPack[4];                 /* Number of cadre in given route*/
66     int good_pack=0;
67     long long int good_number=0;
68     int good_route=0;
69   int adress[4];   int adress[4];
70   int R_er[100];           /* Rout errors number in given rout*/   int R_er[100];           /* Route errors number in given route*/
71   int R_er_interf1[100];   /* Interface Rout error at begin of rout*/   int R_er_interf1[100];   /* Interface Route error at the begining of the route*/
72   int R_er_interf2[100];   /* Interface Rout error at end of rout*/   int R_er_interf2[100];   /* Interface Route error at end of route*/
73   int Mis_group[1000];     /* Number of missed kadrs group   int Mis_group[1000];     /* Number of missed cadres group
74                               for each rout*/                               for each route*/
75     float first_missed[100];                            
76   int Pack_Er[1000];       /* Number of packets with packet number   int Pack_Er[1000];       /* Number of packets with packet number
77                               error for each rout*/                               error for each route*/
78   int CRC_bad_kadr[3];     /* Current CRC bad kadr number*/   int CRC_bad_kadr[3];     /* Current CRC bad cadre number*/
79   int flag_pack[3];        /* Flag for error in packet number*/   int flag_pack[3];        /* Flag for error in packet number*/
80   int flag_rout[3];        /* Flag for error in rout number*/   int flag_rout[3];        /* Flag for error in route number*/
81   int flag_name[3];        /* Flag for name error*/   int flag_name[3];        /* Flag for name error*/
82   int flag_CRC0[3];        /* Flag for CRC0 error*/   int flag_CRC0[3];        /* Flag for CRC0 error*/
83   int flag_CRC[3];         /* Flag for CRC error*/   int flag_CRC[3];         /* Flag for CRC error*/
# Line 85  void read(char fni[], char fno[], char I Line 95  void read(char fni[], char fno[], char I
95    packet_last[i]=0;    packet_last[i]=0;
96    count_missed[i]=0;    count_missed[i]=0;
97    count_sec_missed[i]=0;    count_sec_missed[i]=0;
98      first_missed[i]=0;
99    count_ident[i]=0;    count_ident[i]=0;
100    count_tot[i]=0;    count_tot[i]=0;
101    R_er[i]=0;    R_er[i]=0;
# Line 98  void read(char fni[], char fno[], char I Line 109  void read(char fni[], char fno[], char I
109    flag_pack[i]=0; flag_rout[i]=0;    flag_pack[i]=0; flag_rout[i]=0;
110    flag_name[i]=0; flag_CRC0[i]=0; flag_tr[i]=0;    flag_name[i]=0; flag_CRC0[i]=0; flag_tr[i]=0;
111    flag_CRC[i]=0; flag_CRC_cor[i]=0;}    flag_CRC[i]=0; flag_CRC_cor[i]=0;}
112    
113     stringstream oss;
114    
115     strcpy(connection,MYSQL);
116     strcat(connection,database);
117     TSQLServer *sqlServer = TSQLServer::Connect(connection, "pamelaprod", "rawprod");
118   // open input file "nnnnnmmm.pam"   // open input file "nnnnnmmm.pam"
119   fin.open(fni, ifstream::binary);   fin.open(fni, ifstream::binary);
120   if(fin.good()==0) {cerr<<"can not open input file"; exit(0);}   if(fin.good()==0) {cerr<<"can not open input file"; exit(0);}
121   //open output file "nnnnnmmm_work.dat"   //open output file "nnnnnmmm_work.dat"
122   fout=fopen(fno, "wb");   fout=fopen(fno, "wb");
123   if(fout==0) {cerr<<"can not open outputfile"; exit(0);}   if(fout==0) {cerr<<"can not open work file"; exit(0);}
124    
125     strcpy(ch5,"info");
126     strcpy(ch4,fno);
127     for(i=0;i<4;i++){ch4[strlen(ch4)-8+i]=ch5[i];}
128     // open input file "nnnnnmmm_info.dat"
129     fo_info.open(ch4, ifstream::binary);
130     if(fo_info.good()==0) {cerr<<"can not open info file"; exit(0);}
131    
132   //open output file "nnnnnmmm_protocol.txt"   //open output file "nnnnnmmm_protocol.txt"
133   fo=fopen(protocol, "wb");   fo=fopen(protocol, "wb");
134   if(fo==0) {cerr<<"can not open protocol file"; exit(0);}   if(fo==0) {cerr<<"can not open protocol file"; exit(0);}
  //open output file "nnnnnmmm_cln1.inf"  
  fo1=fopen(INF, "wb");  
  if(fo1==0) {cerr<<"can not open inf file"; exit(0);}  
135   //open output file "nnnnnmmm_route.log"   //open output file "nnnnnmmm_route.log"
136   strncpy(log_route,INF,strlen(INF)-9);   strncpy(log_route,INF,strlen(INF)-9);
137   log_route[strlen(INF)-9]='\0';   log_route[strlen(INF)-9]='\0';
# Line 121  void read(char fni[], char fno[], char I Line 143  void read(char fni[], char fno[], char I
143   k=0;   k=0;
144   while(fin.eof()==0)   while(fin.eof()==0)
145   {   {
146    counter=0;          //pounter for position of kadr truncation    counter=0;          //pounter for position of cadre truncation
147    N_er=0;              N_er=0;          
   M_er=0;  
148    P_er=0;    P_er=0;
149    CRC0_er=0;    CRC0_er=0;
150    if (N1==0) Rcount=1;    if (N1==0) Rcount=1;
# Line 159  void read(char fni[], char fno[], char I Line 180  void read(char fni[], char fno[], char I
180               ||               ||
181    (ch1!=(char)(CODE1)&&ch2!=(char)(CODE2)&&ch3!=(char)(CODE3)))    (ch1!=(char)(CODE1)&&ch2!=(char)(CODE2)&&ch3!=(char)(CODE3)))
182    {fin.get(ch);    {fin.get(ch);
183     if (fin.eof()) {/*cerr<<"\n"<<"gopa";*/ break;}     if (fin.eof()) {break;}
184     ch1=ch2; ch2=ch3; ch3=Short_Buf[0];     ch1=ch2; ch2=ch3; ch3=Short_Buf[0];
185     for (i=0;i<4;i++) Short_Buf[i]=Short_Buf[i+1]; Short_Buf[4]=ch;     for (i=0;i<4;i++) Short_Buf[i]=Short_Buf[i+1]; Short_Buf[4]=ch;
186     icount_total+=1; count_tot[Rcount-1]+=1;     icount_total+=1; count_tot[Rcount-1]+=1;
# Line 170  void read(char fni[], char fno[], char I Line 191  void read(char fni[], char fno[], char I
191    if ((fin.eof())&&(N1==0))    if ((fin.eof())&&(N1==0))
192    {    {
193     fclose(fout);     fclose(fout);
194       fo_info.close();  
195     remove(fno);     remove(fno);
196     fclose(fo);     fclose(fo);
197     remove(protocol);     remove(protocol);
    fclose(fo1);  
198     remove(INF);     remove(INF);
199     exit(0);     exit(0);
200    }    }
 //      cout<<"icount_total="<<icount_total<<" "<<N1<<"\n";  
 //      printf("%x %x %x \n",(int)(unsigned char)(ch1),(int)(unsigned char)(ch2),(int)(unsigned char)(ch3));  
 //      printf("%x %x %x %x %x\n",(int)(unsigned char)(Short_Buf[0]),(int)(unsigned char)(Short_Buf[1]),  
 //                                    (int)(unsigned char)(Short_Buf[2]),(int)(unsigned char)(Short_Buf[3]),  
 //                                                        (int)(unsigned char)(Short_Buf[4]));  
 //      printf("buffer3[3]= %x\n",(int)(unsigned char)buffer3[3]);  
 //      printf("numP[2]+1= %x tmpP= %x\n",numP[2]+1, tmpP);  
 //      printf("buffer3[1023]= %x\n",(int)(unsigned char)buffer3[1023], tmpP);  
201    
202    /*****************************************************************/    /*****************************************************************/
203    /*********************data reading********************************/    /*********************data reading********************************/
# Line 230  void read(char fni[], char fno[], char I Line 243  void read(char fni[], char fno[], char I
243                     CRC_bad_kadr[2]=0;                     CRC_bad_kadr[2]=0;
244                     flag_pack[2]=0; flag_rout[2]=0;                     flag_pack[2]=0; flag_rout[2]=0;
245                     flag_name[2]=0; flag_CRC0[2]=0; flag_tr[2]=0;                     flag_name[2]=0; flag_CRC0[2]=0; flag_tr[2]=0;
246                     flag_CRC[2]=0; flag_CRC_cor[2]=0;                     flag_CRC[2]=0; flag_CRC_cor[2]=0;      
247             for (i=0; i<2; i++)             for (i=0; i<2; i++)
248                     {numP[i]=numP[2];CRC_bad_kadr[i]=CRC_bad_kadr[2];                     {numP[i]=numP[2];CRC_bad_kadr[i]=CRC_bad_kadr[2];
249                      flag_pack[i]=flag_pack[2]; flag_rout[i]=flag_rout[2];                      flag_pack[i]=flag_pack[2]; flag_rout[i]=flag_rout[2];
# Line 283  void read(char fni[], char fno[], char I Line 296  void read(char fni[], char fno[], char I
296              CRC_1[2]=(int)(unsigned char)(ch);                  CRC_1[2]=(int)(unsigned char)(ch);    
297             }             }
298    }      }  
299    //  if (fmod(N1,10000)==0.) cout<<N1<<"\n";        if (print_flag>0 && fmod(N1,10000)==0.) cout<<N1<<"\n";        
 // Caf mod on 31 Aug 2006  
   if ((N1%10000)==0 && N1>0) cout<<N1<<"\n";        
   // end mod  
300    
301  /******************CRC ERROR ANALYSIS***************************/  /******************CRC ERROR ANALYSIS***************************/
302    CalcCheckSum=CRC(buffer3);    CalcCheckSum=CRC(buffer3);
303    if (CalcCheckSum!=(short int)(buffer3[1023]))    if (CalcCheckSum!=(short int)(buffer3[1023]))
304    {    {    
305     /*if (numP[2]==1026)  //write in file nnnnnmmm_info.dat    
306     {for (i=0;i<8;i++) printf(" %x ",buffer3[i]); printf("\n");     fo_info<<"real CRC = "<<hex<<buffer3[1023]<<" (calculated CRC ="<<CalcCheckSum<<") in cadre "<<dec<<numP[2]<<"\n";
307      for (i=8;i<16;i++) printf(" %x ",buffer3[i]); printf("\n");  //write in file nnnnnmmm_work.dat      
308      for (i=1000;i<1008;i++) printf(" %x ",buffer3[i]); printf("\n");     fprintf(fout,"%i %x (%x)", numP[2], CalcCheckSum, buffer3[1023]);
309          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  
 //      cerr<<"checksum error "<<Rcount<<" Cadre="<<numP[2]<<  
 //                " "<<buffer3[3]<<"\n";  
310     flag_CRC[2]=1;     flag_CRC[2]=1;
311     Num_bad_kadrs[Rcount-1]+=1;     Num_bad_kadrs[Rcount-1]+=1;
312     CRC_bad_kadr[2]=numP[2];     CRC_bad_kadr[2]=numP[2];
313     counter=restrict(buffer3,numP,buffer2[3]);     counter=restrict(buffer3,numP,buffer2[3]);
314     if (counter!=0)     if (counter!=0)
315     {fprintf(fout," truncate byte=%i",counter);     {fo_info<<" truncate byte="<<counter<<" in cadre "<<numP[2]<<"\n";
316      cout<<" truncate counter="<<counter<<" route="<<buffer2[3]<<" "<<buffer3[3]<<" kadr="<<numP[2]<<" "<<numP[1]<<" "<<numP[0]<<"\n";      if (print_flag>0)
317        {cout<<" truncate counter="<<counter<<" route="<<buffer2[3]<<" "<<buffer3[3]<<" kadr="<<numP[2]<<" "<<numP[1]<<" "<<numP[0]<<"\n";}
318      fin.seekg(-(1024-counter),ios::cur);      fin.seekg(-(1024-counter),ios::cur);
319      icount-=1024;      icount-=1024;
320      count_ident[Rcount-1]-=1024;      count_ident[Rcount-1]-=1024;
# Line 356  void read(char fni[], char fno[], char I Line 363  void read(char fni[], char fno[], char I
363    /////////////////////////////////////////////////////////////////////    /////////////////////////////////////////////////////////////////////
364    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 ((numP[2]!=(numP[1]+1))&&(N1!=1)&&(buffer3[3]==buffer2[3])&&(CRC_bad_kadr[2]==numP[2])&&(numP[2]!=0)) flag_pack[2]=1;
365    if (flag_pack[2]!=1) packet_last[Rcount-1]=numP[2];    if (flag_pack[2]!=1) packet_last[Rcount-1]=numP[2];
366    ////////Missing kadr analysis////////////////////////////////////    ////////Missing cadre analysis////////////////////////////////////
367    /**route number: 10 (10 19 10) 10 ******************************/    /**route number:  10 (10 19 10) 10 ******************************/
368    /**kadr number:  21 (22 17 18) 19 ******************************/    /**cadre number:  21 (22 17 18) 19 ******************************/
369    /**missed kadrs: 21 20 19 18 17   ******************************/      /**missed cadres: 21  20 19 18  17 ******************************/  
370    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)&&
371            (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0))            (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0))
372    {Mis_group[Rcount-1]+=1;    {Mis_group[Rcount-1]+=1;
373     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];
374     count_sec_missed[Rcount-1]+=numP[0]-1-numP[1]-1;}       fprintf(foMiss,"secondary transmitted cadres may be %lli - %lli number of cadres=%lli\n",numP[1],numP[0],numP[0]-1-numP[1]);
375      }  
376      
377    /**route number: 10 (10 19 10) 10 ******************************/    /**route number: 10 (10 19 10) 10 ******************************/
378    /**kadr number:  11 (12 15 20) 21 ******************************/    /**cadre number: 11 (12 15 20) 21 ******************************/
379    /**missed kadrs:    16 17 18 19   ******************************/      /**missed cadres:   16 17 18 19   ******************************/  
380    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)&&
381            (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0))            (flag_pack[0]==0)&&(flag_pack[1]==0)&&(flag_pack[2]==0))
382    {Mis_group[Rcount-1]+=1;    {Mis_group[Rcount-1]+=1;
383     fprintf(foMiss,"%9i%9i\n",numP[1]+1,numP[2]-1);     fprintf(foMiss,"1) %9lli%9lli\n",numP[1]+1,numP[2]-1);
384     count_missed[Rcount-1]+=numP[2]-numP[1]-1;}     count_missed[Rcount-1]+=numP[2]-numP[1]-1;}
385    /**route number: (10 10 10) 10 10 ******************************/    /**route number: (10 10 10) 10 10 ******************************/
386    /**kadr number:  (11 12 15) 20 21 ******************************/    /**cadre number: (11 12 15) 20 21 ******************************/
387    /**missed kadrs:     13 14        ******************************/      /**missed kadrs:     13 14        ******************************/  
388    else if ((numP[1]>=numP[0]+1)&&(numP[1]<numP[2]-1)&&    else if ((numP[1]>=numP[0]+1)&&(numP[1]<numP[2]-1)&&
389                 (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&                 (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&
390                     (flag_tr[1]==0)&&(flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))                     (flag_tr[1]==0)&&(flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))
391    {Mis_group[Rcount-1]+=1;    {Mis_group[Rcount-1]+=1;
392     fprintf(foMiss,"%9i%9i\n",numP[0]+2,numP[2]-1);     fprintf(foMiss,"2) %9lli%9lli\n",numP[0]+2,numP[2]-1);
393     count_missed[Rcount-1]+=numP[2]-numP[0]-1-1;}     count_missed[Rcount-1]+=numP[2]-numP[0]-1-1;}
394    /**route number: (10 10 10) 10 10 ******************************/    /**route number: (10 10 10) 10 10 ******************************/
395    /**kadr number:  (11 29 15) 14 15 ******************************/    /**cadre number:  (11 29 15) 14 15 ******************************/
396    /**number of missed kadrs:2       ******************************/      /**number of missed cadres:2       ******************************/  
397    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)&&
398             (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&(flag_tr[1]==0)&&             (buffer1[3]==buffer3[3])&&(buffer1[3]==buffer2[3])&&(flag_tr[1]==0)&&
399             (flag_pack[2]==0)&&(flag_pack[1]==1)&&(flag_pack[0]==0))             (flag_pack[2]==0)&&(flag_pack[1]==1)&&(flag_pack[0]==0))
400    {Mis_group[Rcount-1]+=1;    {Mis_group[Rcount-1]+=1;
401     fprintf(foMiss,"%9i%9i -1\n",numP[0]+1,numP[2]-1);     fprintf(foMiss,"3) %9lli%9lli -1\n",numP[0]+1,numP[2]-1);
402     count_missed[Rcount-1]+=numP[2]-numP[0]-1-1;}     count_missed[Rcount-1]+=numP[2]-numP[0]-1-1;}
403    //truncate kadr    /**route number:(10 10 )10 10 10 10 ******************************/
404      /**cadre number:( 1 11 )12 13 14 15 ******************************/
405      /**number of missed cadres: 10       ******************************/  
406      else if ((N1>1)&&(abs(numP[2]-good_pack-2)>abs(Num_kadrs[Rcount-1]-good_number))&&
407               (buffer2[3]==good_route)&&(flag_tr[2]==0)&&(flag_pack[2]==0))
408      {Mis_group[Rcount-1]+=1;  
409       count_missed[Rcount-1]+=abs(numP[2]-good_pack)-1;
410       fprintf(foMiss,"4) %9i%9lli number of missed cadres=%i\n", good_pack, numP[2], abs(numP[2]-good_pack)-1);
411      }
412      //truncate cadres
413    /**route number: 10 (10 10 10) 10 ******************************/    /**route number: 10 (10 10 10) 10 ******************************/
414    /**kadr number:  11 (12 13 14) 15 ******************************/    /**cadre number:  11 (12 13 14) 15 ******************************/
415    /**missed truncate kadr:13        ******************************/      /**missed truncate cadre:13        ******************************/  
416    else if ((numP[0]<numP[2]-1)&&(buffer1[3]==buffer3[3])&&    else if ((numP[0]<numP[2]-1)&&(buffer1[3]==buffer3[3])&&
417                 (flag_tr[0]==0)&&(flag_tr[1]==1)&&(flag_tr[2]==0)&&                 (flag_tr[0]==0)&&(flag_tr[1]==1)&&(flag_tr[2]==0)&&
418                     (flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))                     (flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))
419    {Mis_group[Rcount-1]+=1;    {Mis_group[Rcount-1]+=1;
420     fprintf(foMiss,"%9i%9i\n",numP[0]+1,numP[2]-1);     fprintf(foMiss,"5) %9lli%9lli\n",numP[0]+1,numP[2]-1);
421     count_missed[Rcount-1]+=numP[2]-numP[0]-1;}     count_missed[Rcount-1]+=numP[2]-numP[0]-1;}
422    /***************************************************************/    /***************************************************************/
423    
424      if (flag_CRC[2]!=1) {good_pack=numP[2]; good_route=buffer3[3]; good_number=Num_kadrs[Rcount-1];}
425    
426    /******************ROUTE NUMBER ANALYSIS**************************/    /******************ROUTE NUMBER ANALYSIS**************************/
427    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])&&
428        (flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))        (flag_pack[2]==0)&&(flag_pack[1]==0)&&(flag_pack[0]==0))
429    {    {
430     /******************Interface route error 1***********************/  
    /**route number:7   7  ( 7   4  8) 8 8 **************************/  
    /**kadr number:197 198 (199 200 1) 2 3 **************************/  
    /*if ((numP[1]==numP[0]+1)&&(numP[2]<numP[0]))        
    {Rcount+=1; flag_rout[1]=1;  
     packet_last[Rcount-2]=numP[1];  
     packet_first[Rcount-1]=numP[2];  
     Num_kadrs[Rcount-1]+=1;  
     Num_kadrs[Rcount-2]-=1;                        
     R_er[Rcount-2]+=1;  
     R_er_interf2[Rcount-2]+=1;  
     M_er=N1-1;  
     count_ident[Rcount-1]+=1024; count_ident[Rcount-2]-=1024;  
     count_tot[Rcount-1]+=1024; count_tot[Rcount-2]-=1024;  
     NumRout[Rcount-1]=Rcount;  
     cout<<"\n interface1 "<<NumRout[Rcount-1]<<" "<<Rcount;  
     //    cin>>stop1;  
     if ((CRC_bad_kadr[1]==numP[1])&&(numP[1]!=0))  
     {  
     //Correct error route number//////////////////////////////////////////  
          buffer2[3]=buffer1[3];  
          if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;  
      ////////////////////////////////////////////////////////////////////  
          flags=(short int)(pow((float)(2*flag_name[1]),4)+  
                pow((float)(2*flag_rout[1]),1)+  
                pow((float)(2*flag_pack[1]),2)+  
                pow((float)(2*flag_CRC0[1]),3)+  
                    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 Interface route error 1 %i\n",CRC_bad_kadr[1],flags);  
      CRC_bad_kadr[1]=0;  
     }  
     // CREATE FILE NAMES            
         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");  
         strncpy(log_kadr,file_route,strlen(file_route)-6);  
         log_kadr[strlen(file_route)-6]='\0';  
         strcat(log_kadr,"_bin.log");  
         fclose(foRout);                    //Close files  
         fo_cln1.close();                   //for old  
         fclose(foMiss);                    //route  
         fokadr_cln1.close();  
         fclose(fo_log_kadr);  
     //Open files for new route  
         foRout=fopen(file_route, "wb");          //mmm_nnnnn_cln1Er.log  
         foMiss=fopen(miss_file, "wb");           //mmm_nnnnn_mis_cln1.dat  
         fo_cln1.open(cln1_file, ios::binary);    //mmm_nnnnn_cln1.pam    
         fokadr_cln1.open(cln1_kadr,ios::binary); //mmm_nnnnn_cln1.bin  
         fo_log_kadr=fopen(log_kadr,"wb");        //mmm_nnnnn_cln1_bin.log  
    }*/  
    /******************Interface route error 2***********************/  
    /**route number:7   7  ( 7  4 8) 8 ******************************/  
    /**kadr number:198 199 (200 1 2) 3 ******************************/  
    /*if ((numP[1]==numP[2]-1)&&(numP[1]<numP[0])&&(buffer2[3]!=0))  
    {Rcount+=1; flag_rout[1]=1;  
     packet_last[Rcount-2]=numP[0];  
     packet_first[Rcount-1]=numP[1];  
     Num_kadrs[Rcount-1]+=2;  
     Num_kadrs[Rcount-2]-=2;  
     R_er[Rcount-1]+=1;  
     R_er_interf1[Rcount-1]+=1;  
     M_er=N1-1;  
     count_ident[Rcount-1]+=2048; count_ident[Rcount-2]-=2048;  
     count_tot[Rcount-1]+=2048; count_tot[Rcount-2]-=2048;  
     NumRout[Rcount-1]=Rcount;  
     cout<<"\n interface2 "<<NumRout[Rcount-1]<<" "<<Rcount;  
 //        cin>>stop1;  
     // CREATE FILE NAMES  
         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");  
         strncpy(log_kadr,file_route,strlen(file_route)-6);  
         log_kadr[strlen(file_route)-6]='\0';  
         strcat(log_kadr,"_bin.log");  
         fclose(foRout);                    //Close files  
         fo_cln1.close();                   //for old  
         fclose(foMiss);                    //route  
         fokadr_cln1.close();  
         fclose(fo_log_kadr);  
     //Open files for new route  
         foRout=fopen(file_route, "wb");          //mmm_nnnnn_cln1Er.log  
         foMiss=fopen(miss_file, "wb");           //mmm_nnnnn_mis_cln1.dat  
         fo_cln1.open(cln1_file, ios::binary);    //mmm_nnnnn_cln1.pam  
         fokadr_cln1.open(cln1_kadr,ios::binary); //mmm_nnnnn_cln1.bin  
         fo_log_kadr=fopen(log_kadr,"wb");        //mmm_nnnnn_cln1_bin.log  
         if ((CRC_bad_kadr[1]==numP[1])&&(numP[1]!=0))  
         {Num_bad_kadrs[Rcount-1]+=1; Num_bad_kadrs[Rcount-2]-=1;  
      //Correct error rout number//////////////////////////////////////////  
          buffer2[3]=buffer3[3];  
          if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;  
      ////////////////////////////////////////////////////////////////////  
          flags=(short int)(pow((float)(2*flag_name[1]),4)+  
                 pow((float)(2*flag_rout[1]),1)+  
                 pow((float)(2*flag_pack[1]),2)+  
                 pow((float)(2*flag_CRC0[1]),3)+  
                 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 Interface route error 2 %i\n",CRC_bad_kadr[1],flags);  
          CRC_bad_kadr[1]=0;}  
    }*/  
    /******************Interface route error 3************************/  
    /**route number:7   7  ( 7  8 4) 8 8******************************/  
    /**kadr number:198 199 (200 1 2) 3 4******************************/  
 /*   if ((numP[1]==numP[2]-1)&&(flag_CRC[1]==0)&&(flag_CRC[2]==1))  
    {Rcount+=1;  
     packet_last[Rcount-2]=numP[0];  
     packet_first[Rcount-1]=numP[1];  
     NumRout[Rcount-1]=Rcount;  
     Num_kadrs[Rcount-2]-=2;  
     Num_kadrs[Rcount-1]+=2;  
     count_ident[Rcount-1]+=2048; count_ident[Rcount-2]-=2048;  
     count_tot[Rcount-1]+=2048; count_tot[Rcount-2]-=2048;  
     Num_bad_kadrs[Rcount-1]+=1; Num_bad_kadrs[Rcount-2]-=1;  
     cout<<"\n interface3 "<<NumRout[Rcount-1]<<" "<<Rcount;  
 //        cin>>stop1;  
     // CREATE FILE NAMES  
         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");  
         strncpy(log_kadr,file_route,strlen(file_route)-6);  
         log_kadr[strlen(file_route)-6]='\0';  
         strcat(log_kadr,"_bin.log");  
         fclose(foRout);                    //Close files  
         fo_cln1.close();                   //for old  
         fclose(foMiss);                    //route  
         fokadr_cln1.close();  
         fclose(fo_log_kadr);  
     //Open files for new route  
         foRout=fopen(file_route, "wb");          //mmm_nnnnn_cln1Er.log  
         foMiss=fopen(miss_file, "wb");           //mmm_nnnnn_mis_cln1.dat    
         fo_cln1.open(cln1_file, ios::binary);    //mmm_nnnnn_cln1.pam  
         fokadr_cln1.open(cln1_kadr,ios::binary); //mmm_nnnnn_cln1.bin  
         fo_log_kadr=fopen(log_kadr,"wb");        //mmm_nnnnn_cln1_bin.log    
         flags=(short int)(pow((float)(2*flag_name[1]),4)+  
                pow((float)(2*flag_rout[1]),1)+  
                pow((float)(2*flag_pack[1]),2)+  
                pow((float)(2*flag_CRC0[1]),3)+  
                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 Interface route error 3 %i\n",CRC_bad_kadr[1],flags);  
         CRC_bad_kadr[1]=0;  
    }*/  
431    }    }
   /**************Return to old route*******************************/  
   /**route number: 7 (8 4 8) 8*************************************/  
   /*if ((buffer2[3]!=buffer1[3])&&(buffer2[3]!=buffer3[3])&&  
       (buffer1[3]==buffer3[3])&&(buffer2[3]!=0)&&(N1>2))  
   {                
    R_er[Rcount-1]+=1; flag_rout[1]=1;  
    M_er=N1-1;  
    if ((CRC_bad_kadr[1]==numP[1])&&(numP[1]!=0))  
    {  
     //Correct error route number//////////////////////////////////////////  
     buffer2[3]=buffer1[3];  
     if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;  
     ////////////////////////////////////////////////////////////////////  
     flags=(short int)(pow((float)(2*flag_name[1]),4)+  
                pow((float)(2*flag_rout[1]),1)+  
                pow((float)(2*flag_pack[1]),2)+  
                pow((float)(2*flag_CRC0[1]),3)+  
                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 Old route error %i\n",CRC_bad_kadr[1],flags);  
     CRC_bad_kadr[1]=0;}  
    cout<<"\n"<<Rcount<<" old "<<NumRout[Rcount-1]<<" "<<buffer1[3]<<" "<<buffer2[3]<<" "<<buffer3[3]<<" N1="<<N1-1<<"\n";  
   }*/  
432    /**************First route***************************************/    /**************First route***************************************/
433    /**route number: (0 4) 4*****************************************/    /**route number: (0 4) 4*****************************************/
434    /**kadr number:  (0 1) 2*****************************************/    /**cadre number: (0 1) 2*****************************************/
435    if (N1==2)    if (N1==2)
436    {packet_first[Rcount-1]=numP[1];    {packet_first[Rcount-1]=numP[1];
437     packet_last[Rcount-1]=numP[2];     packet_last[Rcount-1]=numP[2];
438     NumRout[Rcount-1]=1;     NumRout[Rcount-1]=1;
439     // CREATE FILE NAMES     // CREATE FILE NAMES
440     filename(file_route, miss_file, cln1_file, path,     filename(file_route, miss_file, cln1_file, path, NumRout[Rcount-1], name);                    
             NumRout[Rcount-1], name);                      
441     strncpy(cln1_kadr,cln1_file,strlen(cln1_file)-3);     strncpy(cln1_kadr,cln1_file,strlen(cln1_file)-3);
442     cln1_kadr[strlen(cln1_file)-3]='\0';     cln1_kadr[strlen(cln1_file)-3]='\0';
443     strcat(cln1_kadr,"bin");     strcat(cln1_kadr,"bin");
# Line 602  void read(char fni[], char fno[], char I Line 445  void read(char fni[], char fno[], char I
445     log_kadr[strlen(file_route)-6]='\0';     log_kadr[strlen(file_route)-6]='\0';
446     strcat(log_kadr,"_bin.log");     strcat(log_kadr,"_bin.log");
447     //Open files for new route     //Open files for new route
448     foRout=fopen(file_route, "wb");          //mmm_nnnnn_cln1Er.log     foRout=fopen(file_route, "wb");          //nnnnn_mmm_xxx_cln1Er.log
449     foMiss=fopen(miss_file, "wb");           //mmm_nnnnn_mis_cln1.dat     foMiss=fopen(miss_file, "wb");           //nnnnn_mmm_xxx_mis_cln1.dat
450     fo_cln1.open(cln1_file, ios::binary);    //mmm_nnnnn_cln1.pam     if (CLN1_MEPHI>0)
451       {strcpy(ch4,"");
452        strncat(ch4,path,strlen(path)-6);
453        strcat(strcpy(temp1,ch4),"Cln1/tmp.dat");
454        if ((fo_tmp=fopen(temp1,"w"))==NULL)
455        {strcat(strcat(strcpy(command,"mkdir "),ch4),"Cln1");
456         //Creata special directory for cln1 files
457         system(command);
458        }
459        else
460        {
461         fclose(fo_tmp);
462         remove(temp1);
463        }
464        for(i=0;i<22;i++)
465        {name11[i]=cln1_file[i+strlen(cln1_file)-22];}
466        name11[22]='\0';
467        strcat(strcat(ch4,"Cln1/"),name11);
468        strcpy(cln1_file,ch4);
469       }
470       fo_cln1.open(cln1_file, ios::binary);   //nnnnn_mmm_xxx_cln1.pam
471      
472     if (!fo_cln1) cout<<"can not open "<<cln1_file<<"\n";     if (!fo_cln1) cout<<"can not open "<<cln1_file<<"\n";
473     fokadr_cln1.open(cln1_kadr,ios::binary); //mmm_nnnnn_cln1.bin     fokadr_cln1.open(cln1_kadr,ios::binary); //nnnnn_mmm_xxx_cln1.bin
474     if (!fokadr_cln1) cout<<"can not open "<<cln1_kadr<<"\n";     if (!fokadr_cln1) cout<<"can not open "<<cln1_kadr<<"\n";
475     fo_log_kadr=fopen(log_kadr,"wb");        //mmm_nnnnn_cln1_bin.log     fo_log_kadr=fopen(log_kadr,"wb");        //nnnnn_mmm_xxx_cln1_bin.log
    //  cout<<"\n"<<Rcount<<" old "<<NumRout[Rcount-1]<<" "<<numP[1]<<" N1="<<N1-1;  
476    }    }
477    /***************New route***************************************/    /***************New route****************************************/
478    /**route number: 7  ( 7  0 8) 8*********************************/    /**route number:   7 ( 7  0 8) 8*********************************/
479    /**kadr number: 199 (200 0 1) 2*********************************/    /**cadre number: 199 (200 0 1) 2*********************************/
480    /**route number: 7  ( 7  8 8) 8*********************************/    /**route number:   7 ( 7  8 8) 8*********************************/
481    /**kadr number: 199 (200 0 1) 2*********************************/    /**cadre number: 199 (200 0 1) 2*********************************/
482    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))||
483         ((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))||
484         ((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)) )
# Line 623  void read(char fni[], char fno[], char I Line 486  void read(char fni[], char fno[], char I
486     packet_last[Rcount-2]=numP[0];     packet_last[Rcount-2]=numP[0];
487     packet_first[Rcount-1]=numP[1];     packet_first[Rcount-1]=numP[1];
488     NumRout[Rcount-1]=Rcount;     NumRout[Rcount-1]=Rcount;
489     cout<<"New route="<<Rcount<<" "<<buffer1[3]<<" "<<buffer2[3]<<" "<<buffer3[3]<<"\n";     if (print_flag>0)
490     cout<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";cin>>i;     {cout<<"New route="<<Rcount<<" "<<buffer1[3]<<" "<<buffer2[3]<<" "<<buffer3[3]<<"\n";
491        cout<<numP[0]<<" "<<numP[1]<<" "<<numP[2]<<"\n";}
492     if (flag_tr[1]==0)     if (flag_tr[1]==0)
493     {Num_kadrs[Rcount-2]-=2;     {Num_kadrs[Rcount-2]-=2;
494      Num_kadrs[Rcount-1]+=2;      Num_kadrs[Rcount-1]+=2;
# Line 636  void read(char fni[], char fno[], char I Line 500  void read(char fni[], char fno[], char I
500      count_ident[Rcount-1]+=1024; count_ident[Rcount-2]-=1024;      count_ident[Rcount-1]+=1024; count_ident[Rcount-2]-=1024;
501      count_tot[Rcount-1]+=1024; count_tot[Rcount-2]-=1024;}      count_tot[Rcount-1]+=1024; count_tot[Rcount-2]-=1024;}
502     // CREATE FILE NAMES     // CREATE FILE NAMES
503     filename(file_route, miss_file, cln1_file, path,     filename(file_route, miss_file, cln1_file, path, NumRout[Rcount-1], name);                    
             NumRout[Rcount-1], name);                      
504     strncpy(cln1_kadr,cln1_file,strlen(cln1_file)-3);     strncpy(cln1_kadr,cln1_file,strlen(cln1_file)-3);
505     cln1_kadr[strlen(cln1_file)-3]='\0';     cln1_kadr[strlen(cln1_file)-3]='\0';
506     strcat(cln1_kadr,"bin");     strcat(cln1_kadr,"bin");
# Line 650  void read(char fni[], char fno[], char I Line 513  void read(char fni[], char fno[], char I
513     fokadr_cln1.close();     fokadr_cln1.close();
514     fclose(fo_log_kadr);     fclose(fo_log_kadr);
515     //Open files for new route     //Open files for new route
516     foRout=fopen(file_route, "wb");          //mmm_nnnnn_cln1Er.log     foRout=fopen(file_route, "wb");          //nnnnn_mmm_xxx_cln1Er.log
517     foMiss=fopen(miss_file, "wb");           //mmm_nnnnn_mis_cln1.dat         foMiss=fopen(miss_file, "wb");           //nnnnn_mmm_xxx_mis_cln1.dat    
518     fo_cln1.open(cln1_file, ios::binary);    //mmm_nnnnn_cln1.pam     fo_cln1.open(cln1_file, ios::binary);    //nnnnn_mmm_xxx_cln1.pam
519     fokadr_cln1.open(cln1_kadr,ios::binary); //mmm_nnnnn_cln1.bin     fokadr_cln1.open(cln1_kadr,ios::binary); //nnnnn_mmm_xxx_cln1.bin
520     fo_log_kadr=fopen(log_kadr,"wb");        //mmm_nnnnn_cln1_bin.log     fo_log_kadr=fopen(log_kadr,"wb");        //nnnnn_mmm_xxx_cln1_bin.log
521     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[2]==numP[2])&&(numP[2]!=0)) {Num_bad_kadrs[Rcount-1]+=1; Num_bad_kadrs[Rcount-2]-=1;}
522     if ((CRC_bad_kadr[1]==numP[1])&&(numP[1]!=0))     if ((CRC_bad_kadr[1]==numP[1])&&(numP[1]!=0))
523     {Num_bad_kadrs[Rcount-1]+=1; Num_bad_kadrs[Rcount-2]-=1;     {Num_bad_kadrs[Rcount-1]+=1; Num_bad_kadrs[Rcount-2]-=1;
# Line 665  void read(char fni[], char fno[], char I Line 528  void read(char fni[], char fno[], char I
528                 pow((float)(2*flag_CRC[1]),5)+                 pow((float)(2*flag_CRC[1]),5)+
529                 pow((float)(2*flag_CRC_cor[1]),6));                 pow((float)(2*flag_CRC_cor[1]),6));
530      fprintf(foRout,"%1c%1c%1c%1c",(char)(buffer2[4]),(char)(buffer2[5]),(char)(buffer2[6]),(char)(flags));      fprintf(foRout,"%1c%1c%1c%1c",(char)(buffer2[4]),(char)(buffer2[5]),(char)(buffer2[6]),(char)(flags));
531      fprintf(fout,"%i New route %i\n",CRC_bad_kadr[1],flags);      fo_info<<"New route: cadre="<<CRC_bad_kadr[1]<<" error flag="<<flags<<"\n";
532      CRC_bad_kadr[1]=-1;}      CRC_bad_kadr[1]-=1;}
   // cout<<Rcount<<" New "<<NumRout[Rcount-1]<<" "<<numP[1]<<"\n";  
533    }    }
534        
535    if ((numP[1]==0)&&(N1>1))    if (((numP[1]==0)&&(N1>1))||(N1==2))
536    {rr=0;    {rr=0;
537     FindRoute(buffer2, &rr);     FindRoute(buffer2, &rr);
538     if (rr>0)     if (rr>0)
# Line 686  void read(char fni[], char fno[], char I Line 548  void read(char fni[], char fno[], char I
548    if (CRC_bad_kadr[1]==numP[1])    if (CRC_bad_kadr[1]==numP[1])
549    {    {
550     //Correct CRC0 error/////////////////////////////////////////////////     //Correct CRC0 error/////////////////////////////////////////////////
551     buffer2[7]=buffer1[1023];     if (flag_CRC0[1]!=0)
552     if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;     {
553        buffer2[7]=buffer1[1023];
554        if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
555       }
556     //Correct name error/////////////////////////////////////////////////     //Correct name error/////////////////////////////////////////////////
557     buffer2[0]=(int)(CODE1);     if (flag_name[1]!=0)
558     buffer2[1]=(int)(CODE2);     {buffer2[0]=(int)(CODE1);
559     buffer2[2]=(int)(CODE3);      buffer2[1]=(int)(CODE2);
560     if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;      buffer2[2]=(int)(CODE3);
561        if (buffer2[1023]==(int)(CRC(buffer2))) flag_CRC_cor[1]=1;
562       }
563    }    }
564    ////////////////////////////////////////////////////////////////////    ////////////////////////////////////////////////////////////////////
565    /*******************************************************************/    /*******************************************************************/
# Line 706  void read(char fni[], char fno[], char I Line 573  void read(char fni[], char fno[], char I
573                                   pow((float)(2*flag_CRC_cor[1]),6));                                   pow((float)(2*flag_CRC_cor[1]),6));
574     fprintf(foRout,"%1c%1c%1c%1c",(char)(buffer2[4]),     fprintf(foRout,"%1c%1c%1c%1c",(char)(buffer2[4]),
575             (char)(buffer2[5]),(char)(buffer2[6]),(char)(flags));             (char)(buffer2[5]),(char)(buffer2[6]),(char)(flags));
576     fprintf(fout,"%i flags=%i\n",CRC_bad_kadr[1],flags);}     fo_info<<"CRC error in cadre "<<CRC_bad_kadr[1]<<" error flag="<<flags<<"\n";}
577      
578    /***************************************************************/    /***************************************************************/
 //   cout <<"\nCalcChekSum="<<dec <<right<<CalcCheckSum  
 //   <<"\npacket.CheckSum1="<<right<<packet.CheckSum1;  
 //   cout <<"\nPacket number="<<packet.Pckt_Num<<" "<<endl;  
579    //*Packet analysis**********************************************    //*Packet analysis**********************************************
580    if (N_er!=0) fprintf(fout,"Name error %i in kadr %i\n",numP[2],N_er);    if (N_er!=0) fo_info<<"Name error "<<N_er<<" in cadre "<<numP[2]<<"\n";
581    if (M_er!=0) fprintf(fout,"Root number error %i in kadr %i\n",numP[1],M_er);    if (P_er!=0) fo_info<<"Number packet error "<<numP[1]<<" in cadre "<<P_er<<"\n";
582    if (P_er!=0) fprintf(fout,"Number packet error %i in kadr %i\n",numP[1], P_er);    if (CRC0_er!=0) fo_info<<"CRC0 error "<<numP[1]<<" in cadre "<<CRC0_er<<"\n";
   if (CRC0_er!=0) fprintf(fout,"CRC0 error %i in kadr %i\n",numP[1], CRC0_er);  
583    
584    /***************PAMELA data extracting for each route***********/    /***************PAMELA data extracting for each route***********/
585  //  if (N1!=1&&(flag_CRC[1]==0||(flag_rout[1]==0&&flag_pack[1]==0)||  //  if (N1!=1&&(flag_CRC[1]==0||(flag_rout[1]==0&&flag_pack[1]==0)||
586  //        (flag_CRC[1]!=0&&flag_CRC_cor[1]==1)))  //        (flag_CRC[1]!=0&&flag_CRC_cor[1]==1)))
587    if (N1!=1)    if ((N1!=1)&&(flag_tr[1]==0))
588    {end=0;    {end=0;
589     if (buffer2[3]!=buffer1[3]||buffer2[3]!=buffer3[3])     if (buffer2[3]!=buffer1[3]||buffer2[3]!=buffer3[3])
590     {for (i=8;i<1020;i++)     {for (i=8;i<1020;i++)
# Line 734  void read(char fni[], char fno[], char I Line 597  void read(char fni[], char fno[], char I
597     }     }
598     else {for(i=8;i<1023;i++) fo_cln1<<(unsigned char)buffer2[i];}     else {for(i=8;i<1023;i++) fo_cln1<<(unsigned char)buffer2[i];}
599     for (i=0;i<1024;i++) {fokadr_cln1<<(unsigned char)buffer2[i];}     for (i=0;i<1024;i++) {fokadr_cln1<<(unsigned char)buffer2[i];}
600     //Write in file mmm_counter_nnnnn_cln1_bin.log     //Write in file nnnnn_mmm_xxx_cln1_bin.log
601     tmp=numP[1];     tmp=numP[1];
602     for (k=0;k<4;k++) numPack[k]=0;     for (k=0;k<4;k++) numPack[k]=0;
603     for (k=0;k<4;k++)     for (k=0;k<3;k++)
604     {if (tmp/(long int)(pow(256.,3-k))>=1)     {if (tmp/(long int)(pow(256.,2-k))>=1)
605          {numPack[3-k]=(int)(tmp/(long int)(pow(256.,3-k)));          {numPack[2-k]=(int)(tmp/(long int)(pow(256.,2-k)));
606           tmp=tmp-(long int)pow(256.,3-k)*(int)(tmp/(long int)(pow(256.,3-k)));}}                   tmp=tmp-(long int)pow(256.,2-k)*(int)(tmp/(long int)(pow(256.,2-k)));}}        
607     for (i=0;i<4;i++) fputc((int)(unsigned char)(numPack[i]),fo_log_kadr);     for (i=0;i<3;i++) fputc((int)(unsigned char)(numPack[i]),fo_log_kadr);
608     adr1=icount-1024*2;     adr1=icount-1024*2;
609     for (k=0;k<4;k++) adress[k]=0;     for (k=0;k<4;k++) adress[k]=0;
610     for (k=0;k<4;k++)     for (k=0;k<4;k++)
# Line 767  void read(char fni[], char fno[], char I Line 630  void read(char fni[], char fno[], char I
630       (flag_CRC[2]!=0&&flag_CRC_cor[2]==1))       (flag_CRC[2]!=0&&flag_CRC_cor[2]==1))
631   {end=0;   {end=0;
632    for (i=0;i<1024;i++) fokadr_cln1<<(unsigned char)buffer3[i];    for (i=0;i<1024;i++) fokadr_cln1<<(unsigned char)buffer3[i];
633    //Write in file mmm_counter_nnnnn_cln1_bin.log    //Write in file nnnnn_mmm_xxx_cln1_bin.log
634    tmp=numP[2];    tmp=numP[2];
635    for (k=0;k<4;k++) numPack[k]=0;    for (k=0;k<4;k++) numPack[k]=0;
636    for (k=0;k<4;k++)    for (k=0;k<3;k++)
637     {if (tmp/(long int)(pow(256.,3-k))>=1)     {if (tmp/(long int)(pow(256.,2-k))>=1)
638      {numPack[3-k]=(int)(tmp/(long int)(pow(256.,3-k)));      {numPack[2-k]=(int)(tmp/(long int)(pow(256.,2-k)));
639       tmp=tmp-(long int)pow(256.,3-k)*       tmp=tmp-(long int)pow(256.,2-k)*
640       (int)(tmp/(long int)(pow(256.,3-k)));}}           (int)(tmp/(long int)(pow(256.,2-k)));}}    
641     for (i=0;i<4;i++) fputc((int)(unsigned char)(numPack[i]),fo_log_kadr);     for (i=0;i<3;i++) fputc((int)(unsigned char)(numPack[i]),fo_log_kadr);
642     adr1=icount-1024;     adr1=icount-1024;
643     for (k=0;k<4;k++) adress[k]=0;     for (k=0;k<4;k++) adress[k]=0;
644     for (k=0;k<4;k++)     for (k=0;k<4;k++)
# Line 806  void read(char fni[], char fno[], char I Line 669  void read(char fni[], char fno[], char I
669    
670   icount_int = (int)(icount/1024);   icount_int = (int)(icount/1024);
671   icount_frac= icount%1024;   icount_frac= icount%1024;
 // cout<<"\nfile size="<<icount_total<<" Byte   "<<N1  
 //     <<" identified cadrs   Number of routs="<<Rcount<<"\n";  
 // cout<<" size of identified cadrs="<<icount<<" bytes   "<<count_ident[0]<<"\n";  
672    
673  //***protocol file***********************************************  //***protocol file***********************************************
  strcpy(file_name_prot, path);  
  strcat(file_name_prot, name);  
674   file_length=(int)(strlen(name));   file_length=(int)(strlen(name));
675   fprintf(fo,"file name=%s\n",file_name_prot);   fprintf(fo,"file name=/home/ntsomz/downlinks/%s\n",name);
676   fprintf(fo,"file size= %i,%i kByte \n", icount_int,   fprintf(fo,"file size= %i,%i kByte \n", icount_int, icount_frac);
677                                                         icount_frac);   fprintf(fo,"number identified cadres= %i\n", N1);
678   for (i=0;i<Rcount;i++)   fprintf(fo,"number of routes= %i\n", Rcount);
679   {fprintf(fo,"route=%i;",NumRout[i]);  
680    fprintf(fo," size of identified cadrs =%lli bytes;",count_ident[i]);   for(i=0;i<5;i++) {temp[i]=name[i]; temp[5]='\0';}
681     sscanf(temp,"%5d", &cirquit_no);
682     for(i=0;i<3;i++) {temp[i]=name[i+5]; temp[3]='\0';}
683     sscanf(temp,"%3i", &route);
684      
685     for (i=0;i<Rcount;i++)
686     {
687      for(j=0;j<22;j++) temp[j]=cln1_file[j+strlen(cln1_file)-22];
688      temp[22]='\0';
689      oss.str("");
690      oss << "INSERT INTO RAW_FILES (cirquit_n, route_n, counter, NumCadr, CRCerr, cln1_name, cln1_size)"
691      << " VALUES ("<< cirquit_no << ","<< route << ","<< i+1 <<","<< Num_kadrs[i]<< "," << Num_bad_kadrs[i] << ",'"
692      << temp << "'," << count_ident[i]-Num_kadrs[i]*9 << ")";
693      if (print_flag>0) cout<<oss.str().c_str()<<" table RAW_FILES\n";      
694      TSQLResult *res=sqlServer->Query(oss.str().c_str());
695    
696      Num_good_kadrs[i]=Num_kadrs[i]-Num_bad_kadrs[i];
697      if (packet_first[i]>1)
698      {first_missed[i]=(float)(1.)-float(packet_first[i]-1.)/(float)(packet_last[i]);}
699      else if (packet_first[i]==0)
700      {first_missed[i]=1.;}
701      
702    
703      quality[i]=(float)(1.)-(float)(Num_bad_kadrs[i])/(float)(Num_kadrs[i]);
704      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]);
705      fprintf(fo,"size of identified cadres =%lli bytes;",count_ident[i]);
706    fprintf(fo," total size of route=%lli bytes;\n",count_tot[i]);    fprintf(fo," total size of route=%lli bytes;\n",count_tot[i]);
707    fprintf(fo,"number of errors=%i;\n",Num_bad_kadrs[i]);    fprintf(fo,"first cadre number= %lli:", packet_first[i]);
708      fprintf(fo," last cadre number= %lli:", packet_last[i]);
709      fprintf(fo," Number of missed cadres= %lli:\n", count_missed[i]);  
710      fprintf(fo,"number of identified cadres=%lli: quality index=%f: completness=%f \n", Num_kadrs[i], quality[i],
711              first_missed[i]);  
712   }   }
713    
714   //***route's number log file**************************************   //***route's number log file**************************************
715   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]));}
716    
717  //***cln1.inf file*************************************************          
  fprintf(fo1,"file name=");  
  fprintf(fo1,"%s%s",path,name);  
  fprintf(fo1,"\nfile size= %i,%i kByte \n", icount_int,  
                                                        icount_frac);  
  fprintf(fo1,"number identified kadrs= %i\n", N1);  
  fprintf(fo1,"number of routes= %i\n", Rcount);  
   
  for (i=0;i<Rcount;i++)  
  {  
   Num_good_kadrs[i]=Num_kadrs[i]-Num_bad_kadrs[i];  
   quality[i]=(float)(1.)-(float)(Num_bad_kadrs[i])/(float)(Num_kadrs[i]);  
   fprintf(fo1,"\nRoute number=%i: number of normal kadrs= %i:\n",  
                   NumRout[i],Num_good_kadrs[i]);  
   fprintf(fo1,"first kadre number= %lli:", packet_first[i]);  
   fprintf(fo1," last kadre number= %lli:", packet_last[i]);  
   fprintf(fo1," Number of missed kadrs= %i:\n", count_missed[i]);  
 //  fprintf(fo1,"Number of double transmitted kadrs= %i:\n", count_sec_missed[i]);    
   fprintf(fo1,"number of identified kadres=%lli: quality index=%f\n", Num_kadrs[i], quality[i]);  
  }        
718  //***************************************************************  //***************************************************************
719    fin.close();    fin.close();
720    fclose(fout);    fclose(fout);  
721      fo_info.close();
722    fclose(fo);    fclose(fo);
   fclose(fo1);  
723    fclose(foRout);    fclose(foRout);
724    fclose(foMiss);    fclose(foMiss);
725    fo_cln1.close();    fo_cln1.close();
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23