gmaterial.h

Go to the documentation of this file.

#pragma once
 
// gemc
#include <gemc/gbase/gbase.h>
 
// c++
#include <string>
#include <vector>
#include <numeric>

class GMaterial : public GBase<GMaterial> {
public:
    GMaterial(const std::string &system, std::vector<std::string> pars, const std::shared_ptr<GLogger> &logger);

    [[nodiscard]] virtual std::unique_ptr<GMaterial> clone() const {
        return std::make_unique<GMaterial>(*this);
    }

    virtual ~GMaterial() = default;
 
private:
    std::string system; 
    std::string name; 
    std::string description; 

    double density{}; 
    std::vector<std::string> components; 
    std::vector<double> amounts; 

    std::vector<double> photonEnergy; 
    std::vector<double> indexOfRefraction; 
    std::vector<double> absorptionLength; 
    std::vector<double> reflectivity; 
    std::vector<double> efficiency; 

    std::vector<double> fastcomponent; 
    std::vector<double> slowcomponent; 
 
    double scintillationyield{}; 
    double resolutionscale{}; 
    double fasttimeconstant{}; 
    double slowtimeconstant{}; 
    double yieldratio{}; 
    double birksConstant{}; 
 
    bool scintillationyieldSet{}; 
    bool resolutionscaleSet{}; 
    bool fasttimeconstantSet{}; 
    bool slowtimeconstantSet{}; 
    bool yieldratioSet{}; 
    bool birksConstantSet{}; 

    std::vector<double> rayleigh; 

    void setComponentsFromString(const std::string &composition);

    void getMaterialPropertyFromString(const std::string &parameter, const std::string &propertyName);

    friend std::ostream &operator<<(std::ostream &stream, const GMaterial &);

    bool assign_if_set(const std::vector<std::string> &pars, size_t &i, double &out);
 
public:
    [[nodiscard]] std::string getName() const { return name; }
    [[nodiscard]] std::string getDescription() const { return description; }
    [[nodiscard]] double getDensity() const { return density; }
    [[nodiscard]] std::vector<std::string> getComponents() const { return components; }
    [[nodiscard]] std::vector<double> getAmounts() const { return amounts; }

    [[nodiscard]] std::vector<double> getPhotonEnergy() const { return photonEnergy; }
    [[nodiscard]] std::vector<double> getIndexOfRefraction() const { return indexOfRefraction; }
    [[nodiscard]] std::vector<double> getAbsorptionLength() const { return absorptionLength; }
    [[nodiscard]] std::vector<double> getReflectivity() const { return reflectivity; }
    [[nodiscard]] std::vector<double> getEfficiency() const { return efficiency; }

    [[nodiscard]] std::vector<double> getFastComponent() const { return fastcomponent; }
    [[nodiscard]] std::vector<double> getSlowComponent() const { return slowcomponent; }
    [[nodiscard]] double getScintillationYield() const { return scintillationyield; }
    [[nodiscard]] double getResolutionScale() const { return resolutionscale; }
    [[nodiscard]] double getFasttimeConstant() const { return fasttimeconstant; }
    [[nodiscard]] double getSlowtimeConstant() const { return slowtimeconstant; }
    [[nodiscard]] double getYieldratio() const { return yieldratio; }
    [[nodiscard]] double getBirksConstant() const { return birksConstant; }
    [[nodiscard]] bool hasScintillationYield() const { return scintillationyieldSet; }
    [[nodiscard]] bool hasResolutionScale() const { return resolutionscaleSet; }
    [[nodiscard]] bool hasFasttimeConstant() const { return fasttimeconstantSet; }
    [[nodiscard]] bool hasSlowtimeConstant() const { return slowtimeconstantSet; }
    [[nodiscard]] bool hasYieldratio() const { return yieldratioSet; }
    [[nodiscard]] bool hasBirksConstant() const { return birksConstantSet; }

    [[nodiscard]] std::vector<double> getRayleigh() const { return rayleigh; }

    [[nodiscard]] bool isChemicalFormula() const {
        return std::accumulate(amounts.begin(), amounts.end(), 0.0) > 1;
    }
};