Changeset 5143 for trunk/MagicSoft/Mars/mhflux/MHEffectiveOnTime.cc

Timestamp:

09/28/04 12:25:50 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhflux/MHEffectiveOnTime.cc (modified) (11 diffs)

trunk/MagicSoft/Mars/mhflux/MHEffectiveOnTime.cc

@@ -304 +304 @@
     // within the range (yq[0], yq[1]) there must be no empty bin;
     // choose pedestrian approach as long as GetQuantiles is not available
-    Double_t xq[2] = { 0.05, 0.95 };
+    Double_t xq[2] = { 0.6, 0.99 };
     Double_t yq[2];
     h->GetQuantiles(2, yq, xq);
@@ -321 +321 @@
     //func.SetParNames("lambda", "N0", "del");
 
-    func.SetParameter(0, 100);       // Hz
+    func.SetParameter(0, 200);       // Hz
     func.SetParameter(1, Nm);
     func.FixParameter(2, Nmdel/Nm);
@@ -335 +335 @@
 
     // was fit successful ?
-    const Bool_t ok = NDF>0 && chi2<2.5*NDF;
+    const Bool_t ok = NDF>0 && chi2<3*NDF;
 
     if (paint)
@@ -343 +343 @@
         func.Paint("same");
     }
-
-    if (!ok)
-        return kFALSE;
 
     const Double_t lambda = func.GetParameter(0);
@@ -385 +382 @@
     res[5] = NDF;
 
+    // Chi2
+    res[6] = chi2;
+
     // Rdead (from fit) is the fraction from real time lost by the dead time
     //fHRdead.SetBinContent(i, Rdead);
     //fHRdead.SetBinError  (i,dRdead);
 
-    return kTRUE;
+    return ok;
 }
 
@@ -415 +415 @@
         h = fH2DeltaT.ProjectionX(name, i, i, "E");
 
-        Double_t res[6];
+        Double_t res[7];
         if (!FitH(h, res))
             continue;
@@ -457 +457 @@
     // Fit histogram
     //
-    Double_t res[6];
+    Double_t res[7];
     if (!FitH(&fH1DeltaT, res))
         return;
@@ -551 +551 @@
 
     //
-    // If we reached the event number limit for the time-bins fit the histogram
-    //
-    if (fH1DeltaT.GetEntries()>=fNumEvents)
+    // If we reached the event number limit for the time-bins fit the
+    // histogram - if it fails try again when 1.6% more events available
+    //
+    const Int_t n = (Int_t)fH1DeltaT.GetEntries();
+    if (n>=fNumEvents && n%(fNumEvents/60)==0)
         FitTimeBin();
 
@@ -592 +594 @@
 void MHEffectiveOnTime::PaintText(Double_t *res) const
 {
-    TLatex text(0.25, 0.94, Form("T_{eff}=%.1fs\\pm%.1fs  \\labda=%.1f\\pm%.1f  p=%.1f%%  NDF=%d",
-                                 res[0], res[1], res[3], res[4], res[2], res[5]));
+    TLatex text(0.27, 0.94, Form("T_{eff}=%.1fs\\pm%.1fs  \\lambda=%.1f\\pm%.1fHz  p=%.1f%%  \\chi^{2}/%d=%.1f",
+                                 res[0], res[1], res[3], res[4], res[2], (int)res[5], res[6]/res[5]));
     text.SetBit(TLatex::kTextNDC);
     text.SetTextSize(0.04);
@@ -613 +615 @@
         return;
 
-    TLatex text(0.45, 0.94, Form("\\bar{p} = %.1f%%  (n=%d)", sum/n, n));
+    TLatex text(0.47, 0.94, Form("\\bar{p} = %.1f%%", sum/n));
     text.SetBit(TLatex::kTextNDC);
     text.SetTextSize(0.04);
@@ -692 +694 @@
         if ((h = (TH1D*)gPad->FindObject(fNameProjDeltaT)))
         {
-            Double_t res[6];
+            Double_t res[7];
             FitH(h, res, kTRUE);
             PaintText(res);