Changeset 11810 for trunk

Timestamp:

08/06/11 16:05:30 (13 years ago)

Author:

tbretz

Message:

Added reading of FITS files.

File:

• trunk/FACT++/src/DataCalib.cc (modified) (2 diffs)

trunk/FACT++/src/DataCalib.cc

@@ -4 +4 @@
 #include "FitsFile.h"
 #include "DimDescriptionService.h"
+
+#include "externals/fits.h"
 
 using namespace std;
@@ -108 +110 @@
 
     return true;
+}
+
+bool DataCalib::ReadFits(const string &str)
+{
+    try
+    {
+        fits file(str);
+        if (!file)
+            return false;
+
+        if (file.GetStr("TELESCOP")!="FACT")
+        {
+            Error("Reading "+str+" failed: Not a valid FACT file (TELESCOP not FACT in header)");
+            return false;
+        }
+
+        if (!file.HasKey("STEP"))
+        {
+            Error("Reading "+str+" failed: Is not a DRS calib file (STEP not found in header)");
+            return false;
+        }
+
+        if (file.GetNumRows()!=1)
+        {
+            Error("Reading "+str+" failed: Number of rows in table is not 1.");
+            return false;
+        }
+
+        vector<float> data(1440*1024*6+3);
+
+        float *base = data.data();
+
+        file.SetPtrAddress("RunNumberBaseline",      base,   sizeof(float));
+        file.SetPtrAddress("RunNumberGain",          base+1, sizeof(float));
+        file.SetPtrAddress("RunNumberTriggerOffset", base+2, sizeof(float));
+        file.SetPtrAddress("BaselineMean",           base+0*1024*1440+3, sizeof(float)*1024*1440);
+        file.SetPtrAddress("BaselineRms",            base+1*1024*1440+3, sizeof(float)*1024*1440);
+        file.SetPtrAddress("GainMean",               base+2*1024*1440+3, sizeof(float)*1024*1440);
+        file.SetPtrAddress("GainRms",                base+3*1024*1440+3, sizeof(float)*1024*1440);
+        file.SetPtrAddress("TriggerOffsetMean",      base+4*1024*1440+3, sizeof(float)*1024*1440);
+        file.SetPtrAddress("TriggerOffsetRms",       base+5*1024*1440+3, sizeof(float)*1024*1440);
+
+        if (!file.GetNextRow())
+        {
+            Error("Reading data from "+str+" failed.");
+            return false;
+        }
+
+        fStep = file.GetUInt("STEP");
+
+        memcpy(fStats.data(), base, fStats.size()*sizeof(float));
+
+        return true;
+    }
+    catch (const runtime_error &e)
+    {
+        Error("Exception reading "+str+": "+e.what());
+        return false;
+    }
 }