Changeset 8281 for trunk/MagicSoft/Mars/mhist/MHCamera.cc

Timestamp:

02/01/07 11:24:19 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHCamera.cc (modified) (7 diffs)

trunk/MagicSoft/Mars/mhist/MHCamera.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MHCamera.cc,v 1.97 2007-01-30 14:20:47 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MHCamera.cc,v 1.98 2007-02-01 11:24:18 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -319 +319 @@
 // ------------------------------------------------------------------------
 //
-// Return the median value on the y-axis (profile option is correctly taken
-// into account)
-//
-Stat_t MHCamera::GetMedian() const
-{
-    // Just for speed reasons
-    if (!TestBit(kProfile))
-        return TMath::Median(GetSize()-2, GetArray()+1);
-
-    // Copy profiled data into new array (FIXME: Should we take errors into account?)
-    TArrayD arr(fNcells-2);
-    for (int i=1; i<fNcells-1; i++)
-        arr[i-1] = GetBinContent(i);
-
-    // return Median of the profile data
-    return TMath::Median(arr.GetSize(), arr.GetArray());
-}
-
-// ------------------------------------------------------------------------
-//
-// Return the median value (divided by MMath::GausProb(1.0)) of the
-// distribution of abs(y[i]-Median). This is my Median equivalent of the RMS
-//
-Stat_t MHCamera::GetMedianDev() const
-{
-    // Just for speed reasons
-    if (!TestBit(kProfile))
-        return MMath::MedianDev(GetSize()-2, GetArray()+1);
-
-    // Copy profiled data into new array (FIXME: Should we take errors into account?)
-    TArrayD arr(fNcells-2);
-    for (int i=1; i<fNcells-1; i++)
-        arr[i-1] = GetBinContent(i);
-
-    // return MedianDev of the profile data
-    return MMath::MedianDev(arr.GetSize(), arr.GetArray());
-}
-
-// ------------------------------------------------------------------------
-//
 // Return the mean value of all entries which are used if all=kFALSE and
 // of all entries if all=kTRUE if sector<0. If sector>=0 only
@@ -388 +348 @@
 // ------------------------------------------------------------------------
 //
-// Return the sqrt variance of all entries which are used if all=kFALSE and
+// Return the median value of all entries which are used if all=kFALSE and
 // of all entries if all=kTRUE if sector<0. If sector>=0 only
 // entries with match the given sector are taken into account.
 //
