Changeset 12507 for trunk/FACT++/src/PixelMap.h

Timestamp:

11/13/11 18:06:05 (13 years ago)

Author:

tbretz

Message:

Removed the GAPD voltage from the PixelMap and implemented a BiasMap class instead wich content are read from a special file containing the nominal voltages and the offsets.

File:

• trunk/FACT++/src/PixelMap.h (modified) (5 diffs)

trunk/FACT++/src/PixelMap.h

@@ -20 +20 @@
     int   cbpx;                /// Hardware index as CBPX
     int   gapd;                /// gAPD index
-    float Vgapd;               /// gAPD Bias voltage
+//    float Vgapd;               /// gAPD Bias voltage
     int   hv_board;            /// Bias suppply board
     int   hv_channel;          /// Bias supply channel
@@ -62 +62 @@
             std::stringstream str(buf);
 
-            int   idummy;
+            int     idummy;
+            float   fdummy;
 
             PixelMapEntry entry;
@@ -71 +72 @@
             str >> idummy;
             str >> entry.gapd;
-            str >> entry.Vgapd;
+            str >> fdummy; //entry.Vgapd;
             str >> entry.hv_board;
             str >> entry.hv_channel;
@@ -140 +141 @@
     }
 
+    /*
     float Vgapd(int board, int channel) const
     {
@@ -180 +182 @@
 
         return avg;
-    }
-};
-
-#endif
+    }*/
+};
+
+struct BiasMapEntry
+{
+    int   hv_board;            /// Bias suppply board
+    int   hv_channel;          /// Bias supply channel
+    float Vnom;                /// Channel bias voltage nominal
+    float Voff;                /// Channel bias voltage offset
+
+    int hv() const { return hv_channel+hv_board*32; }
+};
+
+class BiasMap : public std::vector<BiasMapEntry>
+{
+public:
+    static const BiasMapEntry empty;
+
+    BiasMap() : std::vector<BiasMapEntry>(416)
+    {
+    }
+
+    bool Read(const std::string &fname)
+    {
+        std::ifstream fin(fname);
+
+        int l = 0;
+
+        std::string buf;
+        while (getline(fin, buf, '\n'))
+        {
+            if (l>416)
+                break;
+
+            buf = Tools::Trim(buf);
+            if (buf[0]=='#')
+                continue;
+
+            std::stringstream str(buf);
+
+            BiasMapEntry entry;
+
+            str >> entry.hv_board;
+            str >> entry.hv_channel;
+            str >> entry.Vnom;
+            str >> entry.Voff;
+
+            if (entry.hv_channel+32*entry.hv_board>=416)
+            {
+#ifdef DEBUG
+                cerr << "Invalid board/channel read from " << fname << "." << endl;
+#endif
+                return false;
+            }
+
+            (*this)[l++] = entry;
+        }
+
+        return l==1440;
+    }
+
+    const BiasMapEntry &hv(int board, int channel) const
+    {
+        for (std::vector<BiasMapEntry>::const_iterator it=begin(); it!=end(); it++)
+            if (it->hv_board==board && it->hv_channel==channel)
+                return *it;
+#ifdef DEBUG
+        std::cerr << "PixelMap: hv " << board << "/" << channel << " not found" << std::endl;
+#endif
+        return empty;
+    }
+
+    const BiasMapEntry &hv(int idx) const
+    {
+        return hv(idx/32, idx%32);
+    }
+
+    float Vgapd(int board, int channel) const
+    {
+        const BiasMapEntry &entry = hv(board, channel);
+        return entry.Vnom+entry.Voff;
+    }
+
+    float Vgapd(int idx) const
+    {
+        return Vgapd(idx/32, idx%32);
+    }
+
+    std::vector<float> Vgapd() const
+    {
+        std::vector<float> volt(416);
+
+        for (std::vector<BiasMapEntry>::const_iterator it=begin(); it!=end(); it++)
+        {
+            const int ch = it->hv_board*32 + it->hv_channel;
+            volt[ch] += it->Vnom+it->Voff;
+        }
+
+        return volt;
+    }
+};
+
+#endif