Index: trunk/MagicSoft/Mars/manalysis/MCT1SupercutsCalc.cc
===================================================================
--- trunk/MagicSoft/Mars/manalysis/MCT1SupercutsCalc.cc	(revision 2274)
+++ trunk/MagicSoft/Mars/manalysis/MCT1SupercutsCalc.cc	(revision 2282)
@@ -158,37 +158,37 @@
 
     TArrayD lup(ncutpar,     par + k0);
-    SetLengthUp(lup);
+    fLengthUp = lup;
     k0 += ncutpar;
 
     TArrayD wup(ncutpar,     par + k0);
-    SetWidthUp(wup);
+    fWidthUp = wup;
     k0 += ncutpar;
 
     TArrayD dup(ncutpar,     par + k0);
-    SetDistUp(dup);
+    fDistUp = dup;
     k0 += ncutpar;
 
     TArrayD llo(ncutpar,     par + k0);
-    SetLengthLo(llo);
+    fLengthLo = llo;
     k0 += ncutpar;
 
     TArrayD wlo(ncutpar,     par + k0);
-    SetWidthLo(wlo);
+    fWidthLo = wlo;
     k0 += ncutpar;
 
     TArrayD dlo(ncutpar,     par + k0);
-    SetDistLo(dlo);
+    fDistLo = dlo;
     k0 += ncutpar;
 
     TArrayD aup(ncutpar,     par + k0);
-    SetAsymUp(aup);
+    fAsymUp = aup;
     k0 += ncutpar;
 
     TArrayD alo(ncutpar,     par + k0);
-    SetAsymLo(alo);
+    fAsymLo = alo;
     k0 += ncutpar;
 
     TArrayD alphaup(ncutpar, par + k0);
-    SetAlphaUp(alphaup);
+    fAlphaUp = alphaup;
 }
 
@@ -277,4 +277,6 @@
 
     InitParams();
+
+    fMatrix = NULL;
 }
 
@@ -423,6 +425,7 @@
 
     const Double_t dist2   = meanx*meanx + meany*meany;
-    const Double_t dd2     = dist2*fMm2Deg;
     const Double_t dist    = sqrt(dist2) * fMm2Deg;
+    const Double_t dd2     = dist*dist;
+
 
     const Double_t dmls    = log(size) - kNomLogSize;