Changeset 7392 for trunk/MagicSoft/Mars/mhflux

Timestamp:

11/10/05 12:03:37 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

• MAlphaFitter.cc (modified) (5 diffs)
• MAlphaFitter.h (modified) (1 diff)

trunk/MagicSoft/Mars/mhflux/MAlphaFitter.cc

@@ -116 +116 @@
     Double_t bgmax =fBgMax;
 
-    //*fLog << inf << "Fit: " << sigmax << " " << fSigInt << "  " << bgmin << " " << bgmax << endl;
-
-    //TF1 fFunc("", Form("gaus(0) + pol%d(3)", fPolynomOrder), 0, 90);
-
-    //fFunc->SetParameters(fCoefficients.GetArray());
-
-    fFunc->FixParameter(1, 0);
-    if (fPolynomOrder!=1)
-        fFunc->FixParameter(4, 0);
-    fFunc->SetParLimits(2, 0, 90);
-    fFunc->SetParLimits(3, -1, 1);
-
     const Double_t alpha0 = h.GetBinContent(1);
     const Double_t alphaw = h.GetXaxis()->GetBinWidth(1);
@@ -137 +125 @@
     // First fit a polynom in the off region
     fFunc->FixParameter(0, 0);
+    fFunc->FixParameter(1, 0);
     fFunc->FixParameter(2, 1);
     fFunc->ReleaseParameter(3);
+    if (fPolynomOrder!=1)
+        fFunc->FixParameter(4, 0);
 
     for (int i=5; i<fFunc->GetNpar(); i++)
@@ -145 +136 @@
         else
             fFunc->SetParameter(i, 0);
+
+    fFunc->SetParLimits(2, 0, 90);
+    //    fFunc->SetParLimits(3, -1, 1);
 
     // options : N  do not store the function, do not draw
@@ -171 +165 @@
     //func.ReleaseParameter(1);
     fFunc->ReleaseParameter(2);
-    fFunc->FixParameter(3, fFunc->GetParameter(3));
-    fFunc->FixParameter(4, fFunc->GetParameter(4));
-    for (int i=5; i<fFunc->GetNpar(); i++)
+    for (int i=3; i<fFunc->GetNpar(); i++)
         fFunc->FixParameter(i, fFunc->GetParameter(i));
 
     // Do not allow signals smaller than the background
-    const Double_t A  = alpha0-fFunc->GetParameter(3);
+    const Double_t s  = fSignalFunc==kGauss ? fFunc->GetParameter(3) : TMath::Exp(fFunc->GetParameter(3));
+    const Double_t A  = alpha0-s;
     const Double_t dA = TMath::Abs(A);
     fFunc->SetParLimits(0, -dA*4, dA*4);
@@ -334 +327 @@
     switch (fSignalFunc)
     {
-    case kGauss:       *fLog << "gauss(x)";        break;
-    case kThetaSq:     *fLog << "gauss(sqrt(x))";  break;
+    case kGauss:       *fLog << "gauss(x)/pol" << fPolynomOrder; break;
+    case kThetaSq:     *fLog << "gauss(sqrt(x))/expo";  break;
     }
     *fLog << endl;

trunk/MagicSoft/Mars/mhflux/MAlphaFitter.h

@@ -136 +136 @@
             break;
         case kThetaSq:
-            fFunc=new TF1 ("", Form("[0]*exp(-0.5*((sqrt(x)-[1])/[2])^2) + pol%d(3)", fPolynomOrder));
+//            if (fPolynomOrder==0)
+//                fFunc=new TF1("", "[0]*exp(-0.5*((sqrt(x)-[1])/[2])^2) + pol0(3)");
+//            else
+            //            {
+            if (fPolynomOrder>0)
+                fPolynomOrder = 1;
+            fFunc=new TF1("", "[0]*exp(-0.5*((sqrt(x)-[1])/[2])^2) + expo(3)");
+//            }
             break;
         }