-Stat_t MHCamera::GetRmsSectors(const TArrayI &sector, const TArrayI &aidx, Bool_t ball) const
+Stat_t MHCamera::GetMedianSectors(const TArrayI &sector, const TArrayI &aidx, Bool_t ball) const
 {
     if (fNcells<=1)
-        return -1;
-
+        return 0;
+
+    TArrayD arr(fNcells-2);
     Int_t n=0;
-
-    Stat_t sum = 0;
-    Stat_t sq  = 0;
+ 
     for (int i=0; i<fNcells-2; i++)
     {
@@ -408 +367 @@
                 continue;
 
+            arr[n++] = TestBit(kProfile) ? fArray[i+1]/fBinEntries[i+1] : fArray[i+1];
+        }
+    }
+
+    // return Median of the profile data
+    return TMath::Median(n, arr.GetArray());
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the sqrt variance of all entries which are used if all=kFALSE and
+// of all entries if all=kTRUE if sector<0. If sector>=0 only
+// entries with match the given sector are taken into account.
+//
+Stat_t MHCamera::GetRmsSectors(const TArrayI &sector, const TArrayI &aidx, Bool_t ball) const
+{
+    if (fNcells<=1)
+        return -1;
+
+    Int_t n=0;
+
+    Stat_t sum = 0;
+    Stat_t sq  = 0;
+    for (int i=0; i<fNcells-2; i++)
+    {
+        if ((ball || IsUsed(i)) && MatchSector(i, sector, aidx))
+        {
+            if (TestBit(kProfile) && fBinEntries[i+1]==0)
+                continue;
+
             const Double_t val = TestBit(kProfile) ? fArray[i+1]/fBinEntries[i+1] : fArray[i+1];
 
@@ -423 +412 @@
 
     return TMath::Sqrt(sq-sum*sum);
+}
+
+// ------------------------------------------------------------------------
+//
+// Return the median value (divided by MMath::GausProb(1.0)=68.3%) of the
+// distribution of abs(y[i]-Median). This is my Median equivalent of the RMS.
+// Return the deviation of all entries which are used if all=kFALSE and
+// of all entries if all=kTRUE if sector<0. If sector>=0 only
+// entries with match the given sector are taken into account.
+//
+Stat_t MHCamera::GetDevSectors(const TArrayI &sector, const TArrayI &aidx, Bool_t ball) const
+{
+    if (fNcells<=1)
+        return 0;
+
+    TArrayD arr(fNcells-2);
+    Int_t n=0;
+ 
+    for (int i=0; i<fNcells-2; i++)
+    {
+        if ((ball || IsUsed(i)) && MatchSector(i, sector, aidx))
+        {
+            if (TestBit(kProfile) && fBinEntries[i+1]==0)
+                continue;
+
+            arr[n++] = TestBit(kProfile) ? fArray[i+1]/fBinEntries[i+1] : fArray[i+1];
+        }
+    }
+
+    // return Median of the profile data
+    return MMath::MedianDev(n, arr.GetArray());
 }
 
@@ -1391 +1411 @@
     {
     case 1:
+        // Under-/Overflow bins not handled!
         for (Int_t idx=0; idx<fNcells-2; idx++)
-            Fill(idx, d.GetBinError(idx+1));
+            if (d.IsUsed(idx))
+                Fill(idx, d.GetBinError(idx+1));
+        fEntries++;
         break;
     case 2:
+        // Under-/Overflow bins not handled!
         for (Int_t idx=0; idx<fNcells-2; idx++)
-            if (d.GetBinContent(idx+1)!=0)
+            if (d.GetBinContent(idx+1)!=0 && d.IsUsed(idx))
                 Fill(idx, TMath::Abs(d.GetBinError(idx+1)/d.GetBinContent(idx+1)));
+        fEntries++;
         break;
     default:
+        if (TestBit(kProfile)!=d.TestBit(kProfile))
+            gLog << warn << "WARNING - You have tried to call AddCamContent for two different kind of histograms (kProfile set or not)." << endl;
+
+        // environment
+        fEntries += d.fEntries;
+        fTsumw   += d.fTsumw;
+        fTsumw2  += d.fTsumw2;
+        fTsumwx  += d.fTsumwx;
+        fTsumwx2 += d.fTsumwx2;
+        // Bin contents
+        for (Int_t idx=1; idx<fNcells-1; idx++)
+        {
+            if (!d.IsUsed(idx-1))
+                continue;
+
+            fArray[idx]          += d.fArray[idx];
+            fBinEntries[idx]     += d.fBinEntries[idx];
+            fSumw2.fArray[idx]   += d.fSumw2.fArray[idx];
+        }
+        // Underflow bin
+        fArray[0]                += d.fArray[0];
+        fBinEntries[0]           += d.fBinEntries[0];
+        fSumw2.fArray[0]         += d.fSumw2.fArray[0];
+        // Overflow bin
+        fArray[fNcells-1]        += d.fArray[fNcells-1];
+        fBinEntries[fNcells-1]   += d.fBinEntries[fNcells-1];
+        fSumw2.fArray[fNcells-1] += d.fSumw2.fArray[fNcells-1];
+        break;
+/*    default:
+        if (TestBit(kProfile)!=d.TestBit(kProfile))
+            gLog << warn << "WARNING - You have tried to call AddCamContent for two different kind of histograms (kProfile set or not)." << endl;
+
         for (Int_t idx=0; idx<fNcells-2; idx++)
             Fill(idx, d.GetBinContent(idx+1));
-        break;
+        break;*/
     }
     fEntries++;
@@ -1629 +1686 @@
     TH1::Reset(opt);
 
-    for (Int_t i=0; i<fNcells-2; i++)
-    {
-        fArray[i+1]=0;
-        fBinEntries[i]=0;
-        ResetUsed(i);
-    }
-
-    fArray[0]         = 0;
-    fArray[fNcells-1] = 0;
+    fUsed.Reset();
+    fBinEntries.Reset();
+
+    for (Int_t i=0; i<fNcells; i++)
+        fArray[i] = 0;
 }