Changeset 7287 for trunk/MagicSoft/Mars/mhflux/MHDisp.cc

Timestamp:

08/16/05 15:38:01 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MHDisp.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/mhflux/MHDisp.cc

@@ -389 +389 @@
     const TAxis &axey = *h2.GetYaxis();
 
-    const Double_t rmax = fWobble ? axex.GetXmax()-0.4 : axex.GetXmax();
+    const Double_t rmax = (fWobble ? axex.GetXmax()-0.4 : axex.GetXmax()) - axex.GetBinWidth(1);
 
     for (int x=1; x<=axex.GetNbins(); x++)
@@ -571 +571 @@
         //Profile(*h2, *h1, x0, y0);
 
-        // Replace with MAlphaFitter?
-        TF1 func("fcn", "gaus + [3]*x*x + [4]");
-        func.SetLineWidth(1);
-        func.SetLineColor(kBlue);
-
-        func.SetParLimits(2, h2->GetBinWidth(1), 1.0);
-
-        func.SetParameter(0, h2->GetBinContent(1));
-        func.FixParameter(1, 0);
-        func.SetParameter(2, 0.12);
-        if (fHistOff)
-            func.FixParameter(3, 0);
-        func.SetParameter(4, 0);//h2->GetBinContent(20));
-        h2->Fit(&func, "IQ", "", 0, 1.0);
-
-        h2->SetTitle(Form("P=(%.2f\\circ/%.2f\\circ) \\omega=%.2f\\circ f=%.2f", x0, y0, func.GetParameter(2), TMath::Abs(scale)));
+        if (h2->GetEntries()>0)
+        {
+            // Replace with MAlphaFitter?
+            TF1 func("fcn", "gaus + [3]*x*x + [4]");
+            func.SetLineWidth(1);
+            func.SetLineColor(kBlue);
+
+            func.SetParLimits(2, h2->GetBinWidth(1), 1.0);
+
+            func.SetParameter(0, h2->GetBinContent(1));
+            func.FixParameter(1, 0);
+            func.SetParameter(2, 0.12);
+            if (fHistOff)
+                func.FixParameter(3, 0);
+            func.SetParameter(4, 0);//h2->GetBinContent(20));
+            h2->Fit(&func, "IQ", "", 0, 1.0);
+
+            h2->SetTitle(Form("P=(%.2f\\circ/%.2f\\circ) \\omega=%.2f\\circ f=%.2f", x0, y0, func.GetParameter(2), TMath::Abs(scale)));
+        }
     }
 
@@ -623 +626 @@
 
     if (fHistOff)
-        GetCatalog()->Draw("white mirror same *");
+        GetCatalog()->Draw("mirror same *");
 
     // ----- Pad number 1 -----
@@ -643 +646 @@
 
     if (fHistOff)
-        GetCatalog()->Draw("white mirror same *");
+        GetCatalog()->Draw("mirror same *");
 
     // ----- Pad number 2 -----