Changeset 5349 for trunk/MagicSoft/Mars/mtemp/mmpi

Timestamp:

11/06/04 18:52:16 (20 years ago)

Author:

mazin

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mmpi

Files:

• MSkyPlot.cc (modified) (2 diffs)
• MSkyPlot.h (modified) (1 diff)

trunk/MagicSoft/Mars/mtemp/mmpi/MSkyPlot.cc

@@ -367 +367 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Set the upper and lower limit for the background region in the alpha plot
+// to estimate the number of background events
+//
+void MSkyPlot::SetAlphaBGLimits(Float_t alphalow, Float_t alphaup)
+{
+    Float_t alphatemp;
+    if (alphalow<0)
+        *fLog << warn << "Alpha<0... taking absolute value." << endl;
+
+    if (alphaup<0)
+        *fLog << warn << "Alpha<0... taking absolute value." << endl;
+
+    if (TMath::Abs(alphaup)<TMath::Abs(alphalow)) {
+        *fLog << warn << "AlphaLow > AlphaUp... exchanging limits." << endl;
+        alphatemp = alphaup;
+        alphaup = alphalow;
+        alphalow = alphatemp;
+    }
+
+    fAlphaBgLow = TMath::Abs(alphalow);
+    fAlphaBgUp  = TMath::Abs(alphaup);
+}
+
 // calculate the values for the cuts:
 Double_t MSkyPlot::CalcLimit(Double_t *a, Double_t ls, Double_t ls2, Double_t dd2)
@@ -701 +726 @@
 
          // fit parabola to a background region
-         (*alpha_iterator).Fit(fitbgpar,"NWR");  // NWR OK?????????????????????????
+         (*alpha_iterator).Fit(fitbgpar,"EQRN");  // NWR OK?????????????????????????
          // get Chi2
          chisquarefit = fitbgpar->GetChisquare();

trunk/MagicSoft/Mars/mtemp/mmpi/MSkyPlot.h

@@ -144 +144 @@
 
     void SetAlphaCut(Float_t alpha); 
+    void SetAlphaBGLimits(Float_t alphalow, Float_t alphalup); 
 
     //std::vector <TH1D>::iterator  GetAlphaPlots()      { std::vector <TH1D>::iterator iter; return  iter = fHistAlpha.begin()  ; }