--- tracker/flight/TrkNuclei/src/TrkNuclei.cpp	2009/01/20 10:29:57	1.1
+++ tracker/flight/TrkNuclei/src/TrkNuclei.cpp	2009/03/05 08:52:29	1.2
@@ -63,6 +63,8 @@
     
     Float_t charge  = 1000.;
     Float_t def     = 1./((Float_t)x[0]);
+// WM 05-03-2009 calibration valid only for def>0
+            def = fabs(def);
     Float_t dedxtrk = (Float_t) x[1]; dedxtrk = dedxtrk/1.3;
     
      Float_t ym = dedxtrk;
@@ -72,7 +74,9 @@
 	 yl = B0[jj] + B1[jj]*xb + B2[jj]*xb*xb + B3[jj]*xb*xb*xb + B4[jj]*xb*xb*xb*xb;
 	 yh = B0[jj+1] + B1[jj+1]*xb + B2[jj+1]*xb*xb + B3[jj+1]*xb*xb*xb + B4[jj+1]*xb*xb*xb*xb;
 	 
-	 if ((yl<ym)&&(ym<yh)){
+//	 if ((yl<ym)&&(ym<yh)){
+//WM 05-03-2009  dEdx=0 => charge=0
+	 if ((yl<=ym)&&(ym<yh)){
 	     m = (C0[jj+1]*C0[jj+1] - C0[jj]*C0[jj]) / (yh - yl);
 	     b = (C0[jj]*C0[jj]) - m*yl;
 	     chelp = m*ym + b;