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

Timestamp:

06/01/05 18:13:13 (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

@@ -197 +197 @@
             return kFALSE;
         }
-        const Double_t m3l = fHilExt->GetM3Long()*TMath::Sign(1.0f, hsrc.GetCosDeltaAlpha())*fMm2Deg;
+        const Double_t m3l = fHilExt->GetM3Long()*TMath::Sign(fMm2Deg, hsrc.GetCosDeltaAlpha());
 
         gweight = m3l>fM3lCut ? 1 : 0;
@@ -217 +217 @@
     }
 
-    // Workaround: Number-of-entries
-    if (gweight>0.25)
-        fHist.Fill(pos1.X(), pos1.Y(), 0.0, w*gweight);
-    if (gweight<0.75)
-        fHist.Fill(pos2.X(), pos2.Y(), 0.0, w*(1-gweight));
+    fHist.Fill(pos1.X(), pos1.Y(), 0.0, w*gweight);
+    fHist.Fill(pos2.X(), pos2.Y(), 0.0, w*(1-gweight));
 
     return kTRUE;
@@ -273 +270 @@
         const Double_t x0 = h1->GetXaxis()->GetBinCenter(ix);
         const Double_t y0 = h1->GetYaxis()->GetBinCenter(iy);
+        const Double_t w0 = h1->GetXaxis()->GetBinWidth(1);
 
         for (int x=0; x<h1->GetNbinsX(); x++)
@@ -294 +292 @@
         h2->Fit(&func, "IMQ", "", 0, 1.0);
 
-        const Double_t r0 = 2*func.GetParameter(2);
-        const Double_t e  = func.Integral(0, r0)/h1->GetBinWidth(1);
-        func.SetParameter(0, 0);
-        const Double_t b  = func.Integral(0, r0)/h1->GetBinWidth(1);
+        // No wintegrate the function f(x) per Delta Area
+        // which is f(x)/(pi*delta r*(2*r+delta r))
+        TF1 func2("fcn2", Form("(gaus + [3]*x*x + [4])/(2*x+%.5f)", w0));
+        for (int i=0; i<5; i++)
+            func2.SetParameter(i, func.GetParameter(i));
+
+        const Double_t r0 = 2*func2.GetParameter(2);
+        const Double_t e  = func2.Integral(0, r0)/(w0*TMath::Pi());
+        func2.SetParameter(0, 0);
+        const Double_t b  = func2.Integral(0, r0)/(w0*TMath::Pi());
         const Double_t s  = MMath::SignificanceLiMa(e, b);