/[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.2 by cafagna, Sun Dec 10 12:21:23 2006 UTC revision 1.3 by cafagna, Sun Dec 10 12:52:45 2006 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];
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, *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;   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;   int good_pack=0;
67   int good_number=0;   long long int good_number=0;
68   int good_route=0;   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 route*/                               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 route*/                               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 89  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 102  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);}
# Line 122  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 171  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);
# Line 222  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 275  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";        
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]);
         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  
309    
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 344  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    /**route number:10 (10 10 10) 10 10 ******************************/    /**route number:(10 10 )10 10 10 10 ******************************/
404    /**kadr number:  1 (11 12 13) 14 15 ******************************/    /**cadre number:( 1 11 )12 13 14 15 ******************************/
405    /**number of missed kadrs: 10       ******************************/      /**number of missed cadres: 10       ******************************/  
406    else if ((N1>1)&&(abs(numP[2]-good_pack-2)>abs(Num_kadrs[Rcount-1]-good_number))&&    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))             (buffer2[3]==good_route)&&(flag_tr[2]==0)&&(flag_pack[2]==0))
408    {Mis_group[Rcount-1]+=1;    {Mis_group[Rcount-1]+=1;  
409     fprintf(foMiss,"%9i%9i \n",good_pack+1,numP[2]-1);     count_missed[Rcount-1]+=abs(numP[2]-good_pack)-1;
410     count_missed[Rcount-1]+=numP[2]-good_pack-1;     fprintf(foMiss,"4) %9i%9lli number of missed cadres=%i\n", good_pack, numP[2], abs(numP[2]-good_pack)-1);
411    }    }
412    //truncate kadr    //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    
# Line 410  void read(char fni[], char fno[], char I Line 431  void read(char fni[], char fno[], char I
431    }    }
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 425  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");          //xxx_mmm_nnnnn_cln1Er.log     foRout=fopen(file_route, "wb");          //nnnnn_mmm_xxx_cln1Er.log
449     foMiss=fopen(miss_file, "wb");           //xxx_mmm_nnnnn_mis_cln1.dat     foMiss=fopen(miss_file, "wb");           //nnnnn_mmm_xxx_mis_cln1.dat
450     fo_cln1.open(cln1_file, ios::binary);    //xxx_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); //xxx_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");        //xxx_mmm_nnnnn_cln1_bin.log     fo_log_kadr=fopen(log_kadr,"wb");        //nnnnn_mmm_xxx_cln1_bin.log
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 445  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";     {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 458  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 472  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");          //xxx_mmm_nnnnn_cln1Er.log     foRout=fopen(file_route, "wb");          //nnnnn_mmm_xxx_cln1Er.log
517     foMiss=fopen(miss_file, "wb");           //xxx_mmm_nnnnn_mis_cln1.dat         foMiss=fopen(miss_file, "wb");           //nnnnn_mmm_xxx_mis_cln1.dat    
518     fo_cln1.open(cln1_file, ios::binary);    //xxx_mmm_nnnnn_cln1.pam     fo_cln1.open(cln1_file, ios::binary);    //nnnnn_mmm_xxx_cln1.pam
519     fokadr_cln1.open(cln1_kadr,ios::binary); //xxx_mmm_nnnnn_cln1.bin     fokadr_cln1.open(cln1_kadr,ios::binary); //nnnnn_mmm_xxx_cln1.bin
520     fo_log_kadr=fopen(log_kadr,"wb");        //xxx_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 487  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;}
533    }    }
534        
535    if (((numP[1]==0)&&(N1>1))||(N1==2))    if (((numP[1]==0)&&(N1>1))||(N1==2))
# Line 507  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 527  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    /***************************************************************/    /***************************************************************/
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)||
# Line 552  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<3;k++)     for (k=0;k<3;k++)
# Line 585  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<3;k++)    for (k=0;k<3;k++)
# Line 631  void read(char fni[], char fno[], char I Line 676  void read(char fni[], char fno[], char I
676   fprintf(fo,"file size= %i,%i kByte \n", icount_int, icount_frac);   fprintf(fo,"file size= %i,%i kByte \n", icount_int, icount_frac);
677   fprintf(fo,"number identified cadres= %i\n", N1);   fprintf(fo,"number identified cadres= %i\n", N1);
678   fprintf(fo,"number of routes= %i\n", Rcount);   fprintf(fo,"number of routes= %i\n", Rcount);
679    
680     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++)   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];    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]);    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]);    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]);    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," first cadre number= %lli:", packet_first[i]);    fprintf(fo,"first cadre number= %lli:", packet_first[i]);
708    fprintf(fo," last cadre number= %lli:", packet_last[i]);    fprintf(fo," last cadre number= %lli:", packet_last[i]);
709    fprintf(fo," Number of missed cadres= %i:\n", count_missed[i]);      fprintf(fo," Number of missed cadres= %lli:\n", count_missed[i]);  
710    fprintf(fo,"number of identified cadres=%lli: quality index=%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],
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                    
718  //***************************************************************  //***************************************************************
719    fin.close();    fin.close();
720    fclose(fout);    fclose(fout);  
721      fo_info.close();
722    fclose(fo);    fclose(fo);
723    fclose(foRout);    fclose(foRout);
724    fclose(foMiss);    fclose(foMiss);
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23