Changeset 7438 for trunk/MagicSoft/Mars/mgui

Timestamp:

12/01/05 16:43:46 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mgui

Files:

• MCamEvent.cc (modified) (2 diffs)
• MCamEvent.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mgui/MCamEvent.cc

@@ -62 +62 @@
 Double_t MCamEvent::GetCameraMean(const MGeomCam &cam, Int_t type) const
 {
+    return GetCameraStat(cam, type)[0];
+    /*
     Int_t    num  = 0;
     Double_t mean = 0;
@@ -74 +76 @@
     }
     return num == 0 ? 0 : mean/num;
+    */
 }
+
+// --------------------------------------------------------------------------
+//
+// Return the rms of all corresponding GetPixelContent
+//
+Double_t MCamEvent::GetCameraRMS(const MGeomCam &cam, Int_t type) const
+{
+    return GetCameraStat(cam, type)[1];
+    /*
+    Int_t    num  = 0;
+    Double_t sum = 0;
+    Double_t sq  = 0;
+
+    for (unsigned int i=0; i<cam.GetNumPixels(); i++)
+    {
+        Double_t val;
+        if (!GetPixelContent(val, i, cam, type))
+            continue;
+
+        sum += val;
+        sq  += val;
+        num ++;
+    }
+
+    if (num==0)
+        return 0;
+
+    sum /= num;
+    sq  /= num;
+
+    return TMath::Sqrt(sq-sum*sum);
+    */
+}
+
+// --------------------------------------------------------------------------
+//
+// Return the rms of all corresponding GetPixelContent
+//
+TArrayD MCamEvent::GetCameraStat(const MGeomCam &cam, Int_t type) const
+{
+    Int_t    num = 0;
+    Double_t sum = 0;
+    Double_t sq  = 0;
+
+    for (unsigned int i=0; i<cam.GetNumPixels(); i++)
+    {
+        Double_t val;
+        if (!GetPixelContent(val, i, cam, type))
+            continue;
+
+        sum += val;
+        sq  += val;
+        num ++;
+    }
+
+    if (num==0)
+        return 0;
+
+    sum /= num;
+    sq  /= num;
+
+    TArrayD stat(2);
+    stat[0] = sum;
+    stat[1] = TMath::Sqrt(sq-sum*sum);
+    return stat;
+}

trunk/MagicSoft/Mars/mgui/MCamEvent.h

@@ -2 +2 @@
 #define MARS_MCamEvent
 
-#ifndef ROOT_TROOT
-#include <TROOT.h>
+#ifndef ROOT_TArrayD
+#include <TArrayD.h>
 #endif
 
@@ -18 +18 @@
 
     virtual Double_t GetCameraMean(const MGeomCam &cam, Int_t type=0) const;
+    virtual Double_t GetCameraRMS(const MGeomCam &cam, Int_t type=0) const;
+    virtual TArrayD  GetCameraStat(const MGeomCam &cam, Int_t type=0) const;
 
     ClassDef(MCamEvent, 0) // A camera event