--- DarthVader/CalorimeterLevel2/src/CaloLevel0.cpp	2010/01/29 05:49:24	1.30
+++ DarthVader/CalorimeterLevel2/src/CaloLevel0.cpp	2011/11/29 13:50:11	1.32
@@ -158,6 +158,7 @@
   this->InitDo(dbc,hs,sgnl,l0tree,isdeb,isverb);
   dbc->Close();
   delete dbc;
+  dbc = 0;
 }
 
 
@@ -277,6 +278,7 @@
   Int_t sig = this->ChkParamDo(dbc,runheader,mechal);
   dbc->Close();
   delete dbc;
+  dbc = 0;
   return(sig);
 }
 
@@ -602,6 +604,7 @@
   Int_t sig = CalcCrossTalkCorrDo(dbc,runheader,usetable);
   dbc->Close();
   delete dbc;
+  dbc = 0;
   //
   return(sig);
   //  
@@ -621,6 +624,7 @@
   Int_t sig = CalcCrossTalkCorrDo(dbc,runheader,true);
   dbc->Close();
   delete dbc;
+  dbc = 0;
   //
   return(sig);
   //
@@ -1160,17 +1164,25 @@
 	if ( strip6s < 4 && base[l][m][pre] > (-0.015*qp+sbase[l][m][pre]) && sbase[l][m][pre] > 0. ){
 	  if ( debug ) printf(" Suspicious calculated baseline: base %f sbase-0.02*qp %f strip6s %i \n",base[l][m][pre],(-qp*0.02+sbase[l][m][pre]),(Int_t)strip6s);
 	  base[l][m][pre] = 31000.;	
+	  nst = 0; // 9RED BUG
+	  qp = 0.; // 9RED BUG
 	  for (Int_t e = pre*16; e < (pre+1)*16 ; e++){
 	    dexyc[l][m][e] = dexy[l][m][e];
 	  };	  
 	};
       } else {
+	if ( debug ) printf(" reset baseline here if ! ( (strip6s >=2 && process == 2) || (strip6s >= 9 and process > 2) ) \n");
 	base[l][m][pre] = 31000.;
+	nst = 0; // 9RED BUG
+	qp = 0.; // 9RED BUG
 	for (Int_t e = pre*16; e < (pre+1)*16 ; e++){
 	  dexyc[l][m][e] = dexy[l][m][e];
 	};
       };
     } else {
+      if ( debug ) printf(" reset baseline here if no minimum find\n");
+      nst = 0; // 9RED BUG
+      qp = 0.; // 9RED BUG
       process += 2;
       base[l][m][pre] = 31000.;
       for (Int_t e = pre*16; e < (pre+1)*16 ; e++){
@@ -1178,6 +1190,7 @@
       };
     };
   };
