Changeset 11871 for trunk/Mars/mraw

Timestamp:

08/09/11 23:55:12 (13 years ago)

Author:

tbretz

Message:

added the possibility to read the pixel map from a file

Location:

trunk/Mars/mraw

Files:

• MRawFitsRead.cc (modified) (4 diffs)
• MRawFitsRead.h (modified) (2 diffs)

trunk/Mars/mraw/MRawFitsRead.cc

@@ -38 +38 @@
 #include "MRawFitsRead.h"
 
+#include <fstream>
+
 #include <TClass.h>
 
@@ -64 +66 @@
     : MRawFileRead(fname, name, title)
 {
+}
+
+Bool_t MRawFitsRead::LoadMap(const char *name)
+{
+    fPixelMap.resize(1440);
+
+    ifstream fin(name);
+
+    int l = 0;
+
+    string buf;
+    while (getline(fin, buf, '\n'))
+    {
+        if (l>1439)
+            break;
+
+        const TString bf = TString(buf.c_str()).Strip(TString::kBoth);
+        if (bf[0]=='#')
+            continue;
+
+        stringstream str(bf.Data());
+
+        int index, cbpx;
+
+        str >> index;
+        str >> cbpx;
+
+        const int c =  cbpx/1000;
+        const int b = (cbpx/100)%10;
+        const int p = (cbpx/10)%10;
+        const int x =  cbpx%10;
+
+        const int hw = x+p*9+b*36+c*360;
+        if (hw>1439)
+            break;
+
+        fPixelMap[hw] = index;
+        l++;
+    }
+
+    if (l!=1440)
+    {
+        gLog << err << "ERROR - Problems reading FACTmapV5.txt" << endl;
+        fPixelMap.resize(0);
+        return kFALSE;
+    }
+
+    return kTRUE;
 }
 
@@ -91 +141 @@
     fRawRunHeader->SetObservation("", "FACT");
     fRawRunHeader->SetRunInfo(0, fin.GetUInt("NIGHT"), fin.GetUInt("RUNID"));
-    fRawRunHeader->InitFact(fin.GetUInt("NPIX")/9, 9, fin.GetUInt("NROI"));
+    fRawRunHeader->InitFact(fin.GetUInt("NPIX")/9, 9, fin.GetUInt("NROI"), fPixelMap.size()==0?0:fPixelMap.data());
     fRawRunHeader->SetFormat(0xf172, fin.GetUInt("BLDVER"));
     fRawRunHeader->SetRunType(0/*data*/);
@@ -129 +179 @@
         return kFALSE;
 
-    fRawEvtData1->SetIndices();
+    fRawEvtData1->SetIndices(fRawRunHeader->GetPixAssignment());
 
     return kTRUE;

trunk/Mars/mraw/MRawFitsRead.h

@@ -11 +11 @@
 {
 private:
-    std::vector<UInt_t> fPCTime; //! Buffer
+    std::vector<UInt_t>   fPCTime;   //! Buffer
+    std::vector<UShort_t> fPixelMap; //! 
 
     istream *OpenFile(const char *filename);
@@ -24 +25 @@
     static Bool_t IsFits(const char *name);
 
+    Bool_t LoadMap(const char *name);
+
     ClassDef(MRawFitsRead, 0)   // Task to read the raw data binary file
 };