--- calo/flight/FQLOOK/macros/FCaloQLOOK.cxx	2006/08/04 15:24:01	1.14
+++ calo/flight/FQLOOK/macros/FCaloQLOOK.cxx	2007/07/18 07:57:13	1.21
@@ -1,12 +1,18 @@
 //
 //   Check for possible errors and shows raw distribution of variables coming from the DSP (no calibrations needed) - Emiliano Mocchiutti
 //
-//   FCaloQLOOK.c      version 1.11  (2006-07-17)
+//   FCaloQLOOK.c      version 1.14  (2006-09-29)
 //
 //   The only input needed is the path to the directory created by YODA for the data file you want to analyze.
 //
 //   Changelog:
 //
+//   1.13 - 1.14 (2006-09-29): Flag the last entry in the txt file, count missing section when checking for DSP entries correspondence.
+//
+//   1.12 - 1.13 (2006-09-28): Give wrong number of events in compress mode, fixed.
+//
+//   1.11 - 1.12 (2006-09-22): Check for different number of events, put output in a file and write down errors for each event. Fixed jumpev bug.
+//
 //   1.09 - 1.11 (2006-07-17): Adapted to flight conditions.
 //
 //   1.07 - 1.09 (2006-05-29): Fixed bug in output filename when input is not in the form DW_YYMMDD_NNN. Bug in latchup alarm not coming from the status word, fixed.
@@ -80,12 +86,30 @@
   } else {
     Int_t posiz2 = 0;
     stringcopy(file2,gSystem->BaseName(filename.Data()),posiz2,posiz);
-    TString pdat(".dat");
+    TString pdat(".pam");
     stringappend(file2,pdat);  
   };
   return file2;
 }
 
+void sighandler(UInt_t pkt_num, UInt_t obt, UInt_t entry, UInt_t section, UInt_t totev, Bool_t &evdone){
+  if ( !evdone ){    
+    printf("\n\n=================================================\n");
+    printf("PSCU-Pkt N. %u - OBT %u ms\n",pkt_num,obt);
+    printf("Total events %u\n",totev);
+    printf("(ROOT-tree entry %u )\n",entry);
+    if ( (totev-1) == entry ) printf("NOTICE: this is last entry! \n");
+    printf("=================================================\n");
+    evdone = true;
+  };
+  TString sec = "UNKNOWN";
+  if ( section == 0 ) sec = "XE";
+  if ( section == 1 ) sec = "XO";
+  if ( section == 2 ) sec = "YE";
+  if ( section == 3 ) sec = "YO";
+  printf(" Section %s (DSP n. %i) message: ",sec.Data(),section);
+};
+
 void FCaloQLOOK(TString filename, Int_t fromevent=0, Int_t toevent=0, TString outDir="", TString saveas = "png", Bool_t iactive =false, Bool_t w4i=false){
   gStyle->SetPaperSize(19.,25.);
   const char* startingdir = gSystem->WorkingDirectory();
@@ -116,7 +140,8 @@
   //
   Float_t ctshitthr = 0.65;
   Float_t cbasethr = 0.95;
-  Float_t cdexythr = 0.995;
+  //  Float_t cdexythr = 0.995;
+  Float_t cdexythr = 0.95;
   //Float_t cdexythr = 0.90;
   Float_t cdexycthr = 0.95;
   Float_t h1rth = 0.90;
@@ -167,7 +192,8 @@
   };
   Int_t const size = nevents;
   Double_t iev2[size+1];
-  Int_t ver[4][23], lver[4][23],shit[4][11], rshit[4][11], se,bl, pl, fulldiff[4],cshit, cpre[4],nullsh[4][11], lalarm[4];
+  //aumentata la lunghezza di ver, prima era 23
+  Int_t ver[4][24], lver[4][23],shit[4][11], rshit[4][11], se,bl, pl, fulldiff[4],cshit, cpre[4],nullsh[4][11], lalarm[4];
   se = 0;
   pl = 0;    
   Float_t allbase = 0.;
@@ -269,7 +295,7 @@
       maxevent = fromevent - 1;
     };
   };
