interaction_data.h

Go to the documentation of this file.

 
#ifndef MC_XRAY_TRANSPORT_PHYSICS_ENGINE_DATA_SERVICE_H
#define MC_XRAY_TRANSPORT_PHYSICS_ENGINE_DATA_SERVICE_H
 
#include "data_access_object.h"
#include "interpolators.h"
#include "material.h"
#include "config.h"
#include <iostream>
#include <memory>
#include <map>
 
namespace {
    // anonymous namespace. only accessible in this file
    const std::string MIDSX_DB_DIR_VAR = MIDSX_DB_DIR;
    const std::string MIDSX_DB_NAME = "midsx.db";
    const std::string MIDSX_DB_PATH = MIDSX_DB_DIR_VAR + "/" + MIDSX_DB_NAME;
}
 
class InteractionData {
public:
    explicit InteractionData(const std::vector<std::string>& material_names);
 
    Material& getMaterialFromId(int id) { return material_map_.at(id); }
 
    Eigen::Matrix<double, Eigen::Dynamic, 2> getMaxTotalCrossSectionsMatrix() const { return max_total_cs_matrix_; }
 
    double interpolateMaxTotalCrossSection(double energy) const { return (max_total_cs_interpolator_)(energy); }
 
    std::string getAnyMaterialNameFromId(int id);
 
    uint8_t getAnyMaterialIdFromName(std::string name);
private:
    DataAccessObject dao_;
    std::vector<std::string> material_names_;
    std::map<int, Material> material_map_;
    Eigen::Matrix<double, Eigen::Dynamic, 2> max_total_cs_matrix_;
    Interpolator::LogLogLinear max_total_cs_interpolator_;
 
    void initializeData();
 
    void setMaterialMap();
 
    void setMaxTotalCrossSectionsAndInterpolator();
 
    Eigen::Matrix<double, Eigen::Dynamic, 2> getTotalMaxCrossSectionsMatrixFromInteractionData();
 
    void fillTotalMaxCrossSectionsMatrix(Eigen::MatrixXd& total_max_cross_sections_matrix, const Eigen::MatrixXd& merged_energy_matrix);
};
 
#endif //MC_XRAY_TRANSPORT_PHYSICS_ENGINE_DATA_SERVICE_H