Index: trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc
===================================================================
--- trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc	(revision 7693)
+++ trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc	(revision 7694)
@@ -138,7 +138,4 @@
             fFunc->SetParameter(i, 0);
 
-    fFunc->SetParLimits(2, 0, 90);
-    //    fFunc->SetParLimits(3, -1, 1);
-
     // options : N  do not store the function, do not draw
     //           I  use integral of function in bin rather than value at bin center
@@ -163,6 +160,6 @@
     // ------------------------------------
 
-    fFunc->ReleaseParameter(0);
-    //func.ReleaseParameter(1);
+    fFunc->ReleaseParameter(0);  // It is also released by SetParLimits later on
+    //func.ReleaseParameter(1);  // It is also released by SetParLimits later on
     fFunc->ReleaseParameter(2);
     for (int i=3; i<fFunc->GetNpar(); i++)
@@ -173,6 +170,6 @@
     const Double_t A  = alpha0-s;
     const Double_t dA = TMath::Abs(A);
-    fFunc->SetParLimits(0, -dA*4, dA*4);
-    fFunc->SetParLimits(2, 0, 90);
+    fFunc->SetParLimits(0, -dA*4, dA*4);  // SetParLimits also releases the parameter
+    fFunc->SetParLimits(2, 0, 90);        // SetParLimits also releases the parameter
 
     // Now fit a gaus in the on region on top of the polynom