pdscalibdata/src/CsPadPedestalsV1.cpp

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//      $Id: CsPadPedestalsV1.cpp 11788 2016-04-27 17:04:41Z dubrovin@SLAC.STANFORD.EDU $
//
// Description:
//      Class CsPadPedestalsV1...
//
// Author List:
//      Andrei Salnikov
//
//------------------------------------------------------------------------

//-----------------------
// This Class's Header --
//-----------------------
#include "pdscalibdata/CsPadPedestalsV1.h"

//-----------------
// C/C++ Headers --
//-----------------
#include <algorithm>
#include <stdexcept>
#include <fstream>

//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "MsgLogger/MsgLogger.h"

//-----------------------------------------------------------------------
// Local Macros, Typedefs, Structures, Unions and Forward Declarations --
//-----------------------------------------------------------------------

//              ----------------------------------------
//              -- Public Function Member Definitions --
//              ----------------------------------------

namespace pdscalibdata {

//----------------
// Constructors --
//----------------
CsPadPedestalsV1::CsPadPedestalsV1 ()
{
  // fill all pedestals with zeros
  std::fill_n(m_pedestals, int(Size), pedestal_t(0));
}

CsPadPedestalsV1::CsPadPedestalsV1 (const std::string& fname) 
{ 
  // open file
  std::ifstream in(fname.c_str());
  if (not in.good()) {
    const std::string msg = "Failed to open pedestals file: "+fname;
    MsgLogRoot(error, msg);
    throw std::runtime_error(msg);
  }

  // read all numbers
  pedestal_t* it = m_pedestals;
  size_t count = 0;

  //===========================
  /*
  while(in and count != Size) {
    in >> *it++;
    ++ count;
  }
  */
  //===========================

  pedestal_t  val;
  std::string str; 
  while(getline(in,str)) { 
        // 1. skip lines with comments (# in the 1st position)
        if(str[0] == '#') continue;

        // 2. skip empty lines 
        else if (str.find_first_not_of(" ")==std::string::npos) continue; 

        // 3. parse one line
        else {
               std::stringstream ss(str);
               while (ss >> val and count != Size) { 
                 *it++ = val;
                 ++ count;
               }
        }
  }

  //===========================

  // check that we read whole array
  if (count != Size) {
    const std::string msg = "Pedestals file does not have enough data: "+fname;
    MsgLogRoot(error, msg);
    throw std::runtime_error(msg);
  }

  // and no data left after we finished reading
  float tmp ;
  if ( in >> tmp ) {
    ++ count;
    const std::string msg = "Pedestals file has extra data: "+fname;
    MsgLogRoot(error, msg);
    MsgLogRoot(error, "read " << count << " numbers, expecting " << Size );
    throw std::runtime_error(msg);
  }
}

//--------------
// Destructor --
//--------------
CsPadPedestalsV1::~CsPadPedestalsV1 ()
{
}

} // namespace pdscalibdata

---