Changeset 8106 for trunk/MagicSoft/Mars/mhflux

Timestamp:

10/17/06 18:18:40 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mhflux

Files:

• MHFalseSource.cc (modified) (6 diffs)
• MHThetaSqN.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/mhflux/MHFalseSource.cc

@@ -1 +1 @@
 /* ======================================================================== *\
+! $Name: not supported by cvs2svn $:$Id: MHFalseSource.cc,v 1.20 2006-10-17 17:16:00 tbretz Exp $
+! --------------------------------------------------------------------------
 !
 ! *
@@ -393 +395 @@
 // the same number of bins than for on-data
 //
-void MHFalseSource::ProjectOff(const TH3D &src, TH2D *h2, TH2D *all)
+void MHFalseSource::ProjectOff(const TH3D &src, TH2D *h2, TH2D *hall)
 {
     TAxis &axe = *src.GetZaxis();
@@ -415 +417 @@
     // Move contents from projection to h2
     h2->Reset();
-    h2->Add(p, all->GetMaximum());
-    h2->Divide(all);
+    h2->Add(p, hall->GetMaximum());
+    h2->Divide(hall);
 
     // Delete p
@@ -434 +436 @@
 // range (0, fAlphaCut)
 //
-void MHFalseSource::ProjectOn(const TH3D &src, TH2D *h3, TH2D *all)
+void MHFalseSource::ProjectOn(const TH3D &src, TH2D *h3, TH2D *hall)
 {
     TAxis &axe = *src.GetZaxis();
@@ -452 +454 @@
     // Move contents from projection to h3
     h3->Reset();
-    h3->Add(p, all->GetMaximum());
-    h3->Divide(all);
+    h3->Add(p, hall->GetMaximum());
+    h3->Divide(hall);
 
     // Delete p
@@ -603 +605 @@
             h->SetTitle(Form("Distribution of \\alpha for x=%.2f y=%.2f (S_{max}=%.1f\\sigma)", x, y, s));
 
-            TH1D *h0=0;
-            if ((h0 = (TH1D*)gPad->FindObject("AlphaOff_z")))
+            TH1D *ho=0;
+            if ((ho = (TH1D*)gPad->FindObject("AlphaOff_z")))
             {
                 fHistOff->ProjectionZ("AlphaOff_z", maxx, maxx, maxy, maxy);
 
                 /* ============= local scaling ================ */
-                const Int_t f = h0->GetXaxis()->FindFixBin(fBgMean-1.5*fAlphaCut);
-                const Int_t l = h0->GetXaxis()->FindFixBin(fAlphaCut*3)+f-1;
-                h0->Scale(h1->Integral(f, l)/h0->Integral(f, l));
+                const Int_t f = ho->GetXaxis()->FindFixBin(fBgMean-1.5*fAlphaCut);
+                const Int_t l = ho->GetXaxis()->FindFixBin(fAlphaCut*3)+f-1;
+                ho->Scale(h1->Integral(f, l)/ho->Integral(f, l));
 
 

trunk/MagicSoft/Mars/mhflux/MHThetaSqN.cc

@@ -1 +1 @@
 /* ======================================================================== *\
+! $Name: not supported by cvs2svn $:$Id: MHThetaSqN.cc,v 1.8 2006-10-17 17:16:00 tbretz Exp $
+! --------------------------------------------------------------------------
 !
 ! *
@@ -282 +284 @@
         //  2: 0.346  ~98%   0.260
 
-        const Double_t dist = src0.Mod()*TMath::Sin(rad/2);
-        const Double_t cut  = TMath::Sqrt(fFit.GetSignalIntegralMax());
-        if (dist<cut)
+        const Double_t dis = src0.Mod()*TMath::Sin(rad/2);
+        const Double_t cut = TMath::Sqrt(fFit.GetSignalIntegralMax());
+        if (dis<cut)
         {
             *fLog << warn << "WARNING - Source regions overlap: distance ";
-            *fLog << dist << " less than theta-sq cut " << cut << "!" << endl;
-            if (dist*1.7<cut*1.0)
+            *fLog << dis << " less than theta-sq cut " << cut << "!" << endl;
+            if (dis*1.7<cut*1.0)
                 fCounter[3]++;
             else
-                if (dist*1.7<cut*1.5)
+                if (dis*1.7<cut*1.5)
                     fCounter[2]++;
                 else