+  if ( debug ) printf(" Baseline calculation: baseline for view %i plane %i pre %i is %f nst %i qp %f \n",l,m,pre,base[l][m][pre],nst,qp);
 }
 
 Int_t CaloLevel0::Calibrate(Int_t ei){
@@ -1388,6 +1401,7 @@
 	    // no suitable new baseline, use old ones!
 	    //
 	    if ( !done ){
+	      if ( debug ) printf(" l %i m %i pre %i ip[i] %i base %f base ip[i] %f sbase %f \n",l,m,pre,ip[i],base[l][m][pre],base[l][m][ip[i]],sbase[l][m][pre]);
 	      if ( (base[l][m][pre] == 31000. || base[l][m][pre] == 0.) ){
 		ck[l][m][pre] = 1;
 		if (pre%2 == 0) {
@@ -1454,6 +1468,7 @@
 		//
 		// OK, now in estrip we have the energy deposit in MIP of all the strips for this event (at the end of loops of course)
 		//
+		if ( debug && l==0 && (m==17 || m==18) ) printf(" view %i plane %i strip %i ener %f calped %f base0 %f base1 %f base2 %f mip %f ENERGIA %f \n",l,m,n,ener,calped[l][m][n],base0,base1,base2,mip[l][m][n],clevel1->estrip[n][m][l]);
 		if ( clevel1->estrip[n][m][l] > 0. ) qpre[l][m][pre] += clevel1->estrip[n][m][l];
 		//
 		//
@@ -1498,7 +1513,7 @@
 	//
 	// Cross-talk corrections  
 	//
-	if ( crosst ){
+	if ( crosst ){  
 	  //
 	  // energy on silicon ladders
 	  //
@@ -1515,7 +1530,9 @@
 	    if ( noselfct ){
 	      for (Int_t j = ladder*32 ; j < (ladder+1)*32 ; j++){
 		ipre = j/16 ;	
+		if ( debug ) printf(" CT STEP1 %i %i %i estrip %f ctprecor %f \n",j,m,l,clevel1->estrip[j][m][l],ctprecor[l][m][ipre]);
 		if ( clevel1->estrip[j][m][l] != 0. ) clevel1->estrip[j][m][l] -=  clevel1->estrip[j][m][l] * ctprecor[l][m][ipre];
+		if ( debug ) printf(" CT STEP2 %i %i %i estrip %f ctprecor %f \n",j,m,l,clevel1->estrip[j][m][l],ctprecor[l][m][ipre]);
 	      };
 	    };
 	    //
@@ -1537,12 +1554,14 @@
               //
               if ( debug ) printf(" CK1 Limit for while: 0.07 \n");
 	      for (Int_t j = ipre*16 ; j < (ipre+1)*16 ; j++){
+		if ( debug ) printf(" CT STEP3 %i %i %i estrip %f ctprecor %f \n",j,m,l,clevel1->estrip[j][m][l],ctprecor[l][m][ipp]);
 		if ( !ctground ){
 		  if ( clevel1->estrip[j][m][l] != 0. ) clevel1->estrip[j][m][l] += - qpre[l][m][ipp] * ctprecor[l][m][ipp];
 		} else {
 		  if ( clevel1->estrip[j][m][l] != 0. ) clevel1->estrip[j][m][l] += - qpre[l][m][ipp] * 0.00478;
 		};
 		if ( clevel1->estrip[j][m][l] > 0. ) nqpre += clevel1->estrip[j][m][l] ;
+		if ( debug ) printf(" CT STEP4 %i %i %i estrip %f ctprecor %f \n",j,m,l,clevel1->estrip[j][m][l],ctprecor[l][m][ipp]);
 	      };
 	      qpre[l][m][ipre] = nqpre;
               nqpre = 0.;
@@ -1553,6 +1572,7 @@
 	      while ( it < 10 && deltaqpre > 0.07 ){
 		nqpre = 0.;
 		for (Int_t j = ipre*16 ; j < (ipre+1)*16 ; j++){
+		  if ( debug ) printf(" CT STEP5 %i %i %i estrip %f ctprecor %f \n",j,m,l,clevel1->estrip[j][m][l],ctprecor[l][m][ipre]);
 		  if ( !ctground ){
 		    if ( debug ) printf(" CK1 pre correction: iteration %i deltaqpre %f  ctprecor %f  TOTAL CORRECTION %f \n",it,deltaqpre,ctprecor[l][m][ipre],deltaqpre * ctprecor[l][m][ipre]);
 		    if ( clevel1->estrip[j][m][l] != 0. ) clevel1->estrip[j][m][l] += deltaqpre * ctprecor[l][m][ipre];
@@ -1560,6 +1580,7 @@
 		    if ( clevel1->estrip[j][m][l] != 0. ) clevel1->estrip[j][m][l] += deltaqpre * 0.00478;
 		  };
  	          if ( clevel1->estrip[j][m][l] > 0. ) nqpre += clevel1->estrip[j][m][l] ; 
+		  if ( debug ) printf(" CT STEP6 %i %i %i estrip %f ctprecor %f \n",j,m,l,clevel1->estrip[j][m][l],ctprecor[l][m][ipre]);
 		};
 		if ( ctground ) it = 100;
                 it++;
@@ -1666,6 +1687,7 @@
 		nqpre[0] = 0.;
 		nqpre[1] = 0.;
  	        for (Int_t j = ladder*32 ; j < (ladder+1)*32 ; j++){
+		  if ( debug ) printf(" CT STEP6 %i %i %i estrip %f ctprecor %f \n",j,m,l,clevel1->estrip[j][m][l],ctsicor[l][m][si2]);
                   ipre = 0;
 		  if ( j > (ladder*32)+15 ) ipre = 1;
                   jpre = j/16 ;
@@ -1681,13 +1703,22 @@
 		  };
  	          if ( clevel1->estrip[j][m][l] > 0. ) nqsi += clevel1->estrip[j][m][l] ;
  	          if ( clevel1->estrip[j][m][l] > 0. ) nqpre[ipre] += clevel1->estrip[j][m][l] ; 
+		  if ( debug ) printf(" CT STEP7 %i %i %i estrip %f ctprecor %f \n",j,m,l,clevel1->estrip[j][m][l],ctsicor[l][m][si2]);
                 };	
 	        if ( ctground ) it = 100;
 		deltaqsi = nqsi-snqsi;
-		snqsi = nqsi;
-                it++;
 		deltaqpre[0] = nqpre[0] - qpre[l][m][pre-1];
 		deltaqpre[1] = nqpre[1] - qpre[l][m][pre];
+		//
+		// Check for divergence and stop if it happens! [9RED bug noticed with plane 18X]
+		//
+		if ( deltaqpre[0] > qpre[l][m][pre-1] || deltaqpre[1] > qpre[l][m][pre] || deltaqsi >snqsi ){
+		  if ( debug ) printf(" WARNING!! DIVERGING CORRECTION EXIT IMMEDIATLY FROM THE LOOP!! dqpre0 %f qpre0 %f // dqpre1 %f qpre1 %f // dqsi %f qsi %f \n",deltaqpre[0],qpre[l][m][pre-1],deltaqpre[1],qpre[l][m][pre],deltaqsi,snqsi);
+		  it = 1000;
+		};
+		//
+		snqsi = nqsi;
+                it++;
                 if ( debug ) printf(" BEFORE: qpre 0 %f qpre 1 %f \n",qpre[l][m][pre-1],qpre[l][m][pre]);   
                 qpre[l][m][pre-1] = nqpre[0];
                 qpre[l][m][pre] = nqpre[1];	
@@ -1723,7 +1754,7 @@
 	jjj++;
 	j4++;
 	if ( j < 96 ) ene[j] = clevel1->estrip[j][m][l];
-	if ( crosst ){
+	if ( crosst ){  
 	  //
 	  // "Real" crosstalk effect on the neighbour strips respect to the one which have seen the energy deposit
 	  //
@@ -1773,8 +1804,12 @@
 	    };
 	    if ( debug ) printf(" HIT STRIP: view %i plane %i strip %i energy: %f \n",l,m,j4,clevel1->estrip[j4][m][l]);
 	    //	    if ( dexyc[l][m][j4] == 32767. ){
-	    if ( dexyc[l][m][j4] > 32000. ){
-	      savel1 += 5000.; 
+	    if ( dexyc[l][m][j4] > 32000. || savel1 > 5000.){ // CaloLevel1 bug with plane 18X [9RED 14/04/2010]
+	      if ( savel1 > 5000 ){
+		if ( debug ) printf(" Absurd plane 18X energy... resetting value to 1100 MIP \n");
+		savel1 = 1100.;            // CaloLevel1 bug with plane 18x [9RED 14/04/2010]
+	      };
+	      savel1 += 5000.; 	      
 	      clevel2->nsatstrip += 1.;
 	    };
 	    //
@@ -1817,6 +1852,7 @@
 	    } else {
 	      svstrip[istrip] = -(fbi*1000000000 + plo*10000000 + j4*100000 + cle);
 	    };
+	    if ( debug ) printf(" svstrip[%i] = %i  fbi %i plo %i j4 %i cle %i \n",istrip,svstrip[istrip],fbi,plo,j4,cle);
 	    //
 	    istrip++;
 	  };