-    
+  //azzera i tipi di errore 
   for (Int_t k = 0; k < 4; k++ ){
     for (Int_t m = 0; m < 23 ; m++ ){
       ver[k][m] = 0 ;
@@ -281,13 +307,15 @@
       };
       lalarm[k] = 0;
     };
+    //il 23 lo azzero cosi per non modificare le altre variabili
+    ver[k][23]=0;
   };
   Int_t pdone, bdone;
   pdone = 0;
   bdone = 0;
-  bool isCOMP = 0;
-  bool isFULL = 0;
-  bool isRAW = 0;
+  bool isCOMP = false;
+  bool isFULL = false;
+  bool isRAW = false;
   Int_t alldexy=0;
   Int_t alldexy2=0;
   Int_t planebases=0;
@@ -306,7 +334,8 @@
   fulldiff[1] = 0;
   fulldiff[2] = 0;
   fulldiff[3] = 0;   
-  printf("\n Processed events: \n\n");
+  //  printf("\n Processed events: \n\n");
+  printf("\n Start processing: \n\n");
   unsigned short int calev0=0;
   unsigned short int oldcalev0;
   unsigned short int calev1=0;
@@ -316,8 +345,8 @@
   unsigned short int calev3=0;
   unsigned short int oldcalev3;
   Int_t i = minevent;
-  Float_t headc = 0.;
-  Float_t headco = 0.;
+  UInt_t headc = 0;
+  UInt_t headco = 0;
   Bool_t h1rcheck = false;
   Int_t h1rin = 0;
   Int_t h1rout = 0;
@@ -331,11 +360,11 @@
   Int_t outdexy = 0;
   Int_t indexyc = 0;
   Int_t outdexyc = 0;
-  Int_t obt = 0;
-  Int_t minobt[4];
-  Int_t maxobt[4];
-  Int_t swminobt[4];
-  Int_t swmaxobt[4];
+  UInt_t obt = 0;
+  UInt_t minobt[4];
+  UInt_t maxobt[4];
+  UInt_t swminobt[4];
+  UInt_t swmaxobt[4];
   Bool_t firstobt[4];
   Bool_t swfirstobt[4];
   for ( Int_t i = 0; i<4; i++){
@@ -346,22 +375,24 @@
     firstobt[i] = true;
     swfirstobt[i] = true;
   };
