Changeset 2178 for trunk/MagicSoft/Mars/mhist/MHCurrents.cc

Timestamp:

06/16/03 13:45:01 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhist/MHCurrents.cc (modified) (9 diffs)

trunk/MagicSoft/Mars/mhist/MHCurrents.cc

@@ -36 +36 @@
 
 #include "MParList.h"
+#include "MBinning.h"
 #include "MCurrents.h"
 #include "MCamDisplay.h"
@@ -54 +55 @@
     // FIXME: Implement a clear function with setmem
     for (int i=0; i<577; i++)
+    {
         fSum[i] = 0;
+        fRms[i] = 0;
+    }
 
     fEntries = 0;
@@ -65 +69 @@
 //
 MHCurrents::MHCurrents(const char *name, const char *title)
-    : fSum(577), fCam(NULL), fEvt(NULL), fDispl(NULL)
+    : fSum(577), fRms(577), fCam(NULL), fEvt(NULL), fDispl(NULL)
 {
     //
@@ -74 +78 @@
 
     Clear();
+
+    fHist.SetName("currents");
+    fHist.SetTitle("Avg.Currents [nA]");
+    fHist.SetDirectory(NULL);
+    fHist.Sumw2();
 }
 
@@ -105 +114 @@
     Clear();
 
+    MBinning bins;
+    bins.SetEdges(577, -0.5, 576.5);
+    bins.Apply(fHist);
+
     return kTRUE;
 }
@@ -121 +134 @@
     }
 
+    for (UInt_t idx=0; idx<577; idx++)
+    {
+        Float_t val;
+        if (!evt->GetPixelContent(val, idx))
+            continue;
+
+        fSum[idx] += val;
+        fRms[idx] += val*val;
+    }
+
+    fEntries++;
+
+    return kTRUE;
+}
+
+// --------------------------------------------------------------------------
+//
+// Scale the sum container with the number of entries
+//
+Bool_t MHCurrents::Finalize()
+{
+    if (fEntries<2)
+    {
+        *fLog << warn << "WARNING - " << GetDescriptor() << " doesn't contain enough entries." << endl;
+        return kTRUE;
+    }
+
     for (UInt_t i=0; i<577; i++)
-        fSum[i] += (*evt)[i];
-
-    fEntries++;
-
+    {
+        // calc sdev^2 for pixel index i
+        // var^2 = (sum[xi^2] - sum[xi]^2/n) / (n-1);
+        fRms[i] -= fSum[i]*fSum[i]/fEntries;
+        fRms[i] /= fEntries-1;
+        fRms[i]  = TMath::Sqrt(fRms[i]);
+
+        // calc mean value for pixel index i
+        fSum[i] /= fEntries;
+
+        fHist.SetBinContent(i+1, fSum[i]);
+        fHist.SetBinError(  i+1, fRms[i]);
+    }
     return kTRUE;
 }
@@ -131 +180 @@
 // --------------------------------------------------------------------------
 //
-// Scale the sum container with the number of entries
-//
-Bool_t MHCurrents::Finalize()
-{
-    for (UInt_t i=0; i<577; i++)
-        fSum[i] /= fEntries;
-
-    return kTRUE;
-}
-
-// --------------------------------------------------------------------------
-//
 // Draw the present 'fill status'
 //
-void MHCurrents::Draw(Option_t *)
+void MHCurrents::Draw(Option_t *o)
 {
     if (!fCam)
@@ -156 +193 @@
     pad->SetBorderMode(0);
 
+    SetDrawOption(o);
     AppendPad("");
 }
@@ -174 +212 @@
         fDispl = new MCamDisplay(fCam);
 
-    fDispl->FillCurrents(fSum);
+    TString opt(GetDrawOption());
+
+    fDispl->Fill(opt.Contains("rms", TString::kIgnoreCase) ? fRms : fSum);
     fDispl->Paint();
 }