Changeset 8689 for trunk/MagicSoft/Mars/mpointing/MHSrcPosCam.cc

Timestamp:

08/21/07 23:20:48 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mpointing/MHSrcPosCam.cc (modified) (7 diffs)

trunk/MagicSoft/Mars/mpointing/MHSrcPosCam.cc

@@ -33 +33 @@
 #include "MHSrcPosCam.h"
 
+#include <TVector2.h>
 #include <TCanvas.h>
 #include <TEllipse.h>
@@ -57 +58 @@
 //
 MHSrcPosCam::MHSrcPosCam(Bool_t wobble, const char *name, const char *title)
-    : fTimeEffOn(NULL), fEffOnTime(NULL), fSourcePos(NULL)
+    : fTimeEffOn(NULL), fEffOnTime(NULL), fSourcePos(NULL),
+    fPositions("TVector2", 50000)
 {
     //
@@ -119 +121 @@
 
     fHist.Reset();
-    fXY            = TVector2();
-    fNum           = 0;
     fTimeLastEffOn = MTime();
     fConvMm2Deg    = geom->GetConvMm2Deg();
+    fNum           = 0;
 
     return kTRUE;
@@ -129 +130 @@
 // --------------------------------------------------------------------------
 //
-// 
-// 
+// All source positions are buffered until the time of the effective on
+// time time stamp changes. Then the observation time is split into
+// identical parts and the histogram is filled by these events. The
+// effective on time time stamp is reset and the buffered source positions
+// deleted.
+//
 Bool_t MHSrcPosCam::Fill(const MParContainer *par, const Stat_t w)
 {
@@ -140 +145 @@
     }
 
-//    if (fName=="MHSrcPosCam")
-//    {
-    fXY += cam->GetXY();
-    fNum++;
-
+    // Increase array size if necessary
+    if (fNum==fPositions.GetSize())
+        fPositions.Expand(fNum*2);
+
+    // buffer position into array (could be speed up a little bit more
+    // by using ExpandCreate and memcpy)
+    new (fPositions[fNum++]) TVector2(cam->GetXY()*fConvMm2Deg);
+
+    // Check if there is a new effective on time
     if (fTimeLastEffOn==MTime())
         fTimeLastEffOn=*fTimeEffOn;
@@ -151 +160 @@
         return kTRUE;
 
-    fXY *= fConvMm2Deg/fNum;
-
-    fHist.Fill(fXY.X(), fXY.Y(), fEffOnTime->GetVal());
-//    }
-//    else
-//        fHist.Fill(cam->GetX()*fConvMm2Deg, cam->GetY()*fConvMm2Deg);
-
-    fXY            = TVector2();
-    fNum           = 0;
+    // Split the observation time to all buffered events
+    const Double_t scale = fEffOnTime->GetVal()/fNum;
+
+    // Fill histogram from array
+    for (int i=0; i<fNum; i++)
+    {
+        const TVector2 &v = (TVector2&)*fPositions[i];
+        fHist.Fill(v.X(), v.Y(), scale);
+    }
+
+    // reset time stamp and remove all buffered positions
     fTimeLastEffOn = *fTimeEffOn;
+    fNum = 0;
 
     return kTRUE;
 }
 
+// --------------------------------------------------------------------------
+//
 void MHSrcPosCam::Paint(Option_t *)
 {
@@ -172 +186 @@
 // --------------------------------------------------------------------------
 //
-// 
-// 
 void MHSrcPosCam::Draw(Option_t *)
 {