Changeset 7744 for trunk/MagicSoft/Mars/mpointing

Timestamp:

05/24/06 19:00:45 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mpointing

Files:

• MPointingDevCalc.cc (modified) (10 diffs)
• MPointingDevCalc.h (modified) (5 diffs)

trunk/MagicSoft/Mars/mpointing/MPointingDevCalc.cc

@@ -84 +84 @@
     }
 
-    fNsbSum   = 0;
-    fNsbSq    = 0;
-    fNsbCount = 0;
+    fNsbSum   =  0;
+    fNsbSq    =  0;
+    fNsbCount =  0;
 
     fRunType = run->GetRunType();
@@ -127 +127 @@
     // We use kRTNone here as a placeholder for data runs.
     fRunType  = MRawRunHeader::kRTNone;
+    fLastMjd  = -1;
 
     fSkip.Reset();
 
     return fDeviation ? kTRUE : kFALSE;
+}
+
+// --------------------------------------------------------------------------
+//
+// Increase fSkip[i] by one. If the data in fDeviation is outdated (older
+// than fMaxAge) and the current report should be skipped reset DevZdAz and
+// DevXY and fSkip[5] is increased by one.
+//
+void MPointingDevCalc::Skip(Int_t i)
+{
+    fSkip[i]++;
+
+    const Double_t diff = (fReport->GetMjd()-fLastMjd)*1440; // [min] 1440=24*60
+    if (diff<fMaxAge && fLastMjd>0)
+        return;
+
+    fDeviation->SetDevZdAz(0, 0);
+    fDeviation->SetDevXY(0, 0);
+    fSkip[5]++;
 }
 
@@ -142 +162 @@
         fDeviation->SetDevXY(0, 0);   //?!?
         fSkip[1]++;
+        fLastMjd = fReport->GetMjd();
         return kTRUE;
     }
@@ -167 +188 @@
             if (nsb<sum-rms || nsb>sum+rms)
             {
-                fSkip[2]++;
+                Skip(2);
                 return kTRUE;
             }
@@ -174 +195 @@
         if (fReport->GetNumIdentifiedStars()<fNumMinStars)
         {
-            fSkip[3]++;
+            Skip(3);
             return kTRUE;
         }
@@ -185 +206 @@
     if (dev*60>fMaxAbsDev)
     {
-        fSkip[4]++;
+        Skip(4);
         return kTRUE;
     }
@@ -198 +219 @@
 
     fSkip[0]++;
+    fLastMjd = fReport->GetMjd();
 
     return kTRUE;
@@ -238 +260 @@
     PrintSkipped(fSkip[3], Form("Number of identified stars < %d", fNumMinStars));
     PrintSkipped(fSkip[4], Form("Absolute deviation > %.1farcmin", fMaxAbsDev));
+    PrintSkipped(fSkip[5], Form("Events set to 0 because older than %.1fmin", fMaxAge));
     *fLog << " " << (int)fSkip[0] << " (" << Form("%5.1f", 100.*fSkip[0]/GetNumExecutions()) << "%) Evts survived calculation!" << endl;
     *fLog << endl;
@@ -251 +274 @@
 // MPointingDevCalc.NsbMax:      60
 // MPointingDevCalc.MaxAbsDev:   15
+// MPointingDevCalc.MaxAge:      1.0
 //
 Int_t MPointingDevCalc::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
@@ -290 +314 @@
         rc = kTRUE;
     }
+    if (IsEnvDefined(env, prefix, "MaxAge", print))
+    {
+        fMaxAge = GetEnvValue(env, prefix, "MaxAge", fMaxAge);
+        rc = kTRUE;
+    }
 
     return rc;

trunk/MagicSoft/Mars/mpointing/MPointingDevCalc.h

@@ -26 +26 @@
 
     TArrayI  fSkip;                //! Counter for execution statistics
+    Double_t fLastMjd;             //! Time of last processed report
 
     UInt_t  fNumMinStars;          // Minimum number of identified stars
@@ -32 +33 @@
     Float_t fNsbMax;               // Maximum NSB to calc mean and rms
     Float_t fMaxAbsDev;            // [arcmin] Maximum considered absolute deviation
+    Float_t fMaxAge;               // [min] Maximum age of reports to be used without an update
 
     Float_t fDx;                   // Starguider calibration dx
@@ -38 +40 @@
     // MPointingDevCalc
     Int_t ProcessStarguiderReport();
+    void  Skip(Int_t i);
 
     // MParContainer
@@ -49 +52 @@
 
 public:
-    MPointingDevCalc() : fReport(0), fDeviation(0), fSkip(5), fNumMinStars(8),
-        fNsbLevel(3), fNsbMin(30), fNsbMax(60), fMaxAbsDev(15), fDx(-7), fDy(16)
+    MPointingDevCalc() : fReport(0), fDeviation(0), fSkip(6), fNumMinStars(8),
+        fNsbLevel(3), fNsbMin(30), fNsbMax(60), fMaxAbsDev(15), fMaxAge(1), fDx(-7), fDy(16)
     {
         fName  = "MPointingDevCalc";
@@ -65 +68 @@
     void SetDx(Float_t dx)         { fDx=dx; }
     void SetDy(Float_t dy)         { fDy=dy; }
+    void SetMaxAge(Float_t age)    { fMaxAge=age; }
 
     ClassDef(MPointingDevCalc, 0) //Task calculating the pointing deviation