Changeset 8879 for trunk/MagicSoft/Mars/mhflux/MHAlpha.cc

Timestamp:

03/03/08 14:59:12 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MHAlpha.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/mhflux/MHAlpha.cc

@@ -884 +884 @@
 }
 
-void MHAlpha::DrawNicePlot(Bool_t newc, const char *title, const char *watermark)
+void MHAlpha::DrawNicePlot(Bool_t newc, const char *title, const char *watermark, Int_t binlo, Int_t binhi)
 {
     if (!newc && !fDisplay)
@@ -909 +909 @@
 
     // Get on-data
-    TH1D *hon = (TH1D*)fHist.ProjectionZ("Proj", -1, 9999, -1, 9999, "E");
+    TH1D *hon = (TH1D*)fHist.ProjectionZ("Proj", -1, 9999, binlo, binhi, "E");
     hon->SetDirectory(NULL);
     hon->SetBit(kCanDelete);
@@ -921 +921 @@
     if (fOffData)
     {
-        hoff = (TH1D*)fOffData->ProjectionZ("ProjOff", -1, 9999, -1, 9999, "E");
+        hoff = (TH1D*)fOffData->ProjectionZ("ProjOff", -1, 9999, binlo, binhi, "E");
         hoff->SetDirectory(NULL);
         hoff->SetBit(kCanDelete);
@@ -947 +947 @@
 
     // scale off-data
-    fFit.Scale(*hoff, *hon);
+
+    MAlphaFitter fit(fFit);
+    fit.ScaleAndFit(*hon, hoff);
 
     hon->SetMinimum(0);
@@ -971 +973 @@
     text.SetTextSize(0.07);
     text.SetTextAngle(2.5);
-    if (watermark)
-        text.DrawLatex(0.45, 0.2, watermark);
+
+    TString wm(watermark);
+    if (binlo>=1 || binhi<hon->GetNbinsX())
+    {
+        wm += wm.IsNull() ? "(" : " (";
+        if (binlo>=1)
+            wm += Form("%.1fGeV", fHist.GetYaxis()->GetBinLowEdge(binlo));
+        wm += "-";
+        if (binhi<hon->GetNbinsX())
+            wm += Form("%.1fGeV", fHist.GetYaxis()->GetBinLowEdge(binhi+1));
+        wm += ")";
+    }
+    if (!wm.IsNull())
+        text.DrawLatex(0.45, 0.2, wm);
     //enum { kTextNDC = BIT(14) };
 
@@ -995 +1009 @@
     //        fFit.GetSignificance(), fFit.GetScaleFactor(), fFit.GetGausSigma());
     sprintf(txt, "Significance %.1f\\sigma,  off-scale %.2f",
-            fFit.GetSignificance(), fFit.GetScaleFactor());
+            fit.GetSignificance(), fit.GetScaleFactor());
     ptxt = pave->AddText(txt);
     ptxt->SetTextAlign(23);
 
     sprintf(txt, "%.1f excess events,  %.1f background events",
-            fFit.GetEventsExcess(), fFit.GetEventsBackground());
+            fit.GetEventsExcess(), fit.GetEventsBackground());
     ptxt = pave->AddText(txt);
     ptxt->SetTextAlign(23);