OutputDataReadReflectometry Class Reference

Class for reading reflectometry data from ASCII file. More...

Public Member Functions
OutputData< double > *	readOutputData (std::istream &input_stream)

Detailed Description

Class for reading reflectometry data from ASCII file.

Definition at line 30 of file OutputDataReadReflectometry.h.

Member Function Documentation

readOutputData()

OutputData< double > * OutputDataReadReflectometry::readOutputData

(

std::istream &

input_stream

)

Definition at line 24 of file OutputDataReadReflectometry.cpp.

{
    OutputData<double>* oData = new OutputData<double>();
    std::string line;
    std::vector<std::vector<double>> vecVec;
    std::map<double, double> QvsR;
    std::map<double, double> QvsDR;
    std::map<double, double> QvsDQ;
 
    // Read numbers from file:
    while (std::getline(fin, line)) {
        line = StringUtils::trim(line);
        try {
            // #bamigration +++ this works only if separator is space or tab; it does not
            // work e.g. with comma or semicolon
            std::vector<double> rowVec = DataFormatUtils::parse_doubles(line);
            vecVec.push_back(rowVec);
        } catch (...) { // #bamigration +++ This eats useful errors away...
            continue;
        }
    }
 
    // validate - There is at least one row and at least two columns
    size_t nrows = vecVec.size();
    if (nrows < 1)
        throw std::runtime_error("No numerical values found");
    size_t ncols = vecVec[0].size();
    if (ncols < 2)
        throw std::runtime_error("Minimum 2 columns required");
 
    // Assign Q vs R, dR, dQ:
    for (size_t row = 0; row < nrows; row++) {
        if (vecVec[row].size() != ncols)
            throw std::runtime_error("The number of columns varies among the rows");
        double Q = vecVec[row][0];
        switch (ncols) {
        case 1:
            break;
        case 2:
            QvsR[Q] = vecVec[row][1];
            QvsDR[Q] = 0;
            QvsDQ[Q] = 0;
            break;
        case 3:
            QvsR[Q] = vecVec[row][1];
            QvsDR[Q] = vecVec[row][2];
            QvsDQ[Q] = 0;
            break;
        default:
            QvsR[Q] = vecVec[row][1];
            QvsDR[Q] = vecVec[row][2];
            QvsDQ[Q] = vecVec[row][3];
            break;
        }
    }
 
    std::vector<double> qVec;
    std::vector<double> rVec;
    for (auto it = QvsR.begin(); it != QvsR.end(); ++it) {
        if (it->second <= 0)
            continue;
        qVec.push_back(it->first);
        rVec.push_back(it->second);
    }
 
    oData->addAxis(PointwiseAxis("qVector", qVec));
    oData->setRawDataVector(rVec);
    return oData;
}

References OutputData< T >::addAxis(), DataFormatUtils::parse_doubles(), OutputData< T >::setRawDataVector(), and StringUtils::trim().

Referenced by AutomaticDataLoader1D::processContents(), and IntensityDataIOFactory::readReflectometryData().

Here is the call graph for this function:

---

The documentation for this class was generated from the following files:

• /home/www/ba/Device/InputOutput/OutputDataReadReflectometry.h
• /home/www/ba/Device/InputOutput/OutputDataReadReflectometry.cpp