Changeset 7744 for trunk/MagicSoft/Mars/mpointing/MPointingDevCalc.cc

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mpointing/MPointingDevCalc.cc (modified) (10 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;