Index: unk/FACT++/src/PixelMap.cc =================================================================== --- /trunk/FACT++/src/PixelMap.cc (revision 15212) +++ (revision ) @@ -1,150 +1,0 @@ -#include "PixelMap.h" - -#ifdef __EXCEPTIONS -#include -#endif - -#ifdef __MARS__ -#include "MLog.h" -#include "MLogManip.h" -#endif - -using namespace std; - -const PixelMapEntry PixelMap::empty = { -1, 0, 0, 0, 0 }; -const BiasMapEntry BiasMap::empty = { -1, 0, 0, 0 }; - -bool BiasMap::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; - -#ifdef __EXCEPTIONS - if (entry.hv_channel+32*entry.hv_board>=416) - throw runtime_error("Invalid board/channel read from "+fname+"."); -#endif -#ifdef __MARS__ - if (entry.hv_channel+32*entry.hv_board>=416) - { - gLog << err << "Invalid board/channel read from " << fname << "." << endl; - return false; - } -#endif - - (*this)[entry.hv()] = entry; - - l++; - } - -#ifdef __EXCEPTIONS - if (l!=416) - throw runtime_error("Number of lines read from "+fname+" does not match 416."); - - if (size()!=416) - throw runtime_error("Number of entries read from "+fname+" does not match 416."); -#endif - -#ifdef __MARS__ - if (l!=416) - { - gLog << err << "Number of lines read from " << fname << " does not match 416." << endl; - return false; - } - - if (size()!=416) - { - gLog << "Number of entries read from " << fname << " does not match 416." << endl; - return false; - } -#endif - - return true; -} - -#ifndef __MARS__ -#include -#include - -void BiasMap::Retrieve(const std::string &database) -{ - static const boost::regex expr("(([[:word:].-]+)(:(.+))?@)?([[:word:].-]+)(:([[:digit:]]+))?(/([[:word:].-]+))"); - // 2: user - // 4: pass - // 5: server - // 7: port - // 9: db - - boost::smatch what; - if (!boost::regex_match(database, what, expr, boost::match_extra)) - throw runtime_error("Couldn't parse '"+database+"'."); - - if (what.size()!=10) - throw runtime_error("Error parsing '"+database+"'."); - - const string user = what[2]; - const string passwd = what[4]; - const string server = what[5]; - const string db = what[9]; - const int port = atoi(string(what[7]).c_str()); - - mysqlpp::Connection conn(db.c_str(), server.c_str(), user.c_str(), passwd.c_str(), port); - - const mysqlpp::StoreQueryResult res = - conn.query("SELECT fPatchNumber, AVG(fVoltageNom), fOffset " - " FROM GapdVoltages " - " LEFT JOIN BiasOffsets USING(fPatchNumber) " - " GROUP BY fPatchNumber").store(); - - clear(); - - int l = 0; - for (vector::const_iterator v=res.begin(); v416) - { - ostringstream str; - str << "Invalid channel id " << id << " received from database."; - throw runtime_error(str.str()); - } - - BiasMapEntry entry; - - entry.hv_board = id/32; - entry.hv_channel = id%32; - entry.Vnom = (*v)[1]; - entry.Voff = (*v)[2]; - - (*this)[id] = entry; - - l++; - } - - if (l!=416) - throw runtime_error("Number of rows retrieved from the database does not match 416."); - - if (size()!=416) - throw runtime_error("Number of entries retrived from database does not match 416."); -} -#endif Index: unk/FACT++/src/PixelMap.h =================================================================== --- /trunk/FACT++/src/PixelMap.h (revision 15212) +++ (revision ) @@ -1,265 +1,0 @@ -#ifndef FACT_PixelMap -#define FACT_PixelMap - -#include -#include -#include -#include - -#include "tools.h" - -#ifdef DEBUG -#include // cerr -- to be removed? -#endif - -// FIXME: Replace 416 by BIAS::kNumChannels - -struct PixelMapEntry -{ - int index; /// Software index - int cbpx; /// Hardware index as CBPX - int gapd; /// gAPD index -// float Vgapd; /// gAPD Bias voltage - int hv_board; /// Bias suppply board - int hv_channel; /// Bias supply channel - - int crate() const { return cbpx/1000; } - int board() const { return (cbpx/100)%10; } - int patch() const { return (cbpx/10)%10; } - int pixel() const { return cbpx%10; } - int hw() const { return pixel()+patch()*9+board()*36+crate()*360; } - int group() const { return pixel()>3; } - int hv() const { return hv_channel+hv_board*32; } - - operator bool() const { return index>=0; } -}; - -class PixelMap : public std::vector -{ -public: - static const PixelMapEntry empty; - - PixelMap() : std::vector(1440) - { - } - - bool Read(const std::string &fname) - { - std::ifstream fin(fname); - - int l = 0; - - std::string buf; - while (getline(fin, buf, '\n')) - { - if (l>1439) - break; - - buf = Tools::Trim(buf); - if (buf[0]=='#') - continue; - - std::stringstream str(buf); - - int idummy; - float fdummy; - - PixelMapEntry entry; - - str >> entry.index; - str >> entry.cbpx; - str >> idummy; - str >> idummy; - str >> entry.gapd; - str >> fdummy; //entry.Vgapd; - str >> entry.hv_board; - str >> entry.hv_channel; - //str >> fdummy; - //str >> fdummy; - //str >> fdummy; - - if (entry.hv_channel+32*entry.hv_board>=416/*BIAS::kNumChannels*/) - { -#ifdef DEBUG - cerr << "Invalid board/channel read from FACTmapV5.txt." << endl; -#endif - return false; - } - - (*this)[l++] = entry; - } - - return l==1440; - } - - const PixelMapEntry &index(int idx) const - { - for (std::vector::const_iterator it=begin(); it!=end(); it++) - if (it->index==idx) - return *it; -#ifdef DEBUG - std::cerr << "PixelMap: index " << idx << " not found" << std::endl; -#endif - return empty; - } - - const PixelMapEntry &cbpx(int c) const - { - for (std::vector::const_iterator it=begin(); it!=end(); it++) - if (it->cbpx==c) - return *it; -#ifdef DEBUG - std::cerr << "PixelMap: cbpx " << c << " not found" << std::endl; -#endif - return empty; - } - - const PixelMapEntry &cbpx(int c, int b, int p, int px) const - { - return cbpx(px + p*10 + b*100 + c*1000); - } - - const PixelMapEntry &hw(int idx) const - { - return cbpx(idx/360, (idx/36)%10, (idx/9)%4, idx%9); - } - - const PixelMapEntry &hv(int board, int channel) const - { - for (std::vector::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 PixelMapEntry &hv(int idx) const - { - return hv(idx/32, idx%32); - } - - /* - float Vgapd(int board, int channel) const - { - float avg = 0; - int num = 0; - - for (std::vector::const_iterator it=begin(); it!=end(); it++) - if (it->hv_board==board && it->hv_channel==channel) - { - avg += it->Vgapd; - num ++; - } - - return num==0 ? 0 : avg/num; - } - - float Vgapd(int idx) const - { - return Vgapd(idx/32, idx%32); - } - - std::vector Vgapd() const - { - std::vector avg(416); - std::vector num(416); - - for (std::vector::const_iterator it=begin(); it!=end(); it++) - { - const int ch = it->hv_board*32 + it->hv_channel; - - avg[ch] += it->Vgapd; - num[ch] ++; - } - - for (int ch=0; ch<416; ch++) - { - if (num[ch]) - avg[ch] /= num[ch]; - } - - return avg; - }*/ -}; - -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 -{ -public: - static const BiasMapEntry empty; - - BiasMap() : std::vector(416) - { - } - -#ifndef __MARS__ - void Retrieve(const std::string &database); -#endif - bool Read(const std::string &fname); - - const BiasMapEntry &hv(int board, int channel) const - { - for (std::vector::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; // use this with GAPDmap_20111126.txt - } - - float Vgapd(int idx) const - { - return Vgapd(idx/32, idx%32); - }*/ - - std::vector Vgapd() const - { - std::vector volt(416); - - for (std::vector::const_iterator it=begin(); it!=end(); it++) - { - const int ch = it->hv_board*32 + it->hv_channel; - volt[ch] = it->Vnom; - } - - return volt; - } - std::vector Voffset() const - { - std::vector volt(416); - - for (std::vector::const_iterator it=begin(); it!=end(); it++) - { - const int ch = it->hv_board*32 + it->hv_channel; - volt[ch] = -it->Voff; - } - - return volt; - } -}; - -#endif