LandauTable.cxx

Go to the documentation of this file.

/* Copyright (C) 2023 FIAS Frankfurt Institute for Advanced Studies, Frankfurt / Main
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Felix Weiglhofer [committer] */
#include "algo/data/sts/LandauTable.h"
 
#include "AlgoFairloggerCompat.h"
 
#include <fstream>
#include <sstream>
 
using namespace cbm::algo;
 
sts::LandauTable sts::LandauTable::FromFile(fs::path path)
{
  sts::LandauTable table;
 
  if (!fs::exists(path)) {
    std::stringstream msg;
    msg << "sts::LandauTable: file " << path.string() << " does not exist";
    throw std::runtime_error(msg.str());
  }
 
  std::vector<f32> charge;
  std::vector<f32> prob;
  std::ifstream file(path.string());
 
  while (!file.eof()) {
 
    f32 q, p;
    file >> q >> p;
    charge.push_back(q);
    prob.push_back(p);
 
    L_(trace) << "Reading Landau table " << path << " q=" << q << " p=" << p;
  }
 
  // TODO: check if charge is monotonically increasing, also more than 2 entries
 
  table.stepSize = charge[1] - charge[0];
  table.values   = std::move(prob);
 
  return table;
}