+  Bool_t evdone = false;
+  //
   while ( i < maxevent+1){
+    evdone = false;
     tshit = 0;
     trshit = 0;
     tr->GetEntry(i);
+
     iev2[i] = ce->iev;     
     //
     ph = eh->GetPscuHeader();
     headco = headc;
     headc = ph->GetCounter();
     obt = ph->GetOrbitalTime();
+
     if ( (int)abs((int)(ce->calevnum[0]-ce->calevnum[1])) ) calevn1++; 
     if ( (int)abs((int)(ce->calevnum[2]-ce->calevnum[3])) ) calevn2++;
     if ( (int)abs((int)(ce->calevnum[0]-ce->calevnum[3])) ) calevn3++;
-    //    calevn1 += (int)abs((int)(ce->calevnum[0]-ce->calevnum[1]));
-    //    calevn2 += (int)abs((int)(ce->calevnum[2]-ce->calevnum[3]));
-    //    calevn3 += (int)abs((int)(ce->calevnum[0]-ce->calevnum[3]));
     calev01->Fill(abs((int)(ce->calevnum[0]-ce->calevnum[1])));
     calev23->Fill(abs((int)(ce->calevnum[2]-ce->calevnum[3])));
     calev03->Fill(abs((int)(ce->calevnum[0]-ce->calevnum[3])));
@@ -376,40 +407,48 @@
     calevv2 = (int)ce->calevnum[2];
     oldcalev3 = calev3;
     calev3 = (int)ce->calevnum[3];
-    goto jumpprintout;
-    if ( (headc - headco -1.) == 0. && ((calev0 - oldcalev0 - 1) > 0 || (calev1 - oldcalev1 - 1) > 0 || (calevv2 - oldcalev2 - 1) > 0 || (calev3 - oldcalev3 - 1) > 0) ){
+    //    goto jumpprintout;
+    if ( (headc - headco -1) == 0. && ((calev0 - oldcalev0 - 1) > 0 || (calev1 - oldcalev1 - 1) > 0 || (calevv2 - oldcalev2 - 1) > 0 || (calev3 - oldcalev3 - 1) > 0) ){
       if ( i != minevent ) {
-	printf(" %f 0 Event %i: %i\n",headco,i,oldcalev0);
-	printf(" %f 0 Event %i: %i\n",headc,i+1,calev0);	   
-	printf(" %f 1 Event %i: %i\n",headco,i,oldcalev1);
-	printf(" %f 1 Event %i: %i\n",headc,i+1,calev1);	   
-	printf(" %f 2 Event %i: %i\n",headco,i,oldcalev2);
-	printf(" %f 2 Event %i: %i\n",headc,i+1,calevv2);	   
-	printf(" %f 3 Event %i: %i\n",headco,i,oldcalev3);
-	printf(" %f 3 Event %i: %i\n",headc,i+1,calev3);	
+// 	printf(" %f 0 Event %i: %i\n",headco,i,oldcalev0);
+// 	printf(" %f 0 Event %i: %i\n",headc,i+1,calev0);	   
+// 	printf(" %f 1 Event %i: %i\n",headco,i,oldcalev1);
+// 	printf(" %f 1 Event %i: %i\n",headc,i+1,calev1);	   
+// 	printf(" %f 2 Event %i: %i\n",headco,i,oldcalev2);
+// 	printf(" %f 2 Event %i: %i\n",headc,i+1,calevv2);	   
+// 	printf(" %f 3 Event %i: %i\n",headco,i,oldcalev3);
+// 	printf(" %f 3 Event %i: %i\n",headc,i+1,calev3);	
 	isRAW = 0;
 	if ( ce->stwerr[0] & (1 << 3)) isRAW = 1;	 
 	if ( (calev0 - oldcalev0 - 1) > 0 && !isRAW && ce->perror[0] != 129 && oldcalev0 != 0) {
 	  ver[0][10]++;
+	  sighandler(headc,obt,i,0,nevents,evdone);
+	  printf(" jump in the counter calev is %i oldcalev was %i \n",calev0,oldcalev0);
 	};
 	isRAW = 0;
 	if ( ce->stwerr[0] & (1 << 3)) isRAW = 1;	 
 	if ( (calev1 - oldcalev1 - 1) > 0 && !isRAW && ce->perror[1] != 129 && oldcalev1 != 0 ){
 	  ver[1][10]++;
+	  sighandler(headc,obt,i,1,nevents,evdone);
+	  printf(" jump in the counter calev is %i oldcalev was %i \n",calev1,oldcalev1);
 	};
 	isRAW = 0;
 	if ( ce->stwerr[0] & (1 << 3)) isRAW = 1;	 
 	if ( (calevv2 - oldcalev2 - 1) > 0 && !isRAW && ce->perror[2] != 129 && oldcalev2 != 0 ){
 	  ver[2][10]++;
+	  sighandler(headc,obt,i,2,nevents,evdone);
+	  printf(" jump in the counter calev is %i oldcalev was %i \n",calevv2,oldcalev2);
 	};
 	isRAW = 0;
 	if ( ce->stwerr[0] & (1 << 3)) isRAW = 1;	 
 	if ( (calev3 - oldcalev3 - 1) > 0 && !isRAW && ce->perror[3] != 129 && oldcalev3 != 0 ){
 	  ver[3][10]++;
+	  sighandler(headc,obt,i,3,nevents,evdone);
+	  printf(" jump in the counter calev is %i oldcalev was %i \n",calev3,oldcalev3);
 	};
-      }
+      };
     };
-  jumpprintout:
+    //  jumpprintout:
     memcpy(shit, nullsh, sizeof(nullsh));
     memcpy(rshit, nullsh, sizeof(nullsh));
     for (Int_t l = 0; l < 2; l++ ){
@@ -433,14 +472,13 @@
 	  se = 1;
 	  pl = (ii-1)/2;
 	  pdone = 0;
-	};
-		
-	isCOMP = 0;
-	isFULL = 0;
-	isRAW = 0;
-	if ( ce->stwerr[se] & (1 << 16) ) isCOMP = 1; 
-	if ( ce->stwerr[se] & (1 << 17) ) isFULL = 1; 
-	if ( ce->stwerr[se] & (1 << 3) ) isRAW = 1; 
+	};		
+	isCOMP = false;
+	isFULL = false;
+	isRAW = false;
+	if ( ce->stwerr[se] & (1 << 16) ) isCOMP = true; 
+	if ( ce->stwerr[se] & (1 << 17) ) isFULL = true; 
+	if ( ce->stwerr[se] & (1 << 3) ) isRAW = true; 
 	bl = -1;
 	for (Int_t kk = 0; kk < 96 ; kk++ ){
 	  if ( kk%16 == 0 ){
@@ -482,12 +520,14 @@
 		  if ( obt > maxobt[se] ) maxobt[se] = obt;
 		};
 		lver[se][2]++ ;		
+		sighandler(headc,obt,i,se,nevents,evdone);
+		printf(" latch up (data) in DSP mode alldexy %i planebases %i \n",alldexy,planebases);
 	      };
 	    }
 	    //
 	    if ( !bdone ){ 
 	      Baseline->Fill(ce->base[l][ii][bl]);
-	      if ( ce->base[l][ii][bl] > 2000. && ce->base[l][ii][bl] < 4500. ){
+	      if ( ce->base[l][ii][bl] > 2000. && ce->base[l][ii][bl] < 5500. ){
 		inbase++;
 	      } else { 
 		outbase++;
@@ -525,7 +565,11 @@
 	      indexyc++; 			    
 	    };
 	    //
-	    if (ce->dexyc[l][ii][kk] < 0 ) ver[se][21]++ ;
+	    if (ce->dexyc[l][ii][kk] < 0 ){
+	      ver[se][21]++ ;
+	      sighandler(headc,obt,i,se,nevents,evdone);
+	      printf(" negative ADC values!! ce->dexyc[%i][%i][%i] = %f \n",l,ii,kk,ce->dexyc[l][ii][kk]);
+	    };
 	    //
 	    if ( allbase == 0. || allbase > 32000.) {
 	      fcheck = 0;
@@ -541,7 +585,7 @@
 	  } else {
 	    if ( !pdone ){
 	      if ( ( alldexy2 > 3072000 ||  alldexy2 == 0) && ce->perror[se] == 0 ) {
-		//	      if ( (alldexy2>512000 || alldexy2 == 0) && ce->perror[se] == 0 ) {
+		//	      if ( (alldexy2>512000 || alldexy2 == 0) && ce->perror[se] == 0 ) 
 		pdone = 1;
 		pshit[se][pl]++ ;			 
 		if ( (ce->stwerr[se] & (1 << 4)) == 0 ){
@@ -551,6 +595,8 @@
 		  if ( obt > maxobt[se] ) maxobt[se] = obt;
 		};
 		lver[se][2]++ ;
+		sighandler(headc,obt,i,se,nevents,evdone);
+		printf(" latch up (data) in RAW mode alldexy2 %i \n",alldexy2);
 	      };
 	    };
 	    if ( ce->dexy[l][ii][kk] > 0 && ce->dexy[l][ii][kk] < 32000 ) {
@@ -569,41 +615,123 @@
     };
     cshit = 0;
     for (Int_t k = 0; k < 4 ; k++ ){
-      isCOMP = 0;
-      isFULL = 0;
-      if ( ce->stwerr[se] & (1 << 16) ) isCOMP = 1; 
-      if ( ce->stwerr[se] & (1 << 17) ) isFULL = 1; 
+      isCOMP = false;
+      isFULL = false;
+      cestw=0;// queste righe qui
+      if ( ce->stwerr[k] ) cestw =  ce->stwerr[k] & cmask ;//
+      if ( ce->stwerr[k] & (1 << 16) && !(cestw & (1 << 3)) ) isCOMP = true; 
+      if ( ce->stwerr[k] & (1 << 17) && !(cestw & (1 << 3)) ) isFULL = true; 
+      //
       if ( isCOMP ) ver[k][16]++;
       if ( isFULL ) ver[k][17]++;
       cshit += (int)ce->calstriphit[k];
-      cestw=0;
-      if ( ce->stwerr[k] ) cestw =  ce->stwerr[k] & cmask ;
       if ( cestw ){
-	if ( cestw & (1 << 0) ) ver[k][6]++ ;
-	if ( cestw & (1 << 1) ) ver[k][5]++ ;
-	if ( cestw & (1 << 2) )	ver[k][4]++ ;
-	if ( cestw & (1 << 3) ) ver[k][3]++ ;
+	if ( cestw & (1 << 0) ){
+	  ver[k][6]++ ;
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" CRC error (st. word) \n");
+	};
+	if ( cestw & (1 << 1) ){
+	  ver[k][5]++ ;
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" execution error \n");
+	};
+	if ( cestw & (1 << 2) ){
+	  ver[k][4]++ ;
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" CMD length error \n");
+	};
+	if ( cestw & (1 << 3) ) ver[k][3]++ ; // raw mode
 	if ( cestw & (1 << 4) ){
 	  ver[k][2]++ ;
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" latch up alarm (st. word)\n");
 	  lupstw[k]->Fill(obt);
 	  if ( swfirstobt[k] ) swminobt[k] = obt;
 	  if ( obt > swmaxobt[k] ) swmaxobt[k] = obt;
 	};
-	if ( cestw & (1 << 5) ) ver[k][1]++ ;
-	if ( cestw & (1 << 6) ) ver[k][0]++ ;		
-      };
+	if ( cestw & (1 << 5) ){
+	  ver[k][1]++ ;
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" temperature alarm  \n");
+	};
+	if ( cestw & (1 << 6) ){
+	  ver[k][0]++ ;		
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" DSP ack error \n");
+	};
+      };// controllo se ci sono errori crc se perr=0 ok!
+      if (  ce->stwerr[k] == 0 && ce->perror[k] == 0 ) ver[k][8]++ ; 
       if ( ce->perror[k] != 0. ){
-	if (ce->perror[k] == 128) ver[k][7]++ ; 
-	if (ce->perror[k] == 129) ver[k][8]++ ; 
-	if (ce->perror[k] == 132) ver[k][11]++ ; 
-	if (ce->perror[k] == 133) ver[k][12]++ ; 
-	if (ce->perror[k] == 134) ver[k][13]++ ; 
-	if (ce->perror[k] == 135) ver[k][14]++ ; 
-	if (ce->perror[k] == 136) ver[k][15]++ ; 
-	if (ce->perror[k] == 139) ver[k][18]++ ; 
-	if (ce->perror[k] == 140) ver[k][19]++ ; 
-	if (ce->perror[k] == 141) ver[k][20]++ ; 
-	if (ce->perror[k] == 142) ver[k][22]++ ; 
+	if (ce->perror[k] == 128){
+	  ver[k][7]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" view or command not recognized \n");
+	};
+	if (ce->perror[k] == 129){
+	  ver[k][8]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" missing section \n");
+	};
+	if (ce->perror[k] == 132){
+	  ver[k][11]++ ; 
+	  if ( ver[k][11] < 30 ){
+	    sighandler(headc,obt,i,k,nevents,evdone);
+	    printf(" CRC error (data) \n");
+	  };
+	  if ( ver[k][11] == 30 ){
+	    sighandler(headc,obt,i,k,nevents,evdone);
+	    printf(" CRC error (data) \n");
+	    printf(" WARNING TOO MANY CRC ERRORS ON DATA FOR THIS SECTION \n");
+	    printf(" THEY WILL NOT BE DISPLAYED HERE ANYMORE FOR THIS FILE! \n");
+	  };
+	};
+	if (ce->perror[k] == 133){
+	  ver[k][12]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" length problems in RAW mode \n");
+	};
+	if (ce->perror[k] == 134){
+	  ver[k][13]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" length problems in COMPRESS mode \n");
+	};
+	if (ce->perror[k] == 135){
+	  ver[k][14]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" length problems in FULL mode \n");
+	};
+	if (ce->perror[k] == 136){
+	  ver[k][15]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" acquisition mode problems \n");
+	};
+	if (ce->perror[k] == 139){
+	  ver[k][18]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" problems with coding\n");
+	};
+	if (ce->perror[k] == 140){
+	  ver[k][19]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" pedestal checksum wrong\n");
+	};
+	if (ce->perror[k] == 141){
+	  ver[k][20]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" thresholds checksum wrong\n");
+	};
+	if (ce->perror[k] == 142){
+	  ver[k][22]++ ; 
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" packet length is zero (YODA input error), skipped\n");
+	};
+	//aggiunto questo errore
+	if (ce->perror[k] == 143){
+	  ver[k][23]++ ;
+	  sighandler(headc,obt,i,k,nevents,evdone);
+	  printf(" corrupted packet lenght \n");
+	};
       };
       for (Int_t kk = 0; kk < 11 ; kk++ ){
 	tshit += shit[k][kk];
@@ -613,7 +741,10 @@
     if (isCOMP || isFULL){
       Calstriphit->Fill(cshit);
       //    if ( (cshit > 0 && cshit < 25) || (cshit > 40 && cshit < 80) ){
-      if ( (cshit > 10 && cshit < 100) ){
+      //
+      //cambiata la banda di soglia  ora 5  prima era 10
+      //
+      if ( (cshit > 5 && cshit < 100) ){
 	incshit++;
       } else {
 	outcshit++;
@@ -621,7 +752,8 @@
       //
       if ( tshit>0 ) h1->Fill(tshit);
       //    if ( (tshit > 0 && tshit < 25) || (tshit > 40 && tshit < 80) ){
-      if ( (tshit > 10 && tshit < 100) ){
+      //cambiata la soglia ora 5  prima era 10
+      if ( (tshit > 5 && tshit < 100) ){
 	intshit++;
       } else {
 	outtshit++;
@@ -636,10 +768,11 @@
 	h1rin++;
       };
     };
-    if ( i%1000 == 0 && i > 0 ) printf("%iK\n",i/1000);
+    //    if ( i%1000 == 0 && i > 0 ) printf("%iK\n",i/1000);
     i++;
   };
   printf("\n");
+  printf("\n end processing. \n\n");
   if ( (float)h1rout/((float)h1rin+(float)h1rout) > h1rth ){
     h1rcheck = true;
   };    
@@ -708,7 +841,9 @@
     h1->Draw();
     //
     h1max = h1->GetMaximum()*1.05;
-    Double_t xc[4] = {10.,100.,100.,10.};
+
+    //cambiata banda  prima era 10  ora 5
+    Double_t xc[4] = {5.,100.,100.,5.};
     Double_t yc[4] = {0.,0.,h1max,h1max};
     banda1 = new TPolyLine(4,xc,yc);
     banda1->SetLineColor(5);
@@ -758,7 +893,8 @@
       h1->SetYTitle("Number of events");
       h1->Draw();
       h1max = h1->GetMaximum()*1.05;
-      Double_t xe[4] = {10.,100.,100.,10.};
+      //cambiato banda, prima era 10,ora 5
+      Double_t xe[4] = {5.,100.,100.,5.};
       Double_t ye[4] = {0.,0.,h1max,h1max};
       banda1 = new TPolyLine(4,xe,ye);
       banda1->SetLineColor(5);
@@ -944,7 +1080,7 @@
     Baseline->SetYTitle("Number of events");
     Baseline->Draw();
     h1max = Baseline->GetMaximum()*1.05;
-    Double_t xc[4] = {2000.,4500.,4500.,2000.};
+    Double_t xc[4] = {2000.,5500.,5500.,2000.};
     Double_t yc[4] = {0.,0.,h1max,h1max};
     banda1 = new TPolyLine(4,xc,yc);
     banda1->SetLineColor(5);
@@ -981,7 +1117,9 @@
     Calstriphit->SetYTitle("Number of events");
     Calstriphit->Draw();
     h1max = Calstriphit->GetMaximum()*1.05;
-    Double_t xg[4] = {10.,100.,100.,10.};
+
+    //banda di sogla portata a 5, prima era a 10
+    Double_t xg[4] = {5.,100.,100.,5.};
     Double_t yg[4] = {0.,0.,h1max,h1max};
     banda1 = new TPolyLine(4,xg,yg);
     banda1->SetLineColor(5);
@@ -1188,7 +1326,14 @@
   pad2->Draw();
   pad3->Draw();
   pad4->Draw();
-
+  //
+  Bool_t redevent = true;
+  //
+  if ( ver[0][3]+ver[0][16]+ver[0][17]+ver[0][8] == ver[1][3]+ver[1][16]+ver[1][17]+ver[1][8] && ver[1][3]+ver[1][16]+ver[1][17]+ver[1][8] == ver[2][3]+ver[2][16]+ver[2][17]+ver[2][8] && ver[2][3]+ver[2][16]+ver[2][17]+ver[2][8] == ver[3][3]+ver[3][16]+ver[3][17]+ver[3][8] ){
+    redevent = false; 
+  };
+  if ( redevent ) check = true;
+  //
   for (Int_t i = 0; i < 4; i++){
     if (i == 2) 
       { 
@@ -1216,7 +1361,8 @@
     t->SetTextAlign(12);
     t->DrawLatex(33.,97.,sezione);
     t->SetTextSize(0.05); 
-    for (Int_t j = 0; j < 23; j++){
+    //aggiunto un errore jmax=23
+    for (Int_t j = 0; j < 24; j++){
       if ( ver[i][j] || lver[i][j] ) {
 	t->SetTextColor(50); 
 	if (j == 0) {
@@ -1275,7 +1421,7 @@
 	  errore.str("");
 	  errore << "RAW mode: " << ver[i][j];
 	  errore << " time(s)";
-	  t->SetTextColor(32); 
+	  if ( !redevent ) t->SetTextColor(32); 
 	  t->DrawLatex(2.,90.,errore.str().c_str());
 	}
 	if (j == 4) {
@@ -1318,7 +1464,7 @@
 	  errore << "Calevnum jump: " << ver[i][j];
 	  errore << " time(s)";
 	  t->DrawLatex(2.,42.,errore.str().c_str());
-	  check = true;
+	  if ( ver[i][j] > 3 ) check = true;
 	} 
 	if (j == 11) {
 	  errore.str("");
@@ -1359,14 +1505,14 @@
 	  errore.str("");
 	  errore << "COMPRESS mode: " << ver[i][j];
 	  errore << " time(s)";
-	  t->SetTextColor(32); 
+	  if ( !redevent ) t->SetTextColor(32); 
 	  t->DrawLatex(2.,86.,errore.str().c_str());
 	}
 	if (j == 17) {
 	  errore.str("");
 	  errore << "FULL mode: " << ver[i][j];
 	  errore << " time(s)";
-	  t->SetTextColor(32); 
+	  if ( !redevent ) t->SetTextColor(32); 
 	  t->DrawLatex(2.,82.,errore.str().c_str());
 	}
 	if (j == 18) {
@@ -1404,8 +1550,16 @@
 	  t->DrawLatex(2.,3.,errore.str().c_str());
 	  check = true;
 	};
-      };
-    }; 
+	//aggiunto questo errore
+	if (j == 23) {
+	  errore.str("");
+	  errore << "Corrupted packet lenght " << ver[i][j];
+	  errore << " time(s) ";
+	  t->DrawLatex(2.,2.,errore.str().c_str());
+	  check = true;
+	};
+      }; 
+    };
     t->SetTextColor(50); 
     if ( fulldiff[i] !=0 ) {
       check = true;
@@ -1432,11 +1586,11 @@
   if ( check ){
     t->SetTextColor(50);
     t->DrawLatex(60.,95.," WARNING, CHECK! ");	
-    printf("cdexyc %f cdexy %f ctshit %f cbase %f \n",cdexyc,cdexy,ctshit,cbase);
   } else {
     t->SetTextColor(32);
     t->DrawLatex(60.,95.," OK! ");	
   };
+  printf("cdexyc %f cdexy %f ctshit %f cbase %f \n",cdexyc,cdexy,ctshit,cbase);
   //
   //  const string fil = (const char*)filename;
   //  Int_t posiz = fil.find("dw_");
@@ -1519,4 +1673,3 @@
     };
   };
 }
-