Changeset 5957 for trunk/MagicSoft/Mars/mhflux

Timestamp:

01/24/05 11:08:09 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

• MHFalseSource.cc (modified) (7 diffs)
• MHFalseSource.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mhflux/MHFalseSource.cc

@@ -220 +220 @@
 // --------------------------------------------------------------------------
 //
-// Calculate Significance as
-// significance = (s-b)/sqrt(s+k*k*b) mit k=s/b
-//
-// s: total number of events in signal region
-// b: number of background events in signal region
-// 
-Double_t MHFalseSource::Significance(Double_t s, Double_t b)
-{
-    return MMath::SignificanceSym(s, b);
-}
-
-// --------------------------------------------------------------------------
-//
-//  calculates the significance according to Li & Ma
-//  ApJ 272 (1983) 317, Formula 17
-//
-//  s                    // s: number of on events
-//  b                    // b: number of off events
-//  alpha = t_on/t_off;  // t: observation time
-//
-Double_t MHFalseSource::SignificanceLiMa(Double_t s, Double_t b, Double_t alpha)
-{
-    return MMath::SignificanceLiMaSigned(s, b);
-}
-
-// --------------------------------------------------------------------------
-//
 // Set binnings (takes BinningFalseSource) and prepare filling of the
 // histogram.
@@ -599 +572 @@
                 const Double_t b = h2->GetBinContent(n);
 
-                const Double_t sig = SignificanceLiMa(s, b);
+                const Double_t sig = MMath::SignificanceLiMaSigned(s, b);
 
                 h4->SetBinContent(n, sig);
@@ -825 +798 @@
     gPad->Modified();
     gPad->cd();
+}
+
+// --------------------------------------------------------------------------
+//
+// The following resources are available:
+//
+//    MHFalseSource.DistMin: 0.5
+//    MHFalseSource.DistMax: 1.4
+//    MHFalseSource.DWMin:   0.1
+//    MHFalseSource.DWMax:   0.3
+//
+Int_t MHFalseSource::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
+{
+    Bool_t rc = kFALSE;
+    if (IsEnvDefined(env, prefix, "DistMin", print))
+    {
+        rc = kTRUE;
+        SetMinDist(GetEnvValue(env, prefix, "DistMin", fMinDist));
+    }
+    if (IsEnvDefined(env, prefix, "DistMax", print))
+    {
+        rc = kTRUE;
+        SetMaxDist(GetEnvValue(env, prefix, "DistMax", fMaxDist));
+    }
+
+    if (IsEnvDefined(env, prefix, "DWMin", print))
+    {
+        rc = kTRUE;
+        SetMinDW(GetEnvValue(env, prefix, "DWMin", fMinDist));
+    }
+    if (IsEnvDefined(env, prefix, "DWMax", print))
+    {
+        rc = kTRUE;
+        SetMaxDW(GetEnvValue(env, prefix, "DWMax", fMaxDist));
+    }
+
+    return rc;
 }
 
@@ -965 +975 @@
                 continue;
 
+            // This is a quick hack to correct for the lower acceptance at
+            // the edge of the camera
             h->Scale(entries/h->GetEntries());
 
@@ -970 +982 @@
             {
                 hoff = fHistOff->ProjectionZ("AlphaFitOff", ix, ix, iy, iy);
+                // This is a quick hack to correct for the lower acceptance at
+                // the edge of the camera
                 hoff->Scale(entries/h->GetEntries());
                 scale = fit.Scale(*hoff, *h);
@@ -982 +996 @@
  
             // Fill results into some histograms
-            h0a.Fill(fit.GetGausA()); // gaus-A
+            h0a.Fill(fit.GetGausA());        // gaus-A
             h0b.Fill(fit.GetCoefficient(3)); // 3
-            h1.Fill(fit.GetGausMu());  // mu
-            h2.Fill(fit.GetGausSigma());  // sigma-gaus
+            h1.Fill(fit.GetGausMu());        // mu
+            h2.Fill(fit.GetGausSigma());     // sigma-gaus
             if (polynom>1 && !fHistOff)
                 h3.Fill(fit.GetCoefficient(5));
@@ -1193 +1207 @@
             const Double_t b = f2.Integral(0, (float)i)/w;
 
-            const Double_t sig = SignificanceLiMa(s, b);
+            const Double_t sig = MMath::SignificanceLiMaSigned(s, b);
 
             g->SetPoint(g->GetN(), i, sig);

trunk/MagicSoft/Mars/mhflux/MHFalseSource.h

@@ -59 +59 @@
     MHFalseSource(const char *name=NULL, const char *title=NULL);
 
-    Bool_t SetupFill(const MParList *pList);
-    Bool_t Fill(const MParContainer *par, const Stat_t w=1);
-
     TH1 *GetHistByName(const TString name) { return &fHist; }
 
@@ -90 +87 @@
     }
 
+    // MParContainer
+    Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print=kFALSE);
+
+    // MH
+    Bool_t SetupFill(const MParList *pList);
+    Bool_t Fill(const MParContainer *par, const Stat_t w=1);
+
+    // TObject
     void Paint(Option_t *opt="");
     void Draw(Option_t *option="");
-
-    static Double_t Significance(Double_t s, Double_t b);
-    static Double_t SignificanceLiMa(Double_t s, Double_t b, Double_t alpha=1);
 
     ClassDef(MHFalseSource, 1) //3D-histogram in alpha, x and y