Changeset 17145

Timestamp:

09/07/13 17:28:31 (11 years ago)

Author:

tbretz

Message:

The macro can now handle files with 82 entries in the 'temp' column as well as 160.

File:

• trunk/Mars/fact/processing/drstemp.C (modified) (3 diffs)

trunk/Mars/fact/processing/drstemp.C

@@ -1 @@
-void drstemp(const char *fname, double beg=0, double end=100000)
+void drstemp82(fits &file, double beg, double end)
 {
-    fits file(fname);
-
-    //file.PrintColumns();
-    //file.PrintKeys();
-
-    Double_t time;
-    Float_t temp[82];
-    file.SetPtrAddress("Time", &time);
-    file.SetPtrAddress("temp", temp);
-
     UInt_t offset = file.GetUInt("MJDREF");
     if (beg < 30000)
@@ -16 +6 @@
     if (end < 30000)
         end+=offset;
+
+    double time;
+    file.SetPtrAddress("Time", &time);
+
+    float *temp = file.SetPtrAddress("temp");
 
     double avgmin = 0;
@@ -90 +85 @@
     cout << "result " << avgmin << " " << avgmax << " " << rmsmin << " " << rmsmax << endl;
 }
+
+void drstemp160(fits &file, double beg, double end)
+{
+    UInt_t offset = file.GetUInt("MJDREF");
+    if (beg < 30000)
+        beg+=offset;
+    if (end < 30000)
+        end+=offset;
+
+    double time;
+    file.SetPtrAddress("Time", &time);
+
+    float *temp = file.SetPtrAddress("temp");
+
+    double avg = 0;
+    double rms = 0;
+
+    int cnt = 0;
+
+    double lastavg = -1;
+    double lastrms = -1;
+
+    double diff  = -1;
+
+    while (file.GetNextRow())
+    {
+        time += offset;
+
+        if (time>end)
+            break;
+
+        double mn = 0;
+        double rn = 0;
+        for (int i=0; i<=160; i++)
+        {
+            mn += temp[i];
+            rn += temp[i]*temp[i];
+        }
+        mn /= 160;
+        rn /= 160;
+
+        if (time<beg)
+        {
+            lastavg = mn;
+            lastrms = sqrt(rn-mn*mn);
+            diff    = beg-time;
+            continue;
+        }
+
+        avg += mn;
+        rms += sqrt(rn-mn*mn);
+
+        cnt ++;
+    }
+
+    if (cnt==0)
+    {
+        if (diff<5./24/3600)
+            return;
+
+        cout << "result " << lastavg << " " << lastavg << " " << lastrms << " " << lastrms << endl;
+        return;
+    }
+
+    avg /= cnt;
+    rms /= cnt;
+
+    cout << "result " << avg << " " << avg << " " << rms << " " << rms << endl;
+}
+
+void drstemp(const char *fname, double beg=0, double end=100000)
+{
+    fits file(fname);
+
+    //file.PrintColumns();
+    //file.PrintKeys();
+
+    float temp[160];
+    if (file.SetPtrAddress("temp", temp, 82))
+    {
+        drstemp82(file, beg, end);
+        return;
+    }
+
+    file.SetPtrAddress("temp", temp, 160);
+    drstemp160(file, beg, end);
+}