DataSet.cpp

Go to the documentation of this file.

#include <goofit/DataSet.h>
#include <goofit/Error.h>
#include <goofit/Variable.h>

#include <climits>
#include <cstdlib>
#include <utility>

namespace GooFit {

DataSet::DataSet(const Observable &var, std::string n)
    : name(std::move(n))
    , observables({var}) {
    generateName();
}

DataSet::DataSet(const std::vector<Observable> &vars, std::string n)
    : name(std::move(n))
    , observables(vars) {
    generateName();
}

DataSet::DataSet(const std::set<Observable> &vars, std::string n)
    : name(std::move(n))
    , observables(std::begin(vars), std::end(vars)) {
    generateName();
}

DataSet::DataSet(std::initializer_list<Observable> vars, std::string n)
    : name(std::move(n))
    , observables(vars) {
    generateName();
}

void DataSet::addWeightedEvent(fptype) {
    throw GooFit::GeneralError("AddWeightedEvent not implemented for this type of DataSet");
}

std::vector<fptype> DataSet::getCurrentValues() const {
    std::vector<fptype> values;

    for(const Observable &v : observables) {
        values.push_back(v.getValue());
    }

    return values;
}

const std::vector<Observable> &DataSet::getObservables() const { return observables; }

size_t DataSet::indexOfVariable(const Observable &var) const {
    for(size_t i = 0; i < observables.size(); ++i)
        if(var == observables[i])
            return i;

    throw GooFit::GeneralError("Invalid variable access into dataset!");
}

void DataSet::generateName() {
    // Create default name as list of variables.
    if(name != "")
        return;

    for(const Observable &v : observables) {
        if(v != observables[0])
            name += ", ";
        name += v.getName();
    }
}

void DataSet::checkAllVars() const {
    for(const Observable &v : observables) {
        if(!v)
            throw GooFit::OutOfRange(v.getName(), v.getValue(), v.getLowerLimit(), v.getUpperLimit());
    }
}
} // namespace GooFit