Index: trunk/MagicSoft/Mars/mhist/MHHadroness.cc
===================================================================
--- trunk/MagicSoft/Mars/mhist/MHHadroness.cc	(revision 1353)
+++ trunk/MagicSoft/Mars/mhist/MHHadroness.cc	(revision 1354)
@@ -266,4 +266,5 @@
     *fLog << "Hadroness histograms:" << endl;
     *fLog << "---------------------" << endl;
+    *fLog << "Used " << fGhness->GetEntries() << " Gammas and " << fPhness->GetEntries() << " Hadrons." << endl;
     *fLog << "Acc Gammas at  1% Hadron-acc: " << Form("%3.0f", GetGammaAcceptance(0.01)*100) << "%" << endl;
     *fLog << "Acc Gammas at  2% Hadron-acc: " << Form("%3.0f", GetGammaAcceptance(0.02)*100) << "%" << endl;
Index: trunk/MagicSoft/Mars/mhist/MHMatrix.cc
===================================================================
--- trunk/MagicSoft/Mars/mhist/MHMatrix.cc	(revision 1353)
+++ trunk/MagicSoft/Mars/mhist/MHMatrix.cc	(revision 1354)
@@ -45,6 +45,4 @@
 #include "MHMatrix.h"
 
-#include <math.h>     // fabs on Alpha
-
 #include <TList.h>
 #include <TArrayD.h>
@@ -353,5 +351,5 @@
         d2 *= m;
 
-        dists[i] = fabs(d2*d);
+        dists[i] = d2*d; // square of distance
     }
 
@@ -365,5 +363,5 @@
         res += dists[idx[i]];
 
-    return res/n;
+    return sqrt(res)/n;
 }