Index: trunk/MagicSoft/Mars/mcalib/MCalibrationPedCam.cc
===================================================================
--- trunk/MagicSoft/Mars/mcalib/MCalibrationPedCam.cc	(revision 5484)
+++ trunk/MagicSoft/Mars/mcalib/MCalibrationPedCam.cc	(revision 5485)
@@ -133,5 +133,5 @@
       break;
     case 5:
-      val = (*this)[idx].GetLoGainMean() > 1. ? 0. : (*this)[idx].GetLoGainMean();
+      val = (*this)[idx].GetLoGainMean();
       break;
     case 6:
Index: trunk/MagicSoft/Mars/mhcalib/MHPedestalCam.cc
===================================================================
--- trunk/MagicSoft/Mars/mhcalib/MHPedestalCam.cc	(revision 5484)
+++ trunk/MagicSoft/Mars/mhcalib/MHPedestalCam.cc	(revision 5485)
@@ -485,5 +485,6 @@
       //
       pix.SetLoGainMean    ( fSum[i] / fNumEvents  );
-      pix.SetLoGainSigma   ( TMath::Sqrt( (fSum2[i] - fSum[i]*fSum[i]/fNumEvents) / (fNumEvents-1) ));
+      const Double_t diff = fSum2[i] - fSum[i]*fSum[i]/fNumEvents;
+      pix.SetLoGainSigma   ( diff > 0. ? TMath::Sqrt( diff / (fNumEvents-1) ) : 0.);
       pix.SetLoGainMeanVar ( pix.GetLoGainSigma() * pix.GetLoGainSigma() / fNumEvents  );
       pix.SetLoGainSigmaVar( pix.GetLoGainMeanVar() / 4. );
@@ -532,6 +533,8 @@
       if (aevts <= 1)
         continue;
+
       pix.SetLoGainMean ( fAreaSum[j] / aevts  );
-      pix.SetLoGainSigma( TMath::Sqrt( (fAreaSum2[j] - fAreaSum[j]*fAreaSum[j]/aevts) / (aevts-1) ));
+      const Double_t diff = fAreaSum2[j] - fAreaSum[j]*fAreaSum[j]/aevts ;
+      pix.SetLoGainSigma( diff > 0. ? TMath::Sqrt( diff / (aevts-1) ) : 0.);
     }
   
@@ -574,6 +577,8 @@
       if (sevts <= 1)
         continue;
-      pix.SetLoGainMean ( fAreaSum[j] / sevts  );
-      pix.SetLoGainSigma( TMath::Sqrt( (fAreaSum2[j] - fAreaSum[j]*fAreaSum[j]/sevts) / (sevts-1) ));
+
+      pix.SetLoGainMean ( fSectorSum[j] / sevts  );
+      const Double_t diff = fSectorSum2[j] - fSectorSum[j]*fSectorSum[j]/sevts ;
+      pix.SetLoGainSigma( diff > 0. ? TMath::Sqrt( diff / (sevts-1) ) : 0.);
     }
 }