Changeset 7205 for trunk/MagicSoft/Mars/mpointing

Timestamp:

07/21/05 17:45:25 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mpointing

Files:

• MPointingDevCalc.cc (modified) (7 diffs)
• MPointingDevCalc.h (modified) (3 diffs)
• MPointingPos.h (modified) (1 diff)
• MPointingPosCalc.h (modified) (1 diff)

trunk/MagicSoft/Mars/mpointing/MPointingDevCalc.cc

@@ -27 +27 @@
 // MPointingDevCalc
 //
+// Calculates the pointing deviation from the starguider information.
+//
+// There are some quality parameters to steer which are the valid
+// starguider data points:
+//  * MPointingDevCalc.NumMinStars: 8
+//    Minimum number of identified stars required to accep the data
+//  * MPointingDevCalc.NsbLevel:    3.0
+//    Minimum deviation (in rms) from the the mean allowed for the measured
+//    NSB (noise in the ccd camera)
+//  * MPointingDevCalc.NsbMin:      30
+//    - minimum NSB to be used in mean/rms calculation
+//  * MPointingDevCalc.NsbMax:      60
+//    - maximum NSB to be used in mean/rms calculation
+//  * MPointingDevCalc.MaxAbsDev:   15
+//    - Maximum absolute deviation which is consideres as valid (arcmin)
+//
+// Starguider data which doens't fullfill this requirements is ignored.
+// If the measures NSB==0 (file too old, starguider didn't write down
+// these values) the checks based on NSB and NumMinStar are skipped.
+//
+// The calculation of NSB mean and rms is reset for each file (ReInit)
+//
+//
 // Input Container:
 //   MReportStarguider
@@ -129 +152 @@
     if (nsb>0)
     {
-        if (nsb>30 && nsb<60)
+        if (nsb>fNsbMin && nsb<fNsbMax)
         {
             fNsbSum += nsb;
@@ -141 +164 @@
             const Double_t sq  = fNsbSq /fNsbCount;
 
-            const Double_t rms = TMath::Sqrt(sq - sum*sum);
-
-            if (nsb<sum-3*rms || nsb>sum+3*rms)
+            const Double_t rms = fNsbLevel*TMath::Sqrt(sq - sum*sum);
+
+            if (nsb<sum-rms || nsb>sum+rms)
             {
                 fSkip[2]++;
@@ -150 +173 @@
         }
 
-        if (fReport->GetNumIdentifiedStars()<8)
+        if (fReport->GetNumIdentifiedStars()<fNumMinStars)
         {
             fSkip[3]++;
@@ -166 +189 @@
 
     // Sanity check... larger deviation are strange and ignored
-    if (dev>0.25)
+    if (dev*60>fMaxAbsDev)
     {
         fSkip[4]++;
@@ -213 +236 @@
     *fLog << GetDescriptor() << " execution statistics:" << endl;
     PrintSkipped(fSkip[1], "Starguider deviation not set, is exactly 0/0");
-    PrintSkipped(fSkip[2], "NSB out of 3 sigma range");
-    PrintSkipped(fSkip[3], "Number of identified stars < 8");
-    PrintSkipped(fSkip[4], "Absolute deviation > 0.25deg");
+    PrintSkipped(fSkip[2], Form("NSB out of %.1f sigma range", fNsbLevel));
+    PrintSkipped(fSkip[3], Form("Number of identified stars < %d", fNumMinStars));
+    PrintSkipped(fSkip[4], Form("Absolute deviation > %.1farcmin", fMaxAbsDev));
     *fLog << " " << (int)fSkip[0] << " (" << Form("%5.1f", 100.*fSkip[0]/GetNumExecutions()) << "%) Evts survived calculation!" << endl;
     *fLog << endl;
@@ -221 +244 @@
     return kTRUE;
 }
+
+// --------------------------------------------------------------------------
+//
+// MPointingDevCalc.NumMinStars: 8
+// MPointingDevCalc.NsbLevel:    3.0
+// MPointingDevCalc.NsbMin:      30
+// MPointingDevCalc.NsbMax:      60
+// MPointingDevCalc.MaxAbsDev:   15
+//
+Int_t MPointingDevCalc::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
+{
+    Bool_t rc = kFALSE;
+    if (IsEnvDefined(env, prefix, "NumMinStars", print))
+    {
+        SetNumMInStars(GetEnvValue(env, prefix, "NumMinStars", fNumMinStars));
+        rc = kTRUE;
+    }
+    if (IsEnvDefined(env, prefix, "NsbLevel", print))
+    {
+        SetNsbLevel(GetEnvValue(env, prefix, "NsbLevel", fNsbLevel));
+        rc = kTRUE;
+    }
+    if (IsEnvDefined(env, prefix, "NsbMin", print))
+    {
+        SetNsbMin(GetEnvValue(env, prefix, "NsbMin", fNsbMin));
+        rc = kTRUE;
+    }
+    if (IsEnvDefined(env, prefix, "NsbMax", print))
+    {
+        SetNsbMax(GetEnvValue(env, prefix, "NsbMax", fNsbMax));
+        rc = kTRUE;
+    }
+    if (IsEnvDefined(env, prefix, "MaxAbsDev", print))
+    {
+        SetMaxAbsDev(GetEnvValue(env, prefix, "MaxAbsDev", fMaxAbsDev));
+        rc = kTRUE;
+    }
+
+    return rc;
+}

trunk/MagicSoft/Mars/mpointing/MPointingDevCalc.h

@@ -27 +27 @@
     TArrayI  fSkip;                //! Counter for execution statistics
 
+    UInt_t  fNumMinStars;          // Minimum number of identified stars
+    Float_t fNsbLevel;             // Minimum deviation from mean in sigma
+    Float_t fNsbMin;               // Minimum NSB to calc mean and rms
+    Float_t fNsbMax;               // Maximum NSB to calc mean and rms
+    Float_t fMaxAbsDev;            // [arcmin] Maximum considered absolute deviation
+
     // MPointingDevCalc
     Int_t ProcessStarguiderReport();
+
+    // MParContainer
 
     // MTask
@@ -37 +45 @@
 
 public:
-    MPointingDevCalc() : fReport(0), fDeviation(0), fSkip(5)
+    MPointingDevCalc() : fReport(0), fDeviation(0), fSkip(5), fNumMinStars(8),
+        fNsbLevel(3), fNsbMin(30), fNsbMax(60), fMaxAbsDev(15)
     {
         fName  = "MPointingDevCalc";
@@ -43 +52 @@
     }
 
+    void SetNumMinStars(UInt_t n)  { fNumMinStars=8; }
+    void SetNsbLevel(Float_t lvl)  { fNsbLevel=lvl;  }
+    void SetNsbMin(Float_t nsb)    { fNsbMin=nsb;    }
+    void SetNsbMax(Float_t nsb)    { fNsbMax=nsb;    }
+    void SetMaxAbsDev(Float_t max) { fMaxAbsDev=max; }
+
     ClassDef(MPointingDevCalc, 0) //Task calculating the pointing deviation
 };

trunk/MagicSoft/Mars/mpointing/MPointingPos.h

@@ -6 +6 @@
 #endif
 
-// FIXME: Should not be here...
+// FIXME: Should not be here... (ZdAz)
 #ifndef MARS_MPointing
 #include "MPointing.h"

trunk/MagicSoft/Mars/mpointing/MPointingPosCalc.h

@@ -1 @@
-#ifndef MARS_MPointingPositionCalc
-#define MARS_MPointingPositionCalc
+#ifndef MARS_MPointingPosCalc
+#define MARS_MPointingPosCalc
 
 #ifndef MARS_MTask