Changeset 1540 for trunk/MagicSoft/Mars/manalysis/MPedestalCam.cc

Timestamp:

10/15/02 17:02:46 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MPedestalCam.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/manalysis/MPedestalCam.cc

@@ -31 +31 @@
 /////////////////////////////////////////////////////////////////////////////
 #include "MPedestalCam.h"
+#include "MPedestalPix.h"
+
+#include <TClonesArray.h>
 
 #include "MLog.h"
+#include "MLogManip.h"
+
+#include "MGeomCam.h"
+
 
 ClassImp(MPedestalCam);
@@ -62 +69 @@
 // --------------------------------------------------------------------------
 //
+// Set the size of the camera
+//
+inline void MPedestalCam::InitSize(const UInt_t i)
+{
+    fArray->ExpandCreate(i);
+}
+
+// --------------------------------------------------------------------------
+//
+// Get the size of the MPedestalCam
+//
+inline Int_t MPedestalCam::GetSize() const
+{
+    return fArray->GetEntries();
+}
+
+// --------------------------------------------------------------------------
+//
+// Get i-th pixel (pixel number)
+//
+inline MPedestalPix &MPedestalCam::operator[](Int_t i)
+{
+    return *(MPedestalPix*)fArray->UncheckedAt(i);
+}
+
+// --------------------------------------------------------------------------
+//
+// Get i-th pixel (pixel number)
+//
+inline MPedestalPix &MPedestalCam::operator[](Int_t i) const
+{
+    return *(MPedestalPix*)fArray->UncheckedAt(i);
+}
+
+// --------------------------------------------------------------------------
+//
 // Check if position i is inside bounds
 //
@@ -69 +112 @@
 } 
 
+void MPedestalCam::Clear(Option_t *o)
+{
+    fArray->ForEach(TObject, Clear)();
+}
 
+void MPedestalCam::Print(Option_t *o="") const
+{
+    *fLog << all << GetDescriptor() << ":" << endl;
+    int id = 0;
+
+    TIter Next(fArray);
+    MPedestalPix *pix;
+    while ((pix=(MPedestalPix*)Next()))
+    {
+        id++;
+
+        if (!pix->IsValid())
+            continue;
+
+        *fLog << id-1 << ": ";
+        *fLog << pix->GetMean() << " " << pix->GetSigma() << " ";
+        *fLog << pix->GetMeanRms() << " " << pix->GetSigmaRms() << endl;
+    }
+}
+
+Float_t MPedestalCam::GetMeanMin(const MGeomCam *geom) const
+{
+    if (fArray->GetEntries() <= 0)
+        return 50.;
+
+    Float_t minval = (*this)[0].GetMean();
+
+    for (Int_t i=1; i<fArray->GetEntries(); i++)
+    {
+        const MPedestalPix &pix = (*this)[i];
+
+        Float_t testval = pix.GetMean();
+
+        if (geom)
+            testval *= geom->GetPixRatio(i);
+
+        if (testval < minval)
+            minval = testval;
+    }
+    return minval;
+}
+
+Float_t MPedestalCam::GetMeanMax(const MGeomCam *geom) const
+{
+    if (fArray->GetEntries() <= 0)
+        return 50.;
+
+    Float_t maxval = (*this)[0].GetMean();
+
+    for (Int_t i=1; i<fArray->GetEntries(); i++)
+    {
+        const MPedestalPix &pix = (*this)[i];
+
+        Float_t testval = pix.GetMean();
+
+        if (geom)
+            testval *= geom->GetPixRatio(i);
+
+        if (testval > maxval)
+            maxval = testval;
+    }
+    return maxval;